scholarly journals Comparing total RNA sequencing and metagenomics pipelines for multi-domain taxonomic profiling: implications for ecological assessments

2021 ◽  
Vol 4 ◽  
Author(s):  
Christopher Hempel ◽  
Julia Harvie ◽  
Jose Hleap Lozano ◽  
Natalie Wright ◽  
Sarah Adamowicz ◽  
...  

Ecological assessments are necessary to evaluate the status of our deteriorating ecosystems, however, assessment methods traditionally omit most microbes because unicellular organisms are challenging to identify. This omission is not ideal, as microbes might be better indicators for changes in environmental conditions than taxa traditionally used. DNA- and RNA-based techniques are increasingly applied for ecological assessments to overcome this challenge but require more testing and optimization. In this study, we compare metagenomics and total RNA sequencing (total RNA-Seq) for their taxonomic profiling performance for microbial communities. We applied both techniques on two sample sets, 1) a commercially available microbial mock community consisting of eight bacterial and two eukaryotic species, and 2) a display tank water sample. We processed the data using 1,532 bioinformatics pipelines and evaluated each workflow, i.e., the combination of sample type (metagenomics or total RNA-Seq) and pipeline, in terms of their accuracy and precision. This talk will showcase preliminary results and highlight differences in workflow performances. A recommended workflow to maximize taxonomic profiling accuracy of microbial communities will also be presented.

2019 ◽  
Author(s):  
Celine Everaert ◽  
Hetty Helsmoortel ◽  
Anneleen Decock ◽  
Eva Hulstaert ◽  
Ruben Van Paemel ◽  
...  

AbstractRNA profiling has emerged as a powerful tool to investigate the biomarker potential of human biofluids. However, despite enormous interest in extracellular nucleic acids, RNA sequencing methods to quantify the total RNA content outside cells are rare. Here, we evaluate the performance of the SMARTer Stranded Total RNA-Seq method in human platelet-rich plasma, platelet-free plasma, urine, conditioned medium, and extracellular vesicles (EVs) from these biofluids. We found the method to be accurate, precise, compatible with low-input volumes and able to quantify a few thousand genes. We picked up distinct classes of RNA molecules, including mRNA, lncRNA, circRNA, miscRNA and pseudogenes. Notably, the read distribution and gene content drastically differ among biofluids. In conclusion, we are the first to show that the SMARTer method can be used for unbiased unraveling of the complete transcriptome of a wide range of biofluids and their extracellular vesicles.


Blood ◽  
2016 ◽  
Vol 128 (22) ◽  
pp. 1760-1760 ◽  
Author(s):  
Daniela S Gerhard ◽  
Bruno Grande ◽  
Nicholas Griner ◽  
Corey Casper ◽  
Sarah E. Gerdts ◽  
...  

Abstract Introduction: Burkitt Lymphoma (BL) is an aggressive B-cell lymphoma with a translocation involving MYC and immunoglobulin(Ig) loci. It is most common in children, but also affects adults, and occurs in sporadic, endemic and HIV-associated forms. The Epstein-Barr virus (EBV)-associated endemic subtype is the most common pediatric cancer in equatorial Africa, but also occurs in other parts of the world, for example in the rain forest of Brazil. Intensive chemotherapy is effective, but the associated toxicity requires supportive care that is not readily available in resource-poor regions. Previously published molecular characterization of small numbers of tumors indicated that the mutation profiles of endemic and sporadic cases are similar, but not identical. One goal of the BLGSP is to conduct comprehensive molecular characterization of BL by sequencing DNA and RNA from a large BL cohort - including endemic, sporadic, pediatric and adult cases - in order to define the genetic and phenotypic features that drive these cancers. These data will be analyzed with an intent toward developing new therapeutic strategies that can be deployed worldwide. Methods: The goal is to collect 160 BL cases, of which 50% will be endemic, 38% sporadic (pediatric and adult) and 12% from HIV+ patients. For the discovery phase, each tumor requires case-matching normal DNA as well as treatment, outcome and other clinical information. The optimal source of tumor DNA and RNA is from frozen tissue with at least 50% tumor nuclei, but FFPE immobilization is also accepted. Accrual locations include Africa, Brazil, Europe and the US. The BLGSP has developed extensive standard operating procedures for tissue collection, pathology review and tissue processing to reduce the variation associated with these parameters in the interpretation of the results (see https://ocg.cancer.gov/programs/cgci/projects/burkitt-lymphoma). The project also established procedures that allow sharing of all clinical and sample information through the National Cancer Institute Genomic Data Commons (https://gdc.cancer.gov). Molecular characterization includes whole genome sequencing of tumor and normal DNA (80X and 40X coverage, respectively), RNA-sequencing (RNA-seq) and micro-RNA sequencing. These data will enable the BLGSP to identify chromosomal rearrangements, chromosomal copy number alternations, somatic mutations (single nucleotide, insertions, deletions), viral insertions, expression signatures, viral expressions and miRNA regulation of transcripts. Results: To date we have accrued 80 cases of BL of which 75% passed diagnostic pathology review. There was an additional 25% attrition at the tissue processing stage, either due to low quality nucleic acids or low percent tumor nuclei. We have completed sequencing for 45 cases, all but one of which have a MYC translocation involving one of the 3 Ig loci; one case has a MYC rearrangement by FISH analysis that is being characterized further. We have identified recurrent mutations in ID3, DDX3X, ARID1A, FOXO1, TP53, SMARCA4 and other genes. Most mutations are supported by the RNA-seq data, which is also useful in defining the pattern of EBV genome transcription. Preliminary unsupervised hierarchical clustering and principal component analysis of gene expression data defined sample clusters that do not correspond to mutation status or EBV infection, warranting further investigation. Some genes accumulated somatic mutations in a BL subtype-specific fashion. Discussion: BLGSP is an ongoing international collaborative project that will provide a comprehensive molecular portrait of BL subtypes when completed, with the potential to suggest new molecular targets for therapy that can eventually lead to effective treatments that are less toxic than the current regimens. Disclosures Casper: Janssen: Consultancy, Research Funding; Roche: Consultancy, Other: Travel, Accommodation, Expenses; TempTime: Consultancy, Other: Travel, Accommodation, Expenses; Up to Date: Patents & Royalties; GSK: Other: Travel, Accommodation, Expenses. Abramson:Kite Pharma: Consultancy; Abbvie: Consultancy; Seattle Genetics: Consultancy; Gilead: Consultancy. Noy:Pharmacyclics, LLC, an AbbVie Company: Other: travel, accommodations, expenses, Research Funding.


2016 ◽  
Author(s):  
Daniel R. Garalde ◽  
Elizabeth A. Snell ◽  
Daniel Jachimowicz ◽  
Andrew J. Heron ◽  
Mark Bruce ◽  
...  

AbstractRibonucleic acid sequencing can allow us to monitor the RNAs present in a sample. This enables us to detect the presence and nucleotide sequence of viruses, or to build a picture of how active transcriptional processes are changing – information that is useful for understanding the status and function of a sample. Oxford Nanopore Technologies’ sequencing technology is capable of electronically analysing a sample’s DNA directly, and in real-time. In this manuscript we demonstrate the ability of an array of nanopores to sequence RNA directly, and we apply it to a range of biological situations. Nanopore technology is the only available sequencing technology that can sequence RNA directly, rather than depending on reverse transcription and PCR. There are several potential advantages of this approach over other RNA-seq strategies, including the absence of amplification and reverse transcription biases, the ability to detect nucleotide analogues and the ability to generate full-length, strand-specific RNA sequences. Direct RNA sequencing is a completely new way of analysing the sequence of RNA samples and it will improve the ease and speed of RNA analysis, while yielding richer biological information.


2018 ◽  
Author(s):  
Verboom Karen ◽  
Everaert Celine ◽  
Bolduc Nathalie ◽  
Livak J. Kenneth ◽  
Yigit Nurten ◽  
...  

AbstractSingle cell RNA sequencing methods have been increasingly used to understand cellular heterogeneity. Nevertheless, most of these methods suffer from one or more limitations, such as focusing only on polyadenylated RNA, sequencing of only the 3’ end of the transcript, an exuberant fraction of reads mapping to ribosomal RNA, and the unstranded nature of the sequencing data. Here, we developed a novel single cell strand-specific total RNA library preparation method addressing all the aforementioned shortcomings. Our method was validated on a microfluidics system using three different cancer cell lines undergoing a chemical or genetic perturbation. We demonstrate that our total RNA-seq method detects an equal or higher number of genes compared to classic polyA[+] RNA-seq, including novel and non-polyadenylated genes. The obtained RNA expression patterns also recapitulate the expected biological signal. Inherent to total RNA-seq, our method is also able to detect circular RNAs. Taken together, SMARTer single cell total RNA sequencing is very well suited for any single cell sequencing experiment in which transcript level information is needed beyond polyadenylated genes.


2020 ◽  
Vol 38 (15_suppl) ◽  
pp. e14006-e14006
Author(s):  
Kevin McDonnell ◽  
Amit Kulkarni ◽  
Melissa Woodhouse ◽  
Sidney A Smith ◽  
Christine Hong ◽  
...  

e14006 Background: Next generation sequencing (NGS) allows for reliable, comprehensive and cost-effective identification of clinically actionable genetic and genomic alterations. The increasing adoption of NGS in clinical oncology has increased our ability to identify germline alterations predisposing to cancer development as well as somatic changes enabling prescription of individualized cancer treatment and enhanced clinical trial participation. Here we summarize implementation of an NGS-based precision medicine initiative involving oncology patients from a single institution cancer center. Methods: IRB-approved NGS matched whole exome (WES) germline and solid tumor somatic tumor sequencing together with somatic tumor RNA sequencing (RNA-seq) were performed using germline DNA extracted from peripheral blood lymphocytes and nucleic acids for tumor DNA and RNA sequencing obtained from formalin-fixed, paraffin-embedded tumor specimens. Results of sequencing and analyses were presented to a multi-disciplinary tumor board to establish recommendations for management of germline pathogenic variation, therapeutic drug matching, clinical trials eligibility and molecularly informed patient prognosis. Results: A total of 1,005 patients completed sequencing. Germline and somatic WES exceeded 100X and 250X mean target coverage, respectively; somatic RNA-seq exceeded 200 million mean reads. Patients ranged in age from 17 to 90 years. The study cohort comprised comparable numbers of female (51%) and male (49%) patients. Ethnicities and races were broadly represented with 22% of participants identifying as Hispanic, 14% as Asian, 4% as Black, 55% as Non-Hispanic White and 5% as other. The most common solid tumor histological classification was colorectal (18%), followed by breast (16%), prostate (7%), head and neck (7%), sarcoma (7%), ovarian (5%), melanoma (4%) and lung (3%). Bioinformatic analyses and precision medicine tumor board review established that 12% of patients harbored a germline pathogenic variant and 43% carried clinically actionable genetic/genomic alterations; a majority of patients met molecular requirements for participation in a clinical trial. Conclusions: This study confirms the feasibility and utility of clinical NGS and precision medicine tumor board review in clinical oncology to identify germline genetic pathology, deliver personalized cancer therapeutics, increase clinical trial enrollment and clarify diagnosis and prognosis.


Plants ◽  
2020 ◽  
Vol 9 (8) ◽  
pp. 1052
Author(s):  
Yeonhwa Jo ◽  
Chang-Gi Back ◽  
Hoseong Choi ◽  
Won Kyong Cho

The dried peach fruits clinging to peach trees or lying on the ground nearby are known as mummified peach fruits. Here, we examined the microbiome communities of three different mummified peach fruits from the nectarine cultivar “Hahong” by DNA- and RNA-sequencing. We found the dominance of Monilinia fructigena followed by Sclerotinia borealis, S. sclerotiorum, and Botrytis cinerea in the mummified peach fruits. Moreover, we found a high number of Proteobacteria, including Frateuria aurantia, Neoasaia chiangmaiensis, Robbsia andropogonis, and Ewingella Americana. Furthermore, we identified several viruses and viroids. Bacteriophages were identified by DNA- and RNA-sequencing, while viruses and viroids with RNA genomes were identified by only RNA-sequencing. Moreover, we identified a novel mycovirus referred to as Monilinia umbra-like virus 1 (MULV1) from M. fructigena. Our results revealed the co-inhabitance of fungi and bacteria in the mummified peach fruits, although dominant microorganisms were present. RNA-sequencing revealed that several fungal and bacterial genes were actively transcribed. Comparative analyses suggested that RNA-sequencing provides more detailed information on microbial communities; however, combining DNA- and RNA-sequencing results increased the diversity of microorganisms, suggesting the importance of databases and analysis tools for microbiome studies. Taken together, our study provides a comprehensive overview of microbial communities in mummified peach fruits by DNA shotgun sequencing and RNA-sequencing.


2019 ◽  
Vol 9 (1) ◽  
Author(s):  
Celine Everaert ◽  
Hetty Helsmoortel ◽  
Anneleen Decock ◽  
Eva Hulstaert ◽  
Ruben Van Paemel ◽  
...  

AbstractRNA profiling has emerged as a powerful tool to investigate the biomarker potential of human biofluids. However, despite enormous interest in extracellular nucleic acids, RNA sequencing methods to quantify the total RNA content outside cells are rare. Here, we evaluate the performance of the SMARTer Stranded Total RNA-Seq method in human platelet-rich plasma, platelet-free plasma, urine, conditioned medium, and extracellular vesicles (EVs) from these biofluids. We found the method to be accurate, precise, compatible with low-input volumes and able to quantify a few thousand genes. We picked up distinct classes of RNA molecules, including mRNA, lncRNA, circRNA, miscRNA and pseudogenes. Notably, the read distribution and gene content drastically differ among biofluids. In conclusion, we are the first to show that the SMARTer method can be used for unbiased unraveling of the complete transcriptome of a wide range of biofluids and their extracellular vesicles.


Author(s):  
Peter H. Culviner ◽  
Chantal K. Guegler ◽  
Michael T. Laub

AbstractThe profiling of gene expression by RNA-sequencing (RNA-seq) has enabled powerful studies of global transcriptional patterns in all organisms, including bacteria. Because the vast majority of RNA in bacteria is ribosomal RNA (rRNA), it is standard practice to deplete the rRNA from a total RNA sample such that the reads in an RNA-seq experiment derive predominantly from mRNA. One of the most commonly used commercial kits for rRNA depletion, the Ribo-Zero kit from Illumina, was recently discontinued. Here, we report the development a simple, cost-effective, and robust method for depleting rRNA that can be easily implemented by any lab or facility. We first developed an algorithm for designing biotinylated oligonucleotides that will hybridize tightly and specifically to the 23S, 16S, and 5S rRNAs from any species of interest. Precipitation of these oligonucleotides bound to rRNA by magnetic streptavidin beads then depletes rRNA from a complex, total RNA sample such that ~75-80% of reads in a typical RNA-seq experiment derive from mRNA. Importantly, we demonstrate a high correlation of RNA abundance or fold-change measurements in RNA-seq experiments between our method and the previously available Ribo-Zero kit. Complete details on the methodology are provided, including open-source software for designing oligonucleotides optimized for any bacterial species or metagenomic sample of interest.ImportanceThe ability to examine global patterns of gene expression in microbes through RNA-sequencing has fundamentally transformed microbiology. However, RNA-seq depends critically on the removal of ribosomal RNA from total RNA samples. Otherwise, rRNA would comprise upwards of 90% of the reads in a typical RNA-seq experiment, limiting the reads coming from messenger RNA or requiring high total read depth. A commonly used, kit for rRNA subtraction from Illumina was recently discontinued. Here, we report the development of a ‘do-it-yourself’ kit for rapid, cost-effective, and robust depletion of rRNA from total RNA. We present an algorithm for designing biotinylated oligonucleotides that will hybridize to the rRNAs from a target set of species. We then demonstrate that the designed oligos enable sufficient rRNA depletion to produce RNA-seq data with 75-80% of reads comming from mRNA. The methodology presented should enable RNA-seq studies on any species or metagenomic sample of interest.


2019 ◽  
Vol 47 (16) ◽  
pp. e93-e93 ◽  
Author(s):  
Karen Verboom ◽  
Celine Everaert ◽  
Nathalie Bolduc ◽  
Kenneth J Livak ◽  
Nurten Yigit ◽  
...  

Abstract Single cell RNA sequencing methods have been increasingly used to understand cellular heterogeneity. Nevertheless, most of these methods suffer from one or more limitations, such as focusing only on polyadenylated RNA, sequencing of only the 3′ end of the transcript, an exuberant fraction of reads mapping to ribosomal RNA, and the unstranded nature of the sequencing data. Here, we developed a novel single cell strand-specific total RNA library preparation method addressing all the aforementioned shortcomings. Our method was validated on a microfluidics system using three different cancer cell lines undergoing a chemical or genetic perturbation and on two other cancer cell lines sorted in microplates. We demonstrate that our total RNA-seq method detects an equal or higher number of genes compared to classic polyA[+] RNA-seq, including novel and non-polyadenylated genes. The obtained RNA expression patterns also recapitulate the expected biological signal. Inherent to total RNA-seq, our method is also able to detect circular RNAs. Taken together, SMARTer single cell total RNA sequencing is very well suited for any single cell sequencing experiment in which transcript level information is needed beyond polyadenylated genes.


2015 ◽  
Author(s):  
Jerod Parsons ◽  
Sarah Munro ◽  
P. Scott Pine ◽  
Jennifer McDaniel ◽  
Michele Mehaffey ◽  
...  

Genome-scale ?-omics? measurements are challenging to benchmark due to the enormous variety of unique biological molecules involved. Mixtures of previously-characterized samples can be used to benchmark repeatability and reproducibility using component proportions as truth for the measurement. We describe and evaluate experiments characterizing the performance of RNA-sequencing (RNA-Seq) measurements, and discuss cases where mixtures can serve as effective process controls. We apply a linear model to total RNA mixture samples in RNA-seq experiments. This model provides a context for performance benchmarking. The parameters of the model fit to experimental results can be evaluated to assess bias and variability of the measurement of a mixture. A linear model describes the behavior of mixture expression measures and provides a context for performance benchmarking. Residuals from fitting the model to experimental data can be used as a metric for evaluating the effect that an individual step in an experimental process has on the linear response function and precision of the underlying measurement while identifying signals affected by interference from other sources. Effective benchmarking requires well-defined mixtures, which for RNA-Seq requires knowledge of the messenger RNA (mRNA) content of the individual total RNA components. We demonstrate and evaluate an experimental method suitable for use in genome-scale process control and lay out a method utilizing spike-in controls to determine mRNA content of total RNA in samples. Genome-scale process controls can be derived from mixtures. These controls relate prior knowledge of individual components to a complex mixture, allowing assessment of measurement performance. The mRNA fraction accounts for differential enrichment of mRNA from varying total RNA samples. Spike-in controls can be utilized to measure this relationship between mRNA content and input total RNA. Our mixture analysis method also enables estimation of the proportions of an unknown mixture, even when component-specific markers are not previously known, whenever pure components are measured alongside the mixture.


Sign in / Sign up

Export Citation Format

Share Document