Two strategies for gene regulation by promoter nucleosomes

The impact of phenotypic heterogeneity on fungal pathogenicity and drug resistance

Phenotypic heterogeneity in genetically clonal populations facilitates cellular adaptation to adverse environmental conditions while enabling a return to the basal physiological state. It also plays a crucial role in pathogenicity and the acquisition of drug resistance in unicellular organisms and cancer cells, yet the exact contributing factors remain elusive. In this review, we outline the current state of understanding concerning the contribution of phenotypic heterogeneity to fungal pathogenesis and antifungal drug resistance.

Revisiting the model for coactivator recruitment: Med15 can select its target sites independent of promoter-bound transcription factors

Activation domains (ADs) within transcription factors (TFs) induce gene expression by recruiting coactivators such as the Mediator complex. Coactivators lack DNA binding domains (DBDs) and are assumed to passively follow their recruiting TFs. This is supported by direct AD-coactivator interactions seen in vitro but has not yet been tested in living cells. To examine that, we targeted two Med15-recruiting ADs to a range of budding yeast promoters through fusion with different DBDs. The DBD-AD fusions localized to hundreds of genomic sites but recruited Med15 and induced transcription in only a subset of bound promoters, characterized by a fuzzy-nucleosome architecture. Direct DBD-Med15 fusions shifted DBD localization towards fuzzy-nucleosome promoters, including promoters devoid of the endogenous Mediator. We propose that Med15, and perhaps other coactivators, possess inherent promoter preference and thus actively contribute to the selection of TF-induced genes.

Living in a noisy world-origins of gene expression noise and its impact on cellular decision-making

The expression level of a gene can vary between genetically identical cells under the same environmental condition-a phenomenon referred to as gene expression noise. Several studies have now elucidated a central role of transcription factors in the generation of expression noise. Transcription factors, as the key components of gene regulatory networks, drive many important cellular decisions in response to cellular and environmental signals. Therefore, a very relevant question is how expression noise impacts gene regulation and influences cellular decision-making. In this Review, we summarize the current understanding of the molecular origins of expression noise, highlighting the role of transcription factors in this process, and discuss the ways in which noise can influence cellular decision-making. As advances in single-cell technologies open new avenues for studying expression noise as well as gene regulatory circuits, a better understanding of the influence of noise on cellular decisions will have important implications for many biological processes.

Contributions of transcriptional noise to leukaemia evolution: KAT2A as a case-study

Transcriptional noise is proposed to participate in cell fate changes, but contributions to mammalian cell differentiation systems, including cancer, remain associative. Cancer evolution is driven by genetic variability, with modulatory or contributory participation of epigenetic variants. Accumulation of epigenetic variants enhances transcriptional noise, which can facilitate cancer cell fate transitions. Acute myeloid leukaemia (AML) is an aggressive cancer with strong epigenetic dependencies, characterized by blocked differentiation. It constitutes an attractive model to probe links between transcriptional noise and malignant cell fate regulation. Gcn5/KAT2A is a classical epigenetic transcriptional noise regulator. Its loss increases transcriptional noise and modifies cell fates in stem and AML cells. By reviewing the analysis of KAT2A-depleted pre-leukaemia and leukaemia models, I discuss that the net result of transcriptional noise is diversification of cell fates secondary to alternative transcriptional programmes. Cellular diversification can enable or hinder AML progression, respectively, by differentiation of cell types responsive to mutations, or by maladaptation of leukaemia stem cells. KAT2A-dependent noise-responsive genes participate in ribosome biogenesis and KAT2A loss destabilizes translational activity. I discuss putative contributions of perturbed translation to AML biology, and propose KAT2A loss as a model for mechanistic integration of transcriptional and translational control of noise and fate decisions. This article is part of a discussion meeting issue 'Causes and consequences of stochastic processes in development and disease'.

Characterizing Different Modes of Interplay Between Rap1 and H3 Using Inducible H3-depletion Yeast

Histones and transcription factors (TFs) are two important DNA-binding proteins that interact, compete, and together regulate transcriptional processes in response to diverse internal and external stimuli. Condition-specific depletion of histones in Saccharomyces cerevisiae using a galactose-inducible H3 promoter provides a suitable framework for examining transcriptional alteration resulting from reduced nucleosome content. However, the effect on DNA binding activities of TFs is yet to be fully explored. In this work, we combine ChIP-seq of H3 with RNA-seq to elucidate the genome-scale relationships between H3 occupancy patterns and transcriptional dynamics before and after global H3 depletion. ChIP-seq of Rap1 is also conducted in the H3-depletion and control treatments, to investigate the interplay between this master regulator TF and nucleosomal H3, and to explore the impact on diverse transcriptional responses of different groups of target genes and functions. Ultimately, we propose a working model and testable hypotheses regarding the impact of global and local H3 depletion on transcriptional changes. We also demonstrate different potential modes of interaction between Rap1 and H3, which sheds light on the potential multifunctional regulatory capabilities of Rap1 and potentially other pioneer factors.

Chromatin analysis of adult pluripotent stem cells reveals a unique stemness maintenance strategy

Many highly regenerative organisms maintain adult pluripotent stem cells throughout their life, but how the long-term maintenance of pluripotency is accomplished is unclear. To decipher the regulatory logic of adult pluripotent stem cells, we analyzed the chromatin organization of stem cell genes in the planarian Schmidtea mediterranea. We identify a special chromatin state of stem cell genes, which is distinct from that of tissue-specific genes and resembles constitutive genes. Where tissue-specific promoters have detectable transcription factor binding sites, the promoters of stem cell-specific genes instead have sequence features that broadly decrease nucleosome binding affinity. This genic organization makes pluripotency-related gene expression the default state in these cells, which is maintained by the activity of chromatin remodelers ISWI and SNF2 in the stem cells.

The architecture of binding cooperativity between densely bound transcription factors

The binding of transcription factors (TFs) along genomes is restricted to a subset of sites containing their preferred motifs. TF-binding specificity is often attributed to the co-binding of interacting TFs; however, apart from specific examples, this model remains untested. Here, we define dependencies among budding yeast TFs that localize to overlapping promoters by profiling the genome-wide consequences of co-depleting multiple TFs. We describe unidirectional interactions, revealing Msn2 as a central factor allowing TF binding at its target promoters. By contrast, no case of mutual cooperation was observed. Particularly, Msn2 retained binding at its preferred promoters upon co-depletion of fourteen similarly bound TFs. Overall, the consequences of TF co-depletions were moderate, limited to a subset of promoters, and failed to explain the role of regions outside the DNA-binding domain in directing TF-binding preferences. Our results call for re-evaluating the role of cooperative interactions in directing TF-binding preferences.

Complementary strategies for directing in vivo transcription factor binding through DNA binding domains and intrinsically disordered regions

DNA binding domains (DBDs) of transcription factors (TFs) recognize DNA sequence motifs that are highly abundant in genomes. Within cells, TFs bind a subset of motif-containing sites as directed by either their DBDs or DBD-external (nonDBD) sequences. To define the relative roles of DBDs and nonDBDs in directing binding preferences, we compared the genome-wide binding of 48 (∼30%) budding yeast TFs with their DBD-only, nonDBD-truncated, and nonDBD-only mutants. With a few exceptions, binding locations differed between DBDs and TFs, resulting from the cumulative action of multiple determinants mapped mostly to disordered nonDBD regions. Furthermore, TFs' preferences for promoters of the fuzzy nucleosome architecture were lost in DBD-only mutants, whose binding spread across promoters, implicating nonDBDs' preferences in this hallmark of budding yeast regulatory design. We conclude that DBDs and nonDBDs employ complementary DNA-targeting strategies, whose balance defines TF binding specificity along genomes.

Nucleosome Remodeling at the Yeast PHO8 and PHO84 Promoters without the Putatively Essential SWI/SNF Remodeler

Chromatin remodeling by ATP-dependent remodeling enzymes is crucial for all genomic processes, like transcription or replication. Eukaryotes harbor many remodeler types, and it is unclear why a given chromatin transition requires more or less stringently one or several remodelers. As a classical example, removal of budding yeast PHO8 and PHO84 promoter nucleosomes upon physiological gene induction by phosphate starvation essentially requires the SWI/SNF remodeling complex. This dependency on SWI/SNF may indicate specificity in remodeler recruitment, in recognition of nucleosomes as remodeling substrate or in remodeling outcome. By in vivo chromatin analyses of wild type and mutant yeast under various PHO regulon induction conditions, we found that overexpression of the remodeler-recruiting transactivator Pho4 allowed removal of PHO8 promoter nucleosomes without SWI/SNF. For PHO84 promoter nucleosome removal in the absence of SWI/SNF, an intranucleosomal Pho4 site, which likely altered the remodeling outcome via factor binding competition, was required in addition to such overexpression. Therefore, an essential remodeler requirement under physiological conditions need not reflect substrate specificity, but may reflect specific recruitment and/or remodeling outcomes.

Cis- and Trans-variations of Stearoyl-CoA Desaturase Provide New Insights into the Mechanisms of Diverged Pattern of Phenotypic Plasticity for Temperature Adaptation in Two Congeneric Oyster Species

The evolution of phenotypic plasticity plays an essential role in adaptive responses to climate change; however, its regulatory mechanisms in marine organisms which exhibit high phenotypic plasticity still remain poorly understood. The temperature-responsive trait oleic acid content and its major gene stearoyl-CoA desaturase (Scd) expression have diverged in two allopatric congeneric oyster species, cold-adapted Crassostrea gigas and warm-adapted Crassostrea angulata. In this study, genetic and molecular methods were used to characterize fatty acid desaturation and membrane fluidity regulated by oyster Scd. Sixteen causative single-nucleotide polymorphisms (SNPs) were identified in the promoter/cis-region of the Scd between wild C. gigas and C. angulata. Further functional experiments showed that an SNP (g.-333C [C. gigas allele] >T [C. angulata allele]) may influence Scd transcription by creating/disrupting the binding motif of the positive trans-factor Y-box factor in C. gigas/C. angulata, which mediates the higher/lower constitutive expression of Scd in C. gigas/C. angulata. Additionally, the positive trans-factor sterol-regulatory element-binding proteins (Srebp) were identified to specifically bind to the promoter of Scd in both species, and were downregulated during cold stress in C. gigas compared to upregulated in C. angulata. This partly explains the relatively lower environmental sensitivity (plasticity) of Scd in C. gigas. This study serves as an experimental case to reveal that both cis- and trans-variations shape the diverged pattern of phenotypic plasticity, which provides new insights into the formation of adaptive traits and the prediction of the adaptive potential of marine organisms to future climate change.

Transcription factor binding process is the primary driver of noise in gene expression

Noise in expression of individual genes gives rise to variations in activity of cellular pathways and generates heterogeneity in cellular phenotypes. Phenotypic heterogeneity has important implications for antibiotic persistence, mutation penetrance, cancer growth and therapy resistance. Specific molecular features such as the presence of the TATA box sequence and the promoter nucleosome occupancy have been associated with noise. However, the relative importance of these features in noise regulation is unclear and how well these features can predict noise has not yet been assessed. Here through an integrated statistical model of gene expression noise in yeast we found that the number of regulating transcription factors (TFs) of a gene was a key predictor of noise, whereas presence of the TATA box and the promoter nucleosome occupancy had poor predictive power. With an increase in the number of regulatory TFs, there was a rise in the number of cooperatively binding TFs. In addition, an increased number of regulatory TFs meant more overlaps in TF binding sites, resulting in competition between TFs for binding to the same region of the promoter. Through modeling of TF binding to promoter and application of stochastic simulations, we demonstrated that competition and cooperation among TFs could increase noise. Thus, our work uncovers a process of noise regulation that arises out of the dynamics of gene regulation and is not dependent on any specific transcription factor or specific promoter sequence.

YMLA: A comparative platform to carry out functional enrichment analysis for multiple gene lists in yeast

Comparative analysis among multiple gene lists on their functional features is now a routine task due to the advancement of high-throughput experiments. Several enrichment analysis tools were developed in the past. However, these tools mainly focus on one gene list and contain only gene ontology or interaction features. What makes it worse, comparative investigation and customized feature set reanalysis are still unavailable. Therefore, we constructed the YMLA (Yeast Multiple List Analyzer) platform in this research. YMLA includes 39 yeast features and facilitates comparative analysis among multiple gene lists via tabular views, heatmaps, and network plots. Moreover, the customized feature set reanalysis function was implemented in YMLA to help form mechanism hypotheses based on a selected enriched feature subset. We demonstrated the biological applicability of YMLA via example lists consisting of genes with top/bottom translation efficiency values. The analysis results provided by YMLA reveal novel facts consistent with previous experiments. YMLA is available at https://cosbi7.ee.ncku.edu.tw/YMLA/.

Accurate prediction of functional states of cis-regulatory modules reveals common epigenetic rules in humans and mice

BACKGROUND

Predicting cis-regulatory modules (CRMs) in a genome and their functional states in various cell/tissue types of the organism are two related challenging computational tasks. Most current methods attempt to simultaneously achieve both using data of multiple epigenetic marks in a cell/tissue type. Though conceptually attractive, they suffer high false discovery rates and limited applications. To fill the gaps, we proposed a two-step strategy to first predict a map of CRMs in the genome, and then predict functional states of all the CRMs in various cell/tissue types of the organism. We have recently developed an algorithm for the first step that was able to more accurately and completely predict CRMs in a genome than existing methods by integrating numerous transcription factor ChIP-seq datasets in the organism. Here, we presented machine-learning methods for the second step.

RESULTS

We showed that functional states in a cell/tissue type of all the CRMs in the genome could be accurately predicted using data of only 1~4 epigenetic marks by a variety of machine-learning classifiers. Our predictions are substantially more accurate than the best achieved so far. Interestingly, a model trained on a cell/tissue type in humans can accurately predict functional states of CRMs in different cell/tissue types of humans as well as of mice, and vice versa. Therefore, epigenetic code that defines functional states of CRMs in various cell/tissue types is universal at least in humans and mice. Moreover, we found that from tens to hundreds of thousands of CRMs were active in a human and mouse cell/tissue type, and up to 99.98% of them were reutilized in different cell/tissue types, while as small as 0.02% of them were unique to a cell/tissue type that might define the cell/tissue type.

CONCLUSIONS

Our two-step approach can accurately predict functional states in any cell/tissue type of all the CRMs in the genome using data of only 1~4 epigenetic marks. Our approach is also more cost-effective than existing methods that typically use data of more epigenetic marks. Our results suggest common epigenetic rules for defining functional states of CRMs in various cell/tissue types in humans and mice.

An integrated SAGA and TFIID PIC assembly pathway selective for poised and induced promoters

Genome-wide, little is understood about how proteins organize at inducible promoters before and after induction and to what extent inducible and constitutive architectures depend on cofactors. We report that sequence-specific transcription factors and their tethered cofactors (e.g., SAGA [Spt-Ada-Gcn5-acetyltransferase], Mediator, TUP, NuA4, SWI/SNF, and RPD3-L) are generally bound to promoters prior to induction ("poised"), rather than recruited upon induction, whereas induction recruits the preinitiation complex (PIC) to DNA. Through depletion and/or deletion experiments, we show that SAGA does not function at constitutive promoters, although a SAGA-independent Gcn5 acetylates +1 nucleosomes there. When inducible promoters are poised, SAGA catalyzes +1 nucleosome acetylation but not PIC assembly. When induced, SAGA catalyzes acetylation, deubiquitylation, and PIC assembly. Surprisingly, SAGA mediates induction by creating a PIC that allows TFIID (transcription factor II-D) to stably associate, rather than creating a completely TFIID-independent PIC, as generally thought. These findings suggest that inducible systems, where present, are integrated with constitutive systems.

Transcriptional regulator Taf14 binds DNA and is required for the function of transcription factor TFIID in the absence of histone H2A.Z

The transcriptional regulator Taf14 is a component of multiple protein complexes involved in transcription initiation and chromatin remodeling in yeast cells. Although Taf14 is not required for cell viability, it becomes essential in conditions where the formation of the transcription preinitiation complex is hampered. The specific role of Taf14 in mediating transcription initiation and preinitiation complex formation is unclear. Here, we explored its role in the general transcription factor IID by mapping Taf14 genetic and proteomic interactions and found that it was needed for the function of the complex if Htz1, the yeast homolog of histone H2A.Z, was absent from chromatin. Dissecting the functional domains of Taf14 revealed that the linker region between the YEATS and ET domains was required for cell viability in the absence of Htz1 protein. We further show that the linker region of Taf14 interacts with DNA. We propose that providing additional DNA binding capacity might be a general role of Taf14 in the recruitment of protein complexes to DNA and chromatin.

Antisense-mediated repression of SAGA-dependent genes involves the HIR histone chaperone

Eukaryotic genomes are pervasively transcribed by RNA polymerase II (RNAPII), and transcription of long non-coding RNAs often overlaps with coding gene promoters. This might lead to coding gene repression in a process named Transcription Interference (TI). In Saccharomyces cerevisiae, TI is mainly driven by antisense non-coding transcription and occurs through re-shaping of promoter Nucleosome-Depleted Regions (NDRs). In this study, we developed a genetic screen to identify new players involved in Antisense-Mediated Transcription Interference (AMTI). Among the candidates, we found the HIR histone chaperone complex known to be involved in de novo histone deposition. Using genome-wide approaches, we reveal that HIR-dependent histone deposition represses the promoters of SAGA-dependent genes via antisense non-coding transcription. However, while antisense transcription is enriched at promoters of SAGA-dependent genes, this feature is not sufficient to define the mode of gene regulation. We further show that the balance between HIR-dependent nucleosome incorporation and transcription factor binding at promoters directs transcription into a SAGA- or TFIID-dependent regulation. This study sheds light on a new connection between antisense non-coding transcription and the nature of coding transcription initiation.

Identification of differentially distributed gene expression and distinct sets of cancer-related genes identified by changes in mean and variability

There is increasing evidence that changes in the variability or overall distribution of gene expression are important both in normal biology and in diseases, particularly cancer. Genes whose expression differs in variability or distribution without a difference in mean are ignored by traditional differential expression-based analyses. Using a Bayesian hierarchical model that provides tests for both differential variability and differential distribution for bulk RNA-seq data, we report here an investigation into differential variability and distribution in cancer. Analysis of eight paired tumour-normal datasets from The Cancer Genome Atlas confirms that differential variability and distribution analyses are able to identify cancer-related genes. We further demonstrate that differential variability identifies cancer-related genes that are missed by differential expression analysis, and that differential expression and differential variability identify functionally distinct sets of potentially cancer-related genes. These results suggest that differential variability analysis may provide insights into genetic aspects of cancer that would not be revealed by differential expression, and that differential distribution analysis may allow for more comprehensive identification of cancer-related genes than analyses based on changes in mean or variability alone.

Control of Gene Expression via the Yeast CWI Pathway

Living cells exposed to stressful environmental situations can elicit cellular responses that guarantee maximal cell survival. Most of these responses are mediated by mitogen-activated protein kinase (MAPK) cascades, which are highly conserved from yeast to humans. Cell wall damage conditions in the yeast Saccharomyces cerevisiae elicit rescue mechanisms mainly associated with reprogramming specific transcriptional responses via the cell wall integrity (CWI) pathway. Regulation of gene expression by this pathway is coordinated by the MAPK Slt2/Mpk1, mainly via Rlm1 and, to a lesser extent, through SBF (Swi4/Swi6) transcription factors. In this review, we summarize the molecular mechanisms controlling gene expression upon cell wall stress and the role of chromatin structure in these processes. Some of these mechanisms are also discussed in the context of other stresses governed by different yeast MAPK pathways. Slt2 regulates both transcriptional initiation and elongation by interacting with chromatin at the promoter and coding regions of CWI-responsive genes but using different mechanisms for Rlm1- and SBF-dependent genes. Since MAPK pathways are very well conserved in eukaryotic cells and are essential for controlling cellular physiology, improving our knowledge regarding how they regulate gene expression could impact the future identification of novel targets for therapeutic intervention.

Genome-wide gene expression noise in Escherichia coli is condition-dependent and determined by propagation of noise through the regulatory network

Although it is well appreciated that gene expression is inherently noisy and that transcriptional noise is encoded in a promoter’s sequence, little is known about the extent to which noise levels of individual promoters vary across growth conditions. Using flow cytometry, we here quantify transcriptional noise in Escherichia coli genome-wide across 8 growth conditions and find that noise levels systematically decrease with growth rate, with a condition-dependent lower bound on noise. Whereas constitutive promoters consistently exhibit low noise in all conditions, regulated promoters are both more noisy on average and more variable in noise across conditions. Moreover, individual promoters show highly distinct variation in noise across conditions. We show that a simple model of noise propagation from regulators to their targets can explain a significant fraction of the variation in relative noise levels and identifies TFs that most contribute to both condition-specific and condition-independent noise propagation. In addition, analysis of the genome-wide correlation structure of various gene properties shows that gene regulation, expression noise, and noise plasticity are all positively correlated genome-wide and vary independently of variations in absolute expression, codon bias, and evolutionary rate. Together, our results show that while absolute expression noise tends to decrease with growth rate, relative noise levels of genes are highly condition-dependent and determined by the propagation of noise through the gene regulatory network.

Genome-wide flow cytometry measurements reveal that gene expression noise in bacteria is highly condition-dependent; while absolute noise levels of all genes decrease with growth-rate, theoretical modeling shows that the relative noise levels of different genes are determined by the propagation of noise through the gene regulatory network (GRN). Thus GRN structure controls not only mean expression but also noise levels.

Mediator dynamics during heat shock in budding yeast

The Mediator complex is central to transcription by RNA polymerase II (Pol II) in eukaryotes. In budding yeast (Saccharomyces cerevisiae), Mediator is recruited by activators and associates with core promoter regions, where it facilitates preinitiation complex (PIC) assembly, only transiently before Pol II escape. Interruption of the transcription cycle by inactivation or depletion of Kin28 inhibits Pol II escape and stabilizes this association. However, Mediator occupancy and dynamics have not been examined on a genome-wide scale in yeast grown in nonstandard conditions. Here we investigate Mediator occupancy following heat shock or CdCl₂ exposure, with and without depletion of Kin28. We find that Pol II occupancy shows similar dependence on Mediator under normal and heat shock conditions. However, although Mediator association increases at many genes upon Kin28 depletion under standard growth conditions, little or no increase is observed at most genes upon heat shock, indicating a more stable association of Mediator after heat shock. Unexpectedly, Mediator remains associated upstream of the core promoter at genes repressed by heat shock or CdCl₂ exposure whether or not Kin28 is depleted, suggesting that Mediator is recruited by activators but is unable to engage PIC components at these repressed targets. This persistent association is strongest at promoters that bind the HMGB family member Hmo1, and is reduced but not eliminated in hmo1Δ yeast. Finally, we show a reduced dependence on PIC components for Mediator occupancy at promoters after heat shock, further supporting altered dynamics or stronger engagement with activators under these conditions.

Mutational sources of trans-regulatory variation affecting gene expression in Saccharomyces cerevisiae

Heritable variation in a gene’s expression arises from mutations impacting cis- and trans-acting components of its regulatory network. Here, we investigate how trans-regulatory mutations are distributed within the genome and within a gene regulatory network by identifying and characterizing 69 mutations with trans-regulatory effects on expression of the same focal gene in Saccharomyces cerevisiae. Relative to 1766 mutations without effects on expression of this focal gene, we found that these trans-regulatory mutations were enriched in coding sequences of transcription factors previously predicted to regulate expression of the focal gene. However, over 90% of the trans-regulatory mutations identified mapped to other types of genes involved in diverse biological processes including chromatin state, metabolism, and signal transduction. These data show how genetic changes in diverse types of genes can impact a gene’s expression in trans, revealing properties of trans-regulatory mutations that provide the raw material for trans-regulatory variation segregating within natural populations.

Uncovering Divergence in Gene Expression Regulation in the Adaptation of Yeast to Nitrogen Scarcity

Saccharomyces cerevisiae rewires its transcriptional output to survive stressful environments, such as nitrogen scarcity under fermentative conditions. Although divergence in nitrogen metabolism among natural yeast populations has been reported, the impact of regulatory genetic variants modulating gene expression and nitrogen consumption remains to be investigated. Here, we employed an F1 hybrid from two contrasting S. cerevisiae strains, providing a controlled genetic environment to map cis factors involved in the divergence of gene expression regulation in response to nitrogen scarcity. We used a dual approach to obtain genome-wide allele-specific profiles of chromatin accessibility, transcription factor binding, and gene expression through ATAC-seq (assay for transposase accessible chromatin) and RNA-seq (transcriptome sequencing). We observed large variability in allele-specific expression and accessibility between the two genetic backgrounds, with a third of these differences specific to a deficient nitrogen environment. Furthermore, we discovered events of allelic bias in gene expression correlating with allelic bias in transcription factor binding solely under nitrogen scarcity, where the majority of these transcription factors orchestrates the nitrogen catabolite repression regulatory pathway and demonstrates a cis × environment-specific response. Our approach allowed us to find cis variants modulating gene expression, chromatin accessibility, and allelic differences in transcription factor binding in response to low nitrogen culture conditions. IMPORTANCE Historically, coding variants were prioritized when searching for causal mechanisms driving adaptation of natural populations to stressful environments. However, the recent focus on noncoding variants demonstrated their ubiquitous role in adaptation. Here, we performed genome-wide regulatory variation profiles between two divergent yeast strains when facing nitrogen nutritional stress. The open chromatin availability of several regulatory regions changes in response to nitrogen scarcity. Importantly, we describe regulatory events that deviate between strains. Our results demonstrate a widespread variation in gene expression regulation between naturally occurring populations in response to stressful environments.

Genomics pipelines to investigate susceptibility in whole genome and exome sequenced data for variant discovery, annotation, prediction and genotyping

Over the last few decades, genomics is leading toward audacious future, and has been changing our views about conducting biomedical research, studying diseases, and understanding diversity in our society across the human species. The whole genome and exome sequencing (WGS/WES) are two of the most popular next-generation sequencing (NGS) methodologies that are currently being used to detect genetic variations of clinical significance. Investigating WGS/WES data for the variant discovery and genotyping is based on the nexus of different data analytic applications. Although several bioinformatics applications have been developed, and many of those are freely available and published. Timely finding and interpreting genetic variants are still challenging tasks among diagnostic laboratories and clinicians. In this study, we are interested in understanding, evaluating, and reporting the current state of solutions available to process the NGS data of variable lengths and types for the identification of variants, alleles, and haplotypes. Residing within the scope, we consulted high quality peer reviewed literature published in last 10 years. We were focused on the standalone and networked bioinformatics applications proposed to efficiently process WGS and WES data, and support downstream analysis for gene-variant discovery, annotation, prediction, and interpretation. We have discussed our findings in this manuscript, which include but not are limited to the set of operations, workflow, data handling, involved tools, technologies and algorithms and limitations of the assessed applications.

Measurement of histone replacement dynamics with genetically encoded exchange timers in yeast

Histone exchange between histones carrying position-specific marks and histones bearing general marks is important for gene regulation, but understanding of histone exchange remains incomplete. To overcome the poor time resolution of conventional pulse-chase histone labeling, we present a genetically encoded histone exchange timer sensitive to the duration that two tagged histone subunits co-reside at an individual genomic locus. We apply these sensors to map genome-wide patterns of histone exchange in yeast using single samples. Comparing H3 exchange in cycling and G1-arrested cells suggests that replication-independent H3 exchange occurs at several hundred nucleosomes (<1% of all nucleosomes) per minute, with a maximal rate at histone promoters. We observed substantial differences between the two nucleosome core subcomplexes: H2A-H2B subcomplexes undergo rapid transcription-dependent replacement within coding regions, whereas H3-H4 replacement occurs predominantly within promoter nucleosomes, in association with gene activation or repression. Our timers allow the in vivo study of histone exchange dynamics with minute time scale resolution.

FACT interacts with Set3 HDAC and fine-tunes GAL1 transcription in response to environmental stimulation

The histone chaperone facilitates chromatin transactions (FACT) functions in various DNA transactions. How FACT performs these multiple functions remains largely unknown. Here, we found, for the first time, that the N-terminal domain of its Spt16 subunit interacts with the Set3 histone deacetylase complex (Set3C) and that FACT and Set3C function in the same pathway to regulate gene expression in some settings. We observed that Spt16-G132D mutant proteins show defects in binding to Set3C but not other reported FACT interactors. At the permissive temperature, induction of the GAL1 and GAL10 genes is reduced in both spt16-G132D and set3Δ cells, whereas transient upregulation of GAL10 noncoding RNA (ncRNA), which is transcribed from the 3′ end of the GAL10 gene, is elevated. Mutations that inhibit GAL10 ncRNA transcription reverse the GAL1 and GAL10 induction defects in spt16-G132D and set3Δ mutant cells. Mechanistically, set3Δ and FACT (spt16-G132D) mutants show reduced histone acetylation and increased nucleosome occupancy at the GAL1 promoter under inducing conditions and inhibition of GAL10 ncRNA transcription also partially reverses these chromatin changes. These results indicate that FACT interacts with Set3C, which in turn prevents uncontrolled GAL10 ncRNA expression and fine-tunes the expression of GAL genes upon a change in carbon source.

Functional characterization of the cnidarian antiviral immune response reveals ancestral complexity

Animals evolved a broad repertoire of innate immune sensors and downstream effector cascades for defense against RNA viruses. Yet, this system varies greatly among different bilaterian animals, masking its ancestral state. In this study, we aimed to characterize the antiviral immune response of the cnidarian Nematostella vectensis and decipher the function of the retinoic acid-inducible gene I (RIG-I)-like receptors (RLRs) known to detect viral double-stranded RNA (dsRNA) in bilaterians but activate different antiviral pathways in vertebrates and nematodes. We show that polyinosinic:polycytidylic acid (poly(I:C)), a mimic of long viral dsRNA and a primary ligand for the vertebrate RLR melanoma differentiation-associated protein 5 (MDA5), triggers a complex antiviral immune response bearing features distinctive for both vertebrate and invertebrate systems. Importantly, a well-characterized agonist of the vertebrate RIG-I receptor does not induce a significant transcriptomic response that bears signature of the antiviral immune response, which experimentally supports the results of a phylogenetic analysis indicating clustering of the two N. vectensis RLR paralogs (NveRLRa and NveRLRb) with MDA5. Furthermore, the results of affinity assays reveal that NveRLRb binds poly(I:C) and long dsRNA and its knockdown impairs the expression of putative downstream effector genes including RNA interference components. Our study provides for the first time the functional evidence for the conserved role of RLRs in initiating immune response to dsRNA that originated before the cnidarian–bilaterian split and lay a strong foundation for future research on the evolution of the immune responses to RNA viruses.

Interplay among transacting factors around promoter in the initial phases of transcription

The initiation signals are raised around the promoter by one of the general transcription factors, triggering a sequence of events that lead to mRNA transcript formation from target genes. Both specific noncoding DNA regions and transacting, macromolecular assemblies are intricately involved and indispensable. The transition between the subsequent transcriptional stages is accompanied by stage-specific signals and structural changes in the macromolecular assemblies and facilitated by the recruitment/removal of other chromatin and transcription-associated elements. Here, we discuss the choreography of transacting factors around promoter in the establishment and effectuation of the initial phases of transcription such as NDR formation, +1 nucleosome positioning, promoter DNA opening, and RNAPII promoter escape from a structural viewpoint.

Sophisticated Conversations between Chromatin and Chromatin Remodelers, and Dissonances in Cancer

The establishment and maintenance of genome packaging into chromatin contribute to define specific cellular identity and function. Dynamic regulation of chromatin organization and nucleosome positioning are critical to all DNA transactions-in particular, the regulation of gene expression-and involve the cooperative action of sequence-specific DNA-binding factors, histone modifying enzymes, and remodelers. Remodelers are molecular machines that generate various chromatin landscapes, adjust nucleosome positioning, and alter DNA accessibility by using ATP binding and hydrolysis to perform DNA translocation, which is highly regulated through sophisticated structural and functional conversations with nucleosomes. In this review, I first present the functional and structural diversity of remodelers, while emphasizing the basic mechanism of DNA translocation, the common regulatory aspects, and the hand-in-hand progressive increase in complexity of the regulatory conversations between remodelers and nucleosomes that accompanies the increase in challenges of remodeling processes. Next, I examine how, through nucleosome positioning, remodelers guide the regulation of gene expression. Finally, I explore various aspects of how alterations/mutations in remodelers introduce dissonance into the conversations between remodelers and nucleosomes, modify chromatin organization, and contribute to oncogenesis.

Scaling of intrinsic noise in an autocratic reaction network

Biochemical reactions in living cells often produce stochastic trajectories due to the fluctuations of the finite number of the macromolecular species present inside the cell. A significant number of computational and theoretical studies have previously investigated stochasticity in small regulatory networks to understand its origin and regulation. At the systems level regulatory networks have been determined to be hierarchical resembling social networks. In order to determine the stochasticity in networks with hierarchical architecture, here we computationally investigated intrinsic noise in an autocratic reaction network in which only the upstream regulators regulate the downstream regulators. We studied the effects of the qualitative and quantitative nature of regulatory interactions on the stochasticity in the network. We established an unconventional scaling of noise with average abundance in which the noise passes through a minimum indicating that the network can be noisy both in the low and high abundance regimes. We determined that the bursty kinetics of the trajectories are responsible for such scaling. The scaling of noise remains intact for a mixed network that includes democratic subnetworks within the autocratic network.

Interplay between genetic, epigenetic, and gene expression variability: Considering complexity in evolvability

Genetic variability, epigenetic variability, and gene expression variability (noise) are generally considered independently in their relationship with phenotypic variation. However, they appear to be intrinsically interconnected and influence it in combination. The study of the interplay between genetic and epigenetic variability has the longest history. This article rather considers the introduction of gene expression variability in its relationships with the two others and reviews for the first time experimental evidences over the four relationships connected to gene expression noise. They show how introducing this third source of variability complicates the way of thinking evolvability and the emergence of biological novelty. Finally, cancer cells are proposed to be an ideal model to decipher the dynamic interplay between genetic, epigenetic, and gene expression variability when one of them is either experimentally increased or therapeutically targeted. This interplay is also discussed in an evolutionary perspective in the context of cancer cell drug resistance.

TATA and paused promoters active in differentiated tissues have distinct expression characteristics

Core promoter types differ in the extent to which RNA polymerase II (Pol II) pauses after initiation, but how this affects their tissue-specific gene expression characteristics is not well understood. While promoters with Pol II pausing elements are active throughout development, TATA promoters are highly active in differentiated tissues. We therefore used a genomics approach on late-stage Drosophila embryos to analyze the properties of promoter types. Using tissue-specific Pol II ChIP-seq, we found that paused promoters have high levels of paused Pol II throughout the embryo, even in tissues where the gene is not expressed, while TATA promoters only show Pol II occupancy when the gene is active. The promoter types are associated with different chromatin accessibility in ATAC-seq data and have different expression characteristics in single-cell RNA-seq data. The two promoter types may therefore be optimized for different properties: paused promoters show more consistent expression when active, while TATA promoters have lower background expression when inactive. We propose that tissue-specific genes have evolved to use two different strategies for their differential expression across tissues.

Distinct roles of nucleosome sliding and histone modifications in controlling the fidelity of transcription initiation

Regulation of gene expression starts from the transcription initiation. Regulated transcription initiation is critical for generating correct transcripts with proper abundance. The impact of epigenetic control, such as histone modifications and chromatin remodelling, on gene regulation has been extensively investigated, but their specific role in regulating transcription initiation is far from well understood. Here we aimed to better understand the roles of genes involved in histone H3 methylations and chromatin remodelling on the regulation of transcription initiation at a genome-scale using the budding yeast as a study system. We obtained and compared maps of transcription start site (TSS) at single-nucleotide resolution by nAnT-iCAGE for a strain with depletion of MINC (Mot1-Ino80C-Nc2) by Mot1p and Ino80p anchor-away (Mot1&Ino80AA) and a strain with loss of histone methylation (set1Δset2Δdot1Δ) to their wild-type controls. Our study showed that the depletion of MINC stimulated transcription initiation from many new sites flanking the dominant TSS of genes, while the loss of histone methylation generates more TSSs in the coding region. Moreover, the depletion of MINC led to less confined boundaries of TSS clusters (TCs) and resulted in broader core promoters, and such patterns are not present in the ssdΔ mutant. Our data also exhibits that the MINC has distinctive impacts on TATA-containing and TATA-less promoters. In conclusion, our study shows that MINC is required for accurate identification of bona fide TSSs, particularly in TATA-containing promoters, and histone methylation contributes to the repression of transcription initiation in coding regions.

High-Resolution Mapping of Transcription Initiation in the Asexual Stages of Toxoplasma gondii

Toxoplasma gondii is a common parasite of humans and animals, causing life-threatening disease in the immunocompromized, fetal abnormalities when contracted during gestation, and recurrent ocular lesions in some patients. Central to the prevalence and pathogenicity of this protozoan is its ability to adapt to a broad range of environments, and to differentiate between acute and chronic stages. These processes are underpinned by a major rewiring of gene expression, yet the mechanisms that regulate transcription in this parasite are only partially characterized. Deciphering these mechanisms requires a precise and comprehensive map of transcription start sites (TSSs); however, Toxoplasma TSSs have remained incompletely defined. To address this challenge, we used 5'-end RNA sequencing to genomically assess transcription initiation in both acute and chronic stages of Toxoplasma. Here, we report an in-depth analysis of transcription initiation at promoters, and provide empirically-defined TSSs for 7603 (91%) protein-coding genes, of which only 1840 concur with existing gene models. Comparing data from acute and chronic stages, we identified instances of stage-specific alternative TSSs that putatively generate mRNA isoforms with distinct 5' termini. Analysis of the nucleotide content and nucleosome occupancy around TSSs allowed us to examine the determinants of TSS choice, and outline features of Toxoplasma promoter architecture. We also found pervasive divergent transcription at Toxoplasma promoters, clustered within the nucleosomes of highly-symmetrical phased arrays, underscoring chromatin contributions to transcription initiation. Corroborating previous observations, we asserted that Toxoplasma 5' leaders are among the longest of any eukaryote studied thus far, displaying a median length of approximately 800 nucleotides. Further highlighting the utility of a precise TSS map, we pinpointed motifs associated with transcription initiation, including the binding sites of the master regulator of chronic-stage differentiation, BFD1, and a novel motif with a similar positional arrangement present at 44% of Toxoplasma promoters. This work provides a critical resource for functional genomics in Toxoplasma, and lays down a foundation to study the interactions between genomic sequences and the regulatory factors that control transcription in this parasite.

Of numbers and movement - understanding transcription factor pathogenesis by advanced microscopy

Transcription factors (TFs) are life-sustaining and, therefore, the subject of intensive research. By regulating gene expression, TFs control a plethora of developmental and physiological processes, and their abnormal function commonly leads to various developmental defects and diseases in humans. Normal TF function often depends on gene dosage, which can be altered by copy-number variation or loss-of-function mutations. This explains why TF haploinsufficiency (HI) can lead to disease. Since aberrant TF numbers frequently result in pathogenic abnormalities of gene expression, quantitative analyses of TFs are a priority in the field. In vitro single-molecule methodologies have significantly aided the identification of links between TF gene dosage and transcriptional outcomes. Additionally, advances in quantitative microscopy have contributed mechanistic insights into normal and aberrant TF function. However, to understand TF biology, TF-chromatin interactions must be characterised in vivo, in a tissue-specific manner and in the context of both normal and altered TF numbers. Here, we summarise the advanced microscopy methodologies most frequently used to link TF abundance to function and dissect the molecular mechanisms underlying TF HIs. Increased application of advanced single-molecule and super-resolution microscopy modalities will improve our understanding of how TF HIs drive disease.

Molecular and evolutionary processes generating variation in gene expression

Heritable variation in gene expression is common within and between species. This variation arises from mutations that alter the form or function of molecular gene regulatory networks that are then filtered by natural selection. High-throughput methods for introducing mutations and characterizing their cis- and trans-regulatory effects on gene expression (particularly, transcription) are revealing how different molecular mechanisms generate regulatory variation, and studies comparing these mutational effects with variation seen in the wild are teasing apart the role of neutral and non-neutral evolutionary processes. This integration of molecular and evolutionary biology allows us to understand how the variation in gene expression we see today came to be and to predict how it is most likely to evolve in the future.

Deep learning suggests that gene expression is encoded in all parts of a co-evolving interacting gene regulatory structure

Understanding the genetic regulatory code governing gene expression is an important challenge in molecular biology. However, how individual coding and non-coding regions of the gene regulatory structure interact and contribute to mRNA expression levels remains unclear. Here we apply deep learning on over 20,000 mRNA datasets to examine the genetic regulatory code controlling mRNA abundance in 7 model organisms ranging from bacteria to Human. In all organisms, we can predict mRNA abundance directly from DNA sequence, with up to 82% of the variation of transcript levels encoded in the gene regulatory structure. By searching for DNA regulatory motifs across the gene regulatory structure, we discover that motif interactions could explain the whole dynamic range of mRNA levels. Co-evolution across coding and non-coding regions suggests that it is not single motifs or regions, but the entire gene regulatory structure and specific combination of regulatory elements that define gene expression levels.

Does Gene Expression Noise Play a Functional Role in Plants?

Gene expression in individual cells can be surprisingly noisy. In unicellular organisms this noise can be functional; for example, by allowing a subfraction of the population to prepare for environmental stress. The role of gene expression noise in multicellular organisms has, however, remained unclear. In this review, we discuss how new techniques are revealing an unexpected level of variability in gene expression between and within genetically identical plants. We describe recent progress as well as speculate on the function of transcriptional noise as a mechanism for generating functional phenotypic diversity in plants.

Histone exchange is associated with activator function at transcribed promoters and with repression at histone loci

Transcription in eukaryotes correlates with major chromatin changes, including the replacement of old nucleosomal histones by new histones at the promoters of genes. The role of these histone exchange events in transcription remains unclear. In particular, the causal relationship between histone exchange and activator binding, preinitiation complex (PIC) assembly, and/or subsequent transcription remains unclear. Here, we provide evidence that histone exchange at gene promoters is not simply a consequence of PIC assembly or transcription but instead is mediated by activators. We further show that not all activators up-regulate gene expression by inducing histone turnover. Thus, histone exchange does not simply correlate with transcriptional activity, but instead reflects the mode of action of the activator. Last, we show that histone turnover is not only associated with activator function but also plays a role in transcriptional repression at the histone loci.

Independent evolution of transcript abundance and gene regulatory dynamics

Changes in gene expression drive novel phenotypes, raising interest in how gene expression evolves. In contrast to the static genome, cells modulate gene expression in response to changing environments. Previous comparative studies focused on specific conditions, describing interspecies variation in expression levels, but providing limited information about variation across different conditions. To close this gap, we profiled mRNA levels of two related yeast species in hundreds of conditions and used coexpression analysis to distinguish variation in the dynamic pattern of gene expression from variation in expression levels. The majority of genes whose expression varied between the species maintained a conserved dynamic pattern. Cases of diverged dynamic pattern correspond to genes that were induced under distinct subsets of conditions in the two species. Profiling the interspecific hybrid allowed us to distinguish between genes with predominantly cis- or trans-regulatory variation. We find that trans-varying alleles are dominantly inherited, and that cis-variations are often complemented by variations in trans Based on these results, we suggest that gene expression diverges primarily through changes in expression levels, but does not alter the pattern by which these levels are dynamically regulated.

Universal promoter scanning by Pol II during transcription initiation in Saccharomyces cerevisiae

BACKGROUND

The majority of eukaryotic promoters utilize multiple transcription start sites (TSSs). How multiple TSSs are specified at individual promoters across eukaryotes is not understood for most species. In Saccharomyces cerevisiae, a pre-initiation complex (PIC) comprised of Pol II and conserved general transcription factors (GTFs) assembles and opens DNA upstream of TSSs. Evidence from model promoters indicates that the PIC scans from upstream to downstream to identify TSSs. Prior results suggest that TSS distributions at promoters where scanning occurs shift in a polar fashion upon alteration in Pol II catalytic activity or GTF function.

RESULTS

To determine the extent of promoter scanning across promoter classes in S. cerevisiae, we perturb Pol II catalytic activity and GTF function and analyze their effects on TSS usage genome-wide. We find that alterations to Pol II, TFIIB, or TFIIF function widely alter the initiation landscape consistent with promoter scanning operating at all yeast promoters, regardless of promoter class. Promoter architecture, however, can determine the extent of promoter sensitivity to altered Pol II activity in ways that are predicted by a scanning model.

CONCLUSIONS

Our observations coupled with previous data validate key predictions of the scanning model for Pol II initiation in yeast, which we term the shooting gallery. In this model, Pol II catalytic activity and the rate and processivity of Pol II scanning together with promoter sequence determine the distribution of TSSs and their usage.

Yeast chromatin remodeling complexes and their roles in transcription

The nucleosome is a small unit of chromatin, which is dynamic in eukaryotes. Chromatin conformation and post-translational modifications affect nucleosome dynamics under certain conditions, playing an important role in the epigenetic regulation of transcription, replication and reprogramming. The Snf2 remodeling family is one of the crucial remodeling complexes that tightly regulate chromatin structure and affect nucleosome dynamics. This family alters nucleosome positioning, exchanges histone variants, and assembles and disassembles nucleosomes at certain locations. Moreover, the Snf2 family, in conjunction with other co-factors, regulates gene expression in Saccharomyces cerevisiae. Here we first review recent findings on the Snf2 family remodeling complexes and then use some examples to illustrate the cooperation between different members of Snf2 family, and the cooperation between Snf2 family and other co-factors in gene regulation especially during transcription initiation.

Dynamics of Chromatin and Transcription during Transient Depletion of the RSC Chromatin Remodeling Complex

Nucleosome organization has a key role in transcriptional regulation, yet the precise mechanisms establishing nucleosome locations and their effect on transcription are unclear. Here, we use an induced degradation system to screen all yeast ATP-dependent chromatin remodelers. We characterize how rapid clearance of the remodeler affects nucleosome locations. Specifically, depletion of Sth1, the catalytic subunit of the RSC (remodel the structure of chromatin) complex, leads to rapid fill-in of nucleosome-free regions at gene promoters. These changes are reversible upon reintroduction of Sth1 and do not depend on DNA replication. RSC-dependent nucleosome positioning is pivotal in maintaining promoters of lowly expressed genes free from nucleosomes. In contrast, we observe that upon acute stress, the RSC is not necessary for the transcriptional response. Moreover, RSC-dependent nucleosome positions are tightly related to usage of specific transcription start sites. Our results suggest organizational principles that determine nucleosome positions with and without RSC and how these interact with the transcriptional process.

•

Screen of all yeast ATP-dependent remodelers with a conditional degradation system

•

RSC depletion leads to rapid replication-independent NFR fill-in

•

Recovery of RSC fully reverses NFR fill-in and transcriptional changes

•

RSC-dependent nucleosome positioning directly affect transcription start site choice

Allele-specific single-cell RNA sequencing reveals different architectures of intrinsic and extrinsic gene expression noises

Gene expression noise refers to the variation of the expression level of a gene among isogenic cells in the same environment, and has two sources: extrinsic noise arising from the disparity of the cell state and intrinsic noise arising from the stochastic process of gene expression in the same cell state. Due to the low throughput of the existing method for measuring the two noise components, the architectures of intrinsic and extrinsic expression noises remain elusive. Using allele-specific single-cell RNA sequencing, we here estimate the two noise components of 3975 genes in mouse fibroblast cells. Our analyses verify predicted influences of several factors such as the TATA-box and microRNA targeting on intrinsic or extrinsic noises and reveal gene function-associated noise trends implicating the action of natural selection. These findings unravel differential regulations, optimizations, and biological consequences of intrinsic and extrinsic noises and can aid the construction of desired synthetic circuits.

YeasTSS: an integrative web database of yeast transcription start sites

The transcription initiation landscape of eukaryotic genes is complex and highly dynamic. In eukaryotes, genes can generate multiple transcript variants that differ in 5' boundaries due to usages of alternative transcription start sites (TSSs), and the abundance of transcript isoforms are highly variable. Due to a large number and complexity of the TSSs, it is not feasible to depict details of transcript initiation landscape of all genes using text-format genome annotation files. Therefore, it is necessary to provide data visualization of TSSs to represent quantitative TSS maps and the core promoters (CPs). In addition, the selection and activity of TSSs are influenced by various factors, such as transcription factors, chromatin remodeling and histone modifications. Thus, integration and visualization of functional genomic data related to these features could provide a better understanding of the gene promoter architecture and regulatory mechanism of transcription initiation. Yeast species play important roles for the research and human society, yet no database provides visualization and integration of functional genomic data in yeast. Here, we generated quantitative TSS maps for 12 important yeast species, inferred their CPs and built a public database, YeasTSS (www.yeastss.org). YeasTSS was designed as a central portal for visualization and integration of the TSS maps, CPs and functional genomic data related to transcription initiation in yeast. YeasTSS is expected to benefit the research community and public education for improving genome annotation, studies of promoter structure, regulated control of transcription initiation and inferring gene regulatory network.

Multivariable regulation of gene expression plasticity in metazoans

An important capacity of genes is the rapid change of expression levels to cope with the environment, known as expression responsiveness or plasticity. Elucidating the genomic mechanisms determining expression plasticity is critical for understanding the molecular basis of phenotypic plasticity, fitness and adaptation. In this study, we systematically quantified gene expression plasticity in four metazoan species by integrating changes of expression levels under a large number of genetic and environmental conditions. From this, we demonstrated that expression plasticity measures a distinct feature of gene expression that is orthogonal to other well-studied features, including gene expression level and tissue specificity/broadness. Expression plasticity is conserved across species with important physiological implications. The magnitude of expression plasticity is highly correlated with gene function and genes with high plasticity are implicated in disease susceptibility. Genome-wide analysis identified many conserved promoter cis-elements, trans-acting factors (such as CTCF), and gene body histone modifications (H3K36me3, H3K79me2 and H4K20me1) that are significantly associated with expression plasticity. Analysis of expression changes in perturbation experiments further validated a causal role of specific transcription factors and histone modifications. Collectively, this work reveals the general properties, physiological implications and multivariable regulation of gene expression plasticity in metazoans, extending the mechanistic understanding of gene regulation.

Epigenetic changes in fibroblasts drive cancer metabolism and differentiation

Genomic changes that drive cancer initiation and progression contribute to the co-evolution of the adjacent stroma. The nature of the stromal reprogramming involves differential DNA methylation patterns and levels that change in response to the tumor and systemic therapeutic intervention. Epigenetic reprogramming in carcinoma-associated fibroblasts are robust biomarkers for cancer progression and have a transcriptional impact that support cancer epithelial progression in a paracrine manner. For prostate cancer, promoter hypermethylation and silencing of the RasGAP, RASAL3 that resulted in the activation of Ras signaling in carcinoma-associated fibroblasts. Stromal Ras activity initiated a process of macropinocytosis that provided prostate cancer epithelia with abundant glutamine for metabolic conversion to fuel its proliferation and a signal to transdifferentiate into a neuroendocrine phenotype. This epigenetic oncogenic metabolic/signaling axis seemed to be further potentiated by androgen receptor signaling antagonists and contributed to therapeutic resistance. Intervention of stromal signaling may complement conventional therapies targeting the cancer cell.

Molecular Bases of DNA Packaging in Bacteria Revealed by All-Atom Molecular Dynamics Simulations: The Case of Histone-Like Proteins in Borrelia burgdorferi

DNA compaction is essential to ensure the packaging of the genetic material in living cells and also plays a key role in the epigenetic regulation of gene expression. In both humans and bacteria, DNA packaging is achieved by specific well-conserved proteins. Here, by means of all-atom molecular dynamics simulations, including the determination of relevant free-energy profiles, we rationalize the molecular bases for this remarkable process in bacteria, illustrating the crucial role played by positively charged amino acids of a small histone-like protein. We also present compelling evidence that this histone-like protein alone can induce strong bending of a DNA duplex around its core domain, a process that requires overcoming a major free-energy barrier.

Resolving noise-control conflict by gene duplication

Gene duplication promotes adaptive evolution in two main ways: allowing one duplicate to evolve a new function and splitting ancestral functions between the duplicates. The second scenario may resolve adaptive conflicts that can rise when one gene performs different functions. In an apparent departure from both scenarios, low-expressing transcription factor (TF) duplicates commonly bind to the same DNA motifs and act in overlapping conditions. To examine for possible benefits of this apparent redundancy, we examined the Msn2 and Msn4 duplicates in budding yeast. We show that Msn2,4 function as one unit by inducing the same set of target genes in overlapping conditions. Yet, the two-factor composition allows this unit's expression to be both environmentally responsive and with low noise, resolving an adaptive conflict that limits expression of single genes. We propose that duplication can provide adaptive benefit through cooperation rather than functional divergence, allowing two-factor dynamics with beneficial properties that cannot be achieved by a single gene.