Animal traits develop through intricate patterns of gene expression that are regulated at multiple levels. This regulation includes interactions between transcription factors and cis-regulatory element (CRE) DNA sequences, and the dynamic accessibility of CREs due to chromatin modifications and remodeling. Polycomb Group (PcG) and Trithorax Group (TrxG) genes are evolutionarily conserved regulators of chromatin state. Although the PcG and TrxG genes have well understood roles in developmental gene regulation, the extent to which these complexes contribute to trait evolution remains unclear. Here, we performed a genetic screen to understand how PcG and TrxG genes shape the rapidly evolving gene regulatory network (GRN) responsible for the dimorphic abdomen tergite pigmentation of Drosophila (D.) melanogaster fruit flies. A near-comprehensive screen of TrxG and PcG genes was conducted using RNAi, and numerous genes were identified whose reduced expression caused alterations to tergite pigmentation. For the eight most impactful genes, their roles in regulating this GRN were explored. We assessed their effects on several key transcription factors, and downstream CREs that are responsible for the expression of the GRN's pigmentation enzyme genes. We show that multiple members of the PcG and TrxG complexes are required for distinct tiers in the pigmentation GRN, suggesting potential points where responsive elements for different factors may have evolved. The results set the stage for future studies to identify the direct GRN targets of PcG and TrxG complexes, and how they have participated in the evolution of this GRN.

A major hurdle in understanding the molecular changes responsible for metazoan diversity is the characterization of cis-regulatory elements (CREs) for gene regulatory networks (GRNs). CRE changes are suspected to be commonplace in trait evolution, since such changes circumvent the deleterious effects of pleiotropy. A growing list of genes, though, is known to be regulated by redundant CREs. Such redundant CRE architectures complicate the characterization of GRN evolution, as they compound the effort to characterize each locus, and raise the questions of how and whether genes with redundant architectures evolve expression. Here, we used the evolution of sexually dimorphic abdomen pigmentation of Drosophila (D.) melanogaster as a model to study the function and evolution of CREs. Numerous sequences were evaluated that were previously predicted as potential abdomen CREs. Most of these predictions were validated, including two, four, and ten that, respectively, reside in the homothorax, grainy head, and Eip74EF transcription factor loci. The homothorax CREs were found to be partially redundant for this gene's pigmentation function, and pupal-stage Homothorax expression and the CRE activities were conserved among Drosophila species with the derived dimorphic and ancestral monomorphic phenotypes. Similarly, the Eip74EF CREs were conserved in the monomorphic D. willistoni. Thus, this gene's extensive CRE spatiotemporal redundancy has been conserved for over 30 million years, predating the dimorphic trait. Pigmentation evolution has been connected elsewhere to changes in nonredundant CREs. When these traits evolve, GRN changes may be biased towards the genes with singular nonredundant CREs, while the expression of redundantly regulated genes remains conserved.

Changes in gene expression are a prominent feature of morphological evolution. These changes occur to hierarchical gene regulatory networks (GRNs) of transcription factor genes that regulate the expression of trait-building differentiation genes. While changes in the expression of differentiation genes are essential to phenotypic evolution, they can be caused by mutations within cis-regulatory elements (CREs) that drive their expression (cis-evolution) or within genes for CRE-interacting transcription factors (trans-evolution). Locating these mutations remains a challenge, especially when experiments are limited to one species that possesses the ancestral or derived phenotype. We investigated CREs that control the expression of the differentiation genes tan and yellow, the expression of which evolved during the gain, modification, and loss of dimorphic pigmentation among Sophophora fruit flies. We show these CREs to be necessary components of a pigmentation GRN, as deletion from Drosophila melanogaster (derived dimorphic phenotype) resulted in lost expression and lost male-specific pigmentation. We evaluated the ability of orthologous CRE sequences to drive reporter gene expression in species with modified (Drosophila auraria), secondarily lost (Drosophila ananassae), and ancestrally absent (Drosophila willistoni) pigmentation. We show that the transgene host frequently determines CRE activity, implicating trans-evolution as a significant factor for this trait's diversity. We validated the gain of dimorphic Bab transcription factor expression as a trans-change contributing to the dimorphic trait. Our findings suggest an amenability to change for the landscape of trans-regulators and begs for an explanation as to why this is so common compared to the evolution of differentiation gene CREs.

Animal traits develop through the expression and action of numerous regulatory and realizator genes that comprise a gene regulatory network (GRN). For each GRN, its underlying patterns of gene expression are controlled by cis -regulatory elements (CREs) that bind activating and repressing transcription factors. These interactions drive cell-type and developmental stage-specific transcriptional activation or repression. Most GRNs remain incompletely mapped, and a major barrier to this daunting task is CRE identification. Here, we used an in silico method to identify predicted CREs (pCREs) that comprise the GRN which governs sex-specific pigmentation of Drosophila melanogaster . Through in vivo assays, we demonstrate that many pCREs activate expression in the correct cell-type and developmental stage. We employed genome editing to demonstrate that two CREs control the pupal abdomen expression of trithorax , whose function is required for the dimorphic phenotype. Surprisingly, trithorax had no detectable effect on this GRN’s key trans -regulators, but shapes the sex-specific expression of two realizator genes. Comparison of sequences orthologous to these CREs supports an evolutionary scenario where these trithorax CREs predated the origin of the dimorphic trait. Collectively, this study demonstrates how in silico approaches can shed novel insights on the GRN basis for a trait’s development and evolution.

Traits that appear discontinuously along phylogenies may be explained by independent origins (homoplasy) or repeated loss (homology). While discriminating between these models is difficult, the dissection of gene regulatory networks (GRNs) which drive the development of such repeatedly occurring traits can offer a mechanistic window on this fundamental problem. The GRN responsible for the male-specific pattern of Drosophila (D.) melanogaster melanic tergite pigmentation has received considerable attention. In this system, a metabolic pathway of pigmentation enzyme genes is expressed in spatial and sex-specific (i.e., dimorphic) patterns. The dimorphic expression of several genes is regulated by the Bab transcription factors, which suppress pigmentation enzyme expression in females, by virtue of their high expression in this sex. Here, we analyzed the phylogenetic distribution of species with male-specific pigmentation and show that this dimorphism is phylogenetically widespread among fruit flies. The analysis of pigmentation enzyme gene expression in distantly related dimorphic and monomorphic species shows that dimorphism is driven by the similar deployment of a conserved metabolic pathway. However, sexually dimorphic Bab expression was found only in D. melanogaster and its close relatives. These results suggest that dimorphism evolved by parallel deployment of differentiation genes, but was derived through distinct architectures at the level of regulatory genes. This work demonstrates the interplay of constraint and flexibility within evolving GRNs, findings that may foretell the mechanisms of homoplasy more broadly.

A challenge for evolutionary research is to uncover how new morphological traits evolve the coordinated spatial and temporal expression patterns of genes that govern their formation during development. Detailed studies are often limited to characterizing how one or a few genes contributed to a trait's emergence, and thus our knowledge of how entire GRNs evolve their coordinated expression of each gene remains unresolved. The melanic color patterns decorating the male abdominal tergites of Drosophila (D.) melanogaster evolved in part by novel expression patterns for genes acting at the terminus of a pigment metabolic pathway, driven by cis-regulatory elements (CREs) with distinct mechanisms of Hox regulation. Here, we examined the expression and evolutionary histories of two important enzymes in this pathway, encoded by the pale and Ddc genes. We found that while both genes exhibit dynamic patterns of expression, a robust pattern of Ddc expression specifically evolved in the lineage of fruit flies with pronounced melanic abdomens. Derived Ddc expression requires the activity of a CRE previously shown to activate expression in response to epidermal wounding. We show that a binding site for the Grainy head transcription factor that promotes the ancestral wound healing function of this CRE is also required for abdominal activity. Together with previous findings in this system, our work shows how the GRN for a novel trait emerged by assembling unique yet similarly functioning CREs from heterogeneous starting points.

Posterior prevalence is the general property attributed to HOX proteins describing the dominant effect of more posterior HOX proteins over the function of anterior orthologs in common areas of expression. To explore the HOX group 13 protein domains required for this property, we used the mouse Prx-1 promoter to drive transgenic expression of Hox constructs throughout the entire limb bud during development. This system allowed us to conclusively demonstrate a hierarchy of Hox function in developing limbs. Furthermore, by substituting the HOXD11 or HOXA9 homeodomain for that of HOXD13, we show that a HOXD13 homeodomain is not necessary for posterior prevalence. Proximal expression of these chimeric proteins unexpectedly caused defects consistent with wild-type HOXD13 mediated posterior prevalence. Moreover, group 13 non-homeodomain residues appear to confer the property as proximal expression of HOXA9 containing the HOXD13 homeodomain did not result in limb reductions characteristic of HOXD13. These data are most compatible with models of posterior prevalence based on protein–protein interactions and support examination of the N-terminal non-homeodomain regions of Hox group 13 proteins as necessary agents for posterior prevalence.

Hox genes encode transcription factors that regulate the morphogenesis of developing embryos. In mammals, knowledge of the genetic pathways, including the possible direct or indirect targets, regulated by HOX proteins is extremely limited. To identify the downstream genes regulated by posterior HOX proteins, we expressed HOXA13 in mouse embryonic fibroblasts lacking paralog group 13 expression using a bicistronic HOXA13/EGFP retroviral vector. Microarray analysis identified 68 genes with significant, reproducible RNA expression changes (50 activated; 18 repressed) in stable HOXA13-expressing cells. Genes with the GO annotation terms "extracellular matrix" and "basement membrane" were greatly overrepresented, and several were shown to be regulated by HOX proteins in other studies. Among the genes strongly activated by HOXA13 were Enpp2, a bifunctional enzyme known to modulate tumor and normal cell motility and which is expressed in precartilaginous condensations; Fhl1, a transcription factor implicated in muscle cell differentiation and development; and M32486, a putative integral membrane molecule expressed in the female reproductive tract. Expression differences in the HOXA13-expressing cells were confirmed for selected downstream genes using semi-quantitative RT-PCR, and in vivo coexpression with Hoxa13 in the limb interdigital mesenchyme was demonstrated for many. For two candidates, Igfbp4 and Fstl, interdigital limb bud expression was reduced in Hoxa13 mutants. To explore whether paralogous and nonparalogous HOX proteins could regulate the same genes, we created new HOX cell lines and examined the expression of selected genes identified by the HOXA13 screen. HOXD13 similarly activated/repressed 6 tested candidates, demonstrating that multiple downstream genetic pathways may be regulated by paralog HOX proteins. In contrast, HOXA9 was only able to repress expression of some gene targets. A HOXD13 mutant, HOXD13IQN > AAA, incapable of monomeric DNA-binding, activated the expression of 5 HOXA13-upregulated genes; but was incapable of repressing the expression of Ngef and Casp8ap2. Our results suggest that HOX protein–protein interactions without direct HOX DNA-binding may play a larger role in HOX transcriptional regulation than generally assumed, and DNA-binding appears critical for repression.