Validation of Two Ribosomal RNA Removal Methods for Microbial Metatranscriptomics

Total Page:16

File Type:pdf, Size:1020Kb

Validation of Two Ribosomal RNA Removal Methods for Microbial Metatranscriptomics ANaLYSIS Validation of two ribosomal RNA removal methods for microbial metatranscriptomics Shaomei He1,2,5, Omri Wurtzel3,5, Kanwar Singh1, Jeff L Froula1, Suzan Yilmaz1, Susannah G Tringe1, Zhong Wang1, Feng Chen1, Erika A Lindquist1, Rotem Sorek3 & Philip Hugenholtz1,2,4 The predominance of rRNAs in the transcriptome is a major Among these methods, subtractive hybridization and exo- technical challenge in sequence-based analysis of cDNAs from nuclease digestion have become the most popular owing to the microbial isolates and communities. Several approaches have availability of commercial kits from Ambion (MICROBExpress been applied to deplete rRNAs from (meta)transcriptomes, but Bacterial mRNA Enrichment kit) and Epicentre (mRNA-ONLY no systematic investigation of potential biases introduced by Prokaryotic mRNA Isolation kit). The former kit uses a subtrac- any of these approaches has been reported. Here we validated tive hybridization with capture oligonucleotides specific to 16S the effectiveness and fidelity of the two most commonly and 23S rRNAs. It has been applied to both bacterial isolates used approaches, subtractive hybridization and exonuclease and environmental samples, in one or two rounds6,13–18. The digestion, as well as combinations of these treatments, on Epicentre kit uses exonuclease to preferentially degrade processed two synthetic five-microorganism metatranscriptomes using RNAs with 5′ monophosphate (the majority of which are believed massively parallel sequencing. We found that the effectiveness to be rRNAs)19,20. In some instances, these methods have been of rRNA removal was a function of community composition and used in combination to improve rRNA removal21–23. There is no RNA integrity for these treatments. Subtractive hybridization consensus, however, on the best approach, and existing data are alone introduced the least bias in relative transcript abundance, anecdotal. Here we validated the performance of these kits with whereas exonuclease and in particular combined treatments two synthetic five-member microbial communities using Illumina greatly compromised mRNA abundance fidelity.I llumina sequencing and found that rRNA removal efficiencies were com- sequencing itself also can compromise quantitative data munity- and RNA integrity–dependent and that only subtractive analysis by introducing a G+C bias between runs. hybridization adequately preserved relative transcript abundance for quantitative analyses. Rapid technological advances in ultra-high-throughput sequenc- © All rights reserved. 2010 Inc. Nature America, ing are making de novo sequencing of transcriptomes (RNA-seq) RESULTS a viable alternative to microarray analysis of microbial isolates Experimental design and communities1. A major technical challenge for de novo tran- We constructed two five-member synthetic microbial communities scriptome sequencing is the low relative abundance of mRNAs by pooling equimolar amounts of total RNAs extracted independ- in total cellular RNA (1–5%; ref. 2), the bulk of which is rRNAs ently from microbial isolates with sequenced genomes that span and tRNAs3. Unlike eukaryotic mRNAs, which can be selectively a wide phylogenetic, (G+C) content and genome-size range. Two synthesized into cDNA by virtue of their poly(A) tails4, bacte- of the five species were common to both communities (Table 1). rial and archaeal cDNAs are predominantly rRNA sequences5,6. We tested subtractive hybridization (Hyb), exonuclease digestion Therefore, prokaryotic rRNAs are often removed before sequenc- (Exo) and combined treatments in two experiments using the two ing to improve mRNA detection sensitivity. Different methods synthetic communities. In experiment 1, we tested Hyb, Exo and have been used to eliminate prokaryotic rRNA, including sub- Hyb followed by Exo (Hyb + Exo) using synthetic community 1. tractive hybridization with rRNA-specific probes7,8, digestion Then we conducted experiment 2 using synthetic community 2 to with exonuclease that preferentially acts on rRNA, poly(A) tail confirm the results of experiment 1 and expanded the depletion addition to discriminate against rRNA9,10, reverse transcription methods tested to also include two rounds of Hyb (2Hyb) as well with rRNA-specific primers followed by RNase H digestion as Exo followed by Hyb (Exo + Hyb) (Table 2). As a control, we to degrade rRNA:DNA hybrids11, and gel electrophoresis size used total RNA without rRNA removal. We initially assessed rRNA separation and extraction of non-rRNA bands12. removal using RNA electropherograms (Agilent 2100 Bioanalyzer) 1Department of Energy Joint Genome Institute, Walnut Creek, California, USA. 2Energy Biosciences Institute, University of California-Berkeley, Berkeley, California, USA. 3Department of Molecular Genetics, Weizmann Institute of Science, Rehovot, Israel. 4Australian Centre for Ecogenomics, School of Chemistry and Molecular Biosciences, The University of Queensland, Brisbane, Australia. 5These authors contributed equally to this work. Correspondence should be addressed to P.H. ([email protected]). RECEIVED 12 APrIL; aCCEPTED 30 AUGUST; PUBLISHED ONLINE 19 SEPTEMBEr 2010; DOI:10.1038/NMETH.1507 NaTURE METHODS | ADVANCE ONLINE PUBLICATION | ANaLYSIS Table 1 | Details of microbial isolates used in the two synthetic communities Genome G+C content Match Hyb Organism size (Mbp) (%) Phylum target sitesa 23S/16Sb RINb Community Desulfovibrio vulgaris strain Hildenborough 3.7 63 Proteobacteria Yes 0.9 7.4 1 Streptomyces sp. strain LCC 8–10c 71 Actinobacteria Yes 1.4 5.5 1 Lactococcus lactis subspecies lactis IL1403 2.53 35 Firmicutes Yes 1.6 9.8 1 Spirochaeta aurantia subspecies aurantia M1 4.3 65 Spirochaeta Yes 2.1 9.8 1 and 2 Lactobacillus brevis ATCCd 367 2.3 46 Firmicutes Yes 1.4 9.9 1 and 2 Kangiella koreensis DSMd 16069 2.9 43 Proteobacteria Yes 1.4 10 2 Catenulispora acidiphila DSMd 44928 10.5 70 Actinobacteria Yes 1.3 8.6 2 Halorhabdus utahensis DSMd 12940 3.1 63 Euryarchaeota No 1.9 10 2 aInformation about 16S and 23S rRNAs matching to the capture oligos in the MICROBExpress mRNA Enrichment Kit was obtained from Ambion. bThe ratio of 23S/16S rRNAs and RNA integrity number (RIN) were determined using the Agilent 2100 Bioanalyzer. RIN values range from 1 (most degraded) to 10 (most intact). cDraft assembly; genome size was estimated to be eight to ten mega–base pairs (Mbp), which is in the range of the size of genomes of this organism’s close relatives, such as S. coelicolor A3(2) (9,054,847 bp), S. griseus subspecies griseus NITE Biological Resource Center (NBRC) 13350 (8,545,929 bp) and S. avermitilis MA-468 (9,119,895 bp). dATCC, American Type Culture Collection; DSM, Deutsche Sammlung von Mikrorganismen. (Supplementary Figs. 1 and 2 and Supplementary Note 1). We Linear regression slopes for each organism positively correlated then sequenced technical replicates (after RNA pooling) distrib- with their genomic G+C content (Fig. 1c). This indicates that uted within and between four runs (flowcells; runs 1 and 2 in run 3 systematically underrepresented (G+C)-rich sequences experiment 1, and runs 3 and 4 in experiment 2) on an Illumina relative to run 4, which we believe was due to variability in run Genome Analyzer II sequencer to evaluate intra- and interrun quality (Table 2). Two simple normalization strategies based on variation (Table 2). We generated ~10–17 million 76-base-pair the total mRNA counts from the source organism or gene G+C single-end reads for each sample, and for all but one sample we content (Supplementary Note 3) improved the overall correlation mapped 99% of reads to a reference, indicating good read quality between technical replicates (Fig. 1d), confirming that the run- and negligible contamination (Supplementary Tables 1 and 2 and to-run variation was indeed largely associated with G+C content. Supplementary Note 2). Therefore, to accurately assess rRNA depletion efficiency and mRNA fidelity, we restricted comparisons of treatment-control Technical reproducibility pairs to the same run in experiment 2 (Table 2). We took the We first evaluated technical reproducibility by comparing the rela- interrun average after we performed the intrarun comparison. tive transcript abundance of technical replicates calculated as reads per gene normalized by total mapped mRNA reads. The four intra- Efficiency of rRNA depletion run technical replicates in three independent runs (Hyb in run 1, As expected, the relative rRNA content of the controls was and controls in runs 1, 3 and 4; Table 2) were highly reproducible 95–97% of total RNA, typical of bacteria and archaea, with each (Pearson’s product-moment correlation coefficient, r = 0.997 ± community member being approximately equally represented 0.001; Fig. 1a), suggesting that both sample preparation–associated in the controls (Fig. 2a,b). After the various treatments, the ­variation and intrarun sequencing variation were minor. In addi- observed rRNA content decreased by as little as 3.6% (Exo) and tion, the sole interrun technical replicate between runs 1 and 2 up to 19.9% (Hyb + Exo). This decrease in rRNA percentage in © All rights reserved. 2010 Inc. Nature America, (Exo 1 i and Exo 1 ii; Table 2) was also highly reproducible total RNA reflected the redistribution of rRNA relative to non- (Fig. 1b), suggesting negligible technical variation between these rRNA owing to the treatment. For experiment 1, the rank order two runs. By contrast, the interrun technical replicates between of rRNA removal efficiencies was Exo < Hyb < Hyb + Exo;
Recommended publications
  • Comparison of RIN and Rine Algorithms for the Agilent 2100 Bioanalyzer and the Agilent 2200 Tapestation Systems Technical Overview
    Comparison of RIN and RINe algorithms for the Agilent 2100 Bioanalyzer and the Agilent 2200 TapeStation systems Technical Overview Introduction A typical application for the 96-well plate compatible Agilent 2200 TapeStation system is automated, effi cient, and reliable RNA analysis, including RNA charac- terization and quality assessment with Agilent R6K ScreenTape or Agilent High Sensitivity R6K ScreenTape. When separated by electrophoresis, eukaryotic total RNA samples typically display two major peaks representing the 18S and 28S ribosomal RNA (rRNA). Depending on the purifi cation method, additional smaller peaks can also been seen, which correspond to small rRNAs as well as tRNAs. As total RNA degrades, the 18S and 28S rRNA peaks slowly disappear and the degradation products emerge in the region between the 18S and small RNAs, termed as the fast zone. The RNA integrity number (RIN) functionality of the Agilent 2100 Expert software for the Agilent 2100 Bioanalyzer system uses a proprietary neural network trained algorithm to make an assessment of the entire electrophoretic trace for a given total RNA sample, including ribosomal peak ratios, separation, and the presence or absence of degradation products. RIN values from 1 to 10 can be obtained, where 10 indicates the highest possible RNA quality. Due to fundamental differences in the nature of the raw data, a dedicated soft- ware algorithm has been developed for the 2200 TapeStation system and R6K ScreenTape to determine the RNA integrity number equivalent (RINe). The RINe algorithm is based on a Experimental RNA preparation mathematical model that calculates an Total RNA from HEK 293, Hep G2, and objective quantitative measurement of Material Jurkat cells was purifi ed using the RNA degradation1.
    [Show full text]
  • Low RIN Value for RNA-Seq Library Construction from Long-Term Stored Seeds: a Case Study of Barley Seeds
    G C A T T A C G G C A T genes Technical Note Low RIN Value for RNA-Seq Library Construction from Long-Term Stored Seeds: A Case Study of Barley Seeds Marta Puchta , Maja Boczkowska * and Jolanta Groszyk National Centre for Plant Genetic Resources, Plant Breeding and Acclimatization National Research Institute, Radzików, 05-870 Błonie, Poland; [email protected] (M.P.); [email protected] (J.G.) * Correspondence: [email protected] Received: 20 August 2020; Accepted: 8 October 2020; Published: 13 October 2020 Abstract: Seed aging is a complex biological process and its fundamentals and mechanisms have not yet been fully recognized. This is a key issue faced by research teams involved in the collection and storage of plant genetic resources in gene banks every day. Transcriptomic changes associated with seed aging in the dry state have barely been studied. The aim of the study was to develop an efficient protocol for construction of RNA-Seq libraries from long-term stored seeds with very low viability and low RNA integrity number (RIN). Here, barley seeds that have almost completely lost their viability as a result of long-term storage were used. As a control, fully viable seeds obtained in the course of field regeneration were used. The effectiveness of protocols dedicated to RNA samples with high and low RIN values was compared. The experiment concluded that library construction from low viable or long-term stored seeds with degraded RNA (RIN < 3) should be carried out with extraordinary attention due to the possibility of uneven degradation of different RNA fractions.
    [Show full text]
  • PD-WP-017 Issue2 (Interpreting Bioanaylzer).Indd
    Interpreting bioanalyzer results for RNA collected using Oragene®•RNA The Agilent Bioanalyzer is an instrument for assessing however, for several reasons, it is more difficult to the “quality” of RNA used for microarray analysis. obtain a high RIN value. As with blood, saliva has a The sharpness of the ribosomal RNA peak and the lower RNA to DNA ratio than that of cultured cells or low baseline between peaks on the electropherogram tissue therefore similar care must be taken to remove generated by the instrument are used as measures of the contaminating DNA. Additionally, saliva is not a integrity of the ribosomal RNA and, by implication, the sterile source of nucleic acid; some of the RNA integrity of messenger RNA. The bioanalyzer uses an contained in saliva originates from oral bacteria. algorithm to generate an “RNA integrity number” (RIN) The ratio between the 16S and 23S bacterial ribosomal that helps researchers decide if the RNA is of sufficiently RNAs is less than the human 18S:28S ratio. The high integrity to proceed with microarray analysis. algorithm used by the bioanalyzer was not designed for use with samples containing both human and Most high quality RNA samples share the same bacterial RNA so the instrument generates a lower common features: RIN that is not indicative of the integrity of the RNA. • distinct 18S and 28S ribosomal peaks Instead of relying on the RIN, a visual assessment of (or 16S and 23S for prokaryotic samples) the electropherogram can provide valuable insight • the absence of smaller peaks between the into the quality of the RNA.
    [Show full text]
  • Impact of RNA Degradation on Gene Expression Profiling
    Opitz et al. BMC Medical Genomics 2010, 3:36 http://www.biomedcentral.com/1755-8794/3/36 RESEARCH ARTICLE Open Access Impact of RNA degradation on gene expression profiling Lennart Opitz1, Gabriela Salinas-Riester1, Marian Grade2, Klaus Jung3, Peter Jo2, Georg Emons2, B Michael Ghadimi2, Tim Beißbarth3*†, Jochen Gaedcke2† Abstract Background: Gene expression profiling is a highly sensitive technique which is used for profiling tumor samples for medical prognosis. RNA quality and degradation influence the analysis results of gene expression profiles. The impact of this influence on the profiles and its medical impact is not fully understood. As patient samples are very valuable for clinical studies, it is necessary to establish criteria for the RNA quality to be able to use these samples in later analysis. Methods: To investigate the effects of RNA integrity on gene expression profiling, whole genome expression arrays were used. We used tumor biopsies from patients diagnosed with locally advanced rectal cancer. To simulate degradation, the isolated total RNA of all patients was subjected to heat-induced degradation in a time- dependent manner. Expression profiling was then performed and data were analyzed bioinformatically to assess the differences. Results: The differences introduced by RNA degradation were largely outweighed by the biological differences between the patients. Only a relatively small number of probes (275 out of 41,000) show a significant effect due to degradation. The genes that show the strongest effect due to RNA degradation were, especially, those with short mRNAs and probe positions near the 5′ end. Conclusions: Degraded RNA from tumor samples (RIN > 5) can still be used to perform gene expression analysis.
    [Show full text]
  • Assessing Integrity of Insect RNA
    Assessing Integrity of Insect RNA Application Note Nucleic Acid Analysis Authors Abstract Jeffrey A. Fabrick and J. Joe Hull Assessing total RNA integrity is important for the success of downstream RNA Pest Management and Biocontrol applications. The Agilent 2100 Bioanalyzer system with the RNA Integrity Number Research Unit, (RIN) provides a quantitative measure of RNA degradation. Although RINs may not Agricultural Research Service, be ascertained for RNA from all organisms, namely those with unusual or complex United States Department of RNA profiles, this Application Note demonstrates the applicability of total RNA Agriculture, quality assessment of different insect species with the 2100 Bioanalyzer system Maricopa, AZ, USA and its relevance to downstream quantitative RT-PCR experiments. Introduction Experimental RNA extraction Total RNA was extracted from insect In general, the success of downstream Materials RNA applications depends on the tissues using a TissueLyser and TRI quality and integrity of the extracted TissueLyser was obtained from Qiagen Reagent. For RNA degradation, tissues total RNA. Therefore, evaluating levels (Valencia, CA, USA). TRI Reagent, were homogenized in 2-mL round-bottom of degradation in RNA samples is OPTIZYME DNase I and a NanoDrop microfuge tubes with two steel critical. Electrophoretic separation and 1000 spectrophotometer were from beads in PBS buffer (700 µL) or TRI visual inspection of the highly abundant Thermo Fisher Scientific (Waltham, Reagent (positive controls), and held ribosomal RNAs (rRNAs) – 16S and MA, USA). An Agilent 2100 Bioanalyzer at room temperature for 0, 5, 10, 30, or 23S in prokaryotes and 18S and 28S in system (G2940CA or G2943CA) and the 60 minutes.
    [Show full text]
  • Massey Genome Service
    MASSEY GENOME SERVICE RNA Quality Analysis using Agilent 2100 Bioanalyzer Operational since 2011 BIOANALYZER TECHNICAL INFORMATION BULLETIN January 2017 BULLETIN INCLUDES Checking the quality of RNA Assay Results Details of Bioanalyzer Report Summary Details of Bioanalyzer Service QC Report Enquiries regarding sequencing technical information contact: Xiaoxiao Lin – Acting Laboratory Manager Massey Genome Service Private Bag 11222, Palmerston North 4442, New Zealand Lab +64 6 350 5378 Office +64 6 951 9080 http://genome.massey.ac.nz MASSEY GENOME SERVICE (MGS) RNA Quality Analysis using Agilent 2100 Bioanalyzer and the RNA 6000 Nano Assay TECHNICAL BULLETIN Background This document provides you with information on how to check the quality of your RNA 6000 Nano assay data, and includes information on: Checking the quality of the RNA Ladder used for product sizing Checking the quality of the total RNA profiles and mRNA profiles Massey Genome Service (MGS) provides the following two reports to the customer: Bioanalyzer Report Summary: This report is generated by the Agilent Bioanlyzer 2100 Expert software used for running the instrument. The report is a summary of the quality of your RNA samples and includes the RNA ladder, which is run with each LabChip assay and is used for product sizing. Bioanalyzer Service QC Report: This report is generated by the MGS and is a summary of our services analysis of the data quality. Included with these two reports are the Data Files, which can be opened using the Agilent 2100 expert software. This is the same software used by the MGS for running and analyzing the results. This software can be downloaded from the MGS website at http://genome.massey.ac.nz , by going to the ‘Bioanalyzer Service’ link.
    [Show full text]
  • RNA-QC-Protocol-2013.Pdf
    SAMPLE MANAGEMENT STANDARD OPERATING PROCEDURE RNA (tube based) QC SOP (For Collaborators) Version Number: 1 Version Date: 5.03.13 Author(s): Jennifer Chiniquy Reviewed/Revised by: Mansi Chovatia, Chia-lin Wei and Shweta Deshpande Summary Before shipping your RNA sample(s), please be sure to follow the JGI sample preparation and sample submission guidelines located at http://my.jgi.doe.gov/general/gettingstarted.html This protocol describes how to perform quality control of RNA samples to evaluate the quantity (Qubit RNA BR Assay) and quality (Agilent Bioanalyzer, RNA 6000 Nano Kit) of your RNA sample(s). We recommend all RNA samples to be evaluated with this protocol before they are shipped to JGI. Materials & Reagents Materials/Reagents/Equipment Supplier/Vendor Stock Number Disposables Qubit assay tubes Q32856 10 ml Serological pipettes 1367811E 0.2 ml PCR 8-tube flex-free strip (RNase-free) 1402-4700 0.2 ml PCR tubes (individual) Ambion/Life Technologies AM12225 Bioanalyzer Cleaning Chips Agilent 5065-9951 Petri Dishes, 100x15 mm, BD Falcon 25373-100 1.5 ml microcentrifuge tubes (RNase-free) Ambion AM12450 Air canister 23022523 Agilent RNA 6000 Nano Kit and Chips (includes syringe for Priming Station) Agilent 5067-1511 Reagents RNase Zap Ambion AM9780 DNA AWAY Surface Decontaminant Fisher Scientific 21-236-28 70% Isopropanol Fisher Scientific 19-130-3860 Nuclease Free Water (10 × 50 ml) Ambion AM9937 Qubit RNA BR Assay Kit Invitrogen Q10211 Equipment 96-well cold block E&K Scientific EK-76120 RNA QC.V1.DOC 1/11 SAMPLE MANAGEMENT STANDARD OPERATING PROCEDURE Bioanalyzer Chip Priming Station Agilent 5065-4401 Chip Vortexer IKA 3617036 2100 Electrophoresis Bioanalyzer Agilent G2939AA Qubit 2.0 Fluorometer Invitrogen Q32866 Thermocycler Applied Biosystems N8050001 Microcentrifuge Eppendorf 022620100 Serological pipettor Eppendorf 022230204 Speed Vac Thermo Fisher SPD1010 EH&S PPE Requirements: Safety glasses, lab coat, and nitrile gloves should be worn at all times while performing this protocol.
    [Show full text]
  • REVIEW: RNA Integrity and the Effect on the Real-Time Qrt-PCR Performance
    Molecular Aspects of Medicine 27 (2006) 126–139 www.elsevier.com/locate/mam Review RNA integrity and the effect on the real-time qRT-PCR performance Simone Fleige a, Michael W. Pfaffl a,b,* a Physiology Weihenstephan, Center of Life and Food Sciences (ZIEL), Technical University of Munich, 85350 Freising, Germany b TATAA Biocenter Germany, Freising-Weihenstephan, Germany Abstract The assessment of RNA integrity is a critical first step in obtaining meaningful gene expres- sion data. Working with low-quality RNA may strongly compromise the experimental results of downstream applications which are often labour-intensive, time-consuming, and highly expensive. Using intact RNA is a key element for the successful application of modern mole- cular biological methods, like qRT-PCR or micro-array analysis. To verify RNA quality now- adays commercially available automated capillary-electrophoresis systems are available which are on the way to become the standard in RNA quality assessment. Profiles generated yield information on RNA concentration, allow a visual inspection of RNA integrity, and generate approximated ratios between the mass of ribosomal sub-units. In this review, the importance of RNA quality for the qRT-PCR was analyzed by determining the RNA quality of different bovine tissues and cell culture. Independent analysis systems are described and compared (OD measurement, NanoDrop, Bioanalyzer 2100 and Experion). Advantage and disadvantages of RNA quantity and quality assessment are shown in performed applications of various tissues and cell cultures. Further the comparison and correlation between the total RNA integrity on PCR performance as well as on PCR efficiency is described. On the basis of the derived results we can argue that qRT-PCR performance is affected by the RNA integrity and PCR efficiency in general is not affected by the RNA integrity.
    [Show full text]
  • Evaluation for Integrity of Extracted RNA by Reference Material Of
    al Diag ic no ed s ti M c f M Journal of Medical Diagnostic o l e a t h n Takahashi et al., J Med Diagn Meth 2013, 2:3 r o d u s o DOI: 10.4172/2168-9784.1000128 J Methods ISSN: 2168-9784 Research Article Article OpenOpen Access Access Evaluation for Integrity of Extracted RNA by Reference Material of RNA Junko Takahashi1*, Akiko Takatsu2 and Hitoshi Iwahashi3 1Biomedical Research Institute, National Institute of Advanced Industrial Science and Technology Tsukuba, Ibaraki, Japan 2National Metrology Institute of Japan, National Institute of Advanced Industrial Science and Technology, Tsukuba, Ibaraki, Japan 3Faculty of Applied Biological Sciences, Gifu University, Gifu, Japan Abstract Physiological information that identified from RNA expression level in tissues or organs give us benefit in clinical diagnostics. After extracting RNA from blood, tissue, or organs, diagnostics were carried out by using technique for semi-quantitative estimation of specific mRNA, such as real-time PCR or DNA microarrays. Validation of RNA quality need to quantitative estimation of gene expression at the moment of sampling, because using degraded RNA will lead to incorrect results. In this study, we conducted a series of experiment to evaluate RNA integrity using the reference materials of RNA length of 1000 base pairs (refRNA1000), and compared with RIN. There are significant correlation between RIN and the amount of refRNA1000, in varying state of degradation. There are also significant correlations between the amount of refRNA1000 and endogenous gene such as Actb and Rps18 in varying state of degradation in our study. The amount of endogenous gene reflects the extent of degradation more directly than that of RIN.
    [Show full text]
  • RNA-Seq: Impact of RNA Degradation on Transcript Quantification Irene Gallego Romero1, Athma a Pai1,2, Jenny Tung1,3 and Yoav Gilad1*
    Gallego Romero et al. BMC Biology 2014, 12:42 http://www.biomedcentral.com/1741-7007/12/42 RESEARCH ARTICLE Open Access RNA-seq: impact of RNA degradation on transcript quantification Irene Gallego Romero1, Athma A Pai1,2, Jenny Tung1,3 and Yoav Gilad1* Abstract Background: The use of low quality RNA samples in whole-genome gene expression profiling remains controversial. It is unclear if transcript degradation in low quality RNA samples occurs uniformly, in which case the effects of degradation can be corrected via data normalization, or whether different transcripts are degraded at different rates, potentially biasing measurements of expression levels. This concern has rendered the use of low quality RNA samples in whole-genome expression profiling problematic. Yet, low quality samples (for example, samples collected in the course of fieldwork) are at times the sole means of addressing specific questions. Results: We sought to quantify the impact of variation in RNA quality on estimates of gene expression levels basedonRNA-seqdata.Todoso,wecollectedexpression data from tissue samples that were allowed to decay for varying amounts of time prior to RNA extraction. The RNA samples we collected spanned the entire range of RNA Integrity Number (RIN) values (a metric commonly used to assess RNA quality). We observed widespread effects of RNA quality on measurements of gene expression levels, as well as a slight but significant loss of library complexity in more degraded samples. Conclusions: While standard normalizations failed to account for the effects of degradation, we found that by explicitly controlling for the effects of RIN using a linear model framework we can correct for the majority of these effects.
    [Show full text]
  • Qpcr RNA Quality Matters
    kubista_Layout 1 04/12/2012 10:45 Page 1 watchara / Shutterstock watchara qPCR RNA quality matters Mikael Kubista, Jens Björkman, David Svec and Robert Sjöback TATAA Biocenter RNA levels can be measured with very high specificity, sensitivity and accuracy with illumination. In those cases, damage can be techniques such as real-time quantitative PCR (qPCR), microarray analysis and next reduced by degassing the solution removing generation sequencing. This makes messenger (m) RNAs and potentially microRNAs oxygen. RNA is also prone to degradation and other non-coding RNAs popular as biomarkers. But RNA is less stable and more by RNases. dynamic than DNA, and assays are not always specific for RNA, so can we trust There are many types of RNases including 3’ measured expression values? and 5’ exonucleases as well as endonucleases. RNases are incredibly stable enzymes that A biomarker is a biological molecule found in ribozymes, but also makes RNA intrinsically withstand heat as they rapidly refold upon blood, other body fluids or tissues, and is a sign unstable. In aqueous solution, RNA spon - cooling; they are not even fully destroyed of a normal or abnormal process, or of a taneously degrades through self-cleavage by autoclaving. Chaotropic salts such as condition or disease1. The biomarker may be catalysed by metal ions such as Mg2+, high (>9) guanidinium will temporarily inhibit RNases, used to see how well the body responds to a or low (<2) pH, and temperature. EDTA or citrate but irreversible inhibition requires treat- treatment for a disease or condition. Most is therefore typically added to RNA preserving ment with beta mercaptoethanol or DEPC popular and common molecular biomarkers are solutions to chelate Mg2+2.
    [Show full text]
  • The Spatial RNA Integrity Number Assay for in Situ Evaluation of Transcriptome Quality
    ARTICLE https://doi.org/10.1038/s42003-020-01573-1 OPEN The spatial RNA integrity number assay for in situ evaluation of transcriptome quality Linda Kvastad 1, Konstantin Carlberg1,8, Ludvig Larsson1,8, Eva Gracia Villacampa 1,8, Alexander Stuckey1, Linnea Stenbeck 1, Annelie Mollbrink1, Margherita Zamboni 2, Jens Peter Magnusson 2,7, Elisa Basmaci3,4, Alia Shamikh3,4, Gabriela Prochazka3,4, Anna-Lena Schaupp5, Åke Borg 6, Lars Fugger5, Monica Nistér 3,4 & 1234567890():,; ✉ Joakim Lundeberg 1 The RNA integrity number (RIN) is a frequently used quality metric to assess the com- pleteness of rRNA, as a proxy for the corresponding mRNA in a tissue. Current methods operate at bulk resolution and provide a single average estimate for the whole sample. Spatial transcriptomics technologies have emerged and shown their value by placing gene expres- sion into a tissue context, resulting in transcriptional information from all tissue regions. Thus, the ability to estimate RNA quality in situ has become of utmost importance to overcome the limitation with a bulk rRNA measurement. Here we show a new tool, the spatial RNA integrity number (sRIN) assay, to assess the rRNA completeness in a tissue wide manner at cellular resolution. We demonstrate the use of sRIN to identify spatial variation in tissue quality prior to more comprehensive spatial transcriptomics workflows. 1 Science for Life Laboratory, KTH - Royal Institute of Technology (KTH), SE-171 65, Solna, Sweden. 2 Department of Cell and Molecular Biology, Karolinska Institutet, Stockholm, Sweden. 3 Department of Oncology-Pathology, Karolinska Institutet, Stockholm, Sweden. 4 Department of Clinical Pathology and Cytology, Karolinska University Hospital, Stockholm, Sweden.
    [Show full text]