The origination and diversification of morphological characteristics represents a key problem in understanding the evolution of development. Morphological traits result from gene regulatory networks (GRNs) that form a web of transcription factors, which regulate multiple cis-regulatory element (CRE) sequences to control the coordinated expression of differentiation genes. The formation and modification of GRNs must ultimately be understood at the level of individual regulatory linkages (i.e., transcription factor binding sites within CREs) that constitute the network. Here, we investigate how elements within a network originated and diversified to generate a broad range of abdominal pigmentation phenotypes among Sophophora fruit flies. Our data indicates that the coordinated expression of two melanin synthesis enzymes, Yellow and Tan, recently evolved through novel CRE activities that respond to the spatial patterning inputs of Hox proteins and the sex-specific input of Bric-à-brac transcription factors. Once established, it seems that these newly evolved activities were repeatedly modified by evolutionary changes in the network's trans-regulators to generate large-scale changes in pigment pattern. By elucidating how yellow and tan are connected to the web of abdominal trans-regulators, we discovered that the yellow and tan abdominal CREs are composed of distinct regulatory inputs that exhibit contrasting responses to the same Hox proteins and Hox cofactors. These results provide an example in which CRE origination underlies a recently evolved novel trait, and highlights how coordinated expression patterns can evolve in parallel through the generation of unique regulatory linkages.

Transcriptional cis-regulatory modules (CRMs), or enhancers, are responsible for directing gene expression in specific territories and cell types during development. In some instances, the same gene may be served by two or more enhancers with similar specificities. Here we show that the utilization of dual, or "shadow", enhancers is a common feature of genes that are active specifically in neural precursor (NP) cells in Drosophila. By genome-wide computational discovery of statistically significant clusters of binding motifs for both proneural activator (P) proteins and basic helix-loop-helix (bHLH) repressor (R) factors (a "P+R" regulatory code), we have identified NP-specific enhancer modules associated with multiple genes expressed in this cell type. These CRMs are distinct from those previously identified for the corresponding gene, establishing the existence of a dual-enhancer arrangement in which both modules reside close to the gene they serve. Using wild-type and mutant reporter gene constructs in vivo, we show that P sites in these modules mediate activation by proneural factors in "proneural cluster" territories, whereas R sites mediate repression by bHLH repressors, which serves to restrict expression specifically to NP cells. To our knowledge, our results identify the first direct targets of these bHLH repressors. Finally, using genomic rescue constructs for neuralized (neur), we demonstrate that each of the gene's two NP-specific enhancers is sufficient to rescue neur function in the lateral inhibition process by which adult sensory organ precursor (SOP) cells are specified, but that deletion of both enhancers results in failure of this event.

Among developmental control genes, transcription factor-target gene "linkages"--the direct connections between target genes and the factors that control their patterns of expression--can show remarkable evolutionary stability. However, the specific binding sites that mediate and define these regulatory connections are themselves often subject to rapid turnover. Here we describe several instances in which particular transcription factor binding motif combinations have evidently been conserved upstream of orthologous target genes for extraordinarily long evolutionary periods. This occurs against a backdrop in which other binding sites for the same factors are coming and going rapidly. Our examples include a particular Dpp Silencer Element upstream of insect brinker genes, in combination with a novel motif we refer to as the Downstream Element; combinations of a Suppressor of Hairless Paired Site (SPS) and a specific proneural protein binding site associated with arthropod Notch pathway target genes; and a three-motif combination, also including an SPS, upstream of deuterostome Hes repressor genes, which are also Notch targets. We propose that these stable motif architectures have been conserved intact from a deep ancestor, in part because they mediate a special mode of regulation that cannot be supplied by the other, unstable motif instances.

The evolutionary origins of complex morphological structures such as the vertebrate eye or insect wing remain one of the greatest mysteries of biology. Recent comparative studies of gene expression imply that new structures are not built from scratch, but rather form by co-opting preexisting gene networks. A key prediction of this model is that upstream factors within the network will activate their preexisting targets (i.e., enhancers) to form novel anatomies. Here, we show how a recently derived morphological novelty present in the genitalia of D. melanogaster employs an ancestral Hox-regulated network deployed in the embryo to generate the larval posterior spiracle. We demonstrate how transcriptional enhancers and constituent transcription factor binding sites are used in both ancestral and novel contexts. These results illustrate network co-option at the level of individual connections between regulatory genes and highlight how morphological novelty may originate through the co-option of networks controlling seemingly unrelated structures.

Diverse traits often covary between species [1-3]. The possibility that a single mutation could contribute to the evolution of several characters between species [3] is rarely investigated as relatively few cases are dissected at the nucleotide level. Drosophila santomea has evolved additional sex comb sensory teeth on its legs and has lost two sensory bristles on its genitalia. We present evidence that a single nucleotide substitution in an enhancer of the scute gene contributes to both changes. The mutation alters a binding site for the Hox protein Abdominal-B in the developing genitalia, leading to bristle loss, and for another factor in the developing leg, leading to bristle gain. Our study suggests that morphological evolution between species can occur through a single nucleotide change affecting several sexually dimorphic traits. VIDEO ABSTRACT.

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.

Male genitalia are usually extremely divergent between closely related species, but relatively constant within one species. Here we examine the effect of temperature on the shape of the ventral branches, a male genital structure involved in reproductive isolation, in the sister species Drosophila santomea and Drosophila yakuba. We designed a semi-automatic measurement machine learning pipeline that can reliably identify curvatures and landmarks based on manually digitized contours of the ventral branches. With this method, we observed that temperature does not affect ventral branches in D. yakuba but that in D. santomea ventral branches tend to morph into a D. yakuba-like shape at lower temperature. We found that male genitalia structures involved in reproductive isolation can be relatively variable within one species and can resemble the shape of closely related species' genitalia through plasticity to temperature. Our results suggest that reproductive isolation mechanisms can be dependent on the environmental context.

Changes in gene regulation represent an important path to generate developmental differences affecting anatomical traits. Interspecific divergence in gene expression often results from changes in transcription-stimulating enhancer elements. While gene repression is crucial for precise spatiotemporal expression patterns, the relative contribution of repressive transcriptional silencers to regulatory evolution remains to be addressed. Here, we show that the Drosophila pigmentation gene ebony has mainly evolved through changes in the spatial domains of silencers patterning its abdominal expression. By precisely editing the endogenous ebony locus of D. melanogaster, we demonstrate the requirement of two redundant abdominal enhancers and three silencers that repress the redundant enhancers in a patterned manner. We observe a role for changes in these silencers in every case of ebony evolution observed to date. Our findings suggest that negative regulation by silencers likely has an under-appreciated role in gene regulatory evolution.

The model organism Drosophila melanogaster has become a focal system for investigations of rapidly evolving genital morphology as well as the development and functions of insect reproductive structures. To follow up on a previous paper outlining unifying terminology for the structures of the male terminalia in this species, we offer here a detailed description of the female terminalia of D. melanogaster. Informative diagrams and micrographs are presented to provide a comprehensive overview of the external and internal reproductive structures of females. We propose a collection of terms and definitions to standardize the terminology associated with the female terminalia in D. melanogaster and we provide a correspondence table with the terms previously used. Unifying terminology for both males and females in this species will help to facilitate communication between various disciplines, as well as aid in synthesizing research across publications within a discipline that has historically focused principally on male features. Our efforts to refine and standardize the terminology should expand the utility of this important model system for addressing questions related to the development and evolution of animal genitalia, and morphology in general.

Reproductive isolation is an intrinsic aspect of species formation. For that reason, the identification of the precise isolating traits, and the rates at which they evolve, is crucial to understanding how species originate and persist. Previous work has measured the rates of evolution of prezygotic and postzygotic barriers to gene flow, yet no systematic analysis has studied the rates of evolution of postmating-prezygotic (PMPZ) barriers. We measured the magnitude of two barriers to gene flow that act after mating occurs but before fertilization. We also measured the magnitude of a premating barrier (female mating rate in nonchoice experiments) and two postzygotic barriers (hybrid inviability and hybrid sterility) for all pairwise crosses of all nine known extant species within the melanogaster subgroup. Our results indicate that PMPZ isolation evolves faster than hybrid inviability but slower than premating isolation. Next, we partition postzygotic isolation into different components and find that, as expected, hybrid sterility evolves faster than hybrid inviability. These results lend support for the hypothesis that, in Drosophila, reproductive isolation mechanisms (RIMs) that act early in reproduction (or in development) tend to evolve faster than those that act later in the reproductive cycle. Finally, we tested whether there was evidence for reinforcing selection at any RIM. We found no evidence for generalized evolution of reproductive isolation via reinforcement which indicates that there is no pervasive evidence of this evolutionary process. Our results indicate that PMPZ RIMs might have important evolutionary consequences in initiating speciation and in the persistence of new species.

Sex-limited polymorphisms are an intriguing form of sexual dimorphism that offer unique opportunities to reconstruct the evolutionary changes that decouple male and female traits encoded by a shared genome. We investigated the genetic basis of a Mendelian female-limited color dimorphism (FLCD) that segregates in natural populations of more than 20 species of the Drosophila montium subgroup. In these species, females have alternative abdominal color morphs, light and dark, whereas males have only one color morph in each species. A comprehensive molecular phylogeny of the montium subgroup supports multiple origins of FLCD. Despite this, we mapped FLCD to the same locus in four distantly related species-the transcription factor POU domain motif 3 (pdm3), which acts as a repressor of abdominal pigmentation in D. melanogaster. In D. serrata, FLCD maps to a structural variant in the first intron of pdm3; however, this variant is not found in the three other species-D. kikkawai, D. leontia, and D. burlai-and sequence analysis strongly suggests the pdm3 alleles responsible for FLCD originated independently at least three times. We propose that cis-regulatory changes in pdm3 form sexually dimorphic and monomorphic alleles that segregate within species and are preserved, at least in one species, by structural variation. Surprisingly, pdm3 has not been implicated in the evolution of sex-specific pigmentation outside the montium subgroup, suggesting that the genetic paths to sexual dimorphism may be constrained within a clade but variable across clades.

Enhancers activate gene transcription in spatial and temporal patterns by interactions with gene promoters. These elements typically reside distal to their target promoter, with which they must interact selectively. Additional elements may contribute to enhancer-promoter specificity, including remote control element sequences within enhancers, tethering elements near promoters, and insulator/boundary elements that disrupt off-target interactions. However, few of these elements have been mapped, and as a result, the mechanisms by which these elements interact remain poorly understood. One impediment is their method of study, namely reporter transgenes in which enhancers are placed adjacent to a heterologous promoter, which may circumvent mechanisms controlling enhancer-promoter specificity and long-range interactions. Here, we report an optimized dual reporter transgene system in Drosophila melanogaster that allows the simultaneous comparison of an enhancer's ability to activate proximal and distal fluorescent reporter genes. Testing a panel of fluorescent transgenes in vivo, we found a two-protein combination that allows simultaneous measurement with minimal detection interference. We note differences among four tested enhancers in their ability to regulate a distally placed reporter transgene. These results suggest that enhancers differ in their requirements for promoter interaction and raise important practical considerations when studying enhancer function.

Hox genes pattern the anterior-posterior axis of animals and are posited to drive animal body plan evolution, yet their precise role in evolution has been difficult to determine. Here, we identified evolutionary modifications in the Hox gene Abd-B that dramatically altered its expression along the body plan of Drosophila santomea. Abd-B is required for pigmentation in Drosophila yakuba, the sister species of D. santomea, and changes to Abd-B expression would be predicted to make large contributions to the loss of body pigmentation in D. santomea. However, manipulating Abd-B expression in current-day D. santomea does not affect pigmentation. We attribute this epistatic interaction to four other genes within the D. santomea pigmentation network, three of which have evolved expression patterns that do not respond to Abd-B. Our results demonstrate how body plans may evolve through small evolutionary steps distributed throughout Hox-regulated networks. Polygenicity and epistasis may hinder efforts to identify genes and mechanisms underlying macroevolutionary traits.

Genome-scale sequence data have invigorated the study of hybridization and introgression, particularly in animals. However, outside of a few notable cases, we lack systematic tests for introgression at a larger phylogenetic scale across entire clades. Here, we leverage 155 genome assemblies from 149 species to generate a fossil-calibrated phylogeny and conduct multilocus tests for introgression across 9 monophyletic radiations within the genus Drosophila. Using complementary phylogenomic approaches, we identify widespread introgression across the evolutionary history of Drosophila. Mapping gene-tree discordance onto the phylogeny revealed that both ancient and recent introgression has occurred across most of the 9 clades that we examined. Our results provide the first evidence of introgression occurring across the evolutionary history of Drosophila and highlight the need to continue to study the evolutionary consequences of hybridization and introgression in this genus and across the tree of life.

Molecular evolutionary studies usually focus on genes with clear roles in adult fitness or on developmental genes expressed at multiple time points during the life of the organism. Here, we examine the evolutionary dynamics of Drosophila glue genes, a set of eight genes tasked with a singular primary function during a specific developmental stage: the production of glue that allows animal pupa to attach to a substrate for several days during metamorphosis. Using phenotypic assays and available data from transcriptomics, PacBio genomes, and sequence variation from global populations, we explore the selective forces acting on glue genes within the cosmopolitan Drosophila melanogaster species and its five closely related species, D. simulans, D. sechellia, D. mauritiana, D. yakuba, and D. teissieri. We observe a three-fold difference in glue adhesion between the least and the most adhesive D. melanogaster strain, indicating a strong genetic component to phenotypic variation. These eight glue genes are among the most highly expressed genes in salivary glands yet they display no notable codon bias. New copies of Sgs3 and Sgs7 are found in D. yakuba and D. teissieri with the Sgs3 coding sequence evolving rapidly after duplication in the D. yakuba branch. Multiple sites along the various glue genes appear to be constrained. Our population genetics analysis in D. melanogaster suggests signals of local adaptive evolution for Sgs3, Sgs5, and Sgs5bis and traces of selective sweeps for Sgs1, Sgs3, Sgs7, and Sgs8. Our work shows that stage-specific genes can be subjected to various dynamic evolutionary forces.

The vinegar fly Drosophila melanogaster is a pivotal model for invertebrate development, genetics, physiology, neuroscience, and disease. The whole family Drosophilidae, which contains over 4,400 species, offers a plethora of cases for comparative and evolutionary studies. Despite a long history of phylogenetic inference, many relationships remain unresolved among the genera, subgenera, and species groups in the Drosophilidae. To clarify these relationships, we first developed a set of new genomic markers and assembled a multilocus data set of 17 genes from 704 species of Drosophilidae. We then inferred a species tree with highly supported groups for this family. Additionally, we were able to determine the phylogenetic position of some previously unplaced species. These results establish a new framework for investigating the evolution of traits in fruit flies, as well as valuable resources for systematics.