optix drives the repeated convergent evolution of butterfly wing pattern mimicry

Structural color in Junonia butterflies evolves by tuning scale lamina thickness

In diverse organisms, nanostructures that coherently scatter light create structural color, but how such structures are built remains mysterious. We investigate the evolution and genetic regulation of butterfly scale laminae, which are simple photonic nanostructures. In a lineage of buckeye butterflies artificially selected for blue wing color, we found that thickened laminae caused a color shift from brown to blue. Deletion of the optix patterning gene also altered color via lamina thickening, revealing shared regulation of pigments and lamina thickness. Finally, we show how lamina thickness variation contributes to the color diversity that distinguishes sexes and species throughout the genus Junonia. Thus, quantitatively tuning one dimension of scale architecture facilitates both the microevolution and macroevolution of a broad spectrum of hues. Because the lamina is an intrinsic component of typical butterfly scales, our findings suggest that tuning lamina thickness is an available mechanism to create structural color across the Lepidoptera.

From iridescent blues to vibrant purples, many butterflies display dazzling ‘structural colors’ created not by pigments but by microscopic structures that interfere with light. For instance, the scales that coat their wings can contain thin films of chitin, the substance that normally makes the external skeleton of insects. In slim layers, however, chitin can also scatter light to produce color, the way that oil can create iridescence at the surface of water.

The thickness of the film, which is encoded by the genes of the butterfly, determines what color will be produced. Yet, little is known about how common thin films are in butterflies, exactly how genetic information codes for them, and how their thickness and the colors they produce can evolve.

To investigate, Thayer et al. used a technique called Helium Ion Microscopy and examined the wings of ten related species of butterflies, showing that thin film structures were present across this sample. However, the different species have evolved many different structural colors over the past millions of years by changing the thickness of the films.

Next, Thayer et al. showed that this evolution could be reproduced at a faster pace in the laboratory using common buckeye butterflies. These insects mostly have brown wings, but they can have specks of blue created by thin film structures. Individuals with more blue on their wings were mated and over the course of a year, the thickness of the film structures increased by 74%, leading to shiny blue butterflies. Deleting a gene called optix from the insects also led to blue wings. Optix was already known to control the patterns of pigments in butterflies, but it now appears that it controls structural colors as well.

From solar panels to new fabrics, microscopic structures that can scatter light are useful in a variety of industries. Understanding how these elements exist and evolve in organisms may help to better design them for human purposes.

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.

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.

Contrasting genomic and phenotypic outcomes of hybridization between pairs of mimetic butterfly taxa across a suture zone

Hybrid zones, whereby divergent lineages come into contact and eventually hybridize, can provide insights on the mechanisms involved in population differentiation and reproductive isolation, and ultimately speciation. Suture zones offer the opportunity to compare these processes across multiple species. In this paper we use reduced-complexity genomic data to compare the genetic and phenotypic structure and hybridization patterns of two mimetic butterfly species, Ithomia salapia and Oleria onega (Nymphalidae: Ithomiini), each consisting of a pair of lineages differentiated for their wing colour pattern and that come into contact in the Andean foothills of Peru. Despite similarities in their life history, we highlight major differences, both at the genomic and phenotypic level, between the two species. These differences include the presence of hybrids, variations in wing phenotype, and genomic patterns of introgression and differentiation. In I. salapia, the two lineages appear to hybridize only rarely, whereas in O. onega the hybrids are not only more common, but also genetically and phenotypically more variable. We also detected loci statistically associated with wing colour pattern variation, but in both species these loci were not over-represented among the candidate barrier loci, suggesting that traits other than wing colour pattern may be important for reproductive isolation. Our results contrast with the genomic patterns observed between hybridizing lineages in the mimetic Heliconius butterflies, and call for a broader investigation into the genomics of speciation in Ithomiini - the largest radiation of mimetic butterflies.

Butterfly Mimicry Polymorphisms Highlight Phylogenetic Limits of Gene Reuse in the Evolution of Diverse Adaptations

Some genes have repeatedly been found to control diverse adaptations in a wide variety of organisms. Such gene reuse reveals not only the diversity of phenotypes these unique genes control but also the composition of developmental gene networks and the genetic routes available to and taken by organisms during adaptation. However, the causes of gene reuse remain unclear. A small number of large-effect Mendelian loci control a huge diversity of mimetic butterfly wing color patterns, but reasons for their reuse are difficult to identify because the genetic basis of mimicry has primarily been studied in two systems with correlated factors: female-limited Batesian mimicry in Papilio swallowtails (Papilionidae) and non-sex-limited Müllerian mimicry in Heliconius longwings (Nymphalidae). Here, we break the correlation between phylogenetic relationship and sex-limited mimicry by identifying loci controlling female-limited mimicry polymorphism Hypolimnas misippus (Nymphalidae) and non-sex-limited mimicry polymorphism in Papilio clytia (Papilionidae). The Papilio clytia polymorphism is controlled by the genome region containing the gene cortex, the classic P supergene in Heliconius numata, and loci controlling color pattern variation across Lepidoptera. In contrast, female-limited mimicry polymorphism in Hypolimnas misippus is associated with a locus not previously implicated in color patterning. Thus, although many species repeatedly converged on cortex and its neighboring genes over 120 My of evolution of diverse color patterns, female-limited mimicry polymorphisms each evolved using a different gene. Our results support conclusions that gene reuse occurs mainly within ∼10 My and highlight the puzzling diversity of genes controlling seemingly complex female-limited mimicry polymorphisms.

Whole-chromosome hitchhiking driven by a male-killing endosymbiont

Neo-sex chromosomes are found in many taxa, but the forces driving their emergence and spread are poorly understood. The female-specific neo-W chromosome of the African monarch (or queen) butterfly Danaus chrysippus presents an intriguing case study because it is restricted to a single 'contact zone' population, involves a putative colour patterning supergene, and co-occurs with infection by the male-killing endosymbiont Spiroplasma. We investigated the origin and evolution of this system using whole genome sequencing. We first identify the 'BC supergene', a broad region of suppressed recombination across nearly half a chromosome, which links two colour patterning loci. Association analysis suggests that the genes yellow and arrow in this region control the forewing colour pattern differences between D. chrysippus subspecies. We then show that the same chromosome has recently formed a neo-W that has spread through the contact zone within approximately 2,200 years. We also assembled the genome of the male-killing Spiroplasma, and find that it shows perfect genealogical congruence with the neo-W, suggesting that the neo-W has hitchhiked to high frequency as the male-killer has spread through the population. The complete absence of female crossing-over in the Lepidoptera causes whole-chromosome hitchhiking of a single neo-W haplotype, carrying a single allele of the BC supergene and dragging multiple non-synonymous mutations to high frequency. This has created a population of infected females that all carry the same recessive colour patterning allele, making the phenotypes of each successive generation highly dependent on uninfected male immigrants. Our findings show how hitchhiking can occur between the physically unlinked genomes of host and endosymbiont, with dramatic consequences.

Selective sweeps on novel and introgressed variation shape mimicry loci in a butterfly adaptive radiation

Natural selection leaves distinct signatures in the genome that can reveal the targets and history of adaptive evolution. By analysing high-coverage genome sequence data from 4 major colour pattern loci sampled from nearly 600 individuals in 53 populations, we show pervasive selection on wing patterns in the Heliconius adaptive radiation. The strongest signatures correspond to loci with the greatest phenotypic effects, consistent with visual selection by predators, and are found in colour patterns with geographically restricted distributions. These recent sweeps are similar between co-mimics and indicate colour pattern turn-over events despite strong stabilising selection. Using simulations, we compare sweep signatures expected under classic hard sweeps with those resulting from adaptive introgression, an important aspect of mimicry evolution in Heliconius butterflies. Simulated recipient populations show a distinct 'volcano' pattern with peaks of increased genetic diversity around the selected target, characteristic of sweeps of introgressed variation and consistent with diversity patterns found in some populations. Our genomic data reveal a surprisingly dynamic history of colour pattern selection and co-evolution in this adaptive radiation.

Dissection of Larval and Pupal Wings of Bicyclus anynana Butterflies

The colorful wings of butterflies are emerging as model systems for evolutionary and developmental studies. Some of these studies focus on localizing gene transcripts and proteins in wings at the larval and pupal stages using techniques such as immunostaining and in situ hybridization. Other studies quantify mRNA expression levels or identify regions of open chromatin that are bound by proteins at different stages of wing development. All these techniques require dissection of the wings from the animal but a detailed video protocol describing this procedure has not been available until now. Here, we present a written and accompanying video protocol where we describe the tools and the method we use to remove the larval and pupal wings of the African Squinting Bush Brown butterfly Bicyclus anynana. This protocol should be easy to adapt to other species.

The dynamics of cyanide defences in the life cycle of an aposematic butterfly: Biosynthesis versus sequestration

Heliconius butterflies are highly specialized in Passiflora plants, laying eggs and feeding as larvae only on them. Interestingly, both Heliconius butterflies and Passiflora plants contain cyanogenic glucosides (CNglcs). While feeding on specific Passiflora species, Heliconius melpomene larvae are able to sequester simple cyclopentenyl CNglcs, the most common CNglcs in this plant genus. Yet, aromatic, aliphatic, and modified CNglcs have been reported in Passiflora species and they were never tested for sequestration by heliconiine larvae. As other cyanogenic lepidopterans, H. melpomene also biosynthesize the aliphatic CNglcs linamarin and lotaustralin, and their toxicity does not rely exclusively on sequestration. Although the genes encoding the enzymes in the CNglc biosynthesis have not yet been biochemically characterized in butterflies, the cytochromes P450 CYP405A4, CYP405A5, CYP405A6 and CYP332A1 have been hypothesized to be involved in this pathway in H. melpomene. In this study, we determine how the CNglc composition and expression of the putative P450s involved in the biosynthesis of these compounds vary at different developmental stages of Heliconius butterflies. We also establish which kind of CNglcs H. melpomene larvae can sequester from Passiflora. By analysing the chemical composition of the haemolymph from larvae fed with different Passiflora diets, we show that H. melpomene is able to sequestered prunasin, an aromatic CNglcs, from P. platyloba. They are also able to sequester amygdalin, gynocardin, [C¹³/C¹⁴]linamarin and [C¹³/C¹⁴]lotaustralin painted on the plant leaves. The CNglc tetraphyllin B-sulphate from P. caerulea is not detected in the larval haemolymph, suggesting that such modified CNglcs cannot be sequestered by Heliconius. Although pupae and virgin adults contain dihydrogynocardin resulting from larval sequestration, this compound was metabolized during adulthood, and not used as nuptial gift or transferred to the offspring. Thus, we speculate that dihydrogynocardin is catabolized to recycle nitrogen and glucose, and/or to produce fitness signals during courtship. Mature adults have a higher concentration of CNglcs than any other developmental stages due to increased de novo biosynthesis of linamarin and lotaustralin. Accordingly, all CYP405As are expressed in adults, whereas larvae mostly express CYP405A4. Our results shed light on the importance of CNglcs for Heliconius biology and their coevolution with Passiflora.

Molecular mechanisms underlying simplification of venation patterns in holometabolous insects

How mechanisms of pattern formation evolve has remained a central research theme in the field of evolutionary and developmental biology. The mechanism of wing vein differentiation in Drosophila is a classic text-book example of pattern formation using a system of positional-information, yet very little is known about how species with a different number of veins pattern their wings, and how insect venation patterns evolved. Here, we examine the expression pattern of genes previously implicated in vein differentiation in Drosophila in two butterfly species with more complex venation Bicyclus anynana and Pieris canidia. We also test the function of some of these genes in B. anynana. We identify both conserved as well as new domains of decapentaplegic, engrailed, invected, spalt, optix, wingless, armadillo, blistered, and rhomboid gene expression in butterflies, and propose how the simplified venation in Drosophila might have evolved via loss of decapentaplegic, spalt and optix gene expression domains, silencing of vein inducing programs at Spalt-expression boundaries, and changes in gene expression of vein maintenance genes.

Evolutionary genetics in insect phenotypic radiations: the value of a comparative genomic approach

Expanding genome sequencing and transgenic technologies are enabling the discovery of genes driving phenotypic diversity across insect taxa. Limitations in downstream functional genetic approaches, however, have been an obstacle for developing non-model systems for evolutionary genetics. Phenotypically diverse radiations, such as those exhibiting convergence and divergence as a result of mimicry, are ideal for evolutionary genetics as they can lead to insights using comparative genomic approaches alone. The varied and repeated instances of phenotypes in highly polymorphic systems allow assessment of whether similar loci are repeatedly targeted by selection and can inform how alleles sort across lineages. Comparative genomics of these taxa can be used to decipher components of gene regulatory networks, dissect regulatory regions, and validate genes.

Parallel evolution of ancient, pleiotropic enhancers underlies butterfly wing pattern mimicry

Color pattern mimicry in Heliconius butterflies is a classic case study of complex trait adaptation via selection on a few large effect genes. Association studies have linked color pattern variation to a handful of noncoding regions, yet the presumptive cis-regulatory elements (CREs) that control color patterning remain unknown. Here we combine chromatin assays, DNA sequence associations, and genome editing to functionally characterize 5 cis-regulatory elements of the color pattern gene optix We were surprised to find that the cis-regulatory architecture of optix is characterized by pleiotropy and regulatory fragility, where deletion of individual cis-regulatory elements has broad effects on both color pattern and wing vein development. Remarkably, we found orthologous cis-regulatory elements associate with wing pattern convergence of distantly related comimics, suggesting that parallel coevolution of ancestral elements facilitated pattern mimicry. Our results support a model of color pattern evolution in Heliconius where changes to ancient, multifunctional cis-regulatory elements underlie adaptive radiation.

Morphological Murals: The Scaling and Allometry of Butterfly Wing Patterns

The color patterns of butterflies moths are exceptionally diverse, but are very stable within a species, so that most species can be identified on the basis of their color pattern alone. The color pattern is established in the wing imaginal disc during a prolonged period of growth and differentiation, beginning during the last larval instar and ending during the first few days of the pupal stage. During this period, a variety of diffusion and reaction-diffusion signaling mechanisms determine the positions and sizes of the various elements that make up the overall color pattern. The patterning occurs while the wing is growing from a small imaginal disc to a very large pupal wing. One would therefore expect that some or all aspects of the color pattern would be sensitive to the size of the developmental field on which pattern formation takes place. To study this possibility, we analyzed the color patterns of Junonia coenia from animals whose growth patterns were altered by periodic starvation during larval growth, which produced individuals with a large range of variation in body size and wing size. Analyses of the color patterns showed that the positions and size of the pattern elements scaled perfectly isometrically with wing size. This is a puzzling finding and suggests the operation of a homeostatic or robustness mechanism that stabilizes pattern in spite of variation in the growth rate and final size of the wing.

Parallel genomic architecture underlies repeated sexual signal divergence in Hawaiian Laupala crickets

When the same phenotype evolves repeatedly, we can explore the predictability of genetic changes underlying phenotypic evolution. Theory suggests that genetic parallelism is less likely when phenotypic changes are governed by many small-effect loci compared to few of major effect, because different combinations of genetic changes can result in the same quantitative outcome. However, some genetic trajectories might be favoured over others, making a shared genetic basis to repeated polygenic evolution more likely. To examine this, we studied the genetics of parallel male mating song evolution in the Hawaiian cricket Laupala. We compared quantitative trait loci (QTL) underlying song divergence in three species pairs varying in phenotypic distance. We tested whether replicated song divergence between species involves the same QTL and whether the likelihood of QTL sharing is related to QTL effect size. Contrary to theoretical predictions, we find substantial parallelism in polygenic genetic architectures underlying repeated song divergence. QTL overlapped more frequently than expected based on simulated QTL analyses. Interestingly, QTL effect size did not predict QTL sharing, but did correlate with magnitude of phenotypic divergence. We highlight potential mechanisms driving these constraints on cricket song evolution and discuss a scenario that consolidates empirical quantitative genetic observations with micro-mutational theory.

Natural variations of growth thermo-responsiveness determined by SAUR26/27/28 proteins in Arabidopsis thaliana

How diversity in growth thermo-responsiveness is generated for local adaptation is a long-standing biological question. We investigated molecular genetic basis of natural variations in thermo-responsiveness of plant architecture in Arabidopsis thaliana. We measured the extent of rosette architecture at 22°C and 28°C in a set of 69 natural accessions and determined their thermo-responsiveness of plant architecture. A genome-wide association study was performed to identify major loci for variations in thermo-responsiveness. The SAUR26 subfamily, a new subfamily of SAUR genes, was identified as a major locus for the thermo-responsive architecture variations. The expression of SAUR26/27/28 is modulated by temperature and PIF4. Extensive natural polymorphisms in these genes affect their RNA expression levels and protein activities and influence the thermo-responsiveness of plant architecture. In addition, the SAUR26 subfamily genes exhibit a high variation frequency and their variations are associated with the local temperature climate. This study reveals that the SAUR26 subfamily is a key variation for thermo-responsive architecture and suggests a preference for generating diversity for local adaptation through signaling connectors.

Wing scale ultrastructure underlying convergent and divergent iridescent colours in mimetic Heliconius butterflies

Iridescence is an optical phenomenon whereby colour changes with the illumination and viewing angle. It can be produced by thin film interference or diffraction. Iridescent optical structures are fairly common in nature, but relatively little is known about their production or evolution. Here we describe the structures responsible for producing blue-green iridescent colour in Heliconius butterflies. Overall the wing scale structures of iridescent and non-iridescent Heliconius species are very similar, both having longitudinal ridges joined by cross-ribs. However, iridescent scales have ridges composed of layered lamellae, which act as multilayer reflectors. Differences in brightness between species can be explained by the extent of overlap of the lamellae and their curvature as well as the density of ridges on the scale. Heliconius are well known for their Müllerian mimicry. We find that iridescent structural colour is not closely matched between co-mimetic species. Differences appear less pronounced in models of Heliconius vision than models of avian vision, suggesting that they are not driven by selection to avoid heterospecific courtship by co-mimics. Ridge profiles appear to evolve relatively slowly, being similar between closely related taxa, while ridge density evolves faster and is similar between distantly related co-mimics.

Genomic evidence of human selection on Vavilovian mimicry

Vavilovian mimicry is an evolutionary process by which weeds evolve to resemble domesticated crop plants and is thought to be the result of unintentional selection by humans. Unravelling its molecular mechanisms will extend our knowledge of mimicry and contribute to our understanding of the origin and evolution of agricultural weeds, an important component of crop biology. To this end, we compared mimetic and non-mimetic populations of Echinochloa crus-galli from the Yangtze River basin phenotypically and by genome resequencing, and we show that this weed in rice paddies has evolved a small tiller angle, allowing it to phenocopy cultivated rice at the seedling stage. We demonstrate that mimetic lines evolved from the non-mimetic population as recently as 1,000 yr ago and were subject to a genetic bottleneck, and that genomic regions containing 87 putative plant architecture-related genes (including LAZY1, a key gene controlling plant tiller angle) were under selection during the mimicry process. Our data provide genome-level evidence for the action of human selection on Vavilovian mimicry.

Unravelling the genes forming the wing pattern supergene in the polymorphic butterfly Heliconius numata

Background

Unravelling the genetic basis of polymorphic characters is central to our understanding of the origins and diversification of living organisms. Recently, supergenes have been implicated in a wide range of complex polymorphisms, from adaptive colouration in butterflies and fish to reproductive strategies in birds and plants. The concept of a supergene is now a hot topic in biology, and identification of its functional elements is needed to shed light on the evolution of highly divergent adaptive traits. Here, we apply different gene expression analyses to study the supergene P that controls polymorphism of mimetic wing colour patterns in the neotropical butterfly Heliconius numata.

Results

We performed de novo transcriptome assembly and differential expression analyses using high-throughput Illumina RNA sequencing on developing wing discs of different H. numata morphs. Within the P interval, 30 and 17 of the 191 transcripts were expressed differentially in prepupae and day-1 pupae, respectively. Among these is the gene cortex, known to play a role in wing pattern formation in Heliconius and other Lepidoptera. Our in situ hybridization experiments confirmed the relationship between cortex expression and adult wing patterns.

Conclusions

This study found the majority of genes in the P interval to be expressed in the developing wing discs during the critical stages of colour pattern formation, and detect drastic changes in expression patterns in multiple genes associated with structural variants. The patterns of expression of cortex only partially recapitulate the variation in adult phenotype, suggesting that the remaining phenotypic variation could be controlled by other genes within the P interval. Although functional studies on cortex are now needed to determine its exact developmental role, our results are in accordance with the classical supergene hypothesis, whereby several genes inherited together due to tight linkage control a major developmental switch.

Electronic supplementary material

The online version of this article (10.1186/s13227-019-0129-2) contains supplementary material, which is available to authorized users.

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.

Comparative Transcriptome Analyses Reveal the Role of Conserved Function in Electric Organ Convergence Across Electric Fishes

The independent origins of multiple electric organs (EOs) of fish are fascinating examples of convergent evolution. However, comparative transcriptomics of different electric fish lineages are scarce. In this study, we found that the gene expression of EOs and skeletal muscles from three lineages (Mormyroidea, Siluriformes, and Gymnotiformes) tended to cluster together based on the species of origin, irrespective of the organ from which they are derived. A pairwise comparison of differentially expressed genes (DEGs) revealed that no less than half of shared DEGs exhibited parallel expression differentiation, indicating conserved directionality of differential expression either in or between lineages, but only a few shared DEGs were identified across all focal species. Nevertheless, the functional enrichment analysis of DEGs indicated that there were more parallel gene expression changes at the level of pathways and biological functions. Therefore, we may conclude that there is no parallel evolution of the entire transcriptomes of EOs among different lineages. Further, our results support the hypothesis that it is not different genes but conserved biological functions that play a crucial role in the convergence of complex phenotypes. This study provides insight into the genetic basis underlying the EO convergent evolution; however, more studies in different cases will be needed to demonstrate whether this pattern can be extended to other cases to derive a general rule for convergent evolution.

Plasticity of Promoter-Core Sequences Allows Bacteria to Compensate for the Loss of a Key Global Regulatory Gene

Transcription regulatory networks (TRNs) are of central importance for both short-term phenotypic adaptation in response to environmental fluctuations and long-term evolutionary adaptation, with global regulatory genes often being targets of natural selection in laboratory experiments. Here, we combined evolution experiments, whole-genome resequencing, and molecular genetics to investigate the driving forces, genetic constraints, and molecular mechanisms that dictate how bacteria can cope with a drastic perturbation of their TRNs. The crp gene, encoding a major global regulator in Escherichia coli, was deleted in four different genetic backgrounds, all derived from the Long-Term Evolution Experiment (LTEE) but with different TRN architectures. We confirmed that crp deletion had a more deleterious effect on growth rate in the LTEE-adapted genotypes; and we showed that the ptsG gene, which encodes the major glucose-PTS transporter, gained CRP (cyclic AMP receptor protein) dependence over time in the LTEE. We then further evolved the four crp-deleted genotypes in glucose minimal medium, and we found that they all quickly recovered from their growth defects by increasing glucose uptake. We showed that this recovery was specific to the selective environment and consistently relied on mutations in the cis-regulatory region of ptsG, regardless of the initial genotype. These mutations affected the interplay of transcription factors acting at the promoters, changed the intrinsic properties of the existing promoters, or produced new transcription initiation sites. Therefore, the plasticity of even a single promoter region can compensate by three different mechanisms for the loss of a key regulatory hub in the E. coli TRN.

Conservation and flexibility in the gene regulatory landscape of heliconiine butterfly wings

BACKGROUND

Many traits evolve by cis-regulatory modification, by which changes to noncoding sequences affect the binding affinity for available transcription factors and thus modify the expression profile of genes. Multiple examples of cis-regulatory evolution have been described at pattern switch genes responsible for butterfly wing pattern polymorphism, including in the diverse neotropical genus Heliconius, but the identities of the factors that can regulate these switch genes have not been identified.

RESULTS

We investigated the spatial transcriptomic landscape across the wings of three closely related butterfly species, two of which have a convergently evolved co-mimetic pattern and the other having a divergent pattern. We identified candidate factors for regulating the expression of wing patterning genes, including transcription factors with a conserved expression profile in all three species, and others, including both transcription factors and Wnt pathway genes, with markedly different profiles in each of the three species. We verified the conserved expression profile of the transcription factor homothorax by immunofluorescence and showed that its expression profile strongly correlates with that of the selector gene optix in butterflies with the Amazonian forewing pattern element 'dennis.'

CONCLUSION

Here we show that, in addition to factors with conserved expression profiles like homothorax, there are also a variety of transcription factors and signaling pathway components that appear to vary in their expression profiles between closely related butterfly species, highlighting the importance of genome-wide regulatory evolution between species.

The mimetic wing pattern of Papilio polytes butterflies is regulated by a doublesex-orchestrated gene network

The swallowtail butterfly Papilio polytes is sexually dimorphic and exhibits female-limited Batesian mimicry. This species also has two female forms, a non-mimetic form with male-like wing patterns, and a mimetic form resembling an unpalatable model, Pachliopta aristolochiae. The mimicry locus H constitutes a dimorphic Mendelian 'supergene', including a transcription factor gene doublesex (dsx). However, how the mimetic-type dsx (dsx-H) orchestrates the downstream gene network and causes the mimetic traits remains unclear. Here we performed RNA-seq-based gene screening and found that Wnt1 and Wnt6 are up-regulated by dsx-H during the early pupal stage and are involved in the red/white pigmentation and patterning of mimetic female wings. In contrast, a homeobox gene abdominal-A is repressed by dsx-H and involved in the non-mimetic colouration pattern. These findings suggest that dual regulation by dsx-H, induction of mimetic gene networks and repression of non-mimetic gene networks, is essential for the switch from non-mimetic to mimetic pattern in mimetic female wings.

A homeotic shift late in development drives mimetic color variation in a bumble bee

Natural phenotypic radiations, with their high diversity and convergence, are well-suited for informing how genomic changes translate to natural phenotypic variation. New genomic tools enable discovery in such traditionally nonmodel systems. Here, we characterize the genomic basis of color pattern variation in bumble bees (Hymenoptera, Apidae, Bombus), a group that has undergone extensive convergence of setal color patterns as a result of Müllerian mimicry. In western North America, multiple species converge on local mimicry patterns through parallel shifts of midabdominal segments from red to black. Using genome-wide association, we establish that a cis-regulatory locus between the abdominal fate-determining Hox genes, abd-A and Abd-B, controls the red-black color switch in a western species, Bombus melanopygus Gene expression analysis reveals distinct shifts in Abd-B aligned with the duration of setal pigmentation at the pupal-adult transition. This results in atypical anterior Abd-B expression, a late developmental homeotic shift. Changing expression of Hox genes can have widespread effects, given their important role across segmental phenotypes; however, the late timing reduces this pleiotropy, making Hox genes suitable targets. Analysis of this locus across mimics and relatives reveals that other species follow independent genetic routes to obtain the same phenotypes.

Digest: Ecomorphological convergence across the Atlantic

Did the remarkable helicopter damselflies (family Pseudostigmatidae) evolve their unique feeding and oviposition behaviors independently on two continents? In this issue, Toussaint et al. use molecular phylogenetic approaches to provide convincing evidence that these "forest giants" are in fact an example of ecomorphological convergence across the Atlantic Ocean.

Diversity in warning coloration: selective paradox or the norm?

Aposematic theory has historically predicted that predators should select for warning signals to converge on a single form, as a result of frequency-dependent learning. However, widespread variation in warning signals is observed across closely related species, populations and, most problematically for evolutionary biologists, among individuals in the same population. Recent research has yielded an increased awareness of this diversity, challenging the paradigm of signal monomorphy in aposematic animals. Here we provide a comprehensive synthesis of these disparate lines of investigation, identifying within them three broad classes of explanation for variation in aposematic warning signals: genetic mechanisms, differences among predators and predator behaviour, and alternative selection pressures upon the signal. The mechanisms producing warning coloration are also important. Detailed studies of the genetic basis of warning signals in some species, most notably Heliconius butterflies, are beginning to shed light on the genetic architecture facilitating or limiting key processes such as the evolution and maintenance of polymorphisms, hybridisation, and speciation. Work on predator behaviour is changing our perception of the predator community as a single homogenous selective agent, emphasising the dynamic nature of predator-prey interactions. Predator variability in a range of factors (e.g. perceptual abilities, tolerance to chemical defences, and individual motivation), suggests that the role of predators is more complicated than previously appreciated. With complex selection regimes at work, polytypisms and polymorphisms may even occur in Müllerian mimicry systems. Meanwhile, phenotypes are often multifunctional, and thus subject to additional biotic and abiotic selection pressures. Some of these selective pressures, primarily sexual selection and thermoregulation, have received considerable attention, while others, such as disease risk and parental effects, offer promising avenues to explore. As well as reviewing the existing evidence from both empirical studies and theoretical modelling, we highlight hypotheses that could benefit from further investigation in aposematic species. Finally by collating known instances of variation in warning signals, we provide a valuable resource for understanding the taxonomic spread of diversity in aposematic signalling and with which to direct future research. A greater appreciation of the extent of variation in aposematic species, and of the selective pressures and constraints which contribute to this once-paradoxical phenomenon, yields a new perspective for the field of aposematic signalling.

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.

Novel approach to quantitative spatial gene expression uncovers genetic stochasticity in the developing Drosophila eye

Robustness in development allows for the accumulation of genetically based variation in expression. However, this variation is usually examined in response to large perturbations, and examination of this variation has been limited to being spatial, or quantitative, but because of technical restrictions not both. Here we bridge these gaps by investigating replicated quantitative spatial gene expression using rigorous statistical models, in different genotypes, sexes, and species (Drosophila melanogaster and D. simulans). Using this type of quantitative approach with molecular developmental data allows for comparison among conditions, such as different genetic backgrounds. We apply this approach to the morphogenetic furrow, a wave of differentiation that patterns the developing eye disc. Within the morphogenetic furrow, we focus on four genes, hairy, atonal, hedgehog, and Delta. Hybridization chain reaction quantitatively measures spatial gene expression, co-staining for all four genes simultaneously. We find considerable variation in the spatial expression pattern of these genes in the eye between species, genotypes, and sexes. We also find that there has been evolution of the regulatory relationship between these genes, and that their spatial interrelationships have evolved between species. This variation has no phenotypic effect, and could be buffered by network thresholds or compensation from other genes. Both of these mechanisms could potentially be contributing to long term developmental systems drift.

Contingency in the convergent evolution of a regulatory network: Dosage compensation in Drosophila

The repeatability or predictability of evolution is a central question in evolutionary biology and most often addressed in experimental evolution studies. Here, we infer how genetically heterogeneous natural systems acquire the same molecular changes to address how genomic background affects adaptation in natural populations. In particular, we take advantage of independently formed neo-sex chromosomes in Drosophila species that have evolved dosage compensation by co-opting the dosage-compensation male-specific lethal (MSL) complex to study the mutational paths that have led to the acquisition of hundreds of novel binding sites for the MSL complex in different species. This complex recognizes a conserved 21-bp GA-rich sequence motif that is enriched on the X chromosome, and newly formed X chromosomes recruit the MSL complex by de novo acquisition of this binding motif. We identify recently formed sex chromosomes in the D. melanica and D. robusta species groups by genome sequencing and generate genomic occupancy maps of the MSL complex to infer the location of novel binding sites. We find that diverse mutational paths were utilized in each species to evolve hundreds of de novo binding motifs along the neo-X, including expansions of microsatellites and transposable element (TE) insertions. However, the propensity to utilize a particular mutational path differs between independently formed X chromosomes and appears to be contingent on genomic properties of that species, such as simple repeat or TE density. This establishes the "genomic environment" as an important determinant in predicting the outcome of evolutionary adaptations.

Genetic dissection of assortative mating behavior

The evolution of new species is made easier when traits under divergent ecological selection are also mating cues. Such ecological mating cues are now considered more common than previously thought, but we still know little about the genetic changes underlying their evolution or more generally about the genetic basis for assortative mating behaviors. Both tight physical linkage and the existence of large-effect preference loci will strengthen genetic associations between behavioral and ecological barriers, promoting the evolution of assortative mating. The warning patterns of Heliconius melpomene and H. cydno are under disruptive selection due to increased predation of nonmimetic hybrids and are used during mate recognition. We carried out a genome-wide quantitative trait locus (QTL) analysis of preference behaviors between these species and showed that divergent male preference has a simple genetic basis. We identify three QTLs that together explain a large proportion (approximately 60%) of the difference in preference behavior observed between the parental species. One of these QTLs is just 1.2 (0-4.8) centiMorgans (cM) from the major color pattern gene optix, and, individually, all three have a large effect on the preference phenotype. Genomic divergence between H. cydno and H. melpomene is high but broadly heterogenous, and admixture is reduced at the preference-optix color pattern locus but not the other preference QTLs. The simple genetic architecture we reveal will facilitate the evolution and maintenance of new species despite ongoing gene flow by coupling behavioral and ecological aspects of reproductive isolation.

Recombination rate variation shapes barriers to introgression across butterfly genomes

Hybridisation and introgression can dramatically alter the relationships among groups of species, leading to phylogenetic discordance across the genome and between populations. Introgression can also erode species differences over time, but selection against introgression at certain loci acts to maintain postmating species barriers. Theory predicts that species barriers made up of many loci throughout the genome should lead to a broad correlation between introgression and recombination rate, which determines the extent to which selection on deleterious foreign alleles will affect neutral alleles at physically linked loci. Here, we describe the variation in genealogical relationships across the genome among three species of Heliconius butterflies: H. melpomene (mel), H. cydno (cyd), and H. timareta (tim), using whole genomes of 92 individuals, and ask whether this variation can be explained by heterogeneous barriers to introgression. We find that species relationships vary predictably at the chromosomal scale. By quantifying recombination rate and admixture proportions, we then show that rates of introgression are predicted by variation in recombination rate. This implies that species barriers are highly polygenic, with selection acting against introgressed alleles across most of the genome. In addition, long chromosomes, which have lower recombination rates, produce stronger barriers on average than short chromosomes. Finally, we find a consistent difference between two species pairs on either side of the Andes, which suggests differences in the architecture of the species barriers. Our findings illustrate how the combined effects of hybridisation, recombination, and natural selection, acting at multitudes of loci over long periods, can dramatically sculpt the phylogenetic relationships among species.

The cellular and genetic basis of structural colour in plants

While the pathways that produce plant pigments have been well studied for decades, the use by plants of nanoscale structures to produce colour effects has only recently begun to be studied. A variety of plants from across the plant kingdom have been shown to use different mechanism to generate structural colours in tissues as diverse as leaves, flowers and fruits. In this review we explore the cellular mechanisms by which these nanoscale structures are built and discuss the first insights that have been published into the genetic pathways underpinning these traits.

Genomic Insights into the Adaptive Convergent Evolution

Adaptive convergent evolution, which refers to the same or similar phenotypes produced by species from independent lineages under similar selective pressures, has been widely examined for a long time. Accumulating studies on the adaptive convergent evolution have been reported from many different perspectives (cellular, anatomical, morphological, physiological, biochemical, and behavioral). Recent advances in the genomic technologies have demonstrated that adaptive convergence can arise from specific genetic mechanisms in different hierarchies, ranging from the same nucleotide or amino acid substitutions to the biological functions or pathways. Among these genetic mechanisms, the same amino acid changes in protein-coding genes play an important role in adaptive phenotypic convergence. Methods for detecting adaptive convergence at the protein sequence level have been constantly debated and developed. Here, we review recent progress on using genomic approaches to evaluate the genetic mechanisms of adaptive convergent evolution, summarize the research methods for identifying adaptive amino acid convergence, and discuss the future perspectives for researching adaptive convergent evolu-tion.

Searching for the genes driving assortative mating

Animals display an astonishing array of diverse colors and patterns, and animals also exhibit preferences for these diverse, species-specific traits when choosing a mate (i.e., assortative mate preference). It is hypothesized that in order for both preference and trait to be species specific, alleles for a trait and the preference for that trait must be inherited together and hence maintained as linked loci. This linkage could be maintained by three different genetic architectures: (A) the genes responsible for a species-specific preferred trait also directly influence preference for that trait; (B) genes producing preference and the preferred trait are not identical but are instead in close physical proximity in the genome; and (C) genes for preference and the preferred trait are nonadjacent but are inherited together due to selection. Merrill and colleagues test these hypotheses by performing large-scale genetic mapping of mating behavior using hybrids of two sympatric species of Heliconius butterflies, Heliconius melpomene and H. cydno. They identified three small genomic regions highly associated with mate preference, one of which was adjacent to a gene for the preferred trait, and two of which were not. Their findings illustrate that mate preference may be influenced by a small number of genes, while providing support for multiple hypotheses for the genetic architecture of assortative mate preferences.

Aristaless Controls Butterfly Wing Color Variation Used in Mimicry and Mate Choice

Neotropical Heliconius butterflies display a diversity of warningly colored wing patterns, which serve roles in both Müllerian mimicry and mate choice behavior. Wing pattern diversity in Heliconius is controlled by a small number of unlinked, Mendelian "switch" loci [1]. One of these, termed the K locus, switches between yellow and white color patterns, important mimicry signals as well as mating cues [2-4]. Furthermore, mate preference behavior is tightly linked to this locus [4]. K controls the distribution of white versus yellow scales on the wing, with a dominant white allele and a recessive yellow allele. Here, we combine fine-scale genetic mapping, genome-wide association studies, gene expression analyses, population and comparative genomics, and genome editing with CRISPR/Cas9 to characterize the molecular basis of the K locus in Heliconius and to infer its evolutionary history. We show that white versus yellow color variation in Heliconius cydno is due to alternate haplotypes at a putative cis-regulatory element (CRE) downstream of a tandem duplication of the homeodomain transcription factor aristaless. Aristaless1 (al1) and aristaless2 (al2) are differentially regulated between white and yellow wings throughout development with elevated expression of al1 in developing white wings, suggesting a role in repressing pigmentation. Consistent with this, knockout of al1 causes white wings to become yellow. The evolution of wing color in this group has been marked by retention of the ancestral yellow color in many lineages, a single origin of white coloration in H. cydno, and subsequent introgression of white color from H. cydno into H. melpomene.

The Genomic Basis of Color Pattern Polymorphism in the Harlequin Ladybird

Many animal species comprise discrete phenotypic forms. A common example in natural populations of insects is the occurrence of different color patterns, which has motivated a rich body of ecological and genetic research [1-6]. The occurrence of dark, i.e., melanic, forms displaying discrete color patterns is found across multiple taxa, but the underlying genomic basis remains poorly characterized. In numerous ladybird species (Coccinellidae), the spatial arrangement of black and red patches on adult elytra varies wildly within species, forming strikingly different complex color patterns [7, 8]. In the harlequin ladybird, Harmonia axyridis, more than 200 distinct color forms have been described, which classic genetic studies suggest result from allelic variation at a single, unknown, locus [9, 10]. Here, we combined whole-genome sequencing, population-based genome-wide association studies, gene expression, and functional analyses to establish that the transcription factor Pannier controls melanic pattern polymorphism in H. axyridis. We show that pannier is necessary for the formation of melanic elements on the elytra. Allelic variation in pannier leads to protein expression in distinct domains on the elytra and thus determines the distinct color patterns in H. axyridis. Recombination between pannier alleles may be reduced by a highly divergent sequence of ∼170 kb in the cis-regulatory regions of pannier, with a 50 kb inversion between color forms. This most likely helps maintain the distinct alleles found in natural populations. Thus, we propose that highly variable discrete color forms can arise in natural populations through cis-regulatory allelic variation of a single gene.

De Novo assembly and comparative transcriptome analyses of purple and green morphs of Apostichopus japonicus during body wall pigmentation process

Pigmentation processes provide a traceable and relevant trait for understanding key issues in evolutionary biology such as adaptation, speciation and the maintenance of balanced polymorphisms. The sea cucumber Apostichopus japonicus, which has nutritive and medical properties, is considered the most valuable commercial species in many parts of Asia. Compared with the green morph, the purple morph is rare and has great appeal to consumers. However, little is currently known about the molecular mechanism of body color formation in A. japonicus, even in echinoderm. Here, we employ illumina sequencing to examine expression patterns of the gene network underlying body wall development in purple and green morphs of A. japonicus. Overall, the number of down-regulated genes in the green morph was significantly more than in the purple morph during the pigmentation stage. We observed dynamic expression patterns of a large number of pigment, regulation and growth genes from the "Melanogenesis", "Melanoma", "Wnt signaling pathway", "Notch signaling pathway", "epithelium development", "epidermal growth factor receptor binding","growth factor activity" and "growth", including contrasting expression patterns of these genes in green and purple morph. This study provides comprehensive lists of differentially expressed genes during body wall development in the green and purple morphs, revealing potential candidate genes that may be involved in regulating body color formation and polymorphism. These data will provide valuable information for future genetic studies on sea cucumbers elucidating the molecular mechanisms underlying pigmentation, and may support the culture of desirable color morphs.

Different ommochrome pigment mixtures enable sexually dimorphic Batesian mimicry in disjunct populations of the common palmfly butterfly, Elymnias hypermnestra

With varied, brightly patterned wings, butterflies have been the focus of much work on the evolution and development of phenotypic novelty. However, the chemical structures of wing pigments from few butterfly species have been identified. We characterized the orange wing pigments of female Elymnias hypermnestra butterflies (Lepidoptera: Nymphalidae: Satyrinae) from two Southeast Asian populations. This species is a sexually dimorphic Batesian mimic of several model species. Females are polymorphic: in some populations, females are dark, resemble conspecific males, and mimic Euploea spp. In other populations, females differ from males and mimic orange Danaus spp. Using LC-MS/MS, we identified nine ommochrome pigments: six from a population in Chiang Mai, Thailand, and five compounds from a population in Bali, Indonesia. Two ommochromes were found in both populations, and only two of the nine compounds have been previously reported. The sexually dimorphic Thai and Balinese populations are separated spatially by monomorphic populations in peninsular Malaysia, Singapore, and Sumatra, suggesting independent evolution of mimetic female wing pigments in these disjunct populations. These results indicate that other butterfly wing pigments remain to be discovered.

Intrinsic cooperativity potentiates parallel cis-regulatory evolution

Convergent evolutionary events in independent lineages provide an opportunity to understand why evolution favors certain outcomes over others. We studied such a case where a large set of genes-those coding for the ribosomal proteins-gained cis-regulatory sequences for a particular transcription regulator (Mcm1) in independent fungal lineages. We present evidence that these gains occurred because Mcm1 shares a mechanism of transcriptional activation with an ancestral regulator of the ribosomal protein genes, Rap1. Specifically, we show that Mcm1 and Rap1 have the inherent ability to cooperatively activate transcription through contacts with the general transcription factor TFIID. Because the two regulatory proteins share a common interaction partner, the presence of one ancestral cis-regulatory sequence can 'channel' random mutations into functional sites for the second regulator. At a genomic scale, this type of intrinsic cooperativity can account for a pattern of parallel evolution involving the fixation of hundreds of substitutions.

A morphological novelty evolved by co-option of a reduced gene regulatory network and gene recruitment in a beetle

The mechanisms underlying the evolution of morphological novelties have remained enigmatic but co-option of existing gene regulatory networks (GRNs), recruitment of genes and the evolution of orphan genes have all been suggested to contribute. Here, we study a morphological novelty of beetle pupae called gin-trap. By combining the classical candidate gene approach with unbiased screening in the beetle Tribolium castaneum, we find that 70% of the tested components of the wing network were required for gin-trap development. However, many downstream and even upstream components were not included in the co-opted network. Only one gene was recruited from another biological context, but it was essential for the anteroposterior symmetry of the gin-traps, which represents a gin-trap-unique morphological innovation. Our data highlight the importance of co-option and modification of GRNs. The recruitment of single genes may not be frequent in the evolution of morphological novelties, but may be essential for subsequent diversification of the novelties. Finally, after having screened about 28% of annotated genes in the Tribolium genome to identify the genes required for gin-trap development, we found none of them are orphan genes, suggesting that orphan genes may have played only a minor, if any, role in the evolution of gin-traps.

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.

Parallel evolution of gene expression between trophic specialists despite divergent genotypes and morphologies

Parallel evolution of gene expression commonly underlies convergent niche specialization, but parallel changes in expression could also underlie divergent specialization. We investigated divergence in gene expression and whole-genome genetic variation across three sympatric Cyprinodon pupfishes endemic to San Salvador Island, Bahamas. This recent radiation consists of a generalist and two derived specialists adapted to novel niches: a scale-eating and a snail-eating pupfish. We sampled total mRNA from all three species at two early developmental stages and compared gene expression with whole-genome genetic differentiation among all three species in 42 resequenced genomes. Eighty percent of genes that were differentially expressed between snail-eaters and generalists were up or down regulated in the same direction between scale-eaters and generalists; however, there were no fixed variants shared between species underlying these parallel changes in expression. Genes showing parallel evolution of expression were enriched for effects on metabolic processes, whereas genes showing divergent expression were enriched for effects on cranial skeleton development and pigment biosynthesis, reflecting the most divergent phenotypes observed between specialist species. Our findings reveal that even divergent niche specialists may exhibit convergent adaptation to higher trophic levels through shared genetic pathways. This counterintuitive result suggests that parallel evolution in gene expression can accompany divergent ecological speciation during adaptive radiation.

Moving Speciation Genetics Forward: Modern Techniques Build on Foundational Studies in Drosophila

The question of how new species evolve has been examined at every level, from macroevolutionary patterns of diversification to molecular population genetic analyses of specific genomic regions between species pairs. Drosophila has been at the center of many of these research efforts. Though our understanding of the speciation process has grown considerably over the past few decades, very few genes have been identified that contribute to barriers to reproduction. The development of advanced molecular genetic and genomic methods provides promising avenues for the rapid discovery of more genes that contribute to speciation, particularly those involving prezygotic isolation. The continued expansion of tools and resources, especially for species other than Drosophila melanogaster, will be most effective when coupled with comparative approaches that reveal the genetic basis of reproductive isolation across a range of divergence times. Future research programs in Drosophila have high potential to answer long-standing questions in speciation. These include identifying the selective forces that contribute to divergence between populations and the genetic basis of traits that cause reproductive isolation. The latter can be expanded upon to understand how the genetic basis of reproductive isolation changes over time and whether certain pathways and genes are more commonly involved.

Extremely Divergent Haplotypes in Two Toxin Gene Complexes Encode Alternative Venom Types within Rattlesnake Species

Natural selection is generally expected to favor one form of a given trait within a population. The presence of multiple functional variants of traits involved in activities such as feeding, reproduction, or the defense against predators is relatively uncommon within animal species. The genetic architecture and evolutionary mechanisms underlying the origin and maintenance of such polymorphisms are of special interest. Among rattlesnakes, several instances of the production of biochemically distinct neurotoxic or hemorrhagic venom types within the same species are known. Here, we investigated the genetic basis of this phenomenon in three species and found that neurotoxic and hemorrhagic individuals of the same species possess markedly different haplotypes at two toxin gene complexes. For example, neurotoxic and hemorrhagic Crotalus scutulatus individuals differ by 5 genes at the phospholipase A2 (PLA2) toxin gene complex and by 11 genes at the metalloproteinase (MP) gene complex. A similar set of extremely divergent haplotypes also underlies alternate venom types within C. helleri and C. horridus. We further show that the MP and PLA2 haplotypes of neurotoxic C. helleri appear to have been acquired through hybridization with C. scutulatus-a rare example of the horizontal transfer of a potentially highly adaptive suite of genes. These large structural variants appear analogous to immunity gene complexes in host-pathogen arms races and may reflect the impact of balancing selection at the PLA2 and MP complexes for predation on different prey.