Divergent expression of paralogous genes by modification of shared enhancer activity through a promoter-proximal silencer.

The duplication of genes and their associated cis-regulatory elements, or enhancers, is a key contributor to genome evolution and biological complexity. Moreover, many paralogs, particularly tandem duplicates, are fixed for long periods of time under the control of shared enhancers. However, in most cases, the mechanism by which gene expression and function diverge following duplication is not known. Here, we dissect the regulation and function of the paralogous nubbin/pdm2 genes during wing development in Drosophila melanogaster. We show that these paralogs play a redundant role in the wing and that their expression relies on a single shared wing enhancer. However, the two genes differ in their ability to respond to this enhancer, with nub responding in all wing progenitor cells and pdm2 only in a small subset. This divergence is a result of a pdm2-specific silencer element at the pdm2 promoter that receives repressive input from the transcription factor Rotund. Repression through this silencer also depends on nub, allowing pdm2 to fully respond to the wing enhancer when nub expression is perturbed and functional compensation to occur. Thus, expression divergence downstream of a shared enhancer arises as a consequence of silencing the promoter of one paralog.

Transcription factor paralogs orchestrate alternative gene regulatory networks by context-dependent cooperation with multiple cofactors.

In eukaryotes, members of transcription factor families often exhibit similar DNA binding properties in vitro, yet orchestrate paralog-specific gene regulatory networks in vivo. The serially homologous first (T1) and third (T3) thoracic legs of Drosophila, which are specified by the Hox proteins Scr and Ubx, respectively, offer a unique opportunity to address this paradox in vivo. Genome-wide analyses using epitope-tagged alleles of both Hox loci in the T1 and T3 leg imaginal discs, the precursors to the adult legs and ventral body regions, show that ~8% of Hox binding is paralog-specific. Binding specificity is mediated by interactions with distinct cofactors in different domains: the Hox cofactor Exd acts in the proximal domain and is necessary for Scr to bind many of its paralog-specific targets, while in the distal leg domain, the homeodomain protein Distal-less (Dll) enhances Scr binding to a different subset of loci. These findings reveal how Hox paralogs, and perhaps paralogs of other transcription factor families, orchestrate alternative downstream gene regulatory networks with the help of multiple, context-specific cofactors.

Cell-type-specific Hox regulatory strategies orchestrate tissue identity.

Hox proteins are homeodomain transcription factors that diversify serially homologous segments along the animal body axis, as revealed by the classic bithorax phenotype of Drosophila melanogaster, in which mutations in Ultrabithorax (Ubx) transform the third thoracic segment into the likeness of the second thoracic segment. To specify segment identity, we show that Ubx both increases and decreases chromatin accessibility, coinciding with its dual role as both an activator and repressor of transcription. However, the choice of transcriptional activity executed by Ubx is spatially regulated and depends on the availability of cofactors, with Ubx acting as a repressor in some populations and as an activator in others. Ubx-mediated changes to chromatin accessibility positively and negatively affect the binding of Scalloped (Sd), a transcription factor that is required for appendage development in both segments. These findings illustrate how a single Hox protein can modify complex gene regulatory networks to transform the identity of an entire tissue.

Control of tissue morphogenesis by the HOX gene Ultrabithorax.

Mutations in the Ultrabithorax (Ubx) gene cause homeotic transformation of the normally two-winged Drosophila into a four-winged mutant fly. Ubx encodes a HOX family transcription factor that specifies segment identity, including transformation of the second set of wings into rudimentary halteres. Ubx is known to control the expression of many genes that regulate tissue growth and patterning, but how it regulates tissue morphogenesis to reshape the wing into a haltere is still unclear. Here, we show that Ubx acts by repressing the expression of two genes in the haltere, Stubble and Notopleural, both of which encode transmembrane proteases that remodel the apical extracellular matrix to promote wing morphogenesis. In addition, Ubx induces expression of the Tissue inhibitor of metalloproteases in the haltere, which prevents the basal extracellular matrix remodelling necessary for wing morphogenesis. Our results provide a long-awaited explanation for how Ubx controls morphogenetic transformation.

Context-Dependent Gene Regulation by Homeodomain Transcription Factor Complexes Revealed by Shape-Readout Deficient Proteins.

Eukaryotic transcription factors (TFs) form complexes with various partner proteins to recognize their genomic target sites. Yet, how the DNA sequence determines which TF complex forms at any given site is poorly understood. Here, we demonstrate that high-throughput in vitro DNA binding assays coupled with unbiased computational analysis provide unprecedented insight into how different DNA sequences select distinct compositions and configurations of homeodomain TF complexes. Using inferred knowledge about minor groove width readout, we design targeted protein mutations that destabilize homeodomain binding both in vitro and in vivo in a complex-specific manner. By performing parallel systematic evolution of ligands by exponential enrichment sequencing (SELEX-seq), chromatin immunoprecipitation sequencing (ChIP-seq), RNA sequencing (RNA-seq), and Hi-C assays, we not only classify the majority of in vivo binding events in terms of complex composition but also infer complex-specific functions by perturbing the gene regulatory network controlled by a single complex.

Low affinity binding sites in an activating CRM mediate negative autoregulation of the Drosophila Hox gene Ultrabithorax.

Specification of cell identity and the proper functioning of a mature cell depend on precise regulation of gene expression. Both binary ON/OFF regulation of transcription, as well as more fine-tuned control of transcription levels in the ON state, are required to define cell types. The Drosophila melanogaster Hox gene, Ultrabithorax (Ubx), exhibits both of these modes of control during development. While ON/OFF regulation is needed to specify the fate of the developing wing (Ubx OFF) and haltere (Ubx ON), the levels of Ubx within the haltere differ between compartments along the proximal-distal axis. Here, we identify and molecularly dissect the novel contribution of a previously identified Ubx cis-regulatory module (CRM), anterobithorax (abx), to a negative auto-regulatory loop that decreases Ubx expression in the proximal compartment of the haltere as compared to the distal compartment. We find that Ubx, in complex with the known Hox cofactors, Homothorax (Hth) and Extradenticle (Exd), acts through low-affinity Ubx-Exd binding sites to reduce the levels of Ubx transcription in the proximal compartment. Importantly, we also reveal that Ubx-Exd-binding site mutations sufficient to result in de-repression of abx activity in a transgenic context are not sufficient to de-repress Ubx expression when mutated at the endogenous locus, suggesting the presence of multiple mechanisms through which Ubx-mediated repression occurs. Our results underscore the complementary nature of CRM analysis through transgenic reporter assays and genome modification of the endogenous locus; but, they also highlight the increasing need to understand gene regulation within the native context to capture the potential input of multiple genomic elements on gene control.

Accurate and sensitive quantification of protein-DNA binding affinity.

Transcription factors (TFs) control gene expression by binding to genomic DNA in a sequence-specific manner. Mutations in TF binding sites are increasingly found to be associated with human disease, yet we currently lack robust methods to predict these sites. Here, we developed a versatile maximum likelihood framework named No Read Left Behind (NRLB) that infers a biophysical model of protein-DNA recognition across the full affinity range from a library of in vitro selected DNA binding sites. NRLB predicts human Max homodimer binding in near-perfect agreement with existing low-throughput measurements. It can capture the specificity of the p53 tetramer and distinguish multiple binding modes within a single sample. Additionally, we confirm that newly identified low-affinity enhancer binding sites are functional in vivo, and that their contribution to gene expression matches their predicted affinity. Our results establish a powerful paradigm for identifying protein binding sites and interpreting gene regulatory sequences in eukaryotic genomes.

Origins and Specification of the Drosophila Wing.

The insect wing is a key evolutionary innovation that was essential for insect diversification. Yet despite its importance, there is still debate about its evolutionary origins. Two main hypotheses have been proposed: the paranotal hypothesis, which suggests that wings evolved as an extension of the dorsal thorax, and the gill-exite hypothesis, which proposes that wings were derived from a modification of a pre-existing branch at the dorsal base (subcoxa) of the leg. Here, we address this question by studying how wing fates are initially specified during Drosophila embryogenesis, by characterizing a cis-regulatory module (CRM) from the snail (sna) gene, sna-DP (for dorsal primordia). sna-DP specifically marks the early primordia for both the wing and haltere, collectively referred to as the DP. We found that the inputs that activate sna-DP are distinct from those that activate Distalless, a marker for leg fates. Further, in genetic backgrounds in which the leg primordia are absent, the DP are still partially specified. However, lineage-tracing experiments demonstrate that cells from the early leg primordia contribute to both ventral and dorsal appendage fates. Together, these results suggest that the wings of Drosophila have a dual developmental origin: two groups of cells, one ventral and one more dorsal, give rise to the mature wing. We suggest that the dual developmental origins of the wing may be a molecular remnant of the evolutionary history of this appendage, in which cells of the subcoxa of the leg coalesced with dorsal outgrowths to evolve a dorsal appendage with motor control.

Connectivity of vertebrate genomes: Paired-related homeobox (Prrx) genes in spotted gar, basal teleosts, and tetrapods.

Teleost fish are important models for human biology, health, and disease. Because genome duplication in a teleost ancestor (TGD) impacts the evolution of teleost genome structure and gene repertoires, we must discriminate gene functions that are shared and ancestral from those that are lineage-specific in teleosts or tetrapods to accurately apply inferences from teleost disease models to human health. Generalizations must account both for the TGD and for divergent evolution between teleosts and tetrapods after the likely two rounds of genome duplication shared by all vertebrates. Progress in sequencing techniques provides new opportunities to generate genomic and transcriptomic information from a broad range of phylogenetically informative taxa that facilitate detailed understanding of gene family and gene function evolution. We illustrate here the use of new sequence resources from spotted gar (Lepisosteus oculatus), a rayfin fish that diverged from teleosts before the TGD, as well as RNA-Seq data from gar and multiple teleost lineages to reconstruct the evolution of the Paired-related homeobox (Prrx) transcription factor gene family, which is involved in the development of mesoderm and neural crest-derived mesenchyme. We show that for Prrx genes, the spotted gar genome and gene expression patterns mimic mammals better than teleosts do. Analyses force the seemingly paradoxical conclusion that regulatory mechanisms for the limb expression domains of Prrx genes existed before the evolution of paired appendages. Detailed evolutionary analyses like those reported here are required to identify fish species most similar to the human genome to optimally connect fish models to human gene functions in health and disease.