A Comparison of Low Read Depth QuantSeq 3' Sequencing to Total RNA-Seq in FUS Mutant Mice

Varying conjunctival immune response adaptations of house finch populations to a rapidly evolving bacterial pathogen

Pathogen adaptations during host-pathogen co-evolution can cause the host balance between immunity and immunopathology to rapidly shift. However, little is known in natural disease systems about the immunological pathways optimised through the trade-off between immunity and self-damage. The evolutionary interaction between the conjunctival bacterial infection Mycoplasma gallisepticum (MG) and its avian host, the house finch (Haemorhous mexicanus), can provide insights into such adaptations in immune regulation. Here we use experimental infections to reveal immune variation in conjunctival tissue for house finches captured from four distinct populations differing in the length of their co-evolutionary histories with MG and their disease tolerance (defined as disease severity per pathogen load) in controlled infection studies. To differentiate contributions of host versus pathogen evolution, we compared house finch responses to one of two MG isolates: the original VA1994 isolate and a more evolutionarily derived one, VA2013. To identify differential gene expression involved in initiation of the immune response to MG, we performed 3'-end transcriptomic sequencing (QuantSeq) of samples from the infection site, conjunctiva, collected 3-days post-infection. In response to MG, we observed an increase in general pro-inflammatory signalling, as well as T-cell activation and IL17 pathway differentiation, associated with a decrease in the IL12/IL23 pathway signalling. The immune response was stronger in response to the evolutionarily derived MG isolate compared to the original one, consistent with known increases in MG virulence over time. The host populations differed namely in pre-activation immune gene expression, suggesting population-specific adaptations. Compared to other populations, finches from Virginia, which have the longest co-evolutionary history with MG, showed significantly higher expression of anti-inflammatory genes and Th1 mediators. This may explain the evolution of disease tolerance to MG infection in VA birds. We also show a potential modulating role of BCL10, a positive B- and T-cell regulator activating the NFKB signalling. Our results illuminate potential mechanisms of house finch adaptation to MG-induced immunopathology, contributing to understanding of the host evolutionary responses to pathogen-driven shifts in immunity-immunopathology trade-offs.

Peripheral inflammation-induced changes in songbird brain gene expression: 3' mRNA transcriptomic approach

Species-specific neural inflammation can be induced by profound immune signalling from periphery to brain. Recent advances in transcriptomics offer cost-effective approaches to study this regulation. In a population of captive zebra finch (Taeniopygia guttata), we compare the differential gene expression patterns in lipopolysaccharide (LPS)-triggered peripheral inflammation revealed by RNA-seq and QuantSeq. The RNA-seq approach identified more differentially expressed genes but failed to detect any inflammatory markers. In contrast, QuantSeq results identified specific expression changes in the genes regulating inflammation. Next, we adopted QuantSeq to relate peripheral and brain transcriptomes. We identified subtle changes in the brain gene expression during the peripheral inflammation (e.g. up-regulation in AVD-like and ACOD1 expression) and detected co-structure between the peripheral and brain inflammation. Our results suggest benefits of the 3'end transcriptomics for association studies between peripheral and neural inflammation in genetically heterogeneous models and identify potential targets for the future brain research in birds.

3' RNA-seq is superior to standard RNA-seq in cases of sparse data but inferior at identifying toxicity pathways in a model organism

Introduction: The application of RNA-sequencing has led to numerous breakthroughs related to investigating gene expression levels in complex biological systems. Among these are knowledge of how organisms, such as the vertebrate model organism zebrafish (Danio rerio), respond to toxicant exposure. Recently, the development of 3' RNA-seq has allowed for the determination of gene expression levels with a fraction of the required reads compared to standard RNA-seq. While 3' RNA-seq has many advantages, a comparison to standard RNA-seq has not been performed in the context of whole organism toxicity and sparse data. Methods and results: Here, we examined samples from zebrafish exposed to perfluorobutane sulfonamide (FBSA) with either 3' or standard RNA-seq to determine the advantages of each with regards to the identification of functionally enriched pathways. We found that 3' and standard RNA-seq showed specific advantages when focusing on annotated or unannotated regions of the genome. We also found that standard RNA-seq identified more differentially expressed genes (DEGs), but that this advantage disappeared under conditions of sparse data. We also found that standard RNA-seq had a significant advantage in identifying functionally enriched pathways via analysis of DEG lists but that this advantage was minimal when identifying pathways via gene set enrichment analysis of all genes. Conclusions: These results show that each approach has experimental conditions where they may be advantageous. Our observations can help guide others in the choice of 3' RNA-seq vs standard RNA sequencing to query gene expression levels in a range of biological systems.

Technical assessment of different extraction methods and transcriptome profiling of RNA isolated from small volumes of blood

Transcriptome profiling of human whole blood is used to discover biomarkers of diseases and to assess phenotypic traits. Recently, finger-stick blood collection systems have allowed a less invasive and quicker collection of peripheral blood. Such non-invasive sampling of small volumes of blood offers practical advantages. The quality of gene expression data is strictly dependent on the steps used for the sample collection, extraction, preparation and sequencing. Here we have: (i) compared the manual and automated RNA extraction of small volumes of blood using the Tempus Spin RNA isolation kit and the MagMAX for Stabilized Blood RNA Isolation kit , respectively; and (ii) assessed the effect of TURBO DNA Free treatment on the transcriptomic data of RNA isolated from small volumes of blood. We have used the QuantSeq 3' FWD mRNA-Seq Library Prep kit to prepare RNA-seq libraries, which were sequenced on the Illumina NextSeq 500 system. The samples isolated manually displayed a higher variability in the transcriptomic data as compared to the other samples. The TURBO DNA Free treatment affected the RNA samples negatively, decreasing the RNA yield and reducing the quality and reproducibility of the transcriptomic data. We conclude that automated extraction systems should be preferred over manual extraction systems for data consistency, and that the TURBO DNA Free treatment should be avoided when working on RNA samples isolated manually from small volumes of blood.

Influence of RNA-Seq library construction, sampling methods, and tissue harvesting time on gene expression estimation

RNA sequencing (RNA-Seq) is popular for measuring gene expression in non-model organisms, including wild populations. While RNA-Seq can detect gene expression variation among wild-caught individuals and yield important insights into biological function, sampling methods can also affect gene expression estimates. We examined the influence of multiple technical variables on estimated gene expression in a non-model fish, the westslope cutthroat trout (Oncorhynchus clarkii lewisi), using two RNA-Seq library types: 3' RNA-Seq (QuantSeq) and whole mRNA-Seq (NEB). We evaluated effects of dip netting versus electrofishing, and of harvesting tissue immediately versus 5 min after euthanasia on estimated gene expression in blood, gill, and muscle. We found no significant differences in gene expression between sampling methods or tissue collection times with either library type. When library types were compared using the same blood samples, 58% of genes detected by both NEB and QuantSeq showed significantly different expression between library types, and NEB detected 31% more genes than QuantSeq. Although the two library types recovered different numbers of genes and expression levels, results with NEB and QuantSeq were consistent in that neither library type showed differences in gene expression between sampling methods and tissue harvesting times. Our study suggests that researchers can safely rely on different fish sampling strategies in the field. In addition, while QuantSeq is more cost effective, NEB detects more expressed genes. Therefore, when it is crucial to detect as many genes as possible (especially low expressed genes), when alternative splicing is of interest, or when working with an organism lacking good genomic resources, whole mRNA-Seq is more powerful.

Application of the 3' mRNA-Seq using unique molecular identifiers in highly degraded RNA derived from formalin-fixed, paraffin-embedded tissue

Background

Archival formalin-fixed, paraffin-embedded (FFPE) tissue samples with clinical and histological data are a singularly valuable resource for developing new molecular biomarkers. However, transcriptome analysis remains challenging with standard mRNA-seq methods as FFPE derived-RNA samples are often highly modified and fragmented. The recently developed 3′ mRNA-seq method sequences the 3′ region of mRNA using unique molecular identifiers (UMI), thus generating gene expression data with minimal PCR bias. In this study, we evaluated the performance of 3′ mRNA-Seq using Lexogen QuantSeq 3′ mRNA-Seq Library Prep Kit FWD with UMI, comparing with TruSeq Stranded mRNA-Seq and RNA Exome Capture kit. The fresh-frozen (FF) and FFPE tissues yielded nucleotide sizes range from 13 to > 70% of DV200 values; input amounts ranged from 1 ng to 100 ng for validation.

Results

The total mapped reads of QuantSeq 3′ mRNA-Seq to the reference genome ranged from 99 to 74% across all samples. After PCR bias correction, 3 to 56% of total sequenced reads were retained. QuantSeq 3′ mRNA-Seq data showed highly reproducible data across replicates in Universal Human Reference RNA (UHR, R > 0.94) at input amounts from 1 ng to 100 ng, and FF and FFPE paired samples (R = 0.92) at 10 ng. Severely degraded FFPE RNA with ≤30% of DV200 value showed good concordance (R > 0.87) with 100 ng input. A moderate correlation was observed when directly comparing QuantSeq 3′ mRNA-Seq data with TruSeq Stranded mRNA-Seq (R = 0.78) and RNA Exome Capture data (R > 0.67).

Conclusion

In this study, QuantSeq 3′ mRNA-Seq with PCR bias correction using UMI is shown to be a suitable method for gene quantification in both FF and FFPE RNAs. 3′ mRNA-Seq with UMI may be applied to severely degraded RNA from FFPE tissues generating high-quality sequencing data.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-08068-1.