Search

RNA-seq datasets can contain millions of intron reads per library that are typically removed from downstream analysis. Only reads overlapping annotated exons are considered to be informative since mature mRNA is assumed to be the major component sequenced, especially for poly(A) RNA libraries. In this study, we show that intron reads are informative, and through exploratory data analysis of read coverage that intron signal is representative of both pre-mRNAs and intron retention. We demonstrate how intron reads can be utilized in differential expression analysis using our index method where a unique set of differentially expressed genes can be detected using intron counts. In exploring read coverage, we also developed the superintronic software that quickly and robustly calculates user-defined summary statistics for exonic and intronic regions. Across multiple datasets, superintronic enabled us to identify several genes with distinctly retained introns that had similar coverage levels to that of neighbouring exons. The work and ideas presented in this paper is the first of its kind to consider multiple biological sources for intron reads through exploratory data analysis, minimizing bias in discovery and interpretation of results. Our findings open up possibilities for further methods development for intron reads and RNA-seq data in general.

Influenza A virus (IAV) is rapidly detected in the airways by the immune system, with resident parenchymal cells and leukocytes orchestrating viral sensing and the induction of antiviral inflammatory responses. The airways are innervated by heterogeneous populations of vagal sensory neurons which also play an important role in pulmonary defense. How these neurons respond to IAV respiratory infection remains unclear. Here, we use a murine model to provide the first evidence that vagal sensory neurons undergo significant transcriptional changes following a respiratory IAV infection. RNA sequencing on vagal sensory ganglia showed that IAV infection induced the expression of many genes associated with an antiviral and pro-inflammatory response and this was accompanied by a significant increase in inflammatory cell recruitment into the vagal ganglia. Assessment of gene expression in single-vagal sensory neurons confirmed that IAV infection induced a neuronal inflammatory phenotype, which was most prominent in bronchopulmonary neurons, and also evident in some neurons innervating other organs. The altered transcriptome could be mimicked by intranasal treatment with cytokines and the lung homogenates of infected mice, in the absence of infectious virus. These data argue that IAV pulmonary infection and subsequent inflammation induces vagal sensory ganglia neuroinflammation and this may have important implications for IAV-induced morbidity.