Genome-scale chromatin binding dynamics of RNA Polymerase II general transcription machinery components.

A great deal of work has revealed, in structural detail, the components of the preinitiation complex (PIC) machinery required for initiation of mRNA gene transcription by RNA polymerase II (Pol II). However, less-well understood are the in vivo PIC assembly pathways and their kinetics, an understanding of which is vital for determining how rates of in vivo RNA synthesis are established. We used competition ChIP in budding yeast to obtain genome-scale estimates of the residence times for five general transcription factors (GTFs): TBP, TFIIA, TFIIB, TFIIE and TFIIF. While many GTF-chromatin interactions were short-lived ( < 1 min), there were numerous interactions with residence times in the range of several minutes. Sets of genes with a shared function also shared similar patterns of GTF kinetic behavior. TFIIE, a GTF that enters the PIC late in the assembly process, had residence times correlated with RNA synthesis rates. The datasets and results reported here provide kinetic information for most of the Pol II-driven genes in this organism, offering a rich resource for exploring the mechanistic relationships between PIC assembly, gene regulation, and transcription.

Transcriptional Profile of Exercise-Induced Protection Against Relapse to Cocaine Seeking in a Rat Model.

Background: Exercise has shown promise as a treatment for cocaine use disorder; however, the mechanism underlying its efficacy has remained elusive. Methods: We used a rat model of relapse (cue-induced reinstatement) and exercise (wheel running, 2 hours/day) coupled with RNA sequencing to establish transcriptional profiles associated with the protective effects of exercise (during early withdrawal [days 1-7] or throughout withdrawal [days 1-14]) versus noneffective exercise (during late withdrawal [days 8-14]) against cocaine-seeking and sedentary conditions. Results: As expected, cue-induced cocaine seeking was highest in the sedentary and late-withdrawal exercise groups; both groups also showed upregulation of a Grin1-associated transcript and enrichment of Drd1-Nmdar1 complex and glutamate receptor complex terms. Surprisingly, these glutamate markers were also enriched in the early- and throughout-withdrawal exercise groups, despite lower levels of cocaine seeking. However, a closer examination of the Grin1-associated transcript revealed a robust loss of transcripts spanning exons 9 and 10 in the sedentary condition relative to saline controls that was normalized by early- and throughout-withdrawal exercise, but not late-withdrawal exercise, indicating that these exercise conditions may normalize RNA mis-splicing induced by cocaine seeking. Our findings also revealed novel mechanisms by which exercise initiated during early withdrawal may modulate glutamatergic signaling in dorsomedial prefrontal cortex (e.g., via transcripts associated with non-NMDA glutamate receptors or those affecting signaling downstream of NMDA receptors), along with mechanisms outside of glutamatergic signaling such as circadian rhythm regulation and neuronal survival. Conclusions: These findings provide a rich resource for future studies aimed at manipulating these molecular networks to better understand how exercise decreases cocaine seeking.

Globally elevated levels of histone H3 lysine 9 trimethylation in early infancy are associated with poor growth trajectory in Bangladeshi children.

BACKGROUND: Stunting is a global health problem affecting hundreds of millions of children worldwide and contributing to 45% of deaths in children under the age of five. Current therapeutic interventions have limited efficacy. Understanding the epigenetic changes underlying stunting will elucidate molecular mechanisms and likely lead to new therapies. RESULTS: We profiled the repressive mark histone H3 lysine 9 trimethylation (H3K9me3) genome-wide in peripheral blood mononuclear cells (PBMCs) from 18-week-old infants (n = 15) and mothers (n = 14) enrolled in the PROVIDE study established in an urban slum in Bangladesh. We associated H3K9me3 levels within individual loci as well as genome-wide with anthropometric measurements and other biomarkers of stunting and performed functional annotation of differentially affected regions. Despite the relatively small number of samples from this vulnerable population, we observed globally elevated H3K9me3 levels were associated with poor linear growth between birth and one year of age. A large proportion of the differentially methylated genes code for proteins targeting viral mRNA and highly significant regions were enriched in transposon elements with potential regulatory roles in immune system activation and cytokine production. Maternal data show a similar trend with child's anthropometry; however, these trends lack statistical significance to infer an intergenerational relationship. CONCLUSIONS: We speculate that high H3K9me3 levels may result in poor linear growth by repressing genes involved in immune system activation. Importantly, changes to H3K9me3 were detectable before the overt manifestation of stunting and therefore may be valuable as new biomarkers of stunting.

Genome-scale chromatin interaction dynamic measurements for key components of the RNA Pol II general transcription machinery.

Background: A great deal of work has revealed in structural detail the components of the machinery responsible for mRNA gene transcription initiation. These include the general transcription factors (GTFs), which assemble at promoters along with RNA Polymerase II (Pol II) to form a preinitiation complex (PIC) aided by the activities of cofactors and site-specific transcription factors (TFs). However, less well understood are the in vivo PIC assembly pathways and their kinetics, an understanding of which is vital for determining on a mechanistic level how rates of in vivo RNA synthesis are established and how cofactors and TFs impact them. Results: We used competition ChIP to obtain genome-scale estimates of the residence times for five GTFs: TBP, TFIIA, TFIIB, TFIIE and TFIIF in budding yeast. While many GTF-chromatin interactions were short-lived (< 1 min), there were numerous interactions with residence times in the several minutes range. Sets of genes with a shared function also shared similar patterns of GTF kinetic behavior. TFIIE, a GTF that enters the PIC late in the assembly process, had residence times correlated with RNA synthesis rates. Conclusions: The datasets and results reported here provide kinetic information for most of the Pol II-driven genes in this organism and therefore offer a rich resource for exploring the mechanistic relationships between PIC assembly, gene regulation, and transcription. The relationships between gene function and GTF dynamics suggest that shared sets of TFs tune PIC assembly kinetics to ensure appropriate levels of expression.

GenomicDistributions: fast analysis of genomic intervals with Bioconductor.

BACKGROUND: Epigenome analysis relies on defined sets of genomic regions output by widely used assays such as ChIP-seq and ATAC-seq. Statistical analysis and visualization of genomic region sets is essential to answer biological questions in gene regulation. As the epigenomics community continues generating data, there will be an increasing need for software tools that can efficiently deal with more abundant and larger genomic region sets. Here, we introduce GenomicDistributions, an R package for fast and easy summarization and visualization of genomic region data. RESULTS: GenomicDistributions offers a broad selection of functions to calculate properties of genomic region sets, such as feature distances, genomic partition overlaps, and more. GenomicDistributions functions are meticulously optimized for best-in-class speed and generally outperform comparable functions in existing R packages. GenomicDistributions also offers plotting functions that produce editable ggplot objects. All GenomicDistributions functions follow a uniform naming scheme and can handle either single or multiple region set inputs. CONCLUSIONS: GenomicDistributions offers a fast and scalable tool for exploratory genomic region set analysis and visualization. GenomicDistributions excels in user-friendliness, flexibility of outputs, breadth of functions, and computational performance. GenomicDistributions is available from Bioconductor ( https://bioconductor.org/packages/release/bioc/html/GenomicDistributions.html ).

Histone H3 lysine 27 acetylation profile undergoes two global shifts in undernourished children and suggests altered one-carbon metabolism.

BACKGROUND: Stunting is a condition in which a child does not reach their full growth potential due to chronic undernutrition. It arises during the first 2 years of a child's life and is associated with developmental deficiencies and life-long health problems. Current interventions provide some benefit, but new approaches to prevention and treatment grounded in a molecular understanding of stunting are needed. Epigenetic analyses are critical as they can provide insight into how signals from a poor environment lead to changes in cell function. RESULTS: Here we profiled histone H3 acetylation on lysine 27 (H3K27ac) in peripheral blood mononuclear cells (PBMCs) of 18-week-old (n = 14) and 1-year-old children (n = 22) living in an urban slum in Dhaka, Bangladesh. We show that 18-week-old children destined to become stunted have elevated levels of H3K27ac overall, functional analysis of which indicates activation of the immune system and stress response pathways as a primary response to a poor environment with high pathogen load. Conversely, overt stunting at 1-year-of age is associated with globally reduced H3K27ac that is indicative of metabolic rewiring and downregulation of the immune system and DNA repair pathways that are likely secondary responses to chronic exposure to a poor environment with limited nutrients. Among processes altered in 1-year-old children, we identified one-carbon metabolism, the significance of which is supported by integrative analysis with results from histone H3 trimethylation on lysine 4 (H3K4me3). Together, these results suggest altered one-carbon metabolism in this population of stunted children. CONCLUSIONS: The epigenomes of stunted children undergo two global changes in H3K27ac within their first year of life, which are associated with probable initial hyperactive immune responses followed by reduced metabolic capacity. Limitation of one-carbon metabolites may play a key role in the development of stunting. Trial registration ClinicalTrials.gov NCT01375647. Registered 17 June 2011, retrospectively registered, https://clinicaltrials.gov/ct2/show/NCT01375647 .

Histone H3 lysine 4 methylation signature associated with human undernutrition.

Chronically undernourished children become stunted during their first 2 years and thereafter bear burdens of ill health for the rest of their lives. Contributors to stunting include poor nutrition and exposure to pathogens, and parental history may also play a role. However, the epigenetic impact of a poor environment on young children is largely unknown. Here we show the unfolding pattern of histone H3 lysine 4 trimethylation (H3K4me3) in children and mothers living in an urban slum in Dhaka, Bangladesh. A pattern of chromatin modification in blood cells of stunted children emerges over time and involves a global decrease in methylation at canonical locations near gene start sites and increased methylation at ectopic sites throughout the genome. This redistribution occurs at metabolic and immune genes and was specific for H3K4me3, as it was not observed for histone H3 lysine 27 acetylation in the same samples. Methylation changes in stunting globally resemble changes that occur in vitro in response to altered methylation capacity, suggesting that reduced levels of one-carbon nutrients in the diet play a key role in stunting in this population. A network of differentially expressed genes in stunted children reveals effects on chromatin modification machinery, including turnover of H3K4me3, as well as posttranscriptional gene regulation affecting immune response pathways and lipid metabolism. Consistent with these changes, reduced expression of the endocytic receptor gene LDL receptor 1 (LRP1) is a driver of stunting in a mouse model, suggesting a target for intervention.

Transcription profiling of butanol producer Clostridium beijerinckii NRRL B-598 using RNA-Seq.

BACKGROUND: Thinning supplies of natural resources increase attention to sustainable microbial production of bio-based fuels. The strain Clostridium beijerinckii NRRL B-598 is a relatively well-described butanol producer regarding its genotype and phenotype under various conditions. However, a link between these two levels, lying in the description of the gene regulation mechanisms, is missing for this strain, due to the lack of transcriptomic data. RESULTS: In this paper, we present a transcription profile of the strain over the whole fermentation using an RNA-Seq dataset covering six time-points with the current highest dynamic range among solventogenic clostridia. We investigated the accuracy of the genome sequence and particular genome elements, including pseudogenes and prophages. While some pseudogenes were highly expressed, all three identified prophages remained silent. Furthermore, we identified major changes in the transcriptional activity of genes using differential expression analysis between adjacent time-points. We identified functional groups of these significantly regulated genes and together with fermentation and cultivation kinetics captured using liquid chromatography and flow cytometry, we identified basic changes in the metabolism of the strain during fermentation. Interestingly, C. beijerinckii NRRL B-598 demonstrated different behavior in comparison with the closely related strain C. beijerinckii NCIMB 8052 in the latter phases of cultivation. CONCLUSIONS: We provided a complex analysis of the C. beijerinckii NRRL B-598 fermentation profile using several technologies, including RNA-Seq. We described the changes in the global metabolism of the strain and confirmed the uniqueness of its behavior. The whole experiment demonstrated a good reproducibility. Therefore, we will be able to repeat the experiment under selected conditions in order to investigate particular metabolic changes and signaling pathways suitable for following targeted engineering.

Bioinformatics strategies for taxonomy independent binning and visualization of sequences in shotgun metagenomics.

One of main steps in a study of microbial communities is resolving their composition, diversity and function. In the past, these issues were mostly addressed by the use of amplicon sequencing of a target gene because of reasonable price and easier computational postprocessing of the bioinformatic data. With the advancement of sequencing techniques, the main focus shifted to the whole metagenome shotgun sequencing, which allows much more detailed analysis of the metagenomic data, including reconstruction of novel microbial genomes and to gain knowledge about genetic potential and metabolic capacities of whole environments. On the other hand, the output of whole metagenomic shotgun sequencing is mixture of short DNA fragments belonging to various genomes, therefore this approach requires more sophisticated computational algorithms for clustering of related sequences, commonly referred to as sequence binning. There are currently two types of binning methods: taxonomy dependent and taxonomy independent. The first type classifies the DNA fragments by performing a standard homology inference against a reference database, while the latter performs the reference-free binning by applying clustering techniques on features extracted from the sequences. In this review, we describe the strategies within the second approach. Although these strategies do not require prior knowledge, they have higher demands on the length of sequences. Besides their basic principle, an overview of particular methods and tools is provided. Furthermore, the review covers the utilization of the methods in context with the length of sequences and discusses the needs for metagenomic data preprocessing in form of initial assembly prior to binning.