Zea mays RNA-seq estimated transcript abundances are strongly affected by read mapping bias

Abstract

Background Genetic variation for gene expression is a source of phenotypic variation for natural and agricultural species. The common approach to map and to quantify gene expression from genetically distinct individuals is to assign their RNA-seq reads to a single reference genome. However, RNA-seq reads from alleles dissimilar to this reference genome may fail to map correctly, causing transcript levels to be underestimated. Presently, the extent of this mapping problem is not clear, particularly in highly diverse species. We investigated if mapping bias occurred and if chromosomal features associated with mapping bias. Zea mays presents a model species to assess these questions, given it has genotypically distinct and well-studied genetic lines. Results In Zea mays, the inbred B73 genome is the standard reference genome and template for RNA-seq read assignments. In the absence of mapping bias, B73 and a second inbred line, Mo17, would each have an approximately equal number of regulatory alleles that increase gene expression. Remarkably, Mo17 had 2–4 times fewer such positively acting alleles than did B73 when RNA-seq reads were aligned to the B73 reference genome. Reciprocally, over one-half of the B73 alleles that increased gene expression were not detected when reads were aligned to the Mo17 genome template. Genes at dissimilar chromosomal ends were strongly affected by mapping bias, and genes at more similar pericentromeric regions were less affected. Biased transcript estimates were higher in untranslated regions and lower in splice junctions. Bias occurred across software and alignment parameters. Conclusions Mapping bias very strongly affects gene transcript abundance estimates in maize, and bias varies across chromosomal features. Individual genome or transcriptome templates are likely necessary for accurate transcript estimation across genetically variable individuals in maize and other species.

Divergence of cochlear transcriptomics between reference‑based and reference‑free transcriptome analyses among Rhinolophus ferrumequinum populations

“…These incorrect gene identifications would affect subsequent gene annotation, resulting in errors in functional analysis results. Sequence identification deviations would have a great impact on gene transcript abundance, and eventually lead to the underestimation of transcription levels of some important genes [ 44 , 45 ]. Our results showed that although gene expression levels of shared DEGs obtained by the two methods were strongly correlated ( Fig 2C ), significant expression level differences of the same genes between the two methods were shown by paired Mann–Whitney U test ( Fig 2D , S8 Table ).…”

Variant calling in polyploids for population and quantitative genetics

“…. For example, in a Zea mays RNA-seq study, as much as one-half of alleles with increased gene expression were not detected when reads from the inbred line, B73, were mapped to the reference of a second inbred line, Mo17, because Z. mays has high nucleotide diversity and structural variation (Zhan et al, 2021) .…”

Molecular Interactions BetweenLeptosphaeria maculansandBrassicaSpecies

“…There is also the potential for read-mapping bias inherent in RNA-Seq analyses using a reference genome that differs from the actual fungal isolate or host plant line used in the experiment. This can result in the underestimation of transcript abundance for alleles that differ from the reference (150).…”