Repeated morphological evolution through cis-regulatory changes in a pleiotropic gene

Co-evolving wing spots and mating displays are genetically separable traits in Drosophila

The evolution of sexual traits often involves correlated changes in morphology and behavior. For example, in Drosophila, divergent mating displays are often accompanied by divergent pigment patterns. To better understand how such traits co-evolve, we investigated the genetic basis of correlated divergence in wing pigmentation and mating display between the sibling species Drosophila elegans and Drosophila gunungcola. Drosophila elegans males have an area of black pigment on their wings known as a wing spot and appear to display this spot to females by extending their wings laterally during courtship. By contrast, D. gunungcola lost both of these traits. Using Multiplexed Shotgun Genotyping (MSG), we identified a ∼440 kb region on the X chromosome that behaves like a genetic switch controlling the presence or absence of male-specific wing spots. This region includes the candidate gene optomotor-blind (omb), which plays a critical role in patterning the Drosophila wing. The genetic basis of divergent wing display is more complex, with at least two loci on the X chromosome and two loci on autosomes contributing to its evolution. Introgressing the X-linked region affecting wing spot development from D. gunungcola into D. elegans reduced pigmentation in the wing spots but did not affect the wing display, indicating that these are genetically separable traits. Consistent with this observation, broader sampling of wild D. gunungcola populations confirmed that the wing spot and wing display are evolving independently: some D. gunungcola males performed wing displays similar to D. elegans despite lacking wing spots. These data suggest that correlated selection pressures rather than physical linkage or pleiotropy are responsible for the coevolution of these morphological and behavioral traits. They also suggest that the change in morphology evolved prior to the change in behavior.

Identifying branch-specific positive selection throughout the regulatory genome using an appropriate proxy neutral

BACKGROUND

Adaptive changes in cis-regulatory elements are an essential component of evolution by natural selection. Identifying adaptive and functional noncoding DNA elements throughout the genome is therefore crucial for understanding the relationship between phenotype and genotype.

RESULTS

We used ENCODE annotations to identify appropriate proxy neutral sequences and demonstrate that the conservativeness of the test can be modulated during the filtration of reference alignments. We applied the method to noncoding Human Accelerated Elements as well as open chromatin elements previously identified in 125 human tissues and cell lines to demonstrate its utility. Then, we evaluated the impact of query region length, proxy neutral sequence length, and branch count on test sensitivity and specificity. We found that the length of the query alignment can vary between 150 bp and 1 kb without affecting the estimation of selection, while for the reference alignment, we found that a length of 3 kb is adequate for proper testing. We also simulated sequence alignments under different classes of evolution and validated our ability to distinguish positive selection from relaxation of constraint and neutral evolution. Finally, we re-confirmed that a quarter of all non-coding Human Accelerated Elements are evolving by positive selection.

CONCLUSION

Here, we introduce a method we called adaptiPhy, which adds significant improvements to our earlier method that tests for branch-specific directional selection in noncoding sequences. The motivation for these improvements is to provide a more sensitive and better targeted characterization of directional selection and neutral evolution across the genome.

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.

Multicomponent structures in camouflage and mimicry in butterfly wing patterns

Understanding how morphological structures are built is essential for appreciating the morphological complexity and divergence of organisms. One representative case of morphological structures is the camouflage and mimicry of butterfly wing patterns. Some previous studies have questioned whether camouflage and mimicry are truly structures, considering that they rely on coloration. Nevertheless, our recent study revealed that the leaf pattern of Kallima inachus butterfly wings evolved through the combination of changes in several pigment components in a block-wise manner; it remains unclear whether such block-wise structures are common in other cases of camouflage and mimicry in butterflies and how they come about. Previous studies focused solely on a set of homologous components, termed the nymphalid ground plan. In the present study, we extended the scope of the description of components by including not only the nymphalid ground plan but also other common components (i.e., ripple patterns, dependent patterns, and color fields). This extension allowed us to analyze the combinatorial building logic of structures and examine multicomponent structures of camouflage and mimicry in butterfly wing patterns. We investigated various patterns of camouflage and mimicry (e.g., masquerade, crypsis, Müllerian mimicry, Batesian mimicry) in nine species and decomposed them into an assembly of multiple components. These structural component analyses suggested that camouflage and mimicry in butterfly wing patterns are built up by combining multiple types of components. We also investigated associations between components and the kinds of camouflage and mimicry. Several components are statistically more often used to produce specific types of camouflage or mimicry. Thus, our work provides empirical evidence that camouflage and mimicry patterns of butterfly wings are mosaic structures, opening up a new avenue of studying camouflage, and mimicry from a structural perspective.

Characterization of the Genetic Architecture Underlying Eye Size Variation Within Drosophila melanogaster and Drosophila simulans

The compound eyes of insects exhibit striking variation in size, reflecting adaptation to different lifestyles and habitats. However, the genetic and developmental bases of variation in insect eye size is poorly understood, which limits our understanding of how these important morphological differences evolve. To address this, we further explored natural variation in eye size within and between four species of the Drosophila melanogaster species subgroup. We found extensive variation in eye size among these species, and flies with larger eyes generally had a shorter inter-ocular distance and vice versa We then carried out quantitative trait loci (QTL) mapping of intra-specific variation in eye size and inter-ocular distance in both D. melanogaster and D. simulans This revealed that different genomic regions underlie variation in eye size and inter-ocular distance in both species, which we corroborated by introgression mapping in D. simulans This suggests that although there is a trade-off between eye size and inter-ocular distance, variation in these two traits is likely to be caused by different genes and so can be genetically decoupled. Finally, although we detected QTL for intra-specific variation in eye size at similar positions in D. melanogaster and D. simulans, we observed differences in eye fate commitment between strains of these two species. This indicates that different developmental mechanisms and therefore, most likely, different genes contribute to eye size variation in these species. Taken together with the results of previous studies, our findings suggest that the gene regulatory network that specifies eye size has evolved at multiple genetic nodes to give rise to natural variation in this trait within and among species.

The Role of Chromatin Accessibility in cis-Regulatory Evolution

Transcription factor (TF) binding is determined by sequence as well as chromatin accessibility. Although the role of accessibility in shaping TF-binding landscapes is well recorded, its role in evolutionary divergence of TF binding, which in turn can alter cis-regulatory activities, is not well understood. In this work, we studied the evolution of genome-wide binding landscapes of five major TFs in the core network of mesoderm specification, between Drosophila melanogaster and Drosophila virilis, and examined its relationship to accessibility and sequence-level changes. We generated chromatin accessibility data from three important stages of embryogenesis in both Drosophila melanogaster and Drosophila virilis and recorded conservation and divergence patterns. We then used multivariable models to correlate accessibility and sequence changes to TF-binding divergence. We found that accessibility changes can in some cases, for example, for the master regulator Twist and for earlier developmental stages, more accurately predict binding change than is possible using TF-binding motif changes between orthologous enhancers. Accessibility changes also explain a significant portion of the codivergence of TF pairs. We noted that accessibility and motif changes offer complementary views of the evolution of TF binding and developed a combined model that captures the evolutionary data much more accurately than either view alone. Finally, we trained machine learning models to predict enhancer activity from TF binding and used these functional models to argue that motif and accessibility-based predictors of TF-binding change can substitute for experimentally measured binding change, for the purpose of predicting evolutionary changes in enhancer activity.

Evolution of Reproductive Behavior

Behaviors associated with reproduction are major contributors to the evolutionary success of organisms and are subject to many evolutionary forces, including natural and sexual selection, and sexual conflict. Successful reproduction involves a range of behaviors, from finding an appropriate mate, courting, and copulation, to the successful production and (in oviparous animals) deposition of eggs following mating. As a consequence, behaviors and genes associated with reproduction are often under strong selection and evolve rapidly. Courtship rituals in flies follow a multimodal pattern, mediated through visual, chemical, tactile, and auditory signals. Premating behaviors allow males and females to assess the species identity, reproductive state, and condition of their partners. Conflicts between the "interests" of individual males, and/or between the reproductive strategies of males and females, often drive the evolution of reproductive behaviors. For example, seminal proteins transmitted by males often show evidence of rapid evolution, mediated by positive selection. Postmating behaviors, including the selection of oviposition sites, are highly variable and Drosophila species span the spectrum from generalists to obligate specialists. Chemical recognition features prominently in adaptation to host plants for feeding and oviposition. Selection acting on variation in pre-, peri-, and postmating behaviors can lead to reproductive isolation and incipient speciation. Response to selection at the genetic level can include the expansion of gene families, such as those for detecting pheromonal cues for mating, or changes in the expression of genes leading to visual cues such as wing spots that are assessed during mating. Here, we consider the evolution of reproductive behavior in Drosophila at two distinct, yet complementary, scales. Some studies take a microevolutionary approach, identifying genes and networks involved in reproduction, and then dissecting the genetics underlying complex behaviors in D. melanogaster Other studies take a macroevolutionary approach, comparing reproductive behaviors across the genus Drosophila and how these might correlate with environmental cues. A full synthesis of this field will require unification across these levels.

Integration of Anatomy Ontologies and Evo-Devo Using Structured Markov Models Suggests a New Framework for Modeling Discrete Phenotypic Traits

Modeling discrete phenotypic traits for either ancestral character state reconstruction or morphology-based phylogenetic inference suffers from ambiguities of character coding, homology assessment, dependencies, and selection of adequate models. These drawbacks occur because trait evolution is driven by two key processes-hierarchical and hidden-which are not accommodated simultaneously by the available phylogenetic methods. The hierarchical process refers to the dependencies between anatomical body parts, while the hidden process refers to the evolution of gene regulatory networks (GRNs) underlying trait development. Herein, I demonstrate that these processes can be efficiently modeled using structured Markov models (SMM) equipped with hidden states, which resolves the majority of the problems associated with discrete traits. Integration of SMM with anatomy ontologies can adequately incorporate the hierarchical dependencies, while the use of the hidden states accommodates hidden evolution of GRNs and substitution rate heterogeneity. I assess the new models using simulations and theoretical synthesis. The new approach solves the long-standing "tail color problem," in which the trait is scored for species with tails of different colors or no tails. It also presents a previously unknown issue called the "two-scientist paradox," in which the nature of coding the trait and the hidden processes driving the trait's evolution are confounded; failing to account for the hidden process may result in a bias, which can be avoided by using hidden state models. All this provides a clear guideline for coding traits into characters. This article gives practical examples of using the new framework for phylogenetic inference and comparative analysis.

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 genetic architecture of adaptation: convergence and pleiotropy in Heliconius wing pattern evolution

Unravelling the genetic basis of adaptive traits is a major challenge in evolutionary biology. Doing so informs our understanding of evolution towards an adaptive optimum, the distribution of locus effect sizes, and the influence of genetic architecture on the evolvability of a trait. In the Müllerian co-mimics Heliconius melpomene and Heliconius erato some Mendelian loci affecting mimicry shifts are well known. However, several phenotypes in H. melpomene remain to be mapped, and the quantitative genetics of colour pattern variation has rarely been analysed. Here we use quantitative trait loci (QTL) analyses of crosses between H. melpomene races from Peru and Suriname to map, for the first time, the control of the broken band phenotype to WntA and identify a ~100 kb region controlling this variation. Additionally, we map variation in basal forewing red-orange pigmentation to a locus centred around the gene ventral veins lacking (vvl). The locus also appears to affect medial band shape variation as it was previously known to do in H. erato. This adds to the list of homologous regions controlling convergent phenotypes between these two species. Finally we show that Heliconius wing-patterning genes are strikingly pleiotropic among wing pattern traits. Our results demonstrate how genetic architecture can shape, aid and constrain adaptive evolution.

Independent erosion of conserved transcription factor binding sites points to shared hindlimb, vision and external testes loss in different mammals

Genetic variation in cis-regulatory elements is thought to be a major driving force in morphological and physiological changes. However, identifying transcription factor binding events that code for complex traits remains a challenge, motivating novel means of detecting putatively important binding events. Using a curated set of 1154 high-quality transcription factor motifs, we demonstrate that independently eroded binding sites are enriched for independently lost traits in three distinct pairs of placental mammals. We show that these independently eroded events pinpoint the loss of hindlimbs in dolphin and manatee, degradation of vision in naked mole-rat and star-nosed mole, and the loss of external testes in white rhinoceros and Weddell seal. We additionally show that our method may also be utilized with more than two species. Our study exhibits a novel methodology to detect cis-regulatory mutations which help explain a portion of the molecular mechanism underlying complex trait formation and loss.

Repeated Cis-Regulatory Tuning of a Metabolic Bottleneck Gene during Evolution

Repeated evolutionary events imply underlying genetic constraints that can make evolutionary mechanisms predictable. Morphological traits are thought to evolve frequently through cis-regulatory changes because these mechanisms bypass constraints in pleiotropic genes that are reused during development. In contrast, the constraints acting on metabolic traits during evolution are less well studied. Here we show how a metabolic bottleneck gene has repeatedly adopted similar cis-regulatory solutions during evolution, likely due to its pleiotropic role integrating flux from multiple metabolic pathways. Specifically, the genes encoding phosphoglucomutase activity (PGM1/PGM2), which connect GALactose catabolism to glycolysis, have gained and lost direct regulation by the transcription factor Gal4 several times during yeast evolution. Through targeted mutations of predicted Gal4-binding sites in yeast genomes, we show this galactose-mediated regulation of PGM1/2 supports vigorous growth on galactose in multiple yeast species, including Saccharomyces uvarum and Lachancea kluyveri. Furthermore, the addition of galactose-inducible PGM1 alone is sufficient to improve the growth on galactose of multiple species that lack this regulation, including Saccharomyces cerevisiae. The strong association between regulation of PGM1/2 by Gal4 even enables remarkably accurate predictions of galactose growth phenotypes between closely related species. This repeated mode of evolution suggests that this specific cis-regulatory connection is a common way that diverse yeasts can govern flux through the pathway, likely due to the constraints imposed by this pleiotropic bottleneck gene. Since metabolic pathways are highly interconnected, we argue that cis-regulatory evolution might be widespread at pleiotropic genes that control metabolic bottlenecks and intersections.

Quantifying the extent of morphological homoplasy: A phylogenetic analysis of 490 characters in Drosophila

Homoplasy is a fundamental phenomenon in evolutionary biology but an appraisal of its extent at the morphological level is still lacking. Here, we analyzed the evolution of 490 morphological characters conceptualized among 56 drosophilid species. We found that two thirds of morphological changes were homoplastic and that the level of homoplasy depended on the stage of development and the type of the organ, with the adult terminalia being the least homoplastic. In spite of its predominance at the character change level, homoplasy accounts for only ∼13% of between species similarities in pairwise comparisons. These results provide empirical insights on the limits of morphological changes and the frequency of recurrent evolution.

Development and evolution of color patterns in ladybird beetles: A case study in Harmonia axyridis

Many organisms show various geometric color patterns on their bodies, and the developmental, evolutionary, genetic, and ecological bases of these patterns have been intensely studied in various organisms. Ladybird beetles display highly diverse patterns of wing (elytral) color and are one of the most attractive model organisms for studying these characteristics. In this study, we reviewed the genetic history of elytral color patterns in the Asian multicolored ladybird beetle Harmonia axyridis from the classical genetic studies led by the pupils of Thomas Hunt Morgan and Theodosius Dobzhansky to recent genomic studies that revealed that a single GATA transcription factor gene, pannier, regulates the highly diverse elytral color patterns in this species. We also reviewed and discussed the developmental and evolutionary mechanisms driven by the pannier locus in H. axyridis. In the development sections, we focused on the following two topics: (a) how the red (carotenoid) and black (melanin) pigmentation of elytra is regulated by the pannier and pigmentation genes and (b) how the diverse color patterns are formed by integrating regulatory inputs from other genes involved in wing development. In the evolution section, we subsequently focused on the highly diversified DNA sequences within the first intron of pannier that are 56-76 kb long and that were generated through recurrent multiple inversions. Furthermore, we discussed how these recurrent inversions have driven the diversification of color patterns throughout evolution.

Inverse resource allocation between vision and olfaction across the genus Drosophila

Divergent populations across different environments are exposed to critical sensory information related to locating a host or mate, as well as avoiding predators and pathogens. These sensory signals generate evolutionary changes in neuroanatomy and behavior; however, few studies have investigated patterns of neural architecture that occur between sensory systems, or that occur within large groups of closely-related organisms. Here we examine 62 species within the genus Drosophila and describe an inverse resource allocation between vision and olfaction, which we consistently observe at the periphery, within the brain, as well as during larval development. This sensory variation was noted across the entire genus and appears to represent repeated, independent evolutionary events, where one sensory modality is consistently selected for at the expense of the other. Moreover, we provide evidence of a developmental genetic constraint through the sharing of a single larval structure, the eye-antennal imaginal disc. In addition, we examine the ecological implications of visual or olfactory bias, including the potential impact on host-navigation and courtship.

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.

Evolution of DNA Methylation in Papio Baboons

Changes in gene regulation have long been thought to play an important role in primate evolution. However, although a number of studies have compared genome-wide gene expression patterns across primate species, fewer have investigated the gene regulatory mechanisms that underlie such patterns, or the relative contribution of drift versus selection. Here, we profiled genome-scale DNA methylation levels in blood samples from five of the six extant species of the baboon genus Papio (4-14 individuals per species). This radiation presents the opportunity to investigate DNA methylation divergence at both shallow and deeper timescales (0.380-1.4 My). In contrast to studies in human populations, but similar to studies in great apes, DNA methylation profiles clearly mirror genetic and geographic structure. Divergence in DNA methylation proceeds fastest in unannotated regions of the genome and slowest in regions of the genome that are likely more constrained at the sequence level (e.g., gene exons). Both heuristic approaches and Ornstein-Uhlenbeck models suggest that DNA methylation levels at a small set of sites have been affected by positive selection, and that this class is enriched in functionally relevant contexts, including promoters, enhancers, and CpG islands. Our results thus indicate that the rate and distribution of DNA methylation changes across the genome largely mirror genetic structure. However, at some CpG sites, DNA methylation levels themselves may have been a target of positive selection, pointing to loci that could be important in connecting sequence variation to fitness-related traits.

DNA fragility in the parallel evolution of pelvic reduction in stickleback fish

Evolution generates a remarkable breadth of living forms, but many traits evolve repeatedly, by mechanisms that are still poorly understood. A classic example of repeated evolution is the loss of pelvic hindfins in stickleback fish (Gasterosteus aculeatus). Repeated pelvic loss maps to recurrent deletions of a pelvic enhancer of the Pitx1 gene. Here, we identify molecular features contributing to these recurrent deletions. Pitx1 enhancer sequences form alternative DNA structures in vitro and increase double-strand breaks and deletions in vivo. Enhancer mutability depends on DNA replication direction and is caused by TG-dinucleotide repeats. Modeling shows that elevated mutation rates can influence evolution under demographic conditions relevant for sticklebacks and humans. DNA fragility may thus help explain why the same loci are often used repeatedly during parallel adaptive evolution.

Diversification of DNA binding specificities enabled SREBP transcription regulators to expand the repertoire of cellular functions that they govern in fungi

The Sterol Regulatory Element Binding Proteins (SREBPs) are basic-helix-loop-helix transcription regulators that control the expression of sterol biosynthesis genes in higher eukaryotes and some fungi. Surprisingly, SREBPs do not regulate sterol biosynthesis in the ascomycete yeasts (Saccharomycotina) as this role was handed off to an unrelated transcription regulator in this clade. The SREBPs, nonetheless, expanded in fungi such as the ascomycete yeasts Candida spp., raising questions about their role and evolution in these organisms. Here we report that the fungal SREBPs diversified their DNA binding preferences concomitantly with an expansion in function. We establish that several branches of fungal SREBPs preferentially bind non-palindromic DNA sequences, in contrast to the palindromic DNA motifs recognized by most basic-helix-loop-helix proteins (including SREBPs) in higher eukaryotes. Reconstruction and biochemical characterization of the likely ancestor protein suggest that an intrinsic DNA binding promiscuity in the family was resolved by alternative mechanisms in different branches of fungal SREBPs. Furthermore, we show that two SREBPs in the human commensal yeast Candida albicans drive a transcriptional cascade that inhibits a morphological switch under anaerobic conditions. Preventing this morphological transition enhances C. albicans colonization of the mammalian intestine, the fungus' natural niche. Thus, our results illustrate how diversification in DNA binding preferences enabled the functional expansion of a family of eukaryotic transcription regulators.

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.

Phenotype loss is associated with widespread divergence of the gene regulatory landscape in evolution

Detecting the genomic changes underlying phenotypic changes between species is a main goal of evolutionary biology and genomics. Evolutionary theory predicts that changes in cis-regulatory elements are important for morphological changes. We combined genome sequencing, functional genomics and genome-wide comparative analyses to investigate regulatory elements in lineages that lost morphological traits. We first show that limb loss in snakes is associated with widespread divergence of limb regulatory elements. We next show that eye degeneration in subterranean mammals is associated with widespread divergence of eye regulatory elements. In both cases, sequence divergence results in an extensive loss of transcription factor binding sites. Importantly, diverged regulatory elements are associated with genes required for normal limb patterning or normal eye development and function, suggesting that regulatory divergence contributed to the loss of these phenotypes. Together, our results show that genome-wide decay of the phenotype-specific cis-regulatory landscape is a hallmark of lost morphological traits.

Cis-regulatory elements are important factors for morphological changes. Here, the authors show widespread divergence of limb and eye regulatory elements in limb loss in snakes and eye degeneration in subterranean mammals respectively.

Species and cell-type properties of classically defined human and rodent neurons and glia

Determination of the molecular properties of genetically targeted cell types has led to fundamental insights into mouse brain function and dysfunction. Here, we report an efficient strategy for precise exploration of gene expression and epigenetic events in specific cell types in a range of species, including postmortem human brain. We demonstrate that classically defined, homologous neuronal and glial cell types differ between rodent and human by the expression of hundreds of orthologous, cell specific genes. Confirmation that these genes are differentially active was obtained using epigenetic mapping and immunofluorescence localization. Studies of sixteen human postmortem brains revealed gender specific transcriptional differences, cell-specific molecular responses to aging, and the induction of a shared, robust response to an unknown external event evident in three donor samples. Our data establish a comprehensive approach for analysis of molecular events associated with specific circuits and cell types in a wide variety of human conditions.

Ancient mechanisms for the evolution of the bicoid homeodomain's function in fly development

The ancient mechanisms that caused developmental gene regulatory networks to diversify among distantly related taxa are not well understood. Here we use ancestral protein reconstruction, biochemical experiments, and developmental assays of transgenic animals carrying reconstructed ancestral genes to investigate how the transcription factor Bicoid (Bcd) evolved its central role in anterior-posterior patterning in flies. We show that most of Bcd's derived functions are attributable to evolutionary changes within its homeodomain (HD) during a phylogenetic interval >140 million years ago. A single substitution from this period (Q50K) accounts almost entirely for the evolution of Bcd's derived DNA specificity in vitro. In transgenic embryos expressing the reconstructed ancestral HD, however, Q50K confers activation of only a few of Bcd's transcriptional targets and yields a very partial rescue of anterior development. Adding a second historical substitution (M54R) confers regulation of additional Bcd targets and further rescues anterior development. These results indicate that two epistatically interacting mutations played a major role in the evolution of Bcd's controlling regulatory role in early development. They also show how ancestral sequence reconstruction can be combined with in vivo characterization of transgenic animals to illuminate the historical mechanisms of developmental evolution.

Developmental Bias and Evolution: A Regulatory Network Perspective

Phenotypic variation is generated by the processes of development, with some variants arising more readily than others-a phenomenon known as "developmental bias." Developmental bias and natural selection have often been portrayed as alternative explanations, but this is a false dichotomy: developmental bias can evolve through natural selection, and bias and selection jointly influence phenotypic evolution. Here, we briefly review the evidence for developmental bias and illustrate how it is studied empirically. We describe recent theory on regulatory networks that explains why the influence of genetic and environmental perturbation on phenotypes is typically not uniform, and may even be biased toward adaptive phenotypic variation. We show how bias produced by developmental processes constitutes an evolving property able to impose direction on adaptive evolution and influence patterns of taxonomic and phenotypic diversity. Taking these considerations together, we argue that it is not sufficient to accommodate developmental bias into evolutionary theory merely as a constraint on evolutionary adaptation. The influence of natural selection in shaping developmental bias, and conversely, the influence of developmental bias in shaping subsequent opportunities for adaptation, requires mechanistic models of development to be expanded and incorporated into evolutionary theory. A regulatory network perspective on phenotypic evolution thus helps to integrate the generation of phenotypic variation with natural selection, leaving evolutionary biology better placed to explain how organisms adapt and diversify.

Introgression of regulatory alleles and a missense coding mutation drive plumage pattern diversity in the rock pigeon

Birds and other vertebrates display stunning variation in pigmentation patterning, yet the genes controlling this diversity remain largely unknown. Rock pigeons (Columba livia) are fundamentally one of four color pattern phenotypes, in decreasing order of melanism: T-check, checker, bar (ancestral), or barless. Using whole-genome scans, we identified NDP as a candidate gene for this variation. Allele-specific expression differences in NDP indicate cis-regulatory divergence between ancestral and melanistic alleles. Sequence comparisons suggest that derived alleles originated in the speckled pigeon (Columba guinea), providing a striking example of introgression. In contrast, barless rock pigeons have an increased incidence of vision defects and, like human families with hereditary blindness, carry start-codon mutations in NDP. In summary, we find that both coding and regulatory variation in the same gene drives wing pattern diversity, and post-domestication introgression supplied potentially advantageous melanistic alleles to feral populations of this ubiquitous urban bird.

Gene regulatory network architecture in different developmental contexts influences the genetic basis of morphological evolution

Convergent phenotypic evolution is often caused by recurrent changes at particular nodes in the underlying gene regulatory networks (GRNs). The genes at such evolutionary ‘hotspots’ are thought to maximally affect the phenotype with minimal pleiotropic consequences. This has led to the suggestion that if a GRN is understood in sufficient detail, the path of evolution may be predictable. The repeated evolutionary loss of larval trichomes among Drosophila species is caused by the loss of shavenbaby (svb) expression. svb is also required for development of leg trichomes, but the evolutionary gain of trichomes in the ‘naked valley’ on T2 femurs in Drosophila melanogaster is caused by reduced microRNA-92a (miR-92a) expression rather than changes in svb. We compared the expression and function of components between the larval and leg trichome GRNs to investigate why the genetic basis of trichome pattern evolution differs in these developmental contexts. We found key differences between the two networks in both the genes employed, and in the regulation and function of common genes. These differences in the GRNs reveal why mutations in svb are unlikely to contribute to leg trichome evolution and how instead miR-92a represents the key evolutionary switch in this context. Our work shows that variability in GRNs across different developmental contexts, as well as whether a morphological feature is lost versus gained, influence the nodes at which a GRN evolves to cause morphological change. Therefore, our findings have important implications for understanding the pathways and predictability of evolution.

A major goal of biology is to identify the genetic causes of organismal diversity. Convergent evolution of traits is often caused by changes in the same genes–evolutionary ‘hotspots’. shavenbaby is a ‘hotspot’ for larval trichome loss in Drosophila, but microRNA-92a underlies the gain of leg trichomes. To understand this difference in the genetics of phenotypic evolution, we compared the expression and function of genes in the underlying regulatory networks. We found that the pathway of evolution is influenced by differences in gene regulatory network architecture in different developmental contexts, as well as by whether a trait is lost or gained. Therefore, hotspots in one context may not readily evolve in a different context. This has important implications for understanding the genetic basis of phenotypic change and the predictability of evolution.

Revisiting the developmental and cellular role of the pigmentation gene yellow in Drosophila using a tagged allele

Pigmentation is a diverse and ecologically relevant trait in insects. Pigment formation has been studied extensively at the genetic and biochemical levels. The temporality of pigment formation during animal development, however, is more elusive. Here, we examine this temporality, focusing on yellow, a gene involved in the formation of black melanin. We generated a protein-tagged yellow allele in the fruit fly Drosophila melanogaster, which allowed us to precisely describe Yellow expression pattern at the tissue and cellular levels throughout development. We found Yellow expressed in the pupal epidermis in patterns prefiguring black pigmentation. We also found Yellow expressed in a few central neurons from the second larval instar to adult stages, including a subset of neurons adjacent to the clock neurons marked by the gene Pdf. We then specifically examined the dynamics of Yellow expression domain and subcellular localization in relationship to pigment formation. In particular, we showed how a late step of re-internalization is regulated by the large low-density lipoprotein receptor-related protein Megalin. Finally we suggest a new function for Yellow in the establishment of sharp pigmentation pattern boundaries, whereby this protein may assume a structural role, anchoring pigment deposits or pigmentation enzymes in the cuticle.

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.

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.

Enhancer adoption caused by genomic insertion elicits interdigital Shh expression and syndactyly in mouse

In this study, we reexamined an old mouse mutant named Hammer toe (Hm), which arose spontaneously almost a half century ago and exhibits a limb phenotype with webbing. We revealed that a 150-kb noncoding genomic fragment that was originally located in chromosome 14 has been inserted into a genomic region proximal to Sonic hedgehog (Shh), located in chromosome 5. This inserted fragment possesses enhancer activity to induce Shh expression in the interdigital regions in Hm, which in turn down-regulates bone morphogenetic protein signaling and eventually results in syndactyly and web formation. Since the donor fragment residing in chromosome 14 has enhancer activity to induce interdigital gene expression, the Hm mutation appears to be an archetypal case of enhancer adoption.

Acquisition of new cis-regulatory elements (CREs) can cause alteration of developmental gene regulation and may introduce morphological novelty in evolution. Although structural variation in the genome generated by chromosomal rearrangement is one possible source of new CREs, only a few examples are known, except for cases of retrotransposition. In this study, we show the acquisition of novel regulatory sequences as a result of large genomic insertion in the spontaneous mouse mutation Hammer toe (Hm). Hm mice exhibit syndactyly with webbing, due to suppression of interdigital cell death in limb development. We reveal that, in the Hm genome, a 150-kb noncoding DNA fragment from chromosome 14 is inserted into the region upstream of the Sonic hedgehog (Shh) promoter in chromosome 5. Phenotyping of mouse embryos with a series of CRISPR/Cas9-aided partial deletion of the 150-kb insert clearly indicated that two different regions are necessary for the syndactyly phenotype of Hm. We found that each of the two regions contains at least one enhancer for interdigital regulation. These results show that a set of enhancers brought by the large genomic insertion elicits the interdigital Shh expression and the Hm phenotype. Transcriptome analysis indicates that ectopic expression of Shh up-regulates Chordin (Chrd) that antagonizes bone morphogenetic protein signaling in the interdigital region. Indeed, Chrd-overexpressing transgenic mice recapitulated syndactyly with webbing. Thus, the Hm mutation provides an insight into enhancer acquisition as a source of creation of novel gene regulation.

Geometric morphometrics reveals shifts in flower shape symmetry and size following gene knockdown of CYCLOIDEA and ANTHOCYANIDIN SYNTHASE

BACKGROUND

While floral symmetry has traditionally been assessed qualitatively, recent advances in geometric morphometrics have opened up new avenues to specifically quantify flower shape and size using robust multivariate statistical methods. In this study, we examine, for the first time, the ability of geometric morphometrics to detect morphological differences in floral dorsoventral asymmetry following virus-induced gene silencing (VIGS). Using Fedia graciliflora Fisch. & Meyer (Valerianaceae) as a model, corolla shape of untreated flowers was compared using canonical variate analysis to knockdown phenotypes of CYCLOIDEA2A (FgCYC2A), ANTHOCYANIDIN SYNTHASE (FgANS), and empty vector controls.

RESULTS

Untreated flowers and all VIGS treatments were morphologically distinct from each other, suggesting that VIGS may cause subtle shifts in floral shape. Knockdowns of FgCYC2A were the most dramatic, affecting the position of dorsal petals in relation to lateral petals, thereby resulting in more actinomorphic-like flowers. Additionally, FgANS knockdowns developed larger flowers with wider corolla tube openings.

CONCLUSIONS

These results provide a method to quantify the role that specific genes play in the developmental pathway affecting the dorsoventral axis of symmetry in zygomorphic flowers. Additionally, they suggest that ANS may have an unintended effect on floral size and shape.

The unexpected depths of genome-skimming data: A case study examining Goodeniaceae floral symmetry genes

PREMISE OF THE STUDY

The use of genome skimming allows systematists to quickly generate large data sets, particularly of sequences in high abundance (e.g., plastomes); however, researchers may be overlooking data in low abundance that could be used for phylogenetic or evo-devo studies. Here, we present a bioinformatics approach that explores the low-abundance portion of genome-skimming next-generation sequencing libraries in the fan-flowered Goodeniaceae.

METHODS

Twenty-four previously constructed Goodeniaceae genome-skimming Illumina libraries were examined for their utility in mining low-copy nuclear genes involved in floral symmetry, specifically the CYCLOIDEA (CYC)-like genes. De novo assemblies were generated using multiple assemblers, and BLAST searches were performed for CYC1, CYC2, and CYC3 genes.

RESULTS

Overall Trinity, SOAPdenovo-Trans, and SOAPdenovo implementing lower k-mer values uncovered the most data, although no assembler consistently outperformed the others. Using SOAPdenovo-Trans across all 24 data sets, we recovered four CYC-like gene groups (CYC1, CYC2, CYC3A, and CYC3B) from a majority of the species. Alignments of the fragments included the entire coding sequence as well as upstream and downstream regions.

DISCUSSION

Genome-skimming data sets can provide a significant source of low-copy nuclear gene sequence data that may be used for multiple downstream applications.

Divergence of annual and perennial species in the Brassicaceae and the contribution of cis-acting variation at FLC orthologues

Variation in life history contributes to reproductive success in different environments. Divergence of annual and perennial angiosperm species is an extreme example that has occurred frequently. Perennials survive for several years and restrict the duration of reproduction by cycling between vegetative growth and flowering, whereas annuals live for 1 year and flower once. We used the tribe Arabideae (Brassicaceae) to study the divergence of seasonal flowering behaviour among annual and perennial species. In perennial Brassicaceae, orthologues of FLOWERING LOCUS C (FLC), a floral inhibitor in Arabidopsis thaliana, are repressed by winter cold and reactivated in spring conferring seasonal flowering patterns, whereas in annuals, they are stably repressed by cold. We isolated FLC orthologues from three annual and two perennial Arabis species and found that the duplicated structure of the A. alpina locus is not required for perenniality. The expression patterns of the genes differed between annuals and perennials, as observed among Arabidopsis species, suggesting a broad relevance of these patterns within the Brassicaceae. Also analysis of plants derived from an interspecies cross of A. alpina and annual A. montbretiana demonstrated that cis-regulatory changes in FLC orthologues contribute to their different transcriptional patterns. Sequence comparisons of FLC orthologues from annuals and perennials in the tribes Arabideae and Camelineae identified two regulatory regions in the first intron whose sequence variation correlates with divergence of the annual and perennial expression patterns. Thus, we propose that related cis-acting changes in FLC orthologues occur independently in different tribes of the Brassicaceae during life history evolution.

Perspectives on Gene Regulatory Network Evolution

Animal development proceeds through the activity of genes and their cis-regulatory modules (CRMs) working together in sets of gene regulatory networks (GRNs). The emergence of species-specific traits and novel structures results from evolutionary changes in GRNs. Recent work in a wide variety of animal models, and particularly in insects, has started to reveal the modes and mechanisms of GRN evolution. I discuss here various aspects of GRN evolution and argue that developmental system drift (DSD), in which conserved phenotype is nevertheless a result of changed genetic interactions, should regularly be viewed from the perspective of GRN evolution. Advances in methods to discover related CRMs in diverse insect species, a critical requirement for detailed GRN characterization, are also described.

Evolution of reduced co-activator dependence led to target expansion of a starvation response pathway

Although combinatorial regulation is a common feature in gene regulatory networks, how it evolves and affects network structure and function is not well understood. In S. cerevisiae, the phosphate starvation (PHO) responsive transcription factors Pho4 and Pho2 are required for gene induction and survival during phosphate starvation. In the related human commensal C. glabrata, Pho4 is required but Pho2 is dispensable for survival in phosphate starvation and is only partially required for inducing PHO genes. Phylogenetic survey suggests that reduced dependence on Pho2 evolved in C. glabrata and closely related species. In S. cerevisiae, less Pho2-dependent Pho4 orthologs induce more genes. In C. glabrata, its Pho4 binds to more locations and induces three times as many genes as Pho4 in S. cerevisiae does. Our work shows how evolution of combinatorial regulation allows for rapid expansion of a gene regulatory network’s targets, possibly extending its physiological functions.

The diversity of life on Earth has intrigued generations of scientists and nature lovers alike. Research over recent decades has revealed that much of the diversity we can see did not require the invention of new genes. Instead, living forms diversified mostly by using old genes in new ways – for example, by changing when or where an existing gene became active. This kind of change is referred to as “regulatory evolution”.

A class of proteins called transcription factors are hot spots in regulatory evolution. These proteins recognize specific sequences of DNA to control the activity of other genes, and so represent the “readers” of the genetic information. Small changes to how a transcription factor is regulated, or the genes it targets, can lead to dramatic changes in an organism. Before we can understand how life on Earth evolved to be so diverse, scientists must first answer how transcription factors evolve and what consequences this has on their target genes.

So far, most studies of regulatory evolution have focused on networks of transcription factors and genes that control how an organism develops. He et al. have now studied a regulatory network that is behind a different process, namely how an organism responds to stress or starvation. These two types of regulatory networks are structured differently and work in different ways. These differences made He et al. wonder if the networks evolved differently too.

The chemical phosphate is an essential nutrient for all living things, and He et al. compared how two different species of yeast responded to a lack of phosphate. The key difference was how much a major transcription factor known as Pho4 depended on a so-called co-activator protein named Pho2 to carry out its role. Baker’s yeast (Saccharomyces cerevisiae), which is commonly used in laboratory experiments, requires both Pho4 and Pho2 to activate about 20 genes when inorganic phosphate is not available in its environment. However, in a related yeast species called Candida glabrata, Pho4 has evolved to depend less on Pho2. He et al. went on to show that, as well as being less dependent on Pho2, Pho4 in C. glabrata activates more than three times as many genes as Pho4 in S. cerevisiae does in the absence of phosphate. These additional gene targets for Pho4 in C. glabrata are predicted to extend the network’s activities, and allow it to regulate new process including the yeast’s responses to other types of stress and the building of the yeast’s cell wall.

Together these findings show a new way that regulatory networks can evolve, that is, by reducing its dependence on the co-activator, a transcription factor can expand the number of genes it targets. This has not been seen for regulatory networks related to development, suggesting that different networks can indeed evolve in different ways. Lastly, because disease-causing microbes are often stressed inside their hosts and C. glabrata sometimes infects humans, understanding how this yeast’s response to stress has evolved may lead to new ways to prevent and treat this infection.

Diversification of spatiotemporal expression and copy number variation of the echinoid hbox12/pmar1/micro1 multigene family

Changes occurring during evolution in the cis-regulatory landscapes of individual members of multigene families might impart diversification in their spatiotemporal expression and function. The archetypal member of the echinoid hbox12/pmar1/micro1 family is hbox12-a, a homeobox-containing gene expressed exclusively by dorsal blastomeres, where it governs the dorsal/ventral gene regulatory network during embryogenesis of the sea urchin Paracentrotus lividus. Here we describe the inventory of the hbox12/pmar1/micro1 genes in P. lividus, highlighting that gene copy number variation occurs across individual sea urchins of the same species. We show that the various hbox12/pmar1/micro1 genes group into three subfamilies according to their spatiotemporal expression, which ranges from broad transcription throughout development to transient expression in either the animal hemisphere or micromeres of the early embryo. Interestingly, the promoter regions of those genes showing comparable expression patterns are highly similar, while differing from those of the other subfamilies. Strikingly, phylogenetic analysis suggests that the hbox12/pmar1/micro1 genes are species-specific, exhibiting extensive divergence in their noncoding, but not in their coding, sequences across three distinct sea urchin species. In spite of this, two micromere-specific genes of P. lividus possess a TCF/LEF-binding motif in a similar position, and their transcription relies on Wnt/β-catenin signaling, similar to the pmar1 and micro1 genes, which in other sea urchin species are involved in micromere specification. Altogether, our findings suggest that the hbox12/pmar1/micro1 gene family evolved rather rapidly, generating paralogs whose cis-regulatory sequences diverged following multiple rounds of duplication from a common ancestor.

Evolution of Multiple Sensory Systems Drives Novel Egg-Laying Behavior in the Fruit Pest Drosophila suzukii

The rise of a pest species represents a unique opportunity to address how species evolve new behaviors and adapt to novel ecological niches [1]. We address this question by studying the egg-laying behavior of Drosophila suzukii, an invasive agricultural pest species that has spread from Southeast Asia to Europe and North America in the last decade [2]. While most closely related Drosophila species lay their eggs on decaying plant substrates, D. suzukii oviposits on ripening fruit, thereby causing substantial economic losses to the fruit industry [3, 4, 5, 6, 7, 8]. D. suzukii has evolved an enlarged, serrated ovipositor that presumably plays a key role by enabling females to pierce the skin of ripe fruit [9]. Here, we explore how D. suzukii selects oviposition sites, and how this behavior differs from that of closely related species. We have combined behavioral experiments in multiple species with neurogenetics and mutant analysis in D. suzukii to show that this species has evolved a specific preference for oviposition on ripe fruit. Our results also establish that changes in mechanosensation, olfaction, and presumably gustation have contributed to this ecological shift. Our observations support a model in which the emergence of D. suzukii as an agricultural pest is the consequence of the progressive modification of several sensory systems, which collectively underlie a radical change in oviposition behavior.

•

The pest Drosophila suzukii prefers to lay eggs on ripening fruit

•

Closely related Drosophila species prefer to lay eggs on rotten fruit

•

Female flies use chemosensation and mechanosensation to choose an oviposition site

•

Orco-dependent detection of ripe fruit odors elicits oviposition in D. suzukii

Different contributions of local- and distant-regulatory changes to transcriptome divergence between stickleback ecotypes

Differential gene expression can play an important role in phenotypic evolution and divergent adaptation. Although differential gene expression can be caused by both local- and distant-regulatory changes, we know little about their relative contribution to transcriptome evolution in natural populations. Here, we conducted expression quantitative trait loci (eQTL) analysis to investigate the genetic architecture underlying transcriptome divergence between marine and stream ecotypes of threespine sticklebacks (Gasterosteus aculeatus). We identified both local and distant eQTLs, some of which constitute hotspots, regions with a disproportionate number of significant eQTLs relative to the genomic background. The majority of local eQTLs including those in the hotspots caused expression changes consistent with the direction of transcriptomic divergence between ecotypes. Genome scan analysis showed that many local eQTLs overlapped with genomic regions of high differentiation. In contrast, nearly half of the distant eQTLs including those in the hotspots caused opposite expression changes, and few overlapped with regions of high differentiation, indicating that distant eQTLs may act as a constraint of transcriptome evolution. Finally, a comparison between two salinity conditions revealed that nearly half of eQTL hotspots were environment specific, suggesting that analysis of genetic architecture in multiple conditions is essential for predicting response to selection.

The evolutionary origination of a novel expression pattern through an extreme heterochronic shift

The evolutionary origins of morphological structures are thought to often depend upon the redeployment of old genes into new developmental settings. Although many examples of cis-regulatory divergence have shown how pre-existing patterns of gene expression have been altered, only a small number of case studies have traced the origins of cis-regulatory elements that drive new expression domains. Here, we elucidate the evolutionary history of a novel expression pattern of the yellow gene within the Zaprionus genus of fruit flies. We observed a unique pattern of yellow transcript accumulation in the wing disc during the third larval instar, a stage that precedes its typical expression pattern associated with cuticular melanization by about a week. The region of the Zaprionus wing disc that expresses yellow subsequently develops into a portion of the thorax, a tissue for which yellow expression has been reported for several fruit fly species. Tests of GFP reporter transgenes containing the Zaprionus yellow regulatory region revealed that the wing disc pattern arose by changes in the cis-regulatory region of yellow. Moreover, the wing disc enhancer activity of yellow depends upon a short conserved sequence with ancestral thoracic functions, suggesting that the pupal thorax regulatory sequence was genetically reprogrammed to drive expression that commences much earlier during development. These results highlight how novel domains of gene expression may arise by extreme shifts in timing during the origins of novel traits.

The transcriptional architecture of phenotypic dimorphism

The profound differences in gene expression between the sexes are increasingly used to study the molecular basis of sexual dimorphism, sexual selection and sexual conflict. Studies of transcriptional architecture, based on comparisons of gene expression, have also been implemented for a wide variety of other intra-specific polymorphisms. These efforts are based on key assumptions regarding the relationship between transcriptional architecture, phenotypic variation and the target of selection. Some of these assumptions are better supported by available evidence than others. In all cases, the evidence is largely circumstantial, leaving considerable gaps in our understanding of the relationship between transcriptional and phenotypic dimorphism.

Limb patterning genes and heterochronic development of the emu wing bud

Background

The forelimb of the flightless emu is a vestigial structure, with greatly reduced wing elements and digit loss. To explore the molecular and cellular mechanisms associated with the evolution of vestigial wings and loss of flight in the emu, key limb patterning genes were examined in developing embryos.

Methods

Limb development was compared in emu versus chicken embryos. Immunostaining for cell proliferation markers was used to analyze growth of the emu forelimb and hindlimb buds. Expression patterns of limb patterning genes were studied, using whole-mount in situ hybridization (for mRNA localization) and RNA-seq (for mRNA expression levels).

Results

The forelimb of the emu embryo showed heterochronic development compared to that in the chicken, with the forelimb bud being retarded in its development. Early outgrowth of the emu forelimb bud is characterized by a lower level of cell proliferation compared the hindlimb bud, as assessed by PH3 immunostaining. In contrast, there were no obvious differences in apoptosis in forelimb versus hindlimb buds (cleaved caspase 3 staining). Most key patterning genes were expressed in emu forelimb buds similarly to that observed in the chicken, but with smaller expression domains. However, expression of Sonic Hedgehog (Shh) mRNA, which is central to anterior–posterior axis development, was delayed in the emu forelimb bud relative to other patterning genes. Regulators of Shh expression, Gli3 and HoxD13, also showed altered expression levels in the emu forelimb bud.

Conclusions

These data reveal heterochronic but otherwise normal expression of most patterning genes in the emu vestigial forelimb. Delayed Shh expression may be related to the small and vestigial structure of the emu forelimb bud. However, the genetic mechanism driving retarded emu wing development is likely to rest within the forelimb field of the lateral plate mesoderm, predating the expression of patterning genes.

Electronic supplementary material

The online version of this article (doi:10.1186/s13227-016-0063-5) contains supplementary material, which is available to authorized users.

Enhancer modularity and the evolution of new traits

Animals have modular cis-regulatory regions in their genomes, and expression of a single gene is often regulated by multiple enhancers residing in such a region. In the laboratory, and also in natural populations, loss of an enhancer can result in a loss of gene expression. Although only a few examples have been well characterized to date, some studies have suggested that an evolutionary gain of a new enhancer function can establish a new gene expression domain. Our recent study showed that Drosophila guttifera has more enhancers and additional expression domains of the wingless gene during the pupal stage, compared to D. melanogaster, and that these new features appear to have evolved in the ancestral lineage leading to D. guttifera. (1) Gain of a new expression domain of a developmental regulatory gene (toolkit gene), such as wingless, can cause co-option of the expression of its downstream genes to the new domain, resulting in duplication of a preexisting structure at this new body position. Recently, with the advancement of evo-devo studies, we have learned that the developmental regulatory systems are strikingly similar across various animal taxa, in spite of the great diversity of the animals' morphology. Even behind "new" traits, co-options of essential developmental genes from known systems are very common. We previously provided concrete evidence of gains of enhancer activities of a developmental regulatory gene underlying gains of new traits. (1) Broad occurrence of this scenario is testable and should be validated in the future.

Coupled enhancer and coding sequence evolution of a homeobox gene shaped leaf diversity

In this study, Vuolo et al. investigate the mechanisms underlying the genetic basis for morphological diversity in leaf shape. They show that evolution of an enhancer element in the homeobox gene REDUCED COMPLEXITY (RCO) altered leaf shape by changing gene expression from the distal leaf blade to its base.

Here we investigate mechanisms underlying the diversification of biological forms using crucifer leaf shape as an example. We show that evolution of an enhancer element in the homeobox gene REDUCED COMPLEXITY (RCO) altered leaf shape by changing gene expression from the distal leaf blade to its base. A single amino acid substitution evolved together with this regulatory change, which reduced RCO protein stability, preventing pleiotropic effects caused by its altered gene expression. We detected hallmarks of positive selection in these evolved regulatory and coding sequence variants and showed that modulating RCO activity can improve plant physiological performance. Therefore, interplay between enhancer and coding sequence evolution created a potentially adaptive path for morphological evolution.

Potential Direct Regulators of the Drosophila yellow Gene Identified by Yeast One-Hybrid and RNAi Screens

The regulation of gene expression controls development, and changes in this regulation often contribute to phenotypic evolution. Drosophila pigmentation is a model system for studying evolutionary changes in gene regulation, with differences in expression of pigmentation genes such as yellow that correlate with divergent pigment patterns among species shown to be caused by changes in cis- and trans-regulation. Currently, much more is known about the cis-regulatory component of divergent yellow expression than the trans-regulatory component, in part because very few trans-acting regulators of yellow expression have been identified. This study aims to improve our understanding of the trans-acting control of yellow expression by combining yeast-one-hybrid and RNAi screens for transcription factors binding to yellow cis-regulatory sequences and affecting abdominal pigmentation in adults, respectively. Of the 670 transcription factors included in the yeast-one-hybrid screen, 45 showed evidence of binding to one or more sequence fragments tested from the 5′ intergenic and intronic yellow sequences from D. melanogaster, D. pseudoobscura, and D. willistoni, suggesting that they might be direct regulators of yellow expression. Of the 670 transcription factors included in the yeast-one-hybrid screen, plus another TF previously shown to be genetically upstream of yellow, 125 were also tested using RNAi, and 32 showed altered abdominal pigmentation. Nine transcription factors were identified in both screens, including four nuclear receptors related to ecdysone signaling (Hr78, Hr38, Hr46, and Eip78C). This finding suggests that yellow expression might be directly controlled by nuclear receptors influenced by ecdysone during early pupal development when adult pigmentation is forming.

The Functionality and Evolution of Eukaryotic Transcriptional Enhancers

Enhancers regulate precise spatial and temporal patterns of gene expression in eukaryotes and, moreover, evolutionary changes in these modular cis-regulatory elements may represent the predominant genetic basis for phenotypic evolution. Here, we review approaches to identify and functionally analyze enhancers and their transcription factor binding sites, including assay for transposable-accessible chromatin-sequencing (ATAC-Seq) and clustered regularly interspaced short palindromic repeats (CRISPR)/Cas9, respectively. We also explore enhancer functionality, including how transcription factor binding sites combine to regulate transcription, as well as research on shadow and super enhancers, and how enhancers can act over great distances and even in trans. Finally, we discuss recent theoretical and empirical data on how transcription factor binding sites and enhancers evolve. This includes how the function of enhancers is maintained despite the turnover of transcription factor binding sites as well as reviewing studies where mutations in enhancers have been shown to underlie morphological change.

Folded gastrulation and T48 drive the evolution of coordinated mesoderm internalization in flies

Gastrulation constitutes a fundamental yet diverse morphogenetic process of metazoan development. Modes of gastrulation range from stochastic translocation of individual cells to coordinated infolding of an epithelial sheet. How such morphogenetic differences are genetically encoded and whether they have provided specific developmental advantages is unclear. Here we identify two genes, folded gastrulation and t48, which in the evolution of fly gastrulation acted as a likely switch from an ingression of individual cells to the invagination of the blastoderm epithelium. Both genes are expressed and required for mesoderm invagination in the fruit fly Drosophila melanogaster but do not appear during mesoderm ingression of the midge Chironomus riparius. We demonstrate that early expression of either or both of these genes in C.riparius is sufficient to invoke mesoderm invagination similar to D.melanogaster. The possible genetic simplicity and a measurable increase in developmental robustness might explain repeated evolution of similar transitions in animal gastrulation.

DOI:http://dx.doi.org/10.7554/eLife.18318.001

In animals, gastrulation is a period of time in early development during which a sphere of cells is reorganized into an embryo with cells arranged into three distinct layers (called germ layers). The process has changed substantially during the course of evolution and thus provides a great experimental system to investigate the genetic basis for differences in animal form and shape.

As an example, true flies use at least two different mechanisms to make the middle germ layer (the mesoderm). In both cases, the mesoderm is made up of cells that move inwards from the boundary of the outer germ layer. In midges and some other flies, these cells migrate individually, while in others including fruit flies, the cells move together as a sheet. Fruit flies and midges shared their last common ancestor 250 million years ago and although the genes that make the mesoderm in fruit flies are well understood, little is known about how the mesoderm forms in midges.

Urbansky, González Avalos et al. investigated which genes were responsible for the evolutionary transition between the different types of cell migration seen in flies. The experiments identified two genes – called folded gastrulation and t48 – that seem to operate as a simple switch between the two ways that mesoderm cells migrate. Both of these genes are active in fruit fly embryos and are required for the group migration of mesoderm cells. However, the genes do not appear to play a major role in the movement of individual mesoderm cells in midges.

Further experiments demonstrate that switching on these genes in midge embryos is sufficient to invoke group mesoderm cell migrations similar to those seen in fruit flies. These findings show that it is possible to identify genetic changes that underlie substantial differences in animal form and shape over hundred million of years. The ease by which Urbansky, González Avalos et al. were able to switch between the two types of mesoderm migration might explain why similar transitions in gastrulation have evolved repeatedly in animals. The next step is to test this hypothesis in other animals.

DOI:http://dx.doi.org/10.7554/eLife.18318.002