Search

Identifying the molecular origins by which new morphological structures evolve is one of the long-standing problems in evolutionary biology. To date, vanishingly few examples provide a compelling account of how new morphologies were initially formed, thereby limiting our understanding of how diverse forms of life derived their complex features. Here, we provide evidence that the large projections on the Drosophila eugracilis phallus that are implicated in sexual conflict have evolved through the partial co-option of the trichome genetic network. These unicellular apical projections on the phallus postgonal sheath are reminiscent of trichomes that cover the Drosophila body but are up to 20-fold larger in size. During their development, they express the transcription factor Shavenbaby, the master regulator of the trichome network. Consistent with the co-option of the Shavenbaby network during the evolution of the D. eugracilis projections, somatic mosaic CRISPR-Cas9 mutagenesis shows that shavenbaby is necessary for their proper length. Moreover, misexpression of Shavenbaby in the sheath of D. melanogaster, a naive species that lacks these projections, is sufficient to induce small trichomes. These induced projections rely on a genetic network that is shared to a large extent with the D. eugracilis projections, indicating its partial co-option but also some genetic rewiring. Thus, by leveraging a genetically tractable evolutionary novelty, our work shows that the trichome-forming network is flexible enough that it can be partially co-opted in a new context and subsequently refined to produce unique apical projections that are barely recognizable compared with their simpler ancestral beginnings., Competing Interests: Declaration of interests The authors declare no competing interests., (Copyright © 2024 Elsevier Inc. All rights reserved.)

Identifying the molecular origins by which new morphological structures evolve is one of the long standing problems in evolutionary biology. To date, vanishingly few examples provide a compelling account of how new morphologies were initially formed, thereby limiting our understanding of how diverse forms of life derived their complex features. Here, we provide evidence that the large projections on the Drosophila eugracilis phallus that are implicated in sexual conflict have evolved through co-option of the trichome genetic network. These unicellular apical projections on the phallus postgonal sheath are reminiscent of trichomes that cover the Drosophila body but are up to 20-fold larger in size. During their development, they express the transcription factor Shavenbaby, the master regulator of the trichome network. Consistent with the co-option of the Shavenbaby network during the evolution of the D. eugracilis projections, somatic mosaic CRISPR/Cas9 mutagenesis shows that shavenbaby is necessary for their proper length. Moreover, mis-expression of Shavenbaby in the sheath of D. melanogaster , a naïve species that lacks these extensions, is sufficient to induce small trichomes. These induced extensions rely on a genetic network that is shared to a large extent with the D. eugracilis projections, indicating its co-option but also some genetic rewiring. Thus, by leveraging a genetically tractable evolutionarily novelty, our work shows that the trichome-forming network is flexible enough that it can be co-opted in a new context, and subsequently refined to produce unique apical projections that are barely recognizable compared to their simpler ancestral beginnings.

Congenital heart disease often arises from perturbations of transcription factors (TFs) that guide cardiac development. ISLET1 (ISL1) is a TF that influences early cardiac cell fate, as well as differentiation of other cell types including motor neuron progenitors (MNPs) and pancreatic islet cells. While lineage specificity of ISL1 function is likely achieved through combinatorial interactions, its essential cardiac interacting partners are unknown. By assaying ISL1 genomic occupancy in human induced pluripotent stem cell-derived cardiac progenitors (CPs) or MNPs and leveraging the deep learning approach BPNet, we identified motifs of other TFs that predicted ISL1 occupancy in each lineage, with NKX2.5 and GATA motifs being most closely associated to ISL1 in CPs. Experimentally, nearly two-thirds of ISL1-bound loci were co-occupied by NKX2.5 and/or GATA4. Removal of NKX2.5 from CPs led to widespread ISL1 redistribution, and overexpression of NKX2.5 in MNPs led to ISL1 occupancy of CP-specific loci. These results reveal how ISL1 guides lineage choices through a combinatorial code that dictates genomic occupancy and transcription., Competing Interests: Declaration of interests D.S. is a co-founder and member of the board of directors of Tenaya Therapeutics and has equity in Tenaya Therapeutics. K.N.I. is an employee and shareholder of Tenaya Therapeutics. N.J.K. has received research support from Vir Biotechnology, F. Hoffmann-La Roche, and Rezo Therapeutics. N.J.K. has financially compensated consulting agreements with Maze Therapeutics, Interline Therapeutics, Rezo Therapeutics, and GEn1E Lifesciences, Inc. He is on the Board of Directors of Rezo Therapeutics and is a shareholder in Tenaya Therapeutics, Maze Therapeutics, Rezo Therapeutics, and Interline Therapeutics., (Copyright © 2023 The Authors. Published by Elsevier Inc. All rights reserved.)

Chromatin accessibility is integral to the process by which transcription factors (TFs) read out cis-regulatory DNA sequences, but it is difficult to differentiate between TFs that drive accessibility and those that do not. Deep learning models that learn complex sequence rules provide an unprecedented opportunity to dissect this problem. Using zygotic genome activation in Drosophila as a model, we analyzed high-resolution TF binding and chromatin accessibility data with interpretable deep learning and performed genetic validation experiments. We identify a hierarchical relationship between the pioneer TF Zelda and the TFs involved in axis patterning. Zelda consistently pioneers chromatin accessibility proportional to motif affinity, whereas patterning TFs augment chromatin accessibility in sequence contexts where they mediate enhancer activation. We conclude that chromatin accessibility occurs in two tiers: one through pioneering, which makes enhancers accessible but not necessarily active, and the second when the correct combination of TFs leads to enhancer activation., Competing Interests: Declaration of interests J.Z. owns a patent on ChIP-nexus (no. 10287628). A.K. is on the scientific advisory board of PatchBio, SerImmune, AINovo, TensorBio and OpenTargets, was a consultant with Illumina, and owns shares in Illumina, Deep Genomics, Immunai, and Freenome Inc. All other authors declare no competing interests., (Copyright © 2023 The Authors. Published by Elsevier Inc. All rights reserved.)

Short tandem repeats (STRs) are enriched in eukaryotic cis -regulatory elements and alter gene expression, yet how they regulate transcription remains unknown. We found that STRs modulate transcription factor (TF)-DNA affinities and apparent on-rates by about 70-fold by directly binding TF DNA-binding domains, with energetic impacts exceeding many consensus motif mutations. STRs maximize the number of weakly preferred microstates near target sites, thereby increasing TF density, with impacts well predicted by statistical mechanics. Confirming that STRs also affect TF binding in cells, neural networks trained only on in vivo occupancies predicted effects identical to those observed in vitro. Approximately 90% of TFs preferentially bound STRs that need not resemble known motifs, providing a cis -regulatory mechanism to target TFs to genomic sites.

While the accessibility of enhancers is dynamically regulated during development, promoters tend to be constitutively accessible and poised for activation by paused Pol II. By studying Lola-I, a Drosophila zinc finger transcription factor, we show here that the promoter state can also be subject to developmental regulation independently of gene activation. Lola-I is ubiquitously expressed at the end of embryogenesis and causes its target promoters to become accessible and acquire paused Pol II throughout the embryo. This promoter transition is required but not sufficient for tissue-specific target gene activation. Lola-I mediates this function by depleting promoter nucleosomes, similar to the action of pioneer factors at enhancers. These results uncover a level of regulation for promoters that is normally found at enhancers and reveal a mechanism for the de novo establishment of paused Pol II at promoters., (© 2023. Springer Nature Limited.)

Transcription factors (TF) are proteins that bind DNA in a sequence-specific manner to regulate gene transcription. Despite their unique intrinsic sequence preferences, in vivo genomic occupancy profiles of TFs differ across cellular contexts. Hence, deciphering the sequence determinants of TF binding, both intrinsic and context-specific, is essential to understand gene regulation and the impact of regulatory, non-coding genetic variation. Biophysical models trained on in vitro TF binding assays can estimate intrinsic affinity landscapes and predict occupancy based on TF concentration and affinity. However, these models cannot adequately explain context-specific, in vivo binding profiles. Conversely, deep learning models, trained on in vivo TF binding assays, effectively predict and explain genomic occupancy profiles as a function of complex regulatory sequence syntax, albeit without a clear biophysical interpretation. To reconcile these complementary models of in vitro and in vivo TF binding, we developed Affinity Distillation (AD), a method that extracts thermodynamic affinities de-novo from deep learning models of TF chromatin immunoprecipitation (ChIP) experiments by marginalizing away the influence of genomic sequence context. Applied to neural networks modeling diverse classes of yeast and mammalian TFs, AD predicts energetic impacts of sequence variation within and surrounding motifs on TF binding as measured by diverse in vitro assays with superior dynamic range and accuracy compared to motif-based methods. Furthermore, AD can accurately discern affinities of TF paralogs. Our results highlight thermodynamic affinity as a key determinant of in vivo binding, suggest that deep learning models of in vivo binding implicitly learn high-resolution affinity landscapes, and show that these affinities can be successfully distilled using AD. This new biophysical interpretation of deep learning models enables high-throughput in silico experiments to explore the influence of sequence context and variation on both intrinsic affinity and in vivo occupancy., Competing Interests: Competing interests A.K. is on the scientific advisory board of PatchBio, SerImmune, AINovo, TensorBio and OpenTargets, was a consulting Fellow with Illumina and owns shares in DeepGenomics, Immuni and Freenome. All other authors have no competing interests to declare.

What new questions can we ask about transcriptional regulation given recent developments in large-scale approaches?, Competing Interests: Declaration of interests J.Z. is an investigator at the Stowers Institute for Medical Research and also a professor (affiliate track) at the Department of Pathology and Laboratory Medicine at the University of Kansas Medical Center. H.J.B. is a co-founder and shareholder of Metric Biotechnologies, Inc. Columbia University has filed a patent on a technology tangentially related, on which H.J.B. is one of the inventors, to the topic of the Voices piece that H.J.B. contributed., (Copyright © 2023. Published by Elsevier Inc.)

After a COVID-related hiatus, the fifth biennial symposium on Evolution and Core Processes in Gene Regulation met at the Stowers Institute in Kansas City, Missouri July 21 to 24, 2022. This symposium, sponsored by the American Society for Biochemistry and Molecular Biology (ASBMB), featured experts in gene regulation and evolutionary biology. Topic areas covered enhancer evolution, the cis-regulatory code, and regulatory variation, with an overall focus on bringing the power of deep learning (DL) to decipher DNA sequence information. DL is a machine learning method that uses neural networks to learn complex rules that make predictions about diverse types of data. When DL models are trained to predict genomic data from DNA sequence information, their high prediction accuracy allows the identification of impactful genetic variants within and across species. In addition, the learned sequence rules can be extracted from the model and provide important clues about the mechanistic underpinnings of the cis-regulatory code., Competing Interests: Conflict of interest The authors declare that they have no conflicts of interest with the contents of the article., (Copyright © 2022 The Authors. Published by Elsevier Inc. All rights reserved.)

The Spt/Ada-Gcn5 Acetyltransferase (SAGA) coactivator complex has multiple modules with different enzymatic and non-enzymatic functions. How each module contributes to gene expression is not well understood. During Drosophila oogenesis, the enzymatic functions are not equally required, which may indicate that different genes require different enzymatic functions. An analogy for this phenomenon is the handyman principle: while a handyman has many tools, which tool he uses depends on what requires maintenance. Here we analyzed the role of the non-enzymatic core module during Drosophila oogenesis, which interacts with TBP. We show that depletion of SAGA-specific core subunits blocked egg chamber development at earlier stages than depletion of enzymatic subunits. These results, as well as additional genetic analyses, point to an interaction with TBP and suggest a differential role of SAGA modules at different promoter types. However, SAGA subunits co-occupied all promoter types of active genes in ChIP-seq and ChIP-nexus experiments, and the complex was not specifically associated with distinct promoter types in the ovary. The high-resolution genomic binding profiles were congruent with SAGA recruitment by activators upstream of the start site, and retention on chromatin by interactions with modified histones downstream of the start site. Our data illustrate that a distinct genetic requirement for specific components may conceal the fact that the entire complex is physically present and suggests that the biological context defines which module functions are critical., Competing Interests: This work was supported by funding from the Stowers Institute for Medical Research and a grant from the National Institute of General Medical Sciences (R35GM118068) to Jerry L Workman. I have read the journal’s policy and the authors of this manuscript have the following competing interests: J. Zeitlinger owns a patent on ChIP-nexus. Author Susan M. Abmayr was unable to confirm their authorship contributions. On their behalf, the corresponding author has reported their contributions to the best of their knowledge.

The arrangement (syntax) of transcription factor (TF) binding motifs is an important part of the cis-regulatory code, yet remains elusive. We introduce a deep learning model, BPNet, that uses DNA sequence to predict base-resolution chromatin immunoprecipitation (ChIP)-nexus binding profiles of pluripotency TFs. We develop interpretation tools to learn predictive motif representations and identify soft syntax rules for cooperative TF binding interactions. Strikingly, Nanog preferentially binds with helical periodicity, and TFs often cooperate in a directional manner, which we validate using clustered regularly interspaced short palindromic repeat (CRISPR)-induced point mutations. Our model represents a powerful general approach to uncover the motifs and syntax of cis-regulatory sequences in genomics data.

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., (© 2021 The Authors. Published under the terms of the CC BY 4.0 license.)

Genomics data are now being generated at large quantities, of exquisite high resolution and from single cells. They offer a unique opportunity to develop powerful machine learning algorithms, including neural networks, to uncover the rules of the cis-regulatory code. However, current modeling assumptions are often not based on state-of-the-art knowledge of the cis-regulatory code from transcription, developmental genetics, imaging and structural studies. Here I aim to fill this gap by giving a brief historical overview of the field, describing common misconceptions and providing knowledge that might help to guide computational approaches. I will describe the principles and mechanisms involved in the combinatorial requirement of transcription factor binding motifs for enhancer activity, including the role of chromatin accessibility, repressors and low-affinity motifs in the cis-regulatory code. Deciphering the cis-regulatory code would unlock an enormous amount of regulatory information in the genome and would allow us to locate cis-regulatory genetic variants involved in development and disease., Competing Interests: Conflict of interest statement J.Z. owns a patent on ChIP-nexus (Patent No. 10287628).

Mouse embryonic stem cells (ESCs) cultured in defined medium resemble the pre-implantation epiblast in the ground state, with full developmental capacity including the germline. β-Catenin is required to maintain ground state pluripotency in mouse ESCs, but its exact role is controversial. Here, we reveal a Tcf3-independent role of β-catenin in restraining germline and somatic lineage differentiation genes. We show that β-catenin binds target genes with E2F6 and forms a complex with E2F6 and HMGA2 or E2F6 and HP1γ. Our data indicate that these complexes help β-catenin restrain and fine-tune germ cell and neural developmental potential. Overall, our data reveal a previously unappreciated role of β-catenin in preserving lineage differentiation integrity in ground state ESCs., (Copyright © 2020 The Authors. Published by Elsevier Inc. All rights reserved.)

During development, transcription factors and signaling molecules govern gene regulatory networks to direct the formation of unique morphologies. As changes in gene regulatory networks are often implicated in morphological evolution, mapping transcription factor landscapes is important, especially in tissues that undergo rapid evolutionary change. The terminalia (genital and anal structures) of Drosophila melanogaster and its close relatives exhibit dramatic changes in morphology between species. While previous studies have identified network components important for patterning the larval genital disc, the networks governing adult structures during pupal development have remained uncharted. Here, we performed RNA-seq in whole Drosophila melanogaster male terminalia followed by in situ hybridization for 100 highly expressed transcription factors during pupal development. We find that the male terminalia are highly patterned during pupal stages and that specific transcription factors mark separate structures and substructures. Our results are housed online in a searchable database (https://flyterminalia.pitt.edu/) as a resource for the community. This work lays a foundation for future investigations into the gene regulatory networks governing the development and evolution of Drosophila terminalia., (Copyright © 2019 Vincent et al.)

FACT (facilitates chromatin transcription) is an evolutionarily conserved histone chaperone that was initially identified as an activity capable of promoting RNA polymerase II (Pol II) transcription through nucleosomes in vitro. In this report, we describe a global analysis of FACT function in Pol II transcription in Drosophila. We present evidence that loss of FACT has a dramatic impact on Pol II elongation-coupled processes including histone H3 lysine 4 (H3K4) and H3K36 methylation, consistent with a role for FACT in coordinating histone modification and chromatin architecture during Pol II transcription. Importantly, we identify a role for FACT in the maintenance of promoter-proximal Pol II pausing, a key step in transcription activation in higher eukaryotes. These findings bring to light a broader role for FACT in the regulation of Pol II transcription., (Copyright © 2019 The Author(s). Published by Elsevier Inc. All rights reserved.)

RNA polymerase II (Pol II) pausing is a general regulatory step in transcription, yet the stability of paused Pol II varies widely between genes. Although paused Pol II stability correlates with core promoter elements, the contribution of individual sequences remains unclear, in part because no rapid assay is available for measuring the changes in Pol II pausing as a result of altered promoter sequences. Here, we overcome this hurdle by showing that ChIP-nexus captures the endogenous Pol II pausing on transfected plasmids. Using this reporter-ChIP-nexus assay in Drosophila cells, we show that the pausing stability is influenced by downstream promoter sequences, but that the strongest contribution to Pol II pausing comes from the initiator sequence, in which a single nucleotide, a G at the +2 position, is critical for stable Pol II pausing. These results establish reporter-ChIP-nexus as a valuable tool to analyze Pol II pausing., Competing Interests: WS, SA No competing interests declared, JZ J.Z. owns a patent on ChIP-nexus. EP3234199A1, (© 2019, Shao et al.)

The Drosophila genus is a unique group containing a wide range of species that occupy diverse ecosystems. In addition to the most widely studied species, Drosophila melanogaster , many other members in this genus also possess a well-developed set of genetic tools. Indeed, high-quality genomes exist for several species within the genus, facilitating studies of the function and evolution of cis -regulatory regions and proteins by allowing comparisons across at least 50 million years of evolution. Yet, the available genomes still fail to capture much of the substantial genetic diversity within the Drosophila genus. We have therefore tested protocols to rapidly and inexpensively sequence and assemble the genome from any Drosophila species using single-molecule sequencing technology from Oxford Nanopore. Here, we use this technology to present highly contiguous genome assemblies of 15 Drosophila species: 10 of the 12 originally sequenced Drosophila species ( ananassae , erecta , mojavensis , persimilis , pseudoobscura , sechellia , simulans , virilis , willistoni , and yakuba ), four additional species that had previously reported assemblies ( biarmipes , bipectinata , eugracilis , and mauritiana ), and one novel assembly ( triauraria ). Genomes were generated from an average of 29x depth-of-coverage data that after assembly resulted in an average contig N50 of 4.4 Mb. Subsequent alignment of contigs from the published reference genomes demonstrates that our assemblies could be used to close over 60% of the gaps present in the currently published reference genomes. Importantly, the materials and reagents cost for each genome was approximately $1,000 (USD). This study demonstrates the power and cost-effectiveness of long-read sequencing for genome assembly in Drosophila and provides a framework for the affordable sequencing and assembly of additional Drosophila genomes., (Copyright © 2018 Miller et al.)

In the initial published version of this article, there was an inadvertent omission from the Acknowledgements that this work was supported by Stowers Institute for Medical Research (SIMR-1004) and NIH National Cancer Institute grant to University of Kansas Cancer Center (P30 CA168524). This omission does not affect the description of the results or the conclusions of this work.

Transplantation of hematopoietic stem cells (HSCs) from human umbilical cord blood (hUCB) holds great promise for treating a broad spectrum of hematological disorders including cancer. However, the limited number of HSCs in a single hUCB unit restricts its widespread use. Although extensive efforts have led to multiple methods for ex vivo expansion of human HSCs by targeting single molecules or pathways, it remains unknown whether it is possible to simultaneously manipulate the large number of targets essential for stem cell self-renewal. Recent studies indicate that N 6 -methyladenosine (m 6 A) modulates the expression of a group of mRNAs critical for stem cell-fate determination by influencing their stability. Among several m 6 A readers, YTHDF2 is recognized as promoting targeted mRNA decay. However, the physiological functions of YTHDF2 in adult stem cells are unknown. Here we show that following the conditional knockout (KO) of mouse Ythdf2 the numbers of functional HSC were increased without skewing lineage differentiation or leading to hematopoietic malignancies. Furthermore, knockdown (KD) of human YTHDF2 led to more than a 10-fold increase in the ex vivo expansion of hUCB HSCs, a fivefold increase in colony-forming units (CFUs), and more than an eightfold increase in functional hUCB HSCs in the secondary serial of a limiting dilution transplantation assay. Mapping of m 6 A in RNAs from mouse hematopoietic stem and progenitor cells (HSPCs) as well as from hUCB HSCs revealed its enrichment in mRNAs encoding transcription factors critical for stem cell self-renewal. These m 6 A-marked mRNAs were recognized by Ythdf2 and underwent decay. In Ythdf2 KO HSPCs and YTHDF2 KD hUCB HSCs, these mRNAs were stabilized, facilitating HSC expansion. Knocking down one of YTHDF2's key targets, Tal1 mRNA, partially rescued the phenotype. Our study provides the first demonstration of the function of YTHDF2 in adult stem cell maintenance and identifies its important role in regulating HSC ex vivo expansion by regulating the stability of multiple mRNAs critical for HSC self-renewal, thus identifying potential for future clinical applications.

The HMG-box protein Capicua (Cic) is a conserved transcriptional repressor that functions downstream of receptor tyrosine kinase (RTK) signaling pathways in a relatively simple switch: In the absence of signaling, Cic represses RTK-responsive genes by binding to nearly invariant sites in DNA, whereas activation of RTK signaling down-regulates Cic activity, leading to derepression of its targets. This mechanism controls gene expression in both Drosophila and mammals, but whether Cic can also function via other regulatory mechanisms remains unknown. Here, we characterize an RTK-independent role of Cic in regulating spatially restricted expression of Toll/IL-1 signaling targets in Drosophila embryogenesis. We show that Cic represses those targets by binding to suboptimal DNA sites of lower affinity than its known consensus sites. This binding depends on Dorsal/NF-κB, which translocates into the nucleus upon Toll activation and binds next to the Cic sites. As a result, Cic binds to and represses Toll targets only in regions with nuclear Dorsal. These results reveal a mode of Cic regulation unrelated to the well-established RTK/Cic depression axis and implicate cooperative binding in conjunction with low-affinity binding sites as an important mechanism of enhancer regulation. Given that Cic plays a role in many developmental and pathological processes in mammals, our results raise the possibility that some of these Cic functions are independent of RTK regulation and may depend on cofactor-assisted DNA binding., Competing Interests: Conflict of interest statement: J.Z. owns a patent on ChIP-nexus.

Hoxa1 has important functional roles in neural crest specification, hindbrain patterning and heart and ear development, yet the enhancers and genes that are targeted by Hoxa1 are largely unknown. In this study, we performed a comprehensive analysis of Hoxa1 target genes using genome-wide Hoxa1 binding data in mouse ES cells differentiated with retinoic acid (RA) into neural fates in combination with differential gene expression analysis in Hoxa1 gain- and loss-of-function mouse and zebrafish embryos. Our analyses reveal that Hoxa1-bound regions show epigenetic marks of enhancers, occupancy of Hox cofactors and differential expression of nearby genes, suggesting that these regions are enriched for enhancers. In support of this, 80 of them mapped to regions with known reporter activity in transgenic mouse embryos based on the Vista enhancer database. Two additional enhancers in Dok5 and Wls1 were shown to mediate neural expression in developing mouse and zebrafish. Overall, our analysis of the putative target genes indicate that Hoxa1 has input to components of major signaling pathways, including Wnt, TGF-β, Hedgehog and Hippo, and frequently does so by targeting multiple components of a pathway such as secreted inhibitors, ligands, receptors and down-stream components. We also identified genes implicated in heart and ear development, neural crest migration and neuronal patterning and differentiation, which may underlie major Hoxa1 mutant phenotypes. Finally, we found evidence for a high degree of evolutionary conservation of many binding regions and downstream targets of Hoxa1 between mouse and zebrafish. Our genome-wide analyses in ES cells suggests that we have enriched for in vivo relevant target genes and pathways associated with functional roles of Hoxa1 in mouse development., (Copyright © 2017 Elsevier Inc. All rights reserved.)

Hyperactivating mutations in Ras signaling are hallmarks of carcinomas. Ras signaling mediates cell fate decisions as well as proliferation during development. It is not known what dictates whether Ras signaling drives differentiation versus proliferation. Here we show that the Hippo pathway is critical for this decision. Loss of Hippo switches Ras activation from promoting cellular differentiation to aggressive cellular proliferation. Transcriptome analysis combined with genetic tests show that this excessive proliferation depends on the synergistic induction of Ras target genes. Using ChIP-nexus, we find that Hippo signaling keeps Ras targets in check by directly regulating the expression of two key downstream transcription factors of Ras signaling: the ETS-domain transcription factor Pointed and the repressor Capicua. Our results highlight how independent signaling pathways can impinge on each other at the level of transcription factors, thereby providing a safety mechanism to keep proliferation in check under normal developmental conditions., (Copyright © 2017 Elsevier Inc. All rights reserved.)

Hoxa1 has diverse functional roles in differentiation and development. We identify and characterize properties of regions bound by HOXA1 on a genome-wide basis in differentiating mouse ES cells. HOXA1-bound regions are enriched for clusters of consensus binding motifs for HOX, PBX, and MEIS, and many display co-occupancy of PBX and MEIS. PBX and MEIS are members of the TALE family and genome-wide analysis of multiple TALE members (PBX, MEIS, TGIF, PREP1, and PREP2) shows that nearly all HOXA1 targets display occupancy of one or more TALE members. The combinatorial binding patterns of TALE proteins define distinct classes of HOXA1 targets, which may create functional diversity. Transgenic reporter assays in zebrafish confirm enhancer activities for many HOXA1-bound regions and the importance of HOX-PBX and TGIF motifs for their regulation. Proteomic analyses show that HOXA1 physically interacts on chromatin with PBX, MEIS, and PREP family members, but not with TGIF, suggesting that TGIF may have an independent input into HOXA1-bound regions. Therefore, TALE proteins appear to represent a wide repertoire of HOX cofactors, which may coregulate enhancers through distinct mechanisms. We also discover extensive auto- and cross-regulatory interactions among the Hoxa1 and TALE genes, indicating that the specificity of HOXA1 during development may be regulated though a complex cross-regulatory network of HOXA1 and TALE proteins. This study provides new insight into a regulatory network involving combinatorial interactions between HOXA1 and TALE proteins., (© 2017 De Kumar et al.; Published by Cold Spring Harbor Laboratory Press.)

RNA polymerase II (Pol II) pauses downstream of the transcription initiation site before beginning productive elongation. This pause is a key component of metazoan gene expression regulation. Some promoters have a strong disposition for Pol II pausing and often mediate faster, more synchronous changes in expression. This requires multiple rounds of transcription and thus cannot rely solely on pause release. However, it is unclear how pausing affects the initiation of new transcripts during consecutive rounds of transcription. Using our recently developed ChIP-nexus method, we find that Pol II pausing inhibits new initiation. We propose that paused Pol II helps prevent new initiation between transcription bursts, which may reduce noise.

Histone modifications are frequently used as markers for enhancer states, but how to interpret enhancer states in the context of embryonic development is not clear. The poised enhancer signature, involving H3K4me1 and low levels of H3K27ac, has been reported to mark inactive enhancers that are poised for future activation. However, future activation is not always observed, and alternative reasons for the widespread occurrence of this enhancer signature have not been investigated. By analyzing enhancers during dorsal-ventral (DV) axis formation in the Drosophila embryo, we find that the poised enhancer signature is specifically generated during patterning in the tissue where the enhancers are not induced, including at enhancers that are known to be repressed by a transcriptional repressor. These results suggest that, rather than serving exclusively as an intermediate step before future activation, the poised enhancer state may be a mark for spatial regulation during tissue patterning. We discuss the possibility that the poised enhancer state is more generally the result of repression by transcriptional repressors., (© 2017 Koenecke et al.; Published by Cold Spring Harbor Laboratory Press.)

Background: Drosophila dorso-ventral (DV) patterning is one of the best-understood regulatory networks to date, and illustrates the fundamental role of enhancers in controlling patterning, cell fate specification, and morphogenesis during development. Histone acetylation such as H3K27ac is an excellent marker for active enhancers, but it is challenging to obtain precise locations for enhancers as the highest levels of this modification flank the enhancer regions. How to best identify tissue-specific enhancers in a developmental system de novo with a minimal set of data is still unclear., Results: Using DV patterning as a test system, we develop a simple and effective method to identify tissue-specific enhancers de novo. We sample a broad set of candidate enhancer regions using data on CREB-binding protein co-factor binding or ATAC-seq chromatin accessibility, and then identify those regions with significant differences in histone acetylation between tissues. This method identifies hundreds of novel DV enhancers and outperforms ChIP-seq data of relevant transcription factors when benchmarked with mRNA expression data and transgenic reporter assays. These DV enhancers allow the de novo discovery of the relevant transcription factor motifs involved in DV patterning and contain additional motifs that are evolutionarily conserved and for which the corresponding transcription factors are expressed in a DV-biased fashion. Finally, we identify novel target genes of the regulatory network, implicating morphogenesis genes as early targets of DV patterning., Conclusions: Taken together, our approach has expanded our knowledge of the DV patterning network even further and is a general method to identify enhancers in any developmental system, including mammalian development.

The Drosophila genome activator Vielfaltig (Vfl), also known as Zelda (Zld), is thought to prime enhancers for activation by patterning transcription factors (TFs). Such priming is accompanied by increased chromatin accessibility, but the mechanisms by which this occurs are poorly understood. Here, we analyze the effect of Zld on genome-wide nucleosome occupancy and binding of the patterning TF Dorsal (Dl). Our results show that early enhancers are characterized by an intrinsically high nucleosome barrier. Zld tackles this nucleosome barrier through local depletion of nucleosomes with the effect being dependent on the number and position of Zld motifs. Without Zld, Dl binding decreases at enhancers and redistributes to open regions devoid of enhancer activity. We propose that Zld primes enhancers by lowering the high nucleosome barrier just enough to assist TFs in accessing their binding motifs and promoting spatially controlled enhancer activation if the right patterning TFs are present. We envision that genome activators in general will utilize this mechanism to activate the zygotic genome in a robust and precise manner., (© 2015 Sun et al.; Published by Cold Spring Harbor Laboratory Press.)

Understanding how eukaryotic enhancers are bound and regulated by specific combinations of transcription factors is still a major challenge. To better map transcription factor binding genome-wide at nucleotide resolution in vivo, we have developed a robust ChIP-exo protocol called ChIP-nexus (chromatin immunoprecipitation experiments with nucleotide resolution through exonuclease, unique barcode and single ligation), which utilizes an efficient DNA self-circularization step during library preparation. Application of ChIP-nexus to four proteins--human TBP and Drosophila NFkB, Twist and Max--shows that it outperforms existing ChIP protocols in resolution and specificity, pinpoints relevant binding sites within enhancers containing multiple binding motifs, and allows for the analysis of in vivo binding specificities. Notably, we show that Max frequently interacts with DNA sequences next to its motif, and that this binding pattern correlates with local DNA-sequence features such as DNA shape. ChIP-nexus will be broadly applicable to the study of in vivo transcription factor binding specificity and its relationship to cis-regulatory changes in humans and model organisms.

The TCT core promoter element is present in most ribosomal protein (RP) genes in Drosophila and humans. Here we show that TBP (TATA box-binding protein)-related factor TRF2, but not TBP, is required for transcription of the TCT-dependent RP genes. In cells, TCT-dependent transcription, but not TATA-dependent transcription, increases or decreases upon overexpression or depletion of TRF2. In vitro, purified TRF2 activates TCT but not TATA promoters. ChIP-seq (chromatin immunoprecipitation [ChIP] combined with deep sequencing) experiments revealed the preferential localization of TRF2 at TCT versus TATA promoters. Hence, a specialized TRF2-based RNA polymerase II system functions in the synthesis of RPs and complements the RNA polymerase I and III systems., (© 2014 Wang et al.; Published by Cold Spring Harbor Laboratory Press.)

Throughout Metazoa, developmental processes are controlled by a surprisingly limited number of conserved signaling pathways. Precisely how these signaling cassettes were assembled in early animal evolution remains poorly understood, as do the molecular transitions that potentiated the acquisition of their myriad developmental functions. Here we analyze the molecular evolution of the proto-oncogene yes-associated protein (Yap)/Yorkie, a key effector of the Hippo signaling pathway that controls organ size in both Drosophila and mammals. Based on heterologous functional analysis of evolutionarily distant Yap/Yorkie orthologs, we demonstrate that a structurally distinct interaction interface between Yap/Yorkie and its partner TEAD/Scalloped became fixed in the eumetazoan common ancestor. We then combine transcriptional profiling of tissues expressing phylogenetically diverse forms of Yap/Yorkie with ChIP-seq validation to identify a common downstream gene expression program underlying the control of tissue growth in Drosophila. Intriguingly, a subset of the newly identified Yorkie target genes are also induced by Yap in mammalian tissues, thus revealing a conserved Yap-dependent gene expression signature likely to mediate organ size control throughout bilaterian animals. Combined, these experiments provide new mechanistic insights while revealing the ancient evolutionary history of Hippo signaling., (© The Author 2014. Published by Oxford University Press on behalf of the Society for Molecular Biology and Evolution.)

The rapid expansion of genomics methods has enabled developmental biologists to address fundamental questions of developmental gene regulation on a genome-wide scale. These efforts have demonstrated that transcription of developmental control genes by RNA polymerase II (Pol II) is commonly regulated at the transition to productive elongation, resulting in the promoter-proximal accumulation of transcriptionally engaged but paused Pol II prior to gene induction. Here we review the mechanisms and possible functions of Pol II pausing and their implications for development.

Motivation: Chromatin immunoprecipitation coupled to next-generation sequencing (ChIP-seq) is widely used to study the in vivo binding sites of transcription factors (TFs) and their regulatory targets. Recent improvements to ChIP-seq, such as increased resolution, promise deeper insights into transcriptional regulation, yet require novel computational tools to fully leverage their advantages., Results: To this aim, we have developed peakzilla, which can identify closely spaced TF binding sites at high resolution (i.e. resolves individual binding sites even if spaced closely), as we demonstrate using semisynthetic datasets, performing ChIP-seq for the TF Twist in Drosophila embryos with different experimental fragment sizes, and analyzing ChIP-exo datasets. We show that the increased resolution reached by peakzilla is highly relevant, as closely spaced Twist binding sites are strongly enriched in transcriptional enhancers, suggesting a signature to discriminate functional from abundant non-functional or neutral TF binding. Peakzilla is easy to use, as it estimates all the necessary parameters from the data and is freely available., Availability and Implementation: The peakzilla program is available from https://github.com/steinmann/peakzilla or http://www.starklab.org/data/peakzilla/., Contact: stark@starklab.org., Supplementary Information: Supplementary data are available at Bioinformatics online.

Massive zygotic transcription begins in many organisms during the midblastula transition when the cell cycle of the dividing egg slows down. A few genes are transcribed before this stage but how this differential activation is accomplished is still an open question. We have performed ChIP-seq experiments on tightly staged Drosophila embryos and show that massive recruitment of RNA polymerase II (Pol II) with widespread pausing occurs de novo during the midblastula transition. However, ∼100 genes are strongly occupied by Pol II before this timepoint and most of them do not show Pol II pausing, consistent with a requirement for rapid transcription during the fast nuclear cycles. This global change in Pol II pausing correlates with distinct core promoter elements and associates a TATA-enriched promoter with the rapid early transcription. This suggests that promoters are differentially used during the zygotic genome activation, presumably because they have distinct dynamic properties. DOI:http://dx.doi.org/10.7554/eLife.00861.001.

Paused RNA polymerase (Pol II) is a pervasive feature of Drosophila embryos and mammalian stem cells, but its role in development is uncertain. Here, we demonstrate that a spectrum of paused Pol II determines the "time to synchrony"-the time required to achieve coordinated gene expression across the cells of a tissue. To determine whether synchronous patterns of gene activation are significant in development, we manipulated the timing of snail expression, which controls the coordinated invagination of ∼1,000 mesoderm cells during gastrulation. Replacement of the strongly paused snail promoter with moderately paused or nonpaused promoters causes stochastic activation of snail expression and increased variability of mesoderm invagination. Computational modeling of the dorsal-ventral patterning network recapitulates these variable and bistable gastrulation profiles and emphasizes the importance of timing of gene activation in development. We conclude that paused Pol II and transcriptional synchrony are essential for coordinating cell behavior during morphogenesis., (Copyright © 2013 Elsevier Inc. All rights reserved.)

Poised RNA polymerase II (Pol II) is predominantly found at developmental control genes and is thought to allow their rapid and synchronous induction in response to extracellular signals. How the recruitment of poised RNA Pol II is regulated during development is not known. By isolating muscle tissue from Drosophila embryos at five stages of differentiation, we show that the recruitment of poised Pol II occurs at many genes de novo and this makes them permissive for future gene expression. A comparison with other tissues shows that these changes are stage specific and not tissue specific. In contrast, Polycomb group repression is tissue specific, and in combination with Pol II (the balanced state) marks genes with highly dynamic expression. This suggests that poised Pol II is temporally regulated and is held in check in a tissue-specific fashion. We compare our data with findings in mammalian embryonic stem cells and discuss a framework for predicting developmental programs on the basis of the chromatin state., (Copyright © 2012 The Authors. Published by Elsevier Inc. All rights reserved.)

Chromatin immunoprecipitation (ChIP) followed by deep sequencing can now easily be performed across different conditions, time points and even species. However, analyzing such data is not trivial and standard methods are as yet unavailable. Here we present a protocol to systematically compare ChIP-sequencing (ChIP-seq) data across conditions. We first describe technical guidelines for data preprocessing, read mapping, read-density visualization and peak calling. We then describe methods and provide code with specific examples to compare different data sets across species and across conditions, including a threshold-free approach to measure global similarity, a strategy to assess the binary conservation of binding events and measurements for quantitative changes of binding. We discuss how differences in binding can be related to gene functions, gene expression and sequence changes. Once established, this protocol should take about 2 d to complete and be generally applicable to many data sets.

The binding of some transcription factors has been shown to diverge substantially between closely related species. Here we show that the binding of the developmental transcription factor Twist is highly conserved across six Drosophila species, revealing strong functional constraints at its enhancers. Conserved binding correlates with sequence motifs for Twist and its partners, permitting the de novo discovery of their combinatorial binding. It also includes over 10,000 low-occupancy sites near the detection limit, which tend to mark enhancers of later developmental stages. These results suggest that developmental enhancers can be highly evolutionarily constrained, presumably because of their complex combinatorial nature.

Genetic experiments over the last few decades have identified many developmental control genes critical for pattern formation and cell fate specification during the development of multicellular organisms. A large fraction of these genes encode transcription factors and signaling molecules, show highly dynamic expression patterns during development, and are deeply evolutionarily conserved and deregulated in various human diseases such as cancer. Because of their importance in development, evolution, and disease, a fundamental question in biology is how these developmental control genes are regulated in such an extensive and precise fashion. Using genomics methods, it has become clear that developmental control genes are a distinct group of genes with special regulatory characteristics. However, a systematic analysis of these characteristics has not been presented. Here we review how developmental control genes were discovered, evaluate their genome-wide regulation and gene structure, discuss emerging evidence for their mode of regulation, and estimate their overall abundance in the genome. Understanding the global regulation of developmental control genes may provide a new perspective on development in the era genomics., (Copyright 2010 Elsevier Inc. All rights reserved.)