Gain of cis-regulatory activities underlies novel domains of wingless gene expression in Drosophila

Beyond modular enhancers: new questions in cis-regulatory evolution

Our understanding of how cis-regulatory elements work has advanced rapidly, outpacing our evolutionary models. In this review, we consider the implications of new mechanistic findings for evolutionary developmental biology. We focus on three different debates: whether evolutionary innovation occurs more often via the modification of old cis-regulatory elements or the emergence of new ones; the extent to which individual elements are specific and autonomous or multifunctional and interdependent; and how the robustness of cis-regulatory architectures influences the rate of trait evolution. These discussions lead us to propose new questions for the evo-devo of cis-regulation.

The origin and evolution of Wnt signalling

The Wnt signal transduction pathway has essential roles in the formation of the primary body axis during development, cellular differentiation and tissue homeostasis. This animal-specific pathway has been studied extensively in contexts ranging from developmental biology to medicine for more than 40 years. Despite its physiological importance, an understanding of the evolutionary origin and primary function of Wnt signalling has begun to emerge only recently. Recent studies on very basal metazoan species have shown high levels of conservation of components of both canonical and non-canonical Wnt signalling pathways. Furthermore, some pathway proteins have been described also in non-animal species, suggesting that recruitment and functional adaptation of these factors has occurred in metazoans. In this Review, we summarize the current state of research regarding the evolutionary origin of Wnt signalling, its ancestral function and the characteristics of the primal Wnt ligand, with emphasis on the importance of genomic studies in various pre-metazoan and basal metazoan species.

Thermal plasticity of wing size and wing spot size in Drosophila guttifera

Thermal plasticity of melanin pigmentation patterns in Drosophila species has been studied as a model to investigate developmental mechanisms of phenotypic plasticity. The developmental process of melanin pigmentation patterns on wings of Drosophila is divided into two parts, prepattern specification during the pupal period and wing vein-dependent transportation of melanin precursors after eclosion. Which part can be affected by thermal changes? To address this question, we used polka-dotted melanin spots on wings of Drosophila guttifera, whose spot areas are specified by wingless morphogen. In this research, we reared D. guttifera at different temperatures to test whether wing spots show thermal plasticity. We found that wing size becomes larger at lower temperature and that different spots have different reaction norms. Furthermore, we changed the rearing temperature in the middle of the pupal period and found that the most sensitive developmental periods for wing size and spot size are different. The results suggest that the size control mechanisms for the thermal plasticity of wing size and spot size are independent. We also found that the most sensitive stage for spot size was part of the pupal period including stages at which wingless is expressed in the polka-dotted pattern. Therefore, it is suggested that temperature change might affect the prepattern specification process and might not affect transportation through wing veins.

Species differences in hormonally mediated gene expression underlie the evolutionary loss of sexually dimorphic coloration in Sceloporus lizards

Phenotypic sexual dimorphism often involves the hormonal regulation of sex-biased expression for underlying genes. However, it is generally unknown whether the evolution of hormonally mediated sexual dimorphism occurs through upstream changes in tissue sensitivity to hormone signals, downstream changes in responsiveness of target genes, or both. Here, we use comparative transcriptomics to explore these possibilities in 2 species of Sceloporus lizards exhibiting different patterns of sexual dichromatism. Sexually dimorphic S. undulatus develops blue and black ventral coloration in response to testosterone, while sexually monomorphic S. virgatus does not, despite exhibiting similar sex differences in circulating testosterone levels. We administered testosterone implants to juveniles of each species and used RNAseq to quantify gene expression in ventral skin. Transcriptome-wide responses to testosterone were stronger in S. undulatus than in S. virgatus, suggesting species differences in tissue sensitivity to this hormone signal. Species differences in the expression of genes for androgen metabolism and sex hormone-binding globulin were consistent with this idea, but expression of the androgen receptor gene was higher in S. virgatus, complicating this interpretation. Downstream of androgen signaling, we found clear species differences in hormonal responsiveness of genes related to melanin synthesis, which were upregulated by testosterone in S. undulatus, but not in S. virgatus. Collectively, our results indicate that hormonal regulation of melanin synthesis pathways contributes to the development of sexual dimorphism in S. undulatus, and that changes in the hormonal responsiveness of these genes in S. virgatus contribute to the evolutionary loss of ventral coloration.

The Genetics of Fitness Reorganization during the Transition to Multicellularity: The Volvocine regA-like Family as a Model

The evolutionary transition from single-celled to multicellular individuality requires organismal fitness to shift from the cell level to a cell group. This reorganization of fitness occurs by re-allocating the two components of fitness, survival and reproduction, between two specialized cell types in the multicellular group: soma and germ, respectively. How does the genetic basis for such fitness reorganization evolve? One possible mechanism is the co-option of life history genes present in the unicellular ancestors of a multicellular lineage. For instance, single-celled organisms must regulate their investment in survival and reproduction in response to environmental changes, particularly decreasing reproduction to ensure survival under stress. Such stress response life history genes can provide the genetic basis for the evolution of cellular differentiation in multicellular lineages. The regA-like gene family in the volvocine green algal lineage provides an excellent model system to study how this co-option can occur. We discuss the origin and evolution of the volvocine regA-like gene family, including regA-the gene that controls somatic cell development in the model organism Volvox carteri. We hypothesize that the co-option of life history trade-off genes is a general mechanism involved in the transition to multicellular individuality, making volvocine algae and the regA-like family a useful template for similar investigations in other lineages.

Human and African ape myosin heavy chain content and the evolution of hominin skeletal muscle

Humans are unique among terrestrial mammals in our manner of walking and running, reflecting 7 to 8 Ma of musculoskeletal evolution since diverging with the genus Pan. One component of this is a shift in our skeletal muscle biology towards a predominance of myosin heavy chain (MyHC) I isoforms (i.e. slow fibers) across our pelvis and lower limbs, which distinguishes us from chimpanzees. Here, new MyHC data from 35 pelvis and hind limb muscles of a Western gorilla (Gorilla gorilla) are presented. These data are combined with a similar chimpanzee dataset to assess the MyHC I content of humans in comparison to African apes (chimpanzees and gorillas) and other terrestrial mammals. The responsiveness of human skeletal muscle to behavioral interventions is also compared to the human-African ape differential. Humans are distinct from African apes and among a small group of terrestrial mammals whose pelvis and hind/lower limb muscle is slow fiber dominant, on average. Behavioral interventions, including immobilization, bed rest, spaceflight and exercise, can induce modest decreases and increases in human MyHC I content (i.e. -9.3% to 2.3%, n = 2033 subjects), but these shifts are much smaller than the mean human-African ape differential (i.e. 31%). Taken together, these results indicate muscle fiber content is likely an evolvable trait under selection in the hominin lineage. As such, we highlight potential targets of selection in the genome (e.g. regions that regulate MyHC content) that may play an important role in hominin skeletal muscle evolution.

The Genetic Mechanisms Underlying the Concerted Expression of the yellow and tan Genes in Complex Patterns on the Abdomen and Wings of Drosophila guttifera

How complex morphological patterns form is an intriguing question in developmental biology. However, the mechanisms that generate complex patterns remain largely unknown. Here, we sought to identify the genetic mechanisms that regulate the tan (t) gene in a multi-spotted pigmentation pattern on the abdomen and wings of Drosophila guttifera. Previously, we showed that yellow (y) gene expression completely prefigures the abdominal and wing pigment patterns of this species. In the current study, we demonstrate that the t gene is co-expressed with the y gene in nearly identical patterns, both transcripts foreshadowing the adult abdominal and wing melanin spot patterns. We identified cis-regulatory modules (CRMs) of t, one of which drives reporter expression in six longitudinal rows of spots on the developing pupal abdomen, while the second CRM activates the reporter gene in a spotted wing pattern. Comparing the abdominal spot CRMs of y and t, we found a similar composition of putative transcription factor binding sites that are thought to regulate the complex expression patterns of both terminal pigmentation genes y and t. In contrast, the y and t wing spots appear to be regulated by distinct upstream factors. Our results suggest that the D. guttifera abdominal and wing melanin spot patterns have been established through the co-regulation of y and t, shedding light on how complex morphological traits may be regulated through the parallel coordination of downstream target genes.

Enhancer architecture and chromatin accessibility constrain phenotypic space during Drosophila development

Developmental enhancers bind transcription factors and dictate patterns of gene expression during development. Their molecular evolution can underlie phenotypical evolution, but the contributions of the evolutionary pathways involved remain little understood. Here, using mutation libraries in Drosophila melanogaster embryos, we observed that most point mutations in developmental enhancers led to changes in gene expression levels but rarely resulted in novel expression outside of the native pattern. In contrast, random sequences, often acting as developmental enhancers, drove expression across a range of cell types; random sequences including motifs for transcription factors with pioneer activity acted as enhancers even more frequently. Our findings suggest that the phenotypic landscapes of developmental enhancers are constrained by enhancer architecture and chromatin accessibility. We propose that the evolution of existing enhancers is limited in its capacity to generate novel phenotypes, whereas the activity of de novo elements is a primary source of phenotypic novelty.

Lost in translation: Molecular basis of reduced flower coloration in a self-pollinated monkeyflower (Mimulus) species

Phenotypic evolution is usually attributed to changes in protein function or gene transcription. In principle, mutations that affect protein abundance through enhancing or attenuating protein translation also could be an important source for phenotypic evolution. However, these types of mutations remain largely unexplored in the studies of phenotypic variation in nature. Through fine-scale genetic mapping and functional interrogation, we identify a single nucleotide substitution in an anthocyanin-activating R2R3-MYB gene causing flower color variation between a pair of closely related monkeyflower (Mimulus) species, the hummingbird-pollinated Mimulus cardinalis, and self-pollinated Mimulus parishii. This causal mutation is located in the 5' untranslated region and generates an upstream ATG start codon, leading to attenuated protein translation and reduced flower coloration in the self-pollinated species. Together, our results provide empirical support for the role of mutations affecting protein translation, as opposed to protein function or transcript level, in natural phenotypic variation.

Differentially Active and Conserved Neural Enhancers Define Two Forms of Adaptive Noncoding Evolution in Humans

The human and chimpanzee genomes are strikingly similar, but our neural phenotypes are very different. Many of these differences are likely driven by changes in gene expression, and some of those changes may have been adaptive during human evolution. Yet, the relative contributions of positive selection on regulatory regions or other functional regulatory changes are unclear. Where are these changes located throughout the human genome? Are functional regulatory changes near genes or are they in distal enhancer regions? In this study, we experimentally combined both human and chimpanzee cis-regulatory elements (CREs) that showed either (1) signs of accelerated evolution in humans or (2) that have been shown to be active in the human brain. Using a massively parallel reporter assay, we tested the ability of orthologous human and chimpanzee CREs to activate transcription in induced pluripotent stem-cell-derived neural progenitor cells and neurons. With this assay, we identified 179 CREs with differential activity between human and chimpanzee; in contrast, we found 722 CREs with signs of positive selection in humans. Selection and differentially expressed CREs strikingly differ in level of expression, size, and genomic location. We found a subset of 69 CREs in loci with genetic variants associated with neuropsychiatric diseases, which underscores the consequence of regulatory activity in these loci for proper neural development and function. By combining CREs that either experienced recent selection in humans or CREs that are functional brain enhancers, presents a novel way of studying the evolution of noncoding elements that contribute to human neural phenotypes.

Genomic prediction of growth in a commercially, recreationally, and culturally important marine resource, the Australian snapper (Chrysophrys auratus)

Growth is one of the most important traits of an organism. For exploited species, this trait has ecological and evolutionary consequences as well as economical and conservation significance. Rapid changes in growth rate associated with anthropogenic stressors have been reported for several marine fishes, but little is known about the genetic basis of growth traits in teleosts. We used reduced genome representation data and genome-wide association approaches to identify growth-related genetic variation in the commercially, recreationally, and culturally important Australian snapper (Chrysophrys auratus, Sparidae). Based on 17,490 high-quality single-nucleotide polymorphisms and 363 individuals representing extreme growth phenotypes from 15,000 fish of the same age and reared under identical conditions in a sea pen, we identified 100 unique candidates that were annotated to 51 proteins. We documented a complex polygenic nature of growth in the species that included several loci with small effects and a few loci with larger effects. Overall heritability was high (75.7%), reflected in the high accuracy of the genomic prediction for the phenotype (small vs large). Although the single-nucleotide polymorphisms were distributed across the genome, most candidates (60%) clustered on chromosome 16, which also explains the largest proportion of heritability (16.4%). This study demonstrates that reduced genome representation single-nucleotide polymorphisms and the right bioinformatic tools provide a cost-efficient approach to identify growth-related loci and to describe genomic architectures of complex quantitative traits. Our results help to inform captive aquaculture breeding programs and are of relevance to monitor growth-related evolutionary shifts in wild populations in response to anthropogenic pressures.

Sex-specific evolution of a Drosophila sensory system via interacting cis- and trans-regulatory changes

The evolution of gene expression via cis-regulatory changes is well established as a major driver of phenotypic evolution. However, relatively little is known about the influence of enhancer architecture and intergenic interactions on regulatory evolution. We address this question by examining chemosensory system evolution in Drosophila. Drosophila prolongata males show a massively increased number of chemosensory bristles compared to females and males of sibling species. This increase is driven by sex-specific transformation of ancestrally mechanosensory organs. Consistent with this phenotype, the Pox neuro transcription factor (Poxn), which specifies chemosensory bristle identity, shows expanded expression in D. prolongata males. Poxn expression is controlled by nonadditive interactions among widely dispersed enhancers. Although some D. prolongata Poxn enhancers show increased activity, the additive component of this increase is slight, suggesting that most changes in Poxn expression are due to epistatic interactions between Poxn enhancers and trans-regulatory factors. Indeed, the expansion of D. prolongata Poxn enhancer activity is only observed in cells that express doublesex (dsx), the gene that controls sexual differentiation in Drosophila and also shows increased expression in D. prolongata males due to cis-regulatory changes. Although expanded dsx expression may contribute to increased activity of D. prolongata Poxn enhancers, this interaction is not sufficient to explain the full expansion of Poxn expression, suggesting that cis-trans interactions between Poxn, dsx, and additional unknown genes are necessary to produce the derived D. prolongata phenotype. Overall, our results demonstrate the importance of epistatic gene interactions for evolution, particularly when pivotal genes have complex regulatory architecture.

Common Themes and Future Challenges in Understanding Gene Regulatory Network Evolution

A major driving force behind the evolution of species-specific traits and novel structures is alterations in gene regulatory networks (GRNs). Comprehending evolution therefore requires an understanding of the nature of changes in GRN structure and the responsible mechanisms. Here, we review two insect pigmentation GRNs in order to examine common themes in GRN evolution and to reveal some of the challenges associated with investigating changes in GRNs across different evolutionary distances at the molecular level. The pigmentation GRN in Drosophila melanogaster and other drosophilids is a well-defined network for which studies from closely related species illuminate the different ways co-option of regulators can occur. The pigmentation GRN for butterflies of the Heliconius species group is less fully detailed but it is emerging as a useful model for exploring important questions about redundancy and modularity in cis-regulatory systems. Both GRNs serve to highlight the ways in which redeployment of trans-acting factors can lead to GRN rewiring and network co-option. To gain insight into GRN evolution, we discuss the importance of defining GRN architecture at multiple levels both within and between species and of utilizing a range of complementary approaches.

Repression precedes independent evolutionary gains of a highly specific gene expression pattern

Highly specific expression patterns can be caused by the overlapping activities of activator and repressor sequences in enhancers. However, few studies illuminate how these sequences evolve in the origin of new enhancers. Here, we show that expression of the bond gene in the semicircular wall epithelium (swe) of the Drosophila melanogaster male ejaculatory bulb (EB) is controlled by an enhancer consisting of an activator region that requires Abdominal-B driving expression in the entire EB and a repressor region that restricts this expression to the EB swe. Although this expression pattern is independently gained in the distantly related Scaptodrosophila lebanonensis and does not require Abdominal-B, we show that functionally similar repressor sequences are present in Scaptodrosophila and also in species that do not express bond in the EB. We suggest that during enhancer evolution, repressor sequences can precede the evolution of activator sequences and may lead to similar but independently evolved expression patterns.

Pu et al. show that the independent gain of a highly specific expression pattern across distantly related species may be because of the preexistence of repressor sequences that precedes the diversification of these species. This may reflect a general mechanism underlying the evolution of highly specific enhancers.

From Aedes to Zeugodacus: a review of dipteran body coloration studies regarding evolutionary developmental biology, pest control, and species discovery

Over the past two decades, evo-devo (evolution of development) studies have elucidated genetic mechanisms underlying novel dipteran body color patterns. Here we review the most recent developments, which show some departure from the model organism Drosophila melanogaster, leading the field into the investigation of more complex color patterns. We also discuss how the robust application of transgenic techniques has facilitated the study of many non-model pest species. Furthermore, we see that subtle pigmentation differences guide the discovery and description of new dipterans. Therefore, we argue that the existence of new field guides and the prevalence of pigmentation studies in non-model flies will enable scientists to adopt uninvestigated species into the lab, allowing them to study novel morphologies.

Predicting future from past: The genomic basis of recurrent and rapid stickleback evolution

Similar forms often evolve repeatedly in nature, raising long-standing questions about the underlying mechanisms. Here, we use repeated evolution in stickleback to identify a large set of genomic loci that change recurrently during colonization of freshwater habitats by marine fish. The same loci used repeatedly in extant populations also show rapid allele frequency changes when new freshwater populations are experimentally established from marine ancestors. Marked genotypic and phenotypic changes arise within 5 years, facilitated by standing genetic variation and linkage between adaptive regions. Both the speed and location of changes can be predicted using empirical observations of recurrence in natural populations or fundamental genomic features like allelic age, recombination rates, density of divergent loci, and overlap with mapped traits. A composite model trained on these stickleback features can also predict the location of key evolutionary loci in Darwin's finches, suggesting that similar features are important for evolution across diverse taxa.

Porcupine/Wntless-dependent trafficking of the conserved WntA ligand in butterflies

Wnt ligands are key signaling molecules in animals, but little is known about the evolutionary dynamics and mode of action of the WntA orthologs, which are not present in the vertebrates or in Drosophila. Here we show that the WntA subfamily evolved at the base of the Bilateria + Cnidaria clade, and conserved the thumb region and Ser209 acylation site present in most other Wnts, suggesting WntA requires the core Wnt secretory pathway. WntA proteins are distinguishable from other Wnts by a synapomorphic Iso/Val/Ala216 amino-acid residue that replaces the otherwise ubiquitous Thr216 position. WntA embryonic expression is conserved between beetles and butterflies, suggesting functionality, but the WntA gene was lost three times within arthropods, in podoplean copepods, in the cyclorrhaphan fly radiation, and in ensiferan crickets and katydids. Finally, CRISPR mosaic knockouts (KOs) of porcupine and wntless phenocopied the pattern-specific effects of WntA KOs in the wings of Vanessa cardui butterflies. These results highlight the molecular conservation of the WntA protein across invertebrates, and imply it functions as a typical Wnt ligand that is acylated and secreted through the Porcupine/Wntless secretory pathway.

Modeling the evolutionary architectures of transcribed human enhancer sequences reveals distinct origins, functions, and associations with human-trait variation

Despite the importance of gene regulatory enhancers in human biology and evolution, we lack a comprehensive model of enhancer evolution and function. This substantially limits our understanding of the genetic basis of species divergence and our ability to interpret the effects of noncoding variants on human traits.

To explore enhancer sequence evolution and its relationship to regulatory function, we traced the evolutionary origins of transcribed human enhancer sequences with activity across diverse tissues and cellular contexts from the FANTOM5 consortium. The transcribed enhancers are enriched for sequences of a single evolutionary age (“simple” evolutionary architectures) compared with enhancers that are composites of sequences of multiple evolutionary ages (“complex” evolutionary architectures), likely indicating constraint against genomic rearrangements. Complex enhancers are older, more pleiotropic, and more active across species than simple enhancers. Genetic variants within complex enhancers are also less likely to associate with human traits and biochemical activity. Transposable-element-derived sequences (TEDS) have made diverse contributions to enhancers of both architectures; the majority of TEDS are found in enhancers with simple architectures, while a minority have remodeled older sequences to create complex architectures. Finally, we compare the evolutionary architectures of transcribed enhancers with histone-mark-defined enhancers.

Our results reveal that most human transcribed enhancers are ancient sequences of a single age, and thus the evolution of most human enhancers was not driven by increases in evolutionary complexity over time. Our analyses further suggest that considering enhancer evolutionary histories provides context that can aid interpretation of the effects of variants on enhancer function. Based on these results, we propose a framework for analyzing enhancer evolutionary architecture.

The color pattern inducing gene wingless is expressed in specific cell types of campaniform sensilla of a polka-dotted fruit fly, Drosophila guttifera

A polka-dotted fruit fly, Drosophila guttifera, has a unique pigmentation pattern on its wings and is used as a model for evo-devo studies exploring the mechanism of evolutionary gain of novel traits. In this species, a morphogen-encoding gene, wingless, is expressed in species-specific positions and induces a unique pigmentation pattern. To produce some of the pigmentation spots on wing veins, wingless is thought to be expressed in developing campaniform sensillum cells, but it was unknown which of the four cell types there express(es) wingless. Here we show that two of the cell types, dome cells and socket cells, express wingless, as indicated by in situ hybridization together with immunohistochemistry. This is a unique case in which non-neuronal SOP (sensory organ precursor) progeny cells produce Wingless as an inducer of pigmentation pattern formation. Our finding opens a path to clarifying the mechanism of evolutionary gain of a unique wingless expression pattern by analyzing gene regulation in dome cells and socket cells.

RNA In Situ Hybridization for Detecting Gene Expression Patterns in the Abdomens and Wings of Drosophila Species

RNA in situ hybridization (ISH) is used to visualize spatio-temporal gene expression patterns with broad applications in biology and biomedicine. Here we provide a protocol for mRNA ISH in developing pupal wings and abdomens for model and non-model Drosophila species. We describe best practices in pupal staging, tissue preparation, probe design and synthesis, imaging of gene expression patterns, and image-editing techniques. This protocol has been successfully used to investigate the roles of genes underlying the evolution of novel color patterns in non-model Drosophila species.

R2R3-MYB genes control petal pigmentation patterning in Clarkia gracilis ssp. sonomensis (Onagraceae)

Petal pigmentation patterning is widespread in flowering plants. The genetics of these pattern elements has been of great interest for understanding the evolution of phenotypic diversification. Here, we investigate the genetic changes responsible for the evolution of an unpigmented petal element on a colored background. We used transcriptome analysis, gene expression assays, cosegregation in F₂ plants and functional tests to identify the gene(s) involved in petal coloration in Clarkia gracilis ssp. sonomensis. We identified an R2R3-MYB transcription factor (CgsMYB12) responsible for anthocyanin pigmentation of the basal region ('cup') in the petal of C. gracilis ssp. sonomensis. A functional mutation in CgsMYB12 creates a white cup on a pink petal background. Additionally, we found that two R2R3-MYB genes (CgsMYB6 and CgsMYB11) are also involved in petal background pigmentation. Each of these three R2R3-MYB genes exhibits a different spatiotemporal expression pattern. The functionality of these R2R3-MYB genes was confirmed through stable transformation of Arabidopsis. Distinct spatial patterns of R2R3-MYB expression have created the possibility that pigmentation in different sections of the petal can evolve independently. This finding suggests that recent gene duplication has been central to the evolution of petal pigmentation patterning in C. gracilis ssp. sonomensis.

The genetic basis for the evolution of soma: mechanistic evidence for the co-option of a stress-induced gene into a developmental master regulator

In multicellular organisms with specialized cells, the most significant distinction among cell types is between reproductive (germ) cells and non-reproductive/somatic cells (soma). Although soma contributed to the marked increase in complexity of many multicellular lineages, little is known about its evolutionary origins. We have previously suggested that the evolution of genes responsible for the differentiation of somatic cells involved the co-option of life history trade-off genes that in unicellular organisms enhanced survival at a cost to immediate reproduction. In the multicellular green alga, Volvox carteri, cell fate is established early in development by the differential expression of a master regulatory gene known as regA. A closely related RegA-Like Sequence (RLS1) is present in its single-celled relative, Chlamydomonas reinhardtii. RLS1 is expressed in response to stress, and we proposed that an environmentally induced RLS1-like gene was co-opted into a developmental pathway in the lineage leading to V. carteri. However, the exact evolutionary scenario responsible for the postulated co-option event remains to be determined. Here, we show that in addition to being developmentally regulated, regA can also be induced by environmental cues, indicating that regA has maintained its ancestral regulation. We also found that the absence of a functional RegA protein confers increased sensitivity to stress, consistent with RegA having a direct or indirect role in stress responses. Overall, this study (i) provides mechanistic evidence for the co-option of an environmentally induced gene into a major developmental regulator, (ii) supports the view that major morphological innovations can evolve via regulatory changes and (iii) argues for the role of stress in the evolution of multicellular complexity.

The modular expression patterns of three pigmentation genes prefigure unique abdominal morphologies seen among three Drosophila species

To understand how novel animal body colorations emerged, one needs to ask how the development of color patterns differs among closely related species. Here we examine three species of fruit flies - Drosophila guttifera (D. guttifera), D. palustris, and D. subpalustris - displaying a varying number of abdominal spot rows. Through in situ hybridization experiments, we examine the mRNA expression patterns for the pigmentation genes Dopa decarboxylase (Ddc), tan (t), and yellow (y) during pupal development. Our results show that Ddc, t, and y are co-expressed in modular, identical patterns, each foreshadowing the adult abdominal spots in D. guttifera, D. palustris, and D. subpalustris. We suggest that differences in the expression patterns of these three genes partially underlie the morphological diversity of the quinaria species group.

Polymorphism and Divergence of Novel Gene Expression Patterns in Drosophila melanogaster

Transcriptomes may evolve by multiple mechanisms, including the evolution of novel genes, the evolution of transcript abundance, and the evolution of cell, tissue, or organ expression patterns. Here, we focus on the last of these mechanisms in an investigation of tissue and organ shifts in gene expression in Drosophila melanogaster. In contrast to most investigations of expression evolution, we seek to provide a framework for understanding the mechanisms of novel expression patterns on a short population genetic timescale. To do so, we generated population samples of D. melanogaster transcriptomes from five tissues: accessory gland, testis, larval salivary gland, female head, and first-instar larva. We combined these data with comparable data from two outgroups to characterize gains and losses of expression, both polymorphic and fixed, in D. melanogaster We observed a large number of gain- or loss-of-expression phenotypes, most of which were polymorphic within D. melanogaster Several polymorphic, novel expression phenotypes were strongly influenced by segregating cis-acting variants. In support of previous literature on the evolution of novelties functioning in male reproduction, we observed many more novel expression phenotypes in the testis and accessory gland than in other tissues. Additionally, genes showing novel expression phenotypes tend to exhibit greater tissue-specific expression. Finally, in addition to qualitatively novel expression phenotypes, we identified genes exhibiting major quantitative expression divergence in the D. melanogaster lineage.

Wnt gene regulation and function during maxillary palp development in Drosophila melanogaster

Wnt genes encode secreted ligands that play many important roles in the development of metazoans. There are thirteen known Wnt gene subfamilies and seven of these are represented in Drosophila melanogaster. While wingless (wg) is the best understood and most widely studied Wnt gene in Drosophila, the functions of many of the other Drosophila Wnt genes are less well understood. For example, relatively little is known about Wnt6, which is an ancient paralog of wg and they form a conserved Wnt cluster together with Wnt9 (Dwnt4) and Wnt10. Wg and Wnt6 encode similar proteins and exhibit overlapping expression in several tissues during development. Both wg and Wnt6 were previously shown to regulate the development of maxillary palps, important olfactory organs in flies, but it remained unclear how these two ligands may combine to carry out specific functions and how this is regulated. Here, we have further analysed Wnt6 function in the context of maxillary palp development. Surprisingly, we found that Wnt6 does not appear to be necessary for development of maxillary palps. While a deletion of the 5' region of Wnt6 results in very small maxillary palps, we show that this effect is more likely to be a consequence of removing cis-regulatory elements that may regulate wg expression in this tissue rather than through the loss of Wnt6 function. Although, we cannot completely exclude the possibility that Wnt6 may subtly regulate maxillary palp development in combination with wg, our analysis of Wnt6 loss of function mutants suggests this ligand plays a more general role in regulating growth during development. Taken together our results provide new insights into maxillary palp formation and Wnt6 functions in Drosophila, and further evidence for a complex cis-regulatory landscape in the Wnt9-wg-Wnt6-Wnt10 cluster, which may help explain its evolutionary conservation.

Temporal flexibility of gene regulatory network underlies a novel wing pattern in flies

Developmental genes can be coopted to generate evolutionary novelties by changing their spatial regulation. However, developmental genes seldom act independently, but rather work in a gene regulatory network (GRN). How is it possible to recruit a single gene from a whole GRN? What are the properties that allow parallel cooptions of the same genes during evolution? Here, we show that a novel engrailed gene expression underlies a novel wing color pattern in flies. We show that cooption is facilitated 1) because of GRN flexibility over development and 2) because every single gene of the GRN has its own functional time window. We suggest these two temporal properties could explain why the same gene can be independently recruited several times during evolution.

Organisms have evolved endless morphological, physiological, and behavioral novel traits during the course of evolution. Novel traits were proposed to evolve mainly by orchestration of preexisting genes. Over the past two decades, biologists have shown that cooption of gene regulatory networks (GRNs) indeed underlies numerous evolutionary novelties. However, very little is known about the actual GRN properties that allow such redeployment. Here we have investigated the generation and evolution of the complex wing pattern of the fly Samoaia leonensis. We show that the transcription factor Engrailed is recruited independently from the other players of the anterior–posterior specification network to generate a new wing pattern. We argue that partial cooption is made possible because 1) the anterior–posterior specification GRN is flexible over time in the developing wing and 2) this flexibility results from the fact that every single gene of the GRN possesses its own functional time window. We propose that the temporal flexibility of a GRN is a general prerequisite for its possible cooption during the course of evolution.

Transcriptome analysis reveals wingless regulates neural development and signaling genes in the region of wing pigmentation of a polka-dotted fruit fly

How evolutionary novelties have arisen is one of the central questions in evolutionary biology. Preexisting gene regulatory networks or signaling pathways have been shown to be co-opted for building novel traits in several organisms. However, the structure of entire gene regulatory networks and evolutionary events of gene co-option for emergence of a novel trait are poorly understood. In this study, to explore the genetic and molecular bases of the novel wing pigmentation pattern of a polka-dotted fruit fly (Drosophila guttifera), we performed de novo genome sequencing and transcriptome analyses. As a result, we comprehensively identified the genes associated with the pigmentation pattern. Furthermore, we revealed that 151 of these associated genes were positively or negatively regulated by wingless, a master regulator of wing pigmentation. Genes for neural development, Wnt signaling, Dpp signaling, and effectors (such as enzymes) for melanin pigmentation were included among these 151 genes. None of the known regulatory genes that regulate pigmentation pattern formation in other fruit fly species were included. Our results suggest that the novel pigmentation pattern of a polka-dotted fruit fly might have emerged through multistep co-options of multiple gene regulatory networks, signaling pathways, and effector genes, rather than recruitment of one large gene circuit.

The Making of Transgenic Drosophila guttifera

The complex color patterns on the wings and body of Drosophila guttifera (D. guttifera) are emerging as model systems for studying evolutionary and developmental processes. Studies regarding these processes depend on overexpression and downregulation of developmental genes, which ultimately rely upon an effective transgenic system. Methods describing transgenesis in Drosophila melanogaster (D. melanogaster) have been reported in several studies, but they cannot be applied to D. guttifera due to the low egg production rate and the delicacy of the eggs. In this protocol, we describe extensively a comprehensive method used for generating transgenic D. guttifera. Using the protocol described here, we are able to establish transgenic lines, identifiable by the expression of enhanced green fluorescent protein (EGFP) in the eye disks of D. guttifera larvae. The entire procedure, from injection to screening for transgenic larvae, can be completed in approximately 30 days and should be relatively easy to adapt to other non-model Drosophila species, for which no white-eyed mutants exist.

Evolution of wing pigmentation in Drosophila: Diversity, physiological regulation, and cis-regulatory evolution

Fruit flies (Drosophila and its close relatives, or “drosophilids”) are a group that includes an important model organism, Drosophila melanogaster, and also very diverse species distributed worldwide. Many of these species have black or brown pigmentation patterns on their wings, and have been used as material for evo‐devo research. Pigmentation patterns are thought to have evolved rapidly compared with body plans or body shapes; hence they are advantageous model systems for studying evolutionary gains of traits and parallel evolution. Various groups of drosophilids, including genus Idiomyia (Hawaiian Drosophila), have a variety of pigmentations, ranging from simple black pigmentations around crossveins to a single antero‐distal spot and a more complex mottled pattern. Pigmentation patterns are sometimes obviously used for sexual displays; however, in some cases they may have other functions. The process of wing formation in Drosophila, the general mechanism of pigmentation formation, and the transport of substances necessary for pigmentation, including melanin precursors, through wing veins are summarized here. Lastly, the evolution of the expression of genes regulating pigmentation patterns, the role of cis‐regulatory regions, and the conditions required for the evolutionary emergence of pigmentation patterns are discussed. Future prospects for research on the evolution of wing pigmentation pattern formation in drosophilids are presented, particularly from the point of view of how they compare with other studies of the evolution of new traits.

Multiple Loci Control Eyespot Number Variation on the Hindwings of Bicyclus anynana Butterflies

The underlying genetic changes that regulate the appearance and disappearance of repeated traits, or serial homologs, remain poorly understood. One hypothesis is that variation in genomic regions flanking master regulatory genes, also known as input-output genes, controls variation in trait number, making the locus of evolution almost predictable. Another hypothesis implicates genetic variation in up- or downstream loci of master control genes. Here, we use the butterfly Bicyclus anynana, a species that exhibits natural variation in eyespot number on the dorsal hindwing, to test these two hypotheses. We first estimated the heritability of dorsal hindwing eyespot number by breeding multiple butterfly families differing in eyespot number and regressing eyespot numbers of offspring on midparent values. We then estimated the number and identity of independent genetic loci contributing to eyespot number variation by performing a genome-wide association study with restriction site-associated DNA sequencing from multiple individuals varying in number of eyespots sampled across a freely breeding laboratory population. We found that dorsal hindwing eyespot number has a moderately high heritability of ∼0.50 and is characterized by a polygenic architecture. Previously identified genomic regions involved in eyespot development, and novel ones, display high association with dorsal hindwing eyespot number, suggesting that homolog number variation is likely determined by regulatory changes at multiple loci that build the trait, and not by variation at single master regulators or input-output genes.

Red Light/Green Light, a Dual Fluorescent Protein Reporter System To Study Enhancer-Promoter Specificity in Drosophila

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.

Homeobox Gene Duplication and Divergence in Arachnids

Homeobox genes are key toolkit genes that regulate the development of metazoans and changes in their regulation and copy number have contributed to the evolution of phenotypic diversity. We recently identified a whole genome duplication (WGD) event that occurred in an ancestor of spiders and scorpions (Arachnopulmonata), and that many homeobox genes, including two Hox clusters, appear to have been retained in arachnopulmonates. To better understand the consequences of this ancient WGD and the evolution of arachnid homeobox genes, we have characterized and compared the homeobox repertoires in a range of arachnids. We found that many families and clusters of these genes are duplicated in all studied arachnopulmonates (Parasteatoda tepidariorum, Pholcus phalangioides, Centruroides sculpturatus, and Mesobuthus martensii) compared with nonarachnopulmonate arachnids (Phalangium opilio, Neobisium carcinoides, Hesperochernes sp., and Ixodes scapularis). To assess divergence in the roles of homeobox ohnologs, we analyzed the expression of P. tepidariorum homeobox genes during embryogenesis and found pervasive changes in the level and timing of their expression. Furthermore, we compared the spatial expression of a subset of P. tepidariorum ohnologs with their single copy orthologs in P. opilio embryos. We found evidence for likely subfunctionlization and neofunctionalization of these genes in the spider. Overall our results show a high level of retention of homeobox genes in spiders and scorpions post-WGD, which is likely to have made a major contribution to their developmental evolution and diversification through pervasive subfunctionlization and neofunctionalization, and paralleling the outcomes of WGD in vertebrates.

Accelerated evolution at chaperone promoters among Antarctic notothenioid fishes

BACKGROUND

Antarctic fishes of the Notothenioidei suborder constitutively upregulate multiple inducible chaperones, a highly derived adaptation that preserves proteostasis in extreme cold, and represent a system for studying the evolution of gene frontloading. We screened for Hsf1-binding sites, as Hsf1 is a master transcription factor of the heat shock response, and highly-conserved non-coding elements within proximal promoters of chaperone genes across 10 Antarctic notothens, 2 subpolar notothens, and 17 perciform fishes. We employed phylogenetic models of molecular evolution to determine whether (i) changes in motifs associated with Hsf1-binding and/or (ii) relaxed purifying selection or exaptation at ancestral cis-regulatory elements coincided with the evolution of chaperone frontloading in Antarctic notothens.

RESULTS

Antarctic notothens exhibited significantly fewer Hsf1-binding sites per bp at chaperone promoters than subpolar notothens and Serranoidei, the most closely-related suborder to Notothenioidei included in this study. 90% of chaperone promoters exhibited accelerated substitution rates among Antarctic notothens relative to other perciformes. The proportion of bases undergoing accelerated evolution (i) was significantly greater in Antarctic notothens than in subpolar notothens and Perciformes in 70% of chaperone genes and (ii) increased among bases that were more conserved among perciformes. Lastly, we detected evidence of relaxed purifying selection and exaptation acting on ancestrally conserved cis-regulatory elements in the Antarctic notothen lineage and its major branches.

CONCLUSION

A large degree of turnover has occurred in Notothenioidei at chaperone promoter regions that are conserved among perciform fishes following adaptation to the cooling of the Southern Ocean. Additionally, derived reductions in Hsf1-binding site frequency suggest cis-regulatory modifications to the classical heat shock response. Of note, turnover events within chaperone promoters were less frequent in the ancestral node of Antarctic notothens relative to younger Antarctic lineages. This suggests that cis-regulatory divergence at chaperone promoters may be greater between Antarctic notothen lineages than between subpolar and Antarctic clades. These findings demonstrate that strong selective forces have acted upon cis-regulatory elements of chaperone genes among Antarctic notothens.

Comparative analysis of integument transcriptomes identifies genes that participate in marking pattern formation in three allelic mutants of silkworm, Bombyx mori

The diversity markings and pigment patterns in insects are outcomes of adaptive evolution. The elucidation of the molecular mechanism underlying variations in pigment patterns may improve our understanding of the origin and evolution of these spectacular diverse phenotypes. Melanin, ommochrome, and pteridine are the three main types of insect pigments, and the genes that directly participate in pigment biosynthesis have been extensively studied. However, available information on gene interactions and the whole pigment regulatory network is limited. In this study, we performed integument transcriptome sequencing to analyze three larval marking allelic mutants, namely, multi lunar (L), L^C, and L^Ca, which have similar twin-spot markings on the dorsal side of multiple segments. Further analysis identified 336 differentially expressed genes (DEGs) between L and Dazao (wild type which exhibits normal markings), 68 DEGs between L^C/+ and +^LC/+^LC, and 188 DEGs between L^Ca/+ and +^LCa/+^LCa. Gene Ontology (GO) analysis indicated a significant DEG enrichment of the functional terms catalytic activity, binding, metabolic process, and cellular process. Furthermore, three mutants share six common enriched KEGG pathways. We finally identified eight common DEGs among three pairwise comparisons, including Krueppel-like factor, TATA-binding protein, protein patched, UDP-glycosyltransferase, an unknown secreted protein, and three cuticular proteins. Microarray-based gene expression analysis revealed that the eight genes are upregulated during molting, which coincides with marking formation, and are significantly differentially expressed between marking and non-marking regions. The results suggest that the eight common genes are involved in the construction of the multiple twin-spot marking patterns in the three mutants.

Modular cis-regulatory logic of yellow gene expression in silkmoth larvae

Colour patterns in butterflies and moths are crucial traits for adaptation. Previous investigations have highlighted genes responsible for pigmentation (ie yellow and ebony). However, the mechanisms by which these genes are regulated in lepidopteran insects remain poorly understood. To elucidate this, molecular studies involving dipterans have largely analysed the cis-regulatory regions of pigmentation genes and have revealed cis-regulatory modularity. Here, we used well-developed transgenic techniques in Bombyx mori and demonstrated that cis-regulatory modularity controls tissue-specific expression of the yellow gene. We first identified which body parts are regulated by the yellow gene via black pigmentation. We then isolated three discrete regulatory elements driving tissue-specific gene expression in three regions of B. mori larvae. Finally, we found that there is no apparent sequence conservation of cis-regulatory regions between B. mori and Drosophila melanogaster, and no expression driven by the regulatory regions of one species when introduced into the other species. Therefore, the trans-regulatory landscapes of the yellow gene differ significantly between the two taxa. The results of this study confirm that lepidopteran species use cis-regulatory modules to control gene expression related to pigmentation, and represent a powerful cadre of transgenic tools for studying evolutionary developmental mechanisms.

The Divergent Genomes of Teleosts

Boasting nearly 30,000 species, teleosts account for half of all extant vertebrates and approximately 98% of all ray-finned fish species (Actinopterygii). Teleosts are also the largest and most diverse group of vertebrates, exhibiting an astonishing level of morphological, physiological, and behavioral diversity. Previous studies had indicated that the teleost lineage has experienced an additional whole-genome duplication event. Recent comparative genomic analyses of teleosts and other bony vertebrates using spotted gar (a nonteleost ray-finned fish) and elephant shark (a cartilaginous fish) as outgroups have revealed several divergent features of teleost genomes. These include an accelerated evolutionary rate of protein-coding and nucleotide sequences, a higher rate of intron turnover, loss of many potential cis-regulatory elements and shorter conserved syntenic blocks. A combination of these divergent genomic features might have contributed to the evolution of the amazing phenotypic diversity and morphological innovations of teleosts.

Genome-Wide Identification and Expression Profiling of Wnt Family Genes in the Silkworm, Bombyx mori

Wnt is a family of conserved glycoproteins that participate in a variety of important biological processes including embryo development, cell proliferation and differentiation, and tissue regeneration. The Wnt family is a metazoan novelty found in all animal phyla. Studies have revealed that the number of Wnt genes varies among species, presumably due to reproduction and loss of genes during evolution. However, a comprehensive inventory of Wnt genes in Lepidoptera is lacking. In this study, we identified the repertoire of Wnt genes in the silkworm and seven other species of Lepidoptera and obtained eight Wnt genes (Wnt1, Wnt5–Wnt7, Wnt9–Wnt11, and WntA) in each species. Four of these Wnt genes are clustered in two orientations (5′-Wnt9-Wnt1-Wnt6-Wnt10-3′ and 5′-Wnt10-Wnt6-Wnt1-Wnt9-3′) in both moths and butterflies. Transcript analysis of Wnt in silkworm embryonic stages showed that each BmWnt gene had a unique expression pattern during embryological development. Analysis of a larval stage revealed differential expression of Wnt family members in diverse tissues. Our study provides an overview of the Wnt family in Lepidoptera and will inspire further functional study of the Wnt genes in the silkworm.

Redundant and Cryptic Enhancer Activities of the Drosophila yellow Gene

Cis-regulatory sequences known as enhancers play a key role in regulating gene expression. Evolutionary changes in these DNA sequences contribute to phenotypic evolution. The Drosophila yellow gene, which is required for pigmentation, has emerged as a model system for understanding how cis-regulatory sequences evolve, providing some of the most detailed insights available into how activities of orthologous enhancers have diverged between species. Here, we examine the evolution of yellow cis-regulatory sequences on a broader scale, by comparing the distribution and function of yellow enhancer activities throughout the 5' intergenic and intronic sequences of Drosophila melanogaster, D. pseudoobscura, and D. willistoni We find that cis-regulatory sequences driving expression in a particular tissue are not as modular as previously described, but rather have many redundant and cryptic enhancer activities distributed throughout the regions surveyed. Interestingly, cryptic enhancer activities of sequences from one species often drove patterns of expression observed in other species, suggesting that the frequent evolutionary changes in yellow expression observed among Drosophila species may be facilitated by gaining and losing repression of preexisting cis-regulatory sequences.

Regulatory Diversification of INDEHISCENT in the Capsella Genus Directs Variation in Fruit Morphology

Evolution of gene-regulatory sequences is considered the primary driver of morphological variation [1-3]. In animals, the diversity of body plans between distantly related phyla is due to the differential expression patterns of conserved "toolkit" genes [4]. In plants, variation in expression domains similarly underlie most of the reported diversity of organ shape both in natural evolution and in the domestication of crops [5-9]. The heart-shaped fruit from members of the Capsella genus is a morphological novelty that has evolved after Capsella diverged from Arabidopsis ∼8 mya [10]. Comparative studies of fruit growth in Capsella and Arabidopsis revealed that the difference in shape is caused by local control of anisotropic growth [11]. Here, we show that sequence variation in regulatory domains of the fruit-tissue identity gene, INDEHISCENT (IND), is responsible for expansion of its expression domain in the heart-shaped fruits from Capsella rubella. We demonstrate that expression of this CrIND gene in the apical part of the valves in Capsella contributes to the heart-shaped appearance. While studies on morphological diversity have revealed the importance of cis-regulatory sequence evolution, few examples exist where the downstream effects of such variation have been characterized in detail. We describe here how CrIND exerts its function on Capsella fruit shape by binding sequence elements of auxin biosynthesis genes to activate their expression and ensure auxin accumulation into highly localized maxima in the fruit valves. Thus, our data provide a direct link between changes in expression pattern and altered hormone homeostasis in the evolution of morphological novelty.

Assessing the use of wing ornamentation and visual display in female choice sexual selection

Conspicuous sexual dimorphism is often ascribed to sexual selection. When the differences between the sexes are ornamental, this is thought to indicate a role for female choice. In spotted winged Drosophila species courtship, a male positioned in front of a female waves his wings, which have a patch of melanization on the exterior margin. In this study, we examine both female preference for wing spots and the role of vision in mating success in three species of the suzukii group: Drosophila biarmipes, D. suzukii, and D. subpulchrella. To assess female preference for wing spot, we removed the spot with a novel, non-invasive method, and competed spotless males with males with two spots. Phenotype did not affect mating success in any species. To eliminate the potential effect of competitive behavior on male mating success, we also ran a no-choice analysis. Mating frequency and timing was not different between phenotypes within these species. The effect of vision on mating success was assessed by comparing mating success of spotted males between light and dark conditions, both for frequency of mating, as well as timing of multiple courtship parameters. Species varied in the extent that lack of vision negatively affected mating success. Though vision is important for mating success, the spot itself may not be providing the signal that females use to make mating decisions.

Evolutionary emergence of the rac3b/rfng/sgca regulatory cluster refined mechanisms for hindbrain boundaries formation

Developmental programs often rely on parallel morphogenetic mechanisms that guarantee precise tissue architecture. While redundancy constitutes an obvious selective advantage, little is known on how novel morphogenetic mechanisms emerge during evolution. In zebrafish, rhombomeric boundaries behave as an elastic barrier, preventing cell intermingling between adjacent compartments. Here, we identify the fundamental role of the small-GTPase Rac3b in actomyosin cable assembly at hindbrain boundaries. We show that the novel rac3b/rfng/sgca regulatory cluster, which is specifically expressed at the boundaries, emerged in the Ostariophysi superorder by chromosomal rearrangement that generated new cis-regulatory interactions. By combining 4C-seq, ATAC-seq, transgenesis, and CRISPR-induced deletions, we characterized this regulatory domain, identifying hindbrain boundary-specific cis-regulatory elements. Our results suggest that the capacity of boundaries to act as an elastic mesh for segregating rhombomeric cells evolved by cooption of critical genes to a novel regulatory block, refining the mechanisms for hindbrain segmentation.

Pigmentation pattern and developmental constraints: flight muscle attachment sites delimit the thoracic trident of Drosophila melanogaster

In their seminal paper published in 1979, Gould and Lewontin argued that some traits arise as by-products of the development of other structures and not for direct utility in themselves. We show here that this applies to the trident, a pigmentation pattern observed on the thorax of Drosophila melanogaster. Using reporter constructs, we show that the expression domain of several genes encoding pigmentation enzymes follows the trident shape. This domain is complementary to the expression pattern of stripe (sr), which encodes an essential transcription factor specifying flight muscle attachment sites. We demonstrate that sr limits the expression of these pigmentation enzyme genes to the trident by repressing them in its own expression domain, i.e. at the flight muscle attachment sites. We give evidence that repression of not only yellow but also other pigmentation genes, notably tan, is involved in the trident shape. The flight muscle attachment sites and sr expression patterns are remarkably conserved in dipterans reflecting the essential role of sr. Our data suggest that the trident is a by-product of flight muscle attachment site patterning that arose when sr was co-opted for the regulation of pigmentation enzyme coding genes.

Methods for Staging Pupal Periods and Measurement of Wing Pigmentation of Drosophila guttifera

Diversified species of Drosophila (fruit fly) provide opportunities to study mechanisms of development and genetic changes responsible for evolutionary changes. In particular, the adult stage is a rich source of morphological traits for interspecific comparison, including wing pigmentation comparison. To study developmental differences among species, detailed observation and appropriate staging are required for precise comparison. Here we describe protocols for staging of pupal periods and quantification of wing pigmentation in a polka-dotted fruit fly, Drosophila guttifera. First, we describe the method for detailed morphological observation and definition of pupal stages based on morphologies. This method includes a technique for removing the puparium, which is the outer chitinous case of the pupa, to enable detailed observation of pupal morphologies. Second, we describe the method for measuring the duration of defined pupal stages. Finally, we describe the method for quantification of wing pigmentation based on image analysis using digital images and ImageJ software. With these methods, we can establish a solid basis for comparing developmental processes of adult traits during pupal stages.

Functional characterization of adaptive variation within a cis-regulatory element influencing Drosophila melanogaster growth

Gene expression variation is a major contributor to phenotypic diversity within species and is thought to play an important role in adaptation. However, examples of adaptive regulatory polymorphism are rare, especially those that have been characterized at both the molecular genetic level and the organismal level. In this study, we perform a functional analysis of the Drosophila melanogaster CG9509 enhancer, a cis-regulatory element that shows evidence of adaptive evolution in populations outside the species’ ancestral range in sub-Saharan Africa. Using site-directed mutagenesis and transgenic reporter gene assays, we determined that 3 single nucleotide polymorphisms are responsible for the difference in CG9509 expression that is observed between sub-Saharan African and cosmopolitan populations. Interestingly, while 2 of these variants appear to have been the targets of a selective sweep outside of sub-Saharan Africa, the variant with the largest effect on expression remains polymorphic in cosmopolitan populations, suggesting it may be subject to a different mode of selection. To elucidate the function of CG9509, we performed a series of functional and tolerance assays on flies in which CG9509 expression was disrupted. We found that CG9509 plays a role in larval growth and influences adult body and wing size, as well as wing loading. Furthermore, variation in several of these traits was associated with variation within the CG9509 enhancer. The effect on growth appears to result from a modulation of active ecdysone levels and expression of growth factors. Taken together, our findings suggest that selection acted on 3 sites within the CG9509 enhancer to increase CG9509 expression and, as a result, reduce wing loading as D. melanogaster expanded out of sub-Saharan Africa.

Much of the phenotypic variation that is observed within species is thought to be caused by variation in gene expression. Variants within cis-regulatory elements, which affect the expression of nearby genes within the same DNA strand, are thought to be an abundant resource upon which natural selection can act. Understanding the functional consequences of adaptive cis-regulatory changes is important, as it can help elucidate the mechanisms underlying phenotypic evolution in general and provide insight into the development and maintenance of biodiversity. However, functional analyses of these types of changes remain rare. Here we present a functional analysis of an adaptively evolving enhancer element of a D. melanogaster gene called CG9509, of previously unknown function. We show that 3 single nucleotide polymorphisms located within the enhancer of this gene are responsible for an increase in CG9509 expression in cosmopolitan populations (outside of south and central Africa) relative to sub-Saharan populations, which include ancestral populations. We further show that CG9509 is involved in the regulation of growth rate and body size determination and propose that the CG9509 enhancer underwent positive selection to reduce wing loading as the species expanded out of sub-Saharan Africa.

Waiting in the wings: what can we learn about gene co-option from the diversification of butterfly wing patterns?

A major challenge is to understand how conserved gene regulatory networks control the wonderful diversity of form that we see among animals and plants. Butterfly wing patterns are an excellent example of this diversity. Butterfly wings form as imaginal discs in the caterpillar and are constructed by a gene regulatory network, much of which is conserved across the holometabolous insects. Recent work in Heliconius butterflies takes advantage of genomic approaches and offers insights into how the diversification of wing patterns is overlaid onto this conserved network. WntA is a patterning morphogen that alters spatial information in the wing. Optix is a transcription factor that acts later in development to paint specific wing regions red. Both of these loci fit the paradigm of conserved protein-coding loci with diverse regulatory elements and developmental roles that have taken on novel derived functions in patterning wings. These discoveries offer insights into the 'Nymphalid Ground Plan', which offers a unifying hypothesis for pattern formation across nymphalid butterflies. These loci also represent 'hotspots' for morphological change that have been targeted repeatedly during evolution. Both convergent and divergent evolution of a great diversity of patterns is controlled by complex alleles at just a few genes. We suggest that evolutionary change has become focused on one or a few genetic loci for two reasons. First, pre-existing complex cis-regulatory loci that already interact with potentially relevant transcription factors are more likely to acquire novel functions in wing patterning. Second, the shape of wing regulatory networks may constrain evolutionary change to one or a few loci. Overall, genomic approaches that have identified wing patterning loci in these butterflies offer broad insight into how gene regulatory networks evolve to produce diversity.This article is part of the themed issue 'Evo-devo in the genomics era, and the origins of morphological diversity'.

Macroevolutionary shifts of WntA function potentiate butterfly wing-pattern diversity

Butterfly wing patterns provide a rich comparative framework to study how morphological complexity develops and evolves. Here we used CRISPR/Cas9 somatic mutagenesis to test a patterning role for WntA, a signaling ligand gene previously identified as a hotspot of shape-tuning alleles involved in wing mimicry. We show that WntA loss-of-function causes multiple modifications of pattern elements in seven nymphalid butterfly species. In three butterflies with a conserved wing-pattern arrangement, WntA is necessary for the induction of stripe-like patterns known as symmetry systems and acquired a novel eyespot activator role specific to Vanessa forewings. In two Heliconius species, WntA specifies the boundaries between melanic fields and the light-color patterns that they contour. In the passionvine butterfly Agraulis, WntA removal shows opposite effects on adjacent pattern elements, revealing a dual role across the wing field. Finally, WntA acquired a divergent role in the patterning of interveinous patterns in the monarch, a basal nymphalid butterfly that lacks stripe-like symmetry systems. These results identify WntA as an instructive signal for the prepatterning of a biological system of exuberant diversity and illustrate how shifts in the deployment and effects of a single developmental gene underlie morphological change.

Automated tools for comparative sequence analysis of genic regions using the GenePalette application

Comparative sequence analysis methods, such as phylogenetic footprinting, represent one of the most effective ways to decode regulatory sequence functions based upon DNA sequence information alone. The laborious task of assembling orthologous sequences to perform these comparisons is a hurdle to these analyses, which is further aggravated by the relative paucity of tools for visualization of sequence comparisons in large genic regions. Here, we describe a second-generation implementation of the GenePalette DNA sequence analysis software to facilitate comparative studies of gene function and regulation. We have developed an automated module called OrthologGrabber (OG) that performs BLAT searches against the UC Santa Cruz genome database to identify and retrieve segments homologous to a region of interest. Upon acquisition, sequences are compared to identify high-confidence anchor-points, which are graphically displayed. The visualization of anchor-points alongside other DNA features, such as transcription factor binding sites, allows users to precisely examine whether a binding site of interest is conserved, even if the surrounding region exhibits poor sequence identity. This approach also aids in identifying orthologous segments of regulatory DNA, facilitating studies of regulatory sequence evolution. As with previous versions of the software, GenePalette 2.1 takes the form of a platform-independent, single-windowed interface that is simple to use.