Evolution of transcription factor binding through sequence variations and turnover of binding sites

The amniote-conserved DNA-binding domain of CGGBP1 restricts cytosine methylation of transcription factor binding sites in proximal promoters to regulate gene expression

CGGBP1 is a GC-rich DNA-binding protein which is important for genomic integrity, gene expression and epigenome maintenance through regulation of CTCF occupancy and cytosine methylation. It has remained unclear how CGGBP1 integrates multiple diverse functions with its simple architecture of only a DNA-binding domain tethered to a C-terminal tail with low structural rigidity. We have used truncated forms of CGGBP1 with or without the DNA-binding domain (DBD) to assay cytosine methylation and global gene expression. Proximal promoters of CGGBP1-repressed genes, although significantly GC-poor, contain GC-rich transcription factor binding motifs and exhibit base compositions indicative of low C-T transition rates due to prevention of cytosine methylation. Genome-wide analyses of cytosine methylation and binding of CGGBP1 DBD show that CGGBP1 restricts cytosine methylation in a manner that depends on its DBD and its DNA-binding. The CGGBP1-repressed genes show an increase in promoter cytosine methylation alongside a decrease in transcript abundance when the DBD-deficient CGGBP1 is expressed. Our findings suggest that CGGBP1 protects transcription factor binding sites (TFBS) from cytosine methylation-associated loss and thereby regulates gene expression. By analysing orthologous promoter sequences we show that restriction of cytosine methylation is a function of CGGBP1 progressively acquired during vertebrate evolution. A superimposition of our results and evolution of CGGBP1 suggests that mitigation of cytosine methylation is majorly achieved by its N-terminal DBD. Our results position CGGBP1 DNA-binding as a major evolutionarily acquired mechanism through which it keeps cytosine methylation under check and regulates TFBS retention and gene activity.

Transcription factor binding specificities of the oomycete Phytophthora infestans reflect conserved and divergent evolutionary patterns and predict function

BACKGROUND

Identifying the DNA-binding specificities of transcription factors (TF) is central to understanding gene networks that regulate growth and development. Such knowledge is lacking in oomycetes, a microbial eukaryotic lineage within the stramenopile group. Oomycetes include many important plant and animal pathogens such as the potato and tomato blight agent Phytophthora infestans, which is a tractable model for studying life-stage differentiation within the group.

RESULTS

Mining of the P. infestans genome identified 197 genes encoding proteins belonging to 22 TF families. Their chromosomal distribution was consistent with family expansions through unequal crossing-over, which were likely ancient since each family had similar sizes in most oomycetes. Most TFs exhibited dynamic changes in RNA levels through the P. infestans life cycle. The DNA-binding preferences of 123 proteins were assayed using protein-binding oligonucleotide microarrays, which succeeded with 73 proteins from 14 families. Binding sites predicted for representatives of the families were validated by electrophoretic mobility shift or chromatin immunoprecipitation assays. Consistent with the substantial evolutionary distance of oomycetes from traditional model organisms, only a subset of the DNA-binding preferences resembled those of human or plant orthologs. Phylogenetic analyses of the TF families within P. infestans often discriminated clades with canonical and novel DNA targets. Paralogs with similar binding preferences frequently had distinct patterns of expression suggestive of functional divergence. TFs were predicted to either drive life stage-specific expression or serve as general activators based on the representation of their binding sites within total or developmentally-regulated promoters. This projection was confirmed for one TF using synthetic and mutated promoters fused to reporter genes in vivo.

CONCLUSIONS

We established a large dataset of binding specificities for P. infestans TFs, representing the first in the stramenopile group. This resource provides a basis for understanding transcriptional regulation by linking TFs with their targets, which should help delineate the molecular components of processes such as sporulation and host infection. Our work also yielded insight into TF evolution during the eukaryotic radiation, revealing both functional conservation as well as diversification across kingdoms.

Dynamic chromatin accessibility landscapes of osteoblast differentiation and mineralization

Bone acts as a self-healing organ, which undergoes continuous regeneration process that is tightly regulated by the cooperation of osteoclasts with the capability of bone resorption and osteoblasts with the capability of bone formation. Generally, bone marrow derived mesenchymal stem cells (BMSCs) differentiated to final osteoblasts have been considered as critical role in bone remodeling. In this regard, several transcription factors (TFs) whose binding sites are initially hidden deep within accessible chromatin that participate in modulating osteoblast differentiation and bone matrix mineralization. Then, it is necessary to explore further the dynamic changes about the epigenetic transcription machinery during osteoblastogenesis. Here, we performed the chromatin accessibility and transcriptomic landscape of osteoblast differentiation and mineralization by using transposase-accessible chromatin sequencing (ATAC-seq) and RNA sequencing (RNA-Seq). Our data found that global chromatin accessibility during osteoblastogenesis was extensively improved. Above this, it is shown that key target genes including Col6a3, Serpina3n, Ms4a4d, Lyz2, Phf11b and Grin3a were enriched in differential loci RNA-seq and ATAC-Seq peaks with continuous changed tendency during osteoblasts differentiation and mineralization. In addition, Analysis of Motif Enrichment (AME) was used to elucidate TFs which modulated these target genes. In this study, it was shown for the first time that these important TFs including MEF2A, PRRX1, Shox2 and HOXB13 could alter promoter accessibility of target genes during osteoblastogenesis. This helps us understand how TF binding motif accessibility influences osteoblast differentiation. In addition, it also suggests that modulating the chromatin accessibility of osteogenesis could be developed as the promising strategies to regulate bone regeneration.

Cell layer-specific expression of the homeotic MADS-box transcription factor PhDEF contributes to modular petal morphogenesis in petunia

Floral homeotic MADS-box transcription factors ensure the correct morphogenesis of floral organs, which are organized in different cell layers deriving from distinct meristematic layers. How cells from these distinct layers acquire their respective identities and coordinate their growth to ensure normal floral organ morphogenesis is unresolved. Here, we studied petunia (Petunia × hybrida) petals that form a limb and tube through congenital fusion. We identified petunia mutants (periclinal chimeras) expressing the B-class MADS-box gene DEFICIENS in the petal epidermis or in the petal mesophyll, called wico and star, respectively. Strikingly, wico flowers form a strongly reduced tube while their limbs are almost normal, while star flowers form a normal tube but greatly reduced and unpigmented limbs, showing that petunia petal morphogenesis is highly modular. These mutants highlight the layer-specific roles of PhDEF during petal development. We explored the link between PhDEF and petal pigmentation, a well-characterized limb epidermal trait. The anthocyanin biosynthesis pathway was strongly downregulated in star petals, including its major regulator ANTHOCYANIN2 (AN2). We established that PhDEF directly binds to the AN2 terminator in vitro and in vivo, suggesting that PhDEF might regulate AN2 expression and therefore petal epidermis pigmentation. Altogether, we show that cell layer-specific homeotic activity in petunia petals differently impacts tube and limb development, revealing the relative importance of the different cell layers in the modular architecture of petunia petals.

Position-specific evolution in transcription factor binding sites, and a fast likelihood calculation for the F81 model

Transcription factor binding sites (TFBS), like other DNA sequence, evolve via mutation and selection relating to their function. Models of nucleotide evolution describe DNA evolution via single-nucleotide mutation. A stationary vector of such a model is the long-term distribution of nucleotides, unchanging under the model. Neutrally evolving sites may have uniform stationary vectors, but one expects that sites within a TFBS instead have stationary vectors reflective of the fitness of various nucleotides at those positions. We introduce 'position-specific stationary vectors' (PSSVs), the collection of stationary vectors at each site in a TFBS locus, analogous to the position weight matrix (PWM) commonly used to describe TFBS. We infer PSSVs for human TFs using two evolutionary models (Felsenstein 1981 and Hasegawa-Kishino-Yano 1985). We find that PSSVs reflect the nucleotide distribution from PWMs, but with reduced specificity. We infer ancestral nucleotide distributions at individual positions and calculate 'conditional PSSVs' conditioned on specific choices of majority ancestral nucleotide. We find that certain ancestral nucleotides exert a strong evolutionary pressure on neighbouring sequence while others have a negligible effect. Finally, we present a fast likelihood calculation for the F81 model on moderate-sized trees that makes this approach feasible for large-scale studies along these lines.

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.

The molecular grammar of protein disorder guiding genome-binding locations

Intrinsically disordered regions (IDRs) direct transcription factors (TFs) towards selected genomic occurrences of their binding motif, as exemplified by budding yeast's Msn2. However, the sequence basis of IDR-directed TF binding selectivity remains unknown. To reveal this sequence grammar, we analyze the genomic localizations of >100 designed IDR mutants, each carrying up to 122 mutations within this 567-AA region. Our data points at multivalent interactions, carried by hydrophobic-mostly aliphatic-residues dispersed within a disordered environment and independent of linear sequence motifs, as the key determinants of Msn2 genomic localization. The implications of our results for the mechanistic basis of IDR-based TF binding preferences are discussed.

Transcriptomic signatures reveal a shift towards an anti-inflammatory gene expression profile but also the induction of type I and type II interferon signaling networks through aryl hydrocarbon receptor activation in murine macrophages

Introduction

The aryl hydrocarbon receptor (AhR) is a ligand-activated transcription factor that regulates a broad range of target genes involved in the xenobiotic response, cell cycle control and circadian rhythm. AhR is constitutively expressed in macrophages (Mϕ), acting as key regulator of cytokine production. While proinflammatory cytokines, i.e., IL-1β, IL-6, IL-12, are suppressed through AhR activation, anti-inflammatory IL-10 is induced. However, the underlying mechanisms of those effects and the importance of the specific ligand structure are not yet completely understood.

Methods

Therefore, we have compared the global gene expression pattern in activated murine bone marrow-derived macrophages (BMMs) subsequently to exposure with either benzo[a]pyrene (BaP) or indole-3-carbinol (I3C), representing high-affinity vs. low-affinity AhR ligands, respectively, by means of mRNA sequencing. AhR dependency of observed effects was proved using BMMs from AhR-knockout (Ahr^-/-) mice.

Results and discussion

In total, more than 1,000 differentially expressed genes (DEGs) could be mapped, covering a plethora of AhR-modulated effects on basal cellular processes, i.e., transcription and translation, but also immune functions, i.e., antigen presentation, cytokine production, and phagocytosis. Among DEGs were genes that are already known to be regulated by AhR, i.e., Irf1, Ido2, and Cd84. However, we identified DEGs not yet described to be AhR-regulated in Mϕ so far, i.e., Slpi, Il12rb1, and Il21r. All six genes likely contribute to shifting the Mϕ phenotype from proinflammatory to anti-inflammatory. The majority of DEGs induced through BaP were not affected through I3C exposure, probably due to higher AhR affinity of BaP in comparison to I3C. Mapping of known aryl hydrocarbon response element (AHRE) sequence motifs in identified DEGs revealed more than 200 genes not possessing any AHRE, and therefore being not eligible for canonical regulation. Bioinformatic approaches modeled a central role of type I and type II interferons in the regulation of those genes. Additionally, RT-qPCR and ELISA confirmed a AhR-dependent expressional induction and AhR-dependent secretion of IFN-γ in response to BaP exposure, suggesting an auto- or paracrine activation pathway of Mϕ.

Contributions of mutation and selection to regulatory variation: lessons from the Saccharomyces cerevisiae TDH3 gene

Heritable variation in gene expression is common within and among species and contributes to phenotypic diversity. Mutations affecting either cis- or trans-regulatory sequences controlling gene expression give rise to variation in gene expression, and natural selection acting on this variation causes some regulatory variants to persist in a population for longer than others. To understand how mutation and selection interact to produce the patterns of regulatory variation we see within and among species, my colleagues and I have been systematically determining the effects of new mutations on expression of the TDH3 gene in Saccharomyces cerevisiae and comparing them to the effects of polymorphisms segregating within this species. We have also investigated the molecular mechanisms by which regulatory variants act. Over the past decade, this work has revealed properties of cis- and trans-regulatory mutations including their relative frequency, effects, dominance, pleiotropy and fitness consequences. Comparing these mutational effects to the effects of polymorphisms in natural populations, we have inferred selection acting on expression level, expression noise and phenotypic plasticity. Here, I summarize this body of work and synthesize its findings to make inferences not readily discernible from the individual studies alone. This article is part of the theme issue 'Interdisciplinary approaches to predicting evolutionary biology'.

The Memory Problem for Neutral Mutational Models of Evolution

Models for the evolution of various phenotypes are sometimes constructed with an assumption that mutational effects will be symmetrically distributed about a static mean. This model produces a memory effect that over long evolutionary times results in an expectation that randomized sequences underlying the genetic architecture of the trait will on average retain the ancestral phenotype. This expectation is unrealistic and also inconsistent with our current understanding of processes of molecular evolution.

Mechanisms of regulatory evolution in yeast

Studies of regulatory variation in yeast - at the level of new mutations, polymorphisms within a species, and divergence between species - have provided great insight into the molecular and evolutionary processes responsible for the evolution of gene expression in eukaryotes. The increasing ease with which yeast genomes can be manipulated and expression quantified in a high-throughput manner has recently accelerated mechanistic studies of cis- and trans-regulatory variation at multiple evolutionary timescales. These studies have, for example, identified differences in the properties of cis- and trans-acting mutations that affect their evolutionary fate, experimentally characterized the molecular mechanisms through which cis- and trans-regulatory variants act, and illustrated how regulatory networks can diverge between species with or without changes in gene expression.