Technical and biological variance structure in mRNA-Seq data
Ontology highlight
ABSTRACT:
ORGANISM(S): Homo sapiens
SUBMITTER: Brian Bot
PROVIDER: E-GEOD-29022 | biostudies-arrayexpress |
REPOSITORIES: biostudies-arrayexpress
Publications
Technical and biological variance structure in mRNA-Seq data: life in the real world.
Oberg Ann L AL Bot Brian M BM Grill Diane E DE Poland Gregory A GA Therneau Terry M TM
BMC genomics 20120707
<h4>Background</h4>mRNA expression data from next generation sequencing platforms is obtained in the form of counts per gene or exon. Counts have classically been assumed to follow a Poisson distribution in which the variance is equal to the mean. The Negative Binomial distribution which allows for over-dispersion, i.e., for the variance to be greater than the mean, is commonly used to model count data as well.<h4>Results</h4>In mRNA-Seq data from 25 subjects, we found technical variation to gen ...[more]
Similar Datasets
Project description:mRNA expression data from next generation sequencing platforms is obtained in the form of counts per gene or exon. Counts are generally assumed to follow a Poisson distribution in which the variance is equal to the mean. However, it is common in observed count data with biological variation for over-dispersion to be present, i.e., for the variance to be greater than the mean.
2012-07-17 | GSE29022 | GEO
Project description:This SuperSeries is composed of the SubSeries listed below. Refer to individual Series
2014-01-01 | E-GEOD-51748 | biostudies-arrayexpress
Project description:Next-generation sequencing is currently the platform of choice for the discovery and quantification of miRNAs. Despite this, there is no clear consensus on how the data should be pre-processed prior to conducting downstream analyses. Often overlooked, data pre-processing is an essential step in data analysis: the presence of unreliable features and noise can affect the conclusions drawn from downstream analyses. Using a spike-in dilution study, we evaluated the effects of several general-purpose aligners (BWA, Bowtie, Bowtie 2 and Novoalign), and normalization methods (counts-per-million, total count scaling, upper quartile scaling, Trimmed Mean of M, DESeq, linear regression, cyclic loess and quantile) with respect to the final miRNA count data distribution, variance, bias and accuracy of differential expression analysis.
2015-03-21 | GSE67074 | GEO
Project description:ChickenM-BM- ovalbumin upstream promoter-transcription factor II (COUP-TFII; NR2F2) is an orphan nuclear receptor involved in cell-fate specification, organogenesis, angiogenesis and metabolism. Ablation of COUP-TFII in the mouse uterus causes infertility due to defects in embryo attachment and impaired uterine stromal cell decidualization. Although the function of COUP-TFII in uterine decidualization has been described in mice, its role in the human uterus remains unknown. To better elucidate the mechanisms with which COUP-TFII regulates target gene transcription, genome-wide COUP-TFII binding sites in human endometrial stromal cells (HESC) treated with deciduogenic hormones were identified using ChIP-seq. A total of 16,298 intervals (binding regions) for COUP-TFII were identified compared with the input in HESC chromatin with a very low false discovery rate (0.17%) using a stringent cutoff of p =1x10-10. Distribution of intervals showed that more than half (58.6%) of the COUP-TFII binding sites are located within 10 kb of gene boundaries. 7.5% of total intervals reside within the 10 kb promoter region. A total of 6,077 unique genes were identified to have COUP-TFII binding sites within 10 kb of their gene boundaries. Examination of NR2F2 binding in pooled primary human endometrial stromal cells from 6 healthy women upon decidualization with a hormone cocktail of cAMP, E2 and medroxyprogesterone acetate.
2013-11-02 | E-GEOD-52008 | biostudies-arrayexpress
Project description:This SuperSeries is composed of the SubSeries listed below. Refer to individual Series
2014-10-01 | E-GEOD-57280 | biostudies-arrayexpress
Project description:The heterotic hybrid offspring of Arabidopsis accessions C24 and Landsberg erecta have altered methylomes. Changes occur most frequently at loci where parental methylation levels are different. There are context-specific biases in the non-additive methylation patterns with mCG generally increased and mCHH decreased relative to the parents. These changes are a result of two main mechanisms, Trans Chromosomal Methylation (TCM) and Trans Chromosomal deMethylation (TCdM), where the methylation level of one parental allele alters to resemble that of the other parent. Regions of altered methylation are enriched around genic regions and are often correlated with changes in siRNA levels. We identified examples of genes with altered expression likely to be due to methylation changes and suggest that in crosses between the C24 and Ler accessions, epigenetic controls can be important in the generation of altered transcription levels which may contribute to the increased biomass of the hybrids. C24, Ler, and C24 x Ler
2012-02-05 | E-GEOD-35542 | biostudies-arrayexpress
Project description:DNA topoisomerases release supercoils in DNA introduced during replication or transcription. How DNA topoisomerases impact transcription in the context of eukaryotic chromatin is poorly understood. In this study, using a floral stem cell model in Arabidopsis, we uncovered a role of TOP1M-NM-1 in Polycomb Group (PcG) protein-mediated histone lysine 27 trimethylation at, and transcriptional repression of, the stem cell maintenance gene WUSCHEL (WUS). The strong genetic interactions between a top1M-NM-1 mutant and mutations in PcG genes and the overwhelming enrichment of PcG targets among genes affected in expression in the top1a mutant revealed a role of TOP1M-NM-1 in PcG-mediated regulation. Intriguingly, not only the repression of some PcG target genes but also the expression of others requires TOP1M-NM-1. The mechanism that unifies the opposing effects of TOP1M-NM-1 on PcG target genes appears to lie in its role in decreasing nucleosome density. Genome-wide nucleosome mapping shows that TOP1M-NM-1 is required for the depletion of nucleosomes at regulatory regions of genes, which probably allows the binding of factors that either recruit PcG, as we show for AGAMOUS at the WUS locus, or counteract PcG-mediated regulation. This study uncovers a strong and previously unknown connection between TOP1M-NM-1 and PcG. 4 samples in total, two replicates from Ler and top1a-2 plants
2014-10-14 | E-GEOD-53950 | biostudies-arrayexpress
Project description:The counts of total white blood cells (WBCs) and WBC subsets are well-established diagnostic factors for various diseases. It has been shown that variations in WBC counts are significantly controlled by individuals’ genetics in swine. However, despite detection of quantitative trait loci (QTLs) for these phenotypes, little is known on the molecular basis underlying their variations. Our aim was to study gene profiling variations according to variations in WBC counts and to connect results with available QTL mapping. Whole blood transcriptome of animals contrasted for levels of WBC counts were compared. A pig generic microarray enriched in immunity-related genes was used. 378 probes representing 334 genes were found significantly differentially expressed between high- and low-count WBC groups. 65 genes were associated with hematological system development and function. 336 probes could be mapped on all autosomes and the X chromosome, and 59 transcripts fell within 28 QTLs reported to affect the counts of WBC and WBC subsets. By combining probe mapping results and biological functions, 6 genes (CDKN2A, TCIRG1, SIPA1, RGL2, FLT1 and CFLAR) were found as putative relevant positional candidate genes for the WBC traits. Genetic linkage experiments are warranted to validate these candidate genes, and further investigating a possible pleiotropic effect of these genes could contribute to elucidate molecular mechanisms involved in WBC development. differentially expressed genes or transccripts between high- and low-count white blood cells groups The dual channel microarray experiments were carried out using a common reference hybridization design. In this study the reference RNA sample was prepared by pooling of total RNAs from different porcine tissues. The test samples were RNA samples isolated from total porcine blood. According to the total white blood cells count, 18 animals were selected from the edges of the distribution, named high-count group (HC, 9 animals) and low-count group (LC, 9 animals).
2014-10-31 | E-GEOD-50589 | biostudies-arrayexpress
Project description:Here, we examined the ramifications of between-species diversity by documenting the transcriptional response of three marine diatoms - Thalassiosira pseudonana, Fragilariopsis cylindrus, and Pseudo-nitzschia multiseries - to the onset of nitrate limitation of growth, a common limiting nutrient in the ocean. Less than 5% of orthologous genes, shared across the three diatoms, displayed the same transcriptional responses across species when growth was limited by nitrate availability. Orthologs, such as those involved in nitrogen uptake and assimilation, as well as carbon metabolism, were differently expressed across the three species. The two pennate diatoms, F. cylindrus and P. multiseries, shared 3,839 clusters without orthologs in the genome of the centric diatom T. pseudonana. A majority of these pennate-clustered genes, as well as the non-orthologous genes in each species, had minimal annotation information, but were often significantly differentially expressed under nitrate limitation, indicating their potential importance in the response to nitrogen availability. Despite these variations in the specific transcriptional response of each diatom, overall transcriptional patterns suggested that all three diatoms displayed a common physiological response to nitrate limitation that consisted of a general reduction in carbon fixation and carbohydrate and fatty acid metabolism and an increase in nitrogen recycling. Transcriptomes were collected for diatom cultures harvested at the onset of stationary phase in low nitrate media (55 M-NM-<M NaNO3, 212 M-NM-<M Na2SiO3, 72.4 M-NM-<M NaH2PO4) or during mid-exponential growth in nutrient-replete media (882 M-NM-<M NaNO3, 106 M-NM-<M Na2SiO3, 36.2 M-NM-<M NaH2PO4) in artificial seawater, maintaining three biological replicates per condition and per diatom (N=18). The SOLiD sequencer (version 4) was used to generate the transcriptomes and the SEAStAR software package was used to process the SOLiD reads and to calculate gene counts. Pooled counts for the nitrate-limited treatment were normalized to pooled counts for the nutrient-replete M-bM-^@M-^\controlM-bM-^@M-^] treatment to generate log fold changes in gene transcription using the R software package edgeR from Bioconductor.
2014-03-25 | E-GEOD-56132 | biostudies-arrayexpress
Project description:Methods: CaMKIIa-creERT2 (Erdmann et al., 2007) and Dicer1f/f (Harfe et al., 2005) were crossed to produce inducible forebrain-restricted Dicer1 knockout mice (Dicer-ifKO) mice. Hippocampal mRNA profiles of 3-month-old wild-type (WT) and (Dicer-ifKO) mice were generated by deep sequencing, in triplicate, using Illumina HiSeq 2500. Each sample included total RNA isolated from the hippocampus of 3 mice. In total, 12 mice per genotype were used. The sequence reads that passed quality filters were mapped to reference genome (GRCm38/mm10) using Bowtie 2 (2.0.5) and TopHat (2.0.6). SAM/BAM files were further processed with Samtools (0.1.18). Read count quantitations were obtained using Seqmonk (0.26.0). Normalization of read counts and differential expression analysis between genotypes was carried out using DESeq2 R package from Bioconductor (Release 2.13). qRT–PCR validation was performed using SYBR Green assays. Results: We mapped about 13-14 million sequence reads per sample to the mouse genome (build GRCm38/mm10) and quantified 76,938 annotated transcripts. DESeq2 R package was used to normalize the counts and perform the differential expression. Differential analysis output was filtered by FDR threshold (padj < 0.1). This approach led us to identify 641 gene isoforms, corresponding to 314 genes that were differentially regulated in the mouse hippocampus upon Dicer ablation. Conclusions: We extend here the characterization of inducible forebrain-restricted Dicer1 mutants confirming the initial memory improvement. Moreover, we describe several novel phenotypes associated with early Dicer loss in the mature brain including an exacerbated response to seizures, increased CA1 neuron excitability, a pronounced weight gain and enhanced induction of immediate early genes (IEGs) in relevant neuronal nuclei. To identify candidate genes that could explain these phenotypes, we conducted two complementary genomic screens for the miRNAs primarily affected and their targets. Overall, our results explain both the initial and late consequences of Dicer loss in excitatory neurons and indicate that Dicer and the miRNA system play a critical role regulating neuronal homeostasis and responsiveness. Hippocampal mRNA profiles of 3-month-old wild-type (WT) and Dicer-ifKO (3 weeks upon tamoxifen administration) male mice were generated by deep sequencing, in triplicate, using Illumina HiSeq 2500. Each sample included total RNA isolated from the hippocampus of 3 mice. In total, 12 mice per genotype were used.
2015-02-04 | E-GEOD-60261 | biostudies-arrayexpress