Hello to everyone.
In my group, we are using RNA-Seq to study the expression of bacteria under different conditions, isolated or as dual-transcriptome when in association with plants.
When checking the table of total counts for each gene, I can observe an important number of reads mapping against the rrna genes, mostly 16S and 23S. I also found a similar amount of read counts mapping against rrnas when analyzing a dual experiment where no bacteria was present, just plant material.
Illumina reads from HiSeq200 were treated by TopHat and Cufflinks. Library preparations followed a RiboZero Gold treatment + TruSeq Stranded Total RNA kit. Additionally, rrna depletion was verified by running Bioanalyzer.
Have any one an idea where this mapping might come?
Thanks in advance for your help.
•
written
by
buthercup_ch •
2 answers
Bacterial contamination is extremely common in many lab reagents.
https://bmcbiol.biomedcentral.com/articles/10.1186/s12915-014-0087-z
Hello
Thanks for your comments.
I was also thinking about if these sequences matching the bacterial rrna genes could actually be product of high sequence duplication (checking the FASTQC file for the samples showed red warning in this point). They might be an artifact product of PCR enrichment, over-sequencing, or contamination of library preparation (I just read about the highly frequent as common library prep and sequencing reagents gets contaminated by bacterial DNA).
What do you think?
Regards
Cristina
Log in to answer this question.
-
DESeq2 sizeFactor for Bacterial RNASeq
written by iXuan •
-
FastQC duplicates questions
written by Beth •
-
RNAseq analysis of phage infection - should I split the counts mapping to host/phage prior to norma…
written by Sean •
-
High percentage of intronic/intergenic reads in RNA-seq
written by lluc.cabus •
-
Low mapping rate due to rRNA
written by garcesj •
-
Identifying and classifying the most abundant sequences in small RNAseq data
written by Mauro •
-
Can't count features with featurecount
written by avelarbio46 •
-
Too low FPKM values from StringTie
written by SJ Basu •
-
STAR output - high multi-mapping, low alignment (human)
written by anikb •
-
Low gene count assignment for ribo-depleted RNA-Seq data?
written by gabriel.rosser
What fraction of the reads are aligning to rRNA? Experimental rRNA depletion may never be perfect so as long as the % is small it should be ok.