Motivation: Measurement precision determines the power of any analysis to reliably identify significant signals, such as in screens for differential expression, independent of whether the ...experimental design incorporates replicates or not. With the compilation of large-scale RNA-Seq datasets with technical replicate samples, however, we can now, for the first time, perform a systematic analysis of the precision of expression level estimates from massively parallel sequencing technology. This then allows considerations for its improvement by computational or experimental means.
Results: We report on a comprehensive study of target identification and measurement precision, including their dependence on transcript expression levels, read depth and other parameters. In particular, an impressive recall of 84% of the estimated true transcript population could be achieved with 331 million 50 bp reads, with diminishing returns from longer read lengths and even less gains from increased sequencing depths. Most of the measurement power (75%) is spent on only 7% of the known transcriptome, however, making less strongly expressed transcripts harder to measure. Consequently, <30% of all transcripts could be quantified reliably with a relative error <20%. Based on established tools, we then introduce a new approach for mapping and analysing sequencing reads that yields substantially improved performance in gene expression profiling, increasing the number of transcripts that can reliably be quantified to over 40%. Extrapolations to higher sequencing depths highlight the need for efficient complementary steps. In discussion we outline possible experimental and computational strategies for further improvements in quantification precision.
Contact:
rnaseq10@boku.ac.at
Supplementary information:
Supplementary data are available at Bioinformatics online.
We used deep sequencing technology to identify transcriptional adaptation of the euryhaline unicellular cyanobacterium Synechococcus sp. PCC 7002 and the marine facultative aerobe Shewanella ...putrefaciens W3-18-1 to growth in a co-culture and infer the effect of carbon flux distributions on photoautotroph-heterotroph interactions. The overall transcriptome response of both organisms to co-cultivation was shaped by their respective physiologies and growth constraints. Carbon limitation resulted in the expansion of metabolic capacities, which was manifested through the transcriptional upregulation of transport and catabolic pathways. Although growth coupling occurred via lactate oxidation or secretion of photosynthetically fixed carbon, there was evidence of specific metabolic interactions between the two organisms. These hypothesized interactions were inferred from the excretion of specific amino acids (for example, alanine and methionine) by the cyanobacterium, which correlated with the downregulation of the corresponding biosynthetic machinery in Shewanella W3-18-1. In addition, the broad and consistent decrease of mRNA levels for many Fe-regulated Synechococcus 7002 genes during co-cultivation may indicate increased Fe availability as well as more facile and energy-efficient mechanisms for Fe acquisition by the cyanobacterium. Furthermore, evidence pointed at potentially novel interactions between oxygenic photoautotrophs and heterotrophs related to the oxidative stress response as transcriptional patterns suggested that Synechococcus 7002 rather than Shewanella W3-18-1 provided scavenging functions for reactive oxygen species under co-culture conditions. This study provides an initial insight into the complexity of photoautotrophic-heterotrophic interactions and brings new perspectives of their role in the robustness and stability of the association.
Geologic carbon dioxide (CO2) sequestration drives physical and geochemical changes in deep subsurface environments that impact indigenous microbial activities. The combined effects of pressurized ...CO2 on a model sulfate-reducing microorganism, Desulfovibrio vulgaris, have been assessed using a suite of genomic and kinetic measurements. Novel high-pressure NMR time-series measurements using (13)C-lactate were used to track D. vulgaris metabolism. We identified cessation of respiration at CO2 pressures of 10 bar, 25 bar, 50 bar, and 80 bar. Concurrent experiments using N2 as the pressurizing phase had no negative effect on microbial respiration, as inferred from reduction of sulfate to sulfide. Complementary pressurized batch incubations and fluorescence microscopy measurements supported NMR observations, and indicated that non-respiring cells were mostly viable at 50 bar CO2 for at least 4 h, and at 80 bar CO2 for 2 h. The fraction of dead cells increased rapidly after 4 h at 80 bar CO2. Transcriptomic (RNA-Seq) measurements on mRNA transcripts from CO2-incubated biomass indicated that cells up-regulated the production of certain amino acids (leucine, isoleucine) following CO2 exposure at elevated pressures, likely as part of a general stress response. Evidence for other poorly understood stress responses were also identified within RNA-Seq data, suggesting that while pressurized CO2 severely limits the growth and respiration of D. vulgaris cells, biomass retains intact cell membranes at pressures up to 80 bar CO2. Together, these data show that geologic sequestration of CO2 may have significant impacts on rates of sulfate reduction in many deep subsurface environments where this metabolism is a key respiratory process.
To understand how cell physiological state affects mRNA translation, we used Shewanella oneidensis MR-1 grown under steady state conditions at either 20% or 8.5% O2. Using a combination of ...quantitative proteomics and RNA-Seq, we generated high-confidence data on >1000 mRNA and protein pairs. By using a steady state model, we found that differences in protein-mRNA ratios were primarily due to differences in the translational efficiency of specific genes. When oxygen levels were lowered, 28% of the proteins showed at least a 2-fold change in expression. Transcription levels were sp. significantly altered for 26% of the protein changes; translational efficiency was significantly altered for 46% and a combination of both was responsible for the remaining 28%. Changes in translational efficiency were significantly correlated with the codon usage pattern of the genes and measurable tRNA pools changed in response to altered O2 levels. Our results suggest that changes in the translational efficiency of proteins, in part due to altered tRNA pools, is a major determinant of regulated alterations in protein expression levels in bacteria.
The isoprenoid pathway converts pyruvate to isoprene and related isoprenoid compounds in plants and some bacteria. Currently, this pathway is of great interest because of the critical role that ...isoprenoids play in basic cellular processes, as well as the industrial value of metabolites such as isoprene. Although the regulation of several pathway genes has been described, there is a paucity of information regarding system level regulation and control of the pathway. To address these limitations, we examined Bacillus subtilis grown under multiple conditions and determined the relationship between altered isoprene production and gene expression patterns. We found that with respect to the amount of isoprene produced, terpenoid genes fall into two distinct subsets with opposing correlations. The group whose expression levels positively correlated with isoprene production included dxs, which is responsible for the commitment step in the pathway, ispD, and two genes that participate in the mevalonate pathway, yhfS and pksG. The subset of terpenoid genes that inversely correlated with isoprene production included ispH, ispF, hepS, uppS, ispE, and dxr. A genome-wide partial least squares regression model was created to identify other genes or pathways that contribute to isoprene production. These analyses showed that a subset of 213 regulated genes was sufficient to create a predictive model of isoprene production under different conditions and showed correlations at the transcriptional level. We conclude that gene expression levels alone are sufficiently informative about the metabolic state of a cell that produces increased isoprene and can be used to build a model that accurately predicts production of this secondary metabolite across many simulated environmental conditions.
DNA Methylation as a Biomarker for Preeclampsia Anderson, Cindy M.; Ralph, Jody L.; Wright, Michelle L. ...
Biological research for nursing,
10/2014, Letnik:
16, Številka:
4
Journal Article
Recenzirano
Background:
Preeclampsia contributes significantly to pregnancy-associated morbidity and mortality as well as future risk of cardiovascular disease in mother and offspring, and preeclampsia in ...offspring. The lack of reliable methods for early detection limits the opportunities for prevention, diagnosis, and timely treatment.
Purpose:
The purpose of this study was to explore distinct DNA methylation patterns associated with preeclampsia in both maternal cells and fetal-derived tissue that represent potential biomarkers to predict future preeclampsia and inheritance in children.
Method:
A convenience sample of nulliparous women (N = 55) in the first trimester of pregnancy was recruited for this prospective study. Genome-wide DNA methylation was quantified in first-trimester maternal peripheral white blood cells and placental chorionic tissue from normotensive women and those with preeclampsia (n = 6/group).
Results:
Late-onset preeclampsia developed in 12.7% of women. Significant differences in DNA methylation were identified in 207 individual linked cytosine and guanine (CpG) sites in maternal white blood cells collected in the first trimester (132 sites with gain and 75 sites with loss of methylation), which were common to approximately 75% of the differentially methylated CpG sites identified in chorionic tissue of fetal origin.
Conclusion:
This study is the first to identify maternal epigenetic targets and common targets in fetal-derived tissue that represent putative biomarkers for early detection and heritable risk of preeclampsia. Findings may pave the way for diagnosis of preeclampsia prior to its clinical presentation and acute damaging effects, and the potential for prevention of the detrimental long-term sequelae.
With currently available RNA-Seq pipelines, expression estimates for most genes are very noisy. We here introduce MapAl, a tool for RNA-Seq expression profiling that builds on the established ...programs Bowtie and Cufflinks. In the post-processing of RNA-Seq reads, it incorporates gene models already at the stage of read alignment, increasing the number of reliably measured known transcripts consistently by 50%. Adding genes identified de novo then allows a reliable assessment of double the total number of transcripts compared to other available pipelines. This substantial improvement is of general relevance: Measurement precision determines the power of any analysis to reliably identify significant signals, such as in screens for differential expression, independent of whether the experimental design incorporates replicates or not.
To understand how cell physiological state affects mRNA translation, we used
Shewanella oneidensis MR-1
grown under steady state conditions at either 20% or 8.5% O
2
. Using a combination of ...quantitative proteomics and RNA-Seq, we generated high-confidence data on >1000 mRNA and protein pairs. By using a steady state model, we found that differences in protein-mRNA ratios were primarily due to differences in the translational efficiency of specific genes. When oxygen levels were lowered, 28% of the proteins showed at least a 2-fold change in expression. Transcription levels were sp. significantly altered for 26% of the protein changes; translational efficiency was significantly altered for 46% and a combination of both was responsible for the remaining 28%. Changes in translational efficiency were significantly correlated with the codon usage pattern of the genes and measurable tRNA pools changed in response to altered O
2
levels. Our results suggest that changes in the translational efficiency of proteins, in part due to altered tRNA pools, is a major determinant of regulated alterations in protein expression levels in bacteria.
We show that environmentally-induced changes in protein abundance in bacteria is primarily due to their alterations in the translational efficiency of mRNAs.
Inflammatory bowel diseases (IBD), including Crohn's disease (CD) and ulcerative colitis (UC), are complex chronic inflammatory conditions of the gastrointestinal tract that are driven by perturbed ...cytokine pathways. Anti-tumor necrosis factor-α (TNF) antibodies are mainstay therapies for IBD. However, up to 40% of patients are nonresponsive to anti-TNF agents, which makes the identification of alternative therapeutic targets a priority. Here we show that, relative to healthy controls, inflamed intestinal tissues from patients with IBD express high amounts of the cytokine oncostatin M (OSM) and its receptor (OSMR), which correlate closely with histopathological disease severity. The OSMR is expressed in nonhematopoietic, nonepithelial intestinal stromal cells, which respond to OSM by producing various proinflammatory molecules, including interleukin (IL)-6, the leukocyte adhesion factor ICAM1, and chemokines that attract neutrophils, monocytes, and T cells. In an animal model of anti-TNF-resistant intestinal inflammation, genetic deletion or pharmacological blockade of OSM significantly attenuates colitis. Furthermore, according to an analysis of more than 200 patients with IBD, including two cohorts from phase 3 clinical trials of infliximab and golimumab, high pretreatment expression of OSM is strongly associated with failure of anti-TNF therapy. OSM is thus a potential biomarker and therapeutic target for IBD, and has particular relevance for anti-TNF-resistant patients.