Long non-coding RNA as environmental regulator in a non-model fish

The majority of the transcribed genome does not have coding potential but is composed of non-coding transcripts that are involved in transcriptional and post-transcriptional regulation of protein-coding genes. Regulation of gene expression is important in determining the response of organisms to changes in the environment, and therefore their persistence as population or species under global change. However, long non-coding RNAs (lncRNAs) are scarcely studied especially in non-model organisms due to the lack of a reliable pipeline for their accurate identification and annotation. Here, we present a pipeline which uses a combination of alignment-dependent and independent methods for the identification of conserved and species-specific lncRNAs from RNA-Seq data. Validation of this pipeline was performed using existing RNA-Seq data from Acanthochromis polyacanthus brain tissue, identifying a total of 4,728 lncRNAs across the genome, the majority of which (3,272) are intergenic. To investigate the possible implications of these intergenic lncRNAs (lincRNAs), we estimated the expression changes of lincRNAs and coding genes in response to ocean acidification. We found lincRNAs which neighbour or possibly trans-regulate differentially expressed coding genes related to pH regulation, neural signal transduction and ion transport, which are known to be important in the response to ocean acidification in fish. Overall, this pipeline enables the use of existing RNA sequencing data to reveal additional underlying molecular mechanisms involved in the response to environmental changes by integrating the study of lncRNAs with gene expression.