Exon-level estimates improve the detection of differentially expressed genes in RNA-seq studies

Transcriptome analyses reveal key roles of alternative splicing regulation in atlantic salmon during the infectious process of Piscirickettsiosis disease

In the Chilean salmon farming industry, infection by Piscirickettsia salmonis is the primary cause of the main bacterial disease known as Piscirickettsiosis, which has an overwhelming economic impact. Although it has been demonstrated that Piscirickettsiosis modifies the expression of numerous salmonids genes, it is yet unknown how alternative splicing (AS) contributes to salmonids bacterial infection. AS, has the potential to create heterogeneity at the protein and RNA levels and has been associated as a relevant molecular mechanism in the immune response of eukaryotes to several diseases. In this study, we used RNA data to survey P. salmonis-induced modifications in the AS of Atlantic salmon and found that P. salmonis infection promoted a substantial number (158,668) of AS events. Differentially spliced genes (DSG) sensitive to Piscirickettsiosis were predominantly enriched in genes involved in RNA processing, splicing and spliceosome processes (e.g., hnRNPm, hnRPc, SRSF7, SRSF45), whereas among the DSG of resistant and susceptible to Piscirickettsiosis, several metabolic and immune processes were found, most notably associated to the regulation of GTPase, lysosome and telomere organization-maintenance. Furthermore, we found that DSG were mostly not differentially expressed (5-7 %) and were implicated in distinct biological pathways. Therefore, our results underpin AS achieving a significant regulatory performance in the response of salmonids to Piscirickettsiosis.

Transcriptomics for Clinical and Experimental Biology Research: Hang on a Seq

Sequencing the human genome empowers translational medicine, facilitating transcriptome-wide molecular diagnosis, pathway biology, and drug repositioning. Initially, microarrays are used to study the bulk transcriptome; but now short-read RNA sequencing (RNA-seq) predominates. Positioned as a superior technology, that makes the discovery of novel transcripts routine, most RNA-seq analyses are in fact modeled on the known transcriptome. Limitations of the RNA-seq methodology have emerged, while the design of, and the analysis strategies applied to, arrays have matured. An equitable comparison between these technologies is provided, highlighting advantages that modern arrays hold over RNA-seq. Array protocols more accurately quantify constitutively expressed protein coding genes across tissue replicates, and are more reliable for studying lower expressed genes. Arrays reveal long noncoding RNAs (lncRNA) are neither sparsely nor lower expressed than protein coding genes. Heterogeneous coverage of constitutively expressed genes observed with RNA-seq, undermines the validity and reproducibility of pathway analyses. The factors driving these observations, many of which are relevant to long-read or single-cell sequencing are discussed. As proposed herein, a reappreciation of bulk transcriptomic methods is required, including wider use of the modern high-density array data-to urgently revise existing anatomical RNA reference atlases and assist with more accurate study of lncRNAs.