Related Experiment Video
Updated: May 20, 2026
10:10
Three Differential Expression Analysis Methods for RNA Sequencing: limma, EdgeR, DESeq2
Published on: September 18, 2021
Technical and biological variance structure in mRNA-Seq data: life in the real world
Ann L Oberg1, Brian M Bot, Diane E Grill
1Division of Biomedical Statistics and Informatics, Department of Health Sciences Research, Mayo Clinic, 200 1st St SW, Rochester, MN 55905, USA. oberg.ann@mayo.edu
BMC Genomics
|July 10, 2012
Summary
RNA sequencing count data analysis benefits from using Negative Binomial distribution models. These models better capture biological variability and improve the identification of true biological signals in gene expression studies.
Area of Science:
- Genomics
- Bioinformatics
- Statistical Modeling
Background:
- Next-generation sequencing generates mRNA expression data as gene or exon counts.
- Traditionally, count data assumed a Poisson distribution (variance equals mean).
- The Negative Binomial distribution is often used for count data, accommodating over-dispersion (variance > mean).
Purpose of the Study:
- To evaluate the suitability of Poisson and Negative Binomial distributions for modeling mRNA sequencing count data.
- To assess the impact of over-dispersion and experimental effects on statistical modeling.
- To guide the development of analytical strategies for accurate variance modeling in gene expression studies.
Main Methods:
- Analysis of mRNA sequencing data from 25 subjects.
- Comparison of Poisson, over-dispersed Poisson, and Negative Binomial models.
- Assessment of goodness-of-fit and over-fitting for different models.
- Inclusion of experimental effects in modeling.
Main Results:
- Technical variation in mRNA-Seq data generally followed a Poisson distribution.
- Biological variability exhibited over-dispersion relative to the Poisson model.
- A quadratic mean-variance relationship indicated a Negative Binomial distribution.
- Over-dispersed Poisson and Negative Binomial models showed improved goodness-of-fit over the standard Poisson model, with some evidence of over-fitting.
- Modeling experimental effects improved goodness-of-fit for high-variance genes but exacerbated over-fitting.
Conclusions:
- The Negative Binomial distribution provides a better fit for mRNA-Seq count data due to biological over-dispersion.
- Accurate variance modeling is crucial for reliable statistical analysis and sample size determination.
- Improved analytical strategies will enhance the identification of true biological signals for a better understanding of biological systems.
Related Concept Videos
03:21RNA-seq
RNA sequencing, or RNA-Seq, is a high-throughput sequencing technology used to study the transcriptome of a cell. Transcriptomics helps to interpret the functional elements of a genome and identify the molecular constituents of an organism. Additionally, it also helps in understanding the development of an organism and the occurrence of diseases.
Before the discovery of RNA-seq, microarray-based methods and Sanger sequencing were used for transcriptome analysis. However, while microarray-based...
Before the discovery of RNA-seq, microarray-based methods and Sanger sequencing were used for transcriptome analysis. However, while microarray-based...
02:24Ribosome Profiling
Ribosome profiling or ribo-sequencing is a deep sequencing technique that produces a snapshot of active translation in a cell. It selectively sequences the mRNAs protected by ribosomes to get an insight into a cell’s translation landscape at any given point in time.
Applications of ribosome profiling
Ribosome profiling has many applications, including in vivo monitoring of translation inside a particular organ or tissue type and quantifying new protein synthesis levels.
The technique helps...
Applications of ribosome profiling
Ribosome profiling has many applications, including in vivo monitoring of translation inside a particular organ or tissue type and quantifying new protein synthesis levels.
The technique helps...
02:05Gene Evolution - Fast or Slow?
The genomes of eukaryotes are punctuated by long stretches of sequence which do not code for proteins or RNAs. Although some of these regions do contain crucial regulatory sequences, the vast majority of this DNA serves no known function. Typically, these regions of the genome are the ones in which the fastest change, in evolutionary terms, is observed, because there is typically little to no selection pressure acting on these regions to preserve their sequences.
In contrast, regions which code...
In contrast, regions which code...
01:25Genetic Variation
Genetic variation is the diversity in DNA sequences found among individuals of the same species. This diversity is crucial for a species' survival because it helps organisms adapt to environmental changes. Genetic variation begins with fertilization, where an egg and sperm cell merge. Each of these cells carries 23 chromosomes, up to 46 in the fertilized egg. Chromosomes are long DNA strands that contain genes, the basic units of heredity.
Genes exist in different versions called alleles, which...
Genes exist in different versions called alleles, which...
02:54Evolutionary Relationships through Genome Comparisons
Genome comparison is one of the excellent ways to interpret the evolutionary relationships between organisms. The basic principle of genome comparison is that if two species share a common feature, it is likely encoded by the DNA sequence conserved between both species. The advent of genome sequencing technologies in the late 20th century enabled scientists to understand the concept of conservation of domains between species and helped them to deduce evolutionary relationships across diverse...
02:28Leaky Scanning
During most eukaryotic translation processes, the small 40S ribosome subunit scans an mRNA from its 5' end until it encounters the first start AUG codon. The large 60S ribosomal subunit then joins the smaller one to initiate protein synthesis. The location of the translation initiation is largely determined by the nucleotides near the start codon as there may be multiple translation initiation sites present on the mRNA. Marilyn Kozak discovered that the sequence RCCAUGG (where R stands for...