Genes controlling mimetic colour pattern variation in butterflies

Genetic Modification of a Hox Locus Drives Mimetic Color Pattern Variation in a Highly Polymorphic Bumble Bee

Müllerian mimicry provides natural replicates ideal for exploring mechanisms underlying adaptive phenotypic divergence and convergence, yet the genetic mechanisms underlying mimetic variation remain largely unknown. The current study investigates the genetic basis of mimetic color pattern variation in a highly polymorphic bumble bee, Bombus breviceps (Hymenoptera, Apidae). In South Asia, this species and multiple comimetic species converge onto local Müllerian mimicry patterns by shifting the abdominal setal color from orange to black. Genetic crossing between the orange and black phenotypes suggested the color dimorphism being controlled by a single Mendelian locus, with the orange allele being dominant over black. Genome-wide association suggests that a locus at the intergenic region between 2 abdominal fate-determining Hox genes, abd-A and Abd-B, is associated with the color change. This locus is therefore in the same intergenic region but not the same exact locus as found to drive red black midabdominal variation in a distantly related bumble bee species, Bombus melanopygus. Gene expression analysis and RNA interferences suggest that differential expression of an intergenic long noncoding RNA between abd-A and Abd-B at the onset setal color differentiation may drive the orange black color variation by causing a homeotic shift late in development. Analysis of this same color locus in comimetic species reveals no sequence association with the same color shift, suggesting that mimetic convergence is achieved through distinct genetic routes. Our study establishes Hox regions as genomic hotspots for color pattern evolution in bumble bees and demonstrates how pleiotropic developmental loci can drive adaptive radiations in nature.

Functional analysis of ABCG2 gene in pigment transport of Neocaridina denticulata sinensis

Many strains of Neocaridina denticulata sinensis (cherry shrimp) possess vivid body colors.However, the molecular underpinnings of these various body colors are scarcely understood. To study the role of the ABCG2 gene in the pigmentation of cherry shrimp, four strains (red, yellow, blue, and wild strains) were sampled. The sequence and expression pattern of ABCG2in tissues and embryos were analyzed, and the distribution of ABCG2 was also explored via WFISH (whole mount flurescence in situ hybridization). And further, RNA interference (RNAi) was used to explore the role of ABCG2 in body color deposition. The results showed that the ABCG2 sequence contained the conserved motif of Walker A, Walker B, Q-loop, d-loop, and H-loop. In tissues, ABCG2 was highly expressed in the epidermis of the four strains. During development stages, it was first expressed at the cleavage stage, then decreased at the gastrula stage, with the lowest expression at the pre-nauplius stage. From the metanauplius stage, its expression level was significantly upregulated until it reached the highest level at the membrane-zoea stage. WFISH showed that ABCG2 was first expressed at the cleavage stage in all four strains, and its distribution was similar from the cleavage stage to the before-zoea stage. The strongest positive signals were observed at the membrane-zoea stage, which was consistent with the qPCR results. Moreover, at the membrane-zoea stage, the positive signals of the four strains were mainly distributed in the compound eye and appendages. RNAi knockdown of ABCG2 encumbered the development of compound eye pigment cells (CEPCs) and erythrophores. It had effects on the expressions of other genes related to body color. These results suggest that ABCG2 is involved in the development of compound eye and the proliferation of erythrophores. This study provides new ideas for the cultivation of novel body colors in cherry shrimp at the molecular level.

Repeated genetic adaptation to altitude in two tropical butterflies

Repeated evolution can provide insight into the mechanisms that facilitate adaptation to novel or changing environments. Here we study adaptation to altitude in two tropical butterflies, Heliconius erato and H. melpomene, which have repeatedly and independently adapted to montane habitats on either side of the Andes. We sequenced 518 whole genomes from altitudinal transects and found many regions differentiated between highland (~ 1200 m) and lowland (~ 200 m) populations. We show repeated genetic differentiation across replicate populations within species, including allopatric comparisons. In contrast, there is little molecular parallelism between the two species. By sampling five close relatives, we find that a large proportion of divergent regions identified within species have arisen from standing variation and putative adaptive introgression from high-altitude specialist species. Taken together our study supports a role for both standing genetic variation and gene flow from independently adapted species in promoting parallel local adaptation to the environment.

Here, the authors study adaptation to altitude in 518 whole genomes from two species of tropical butterflies. They find repeated genetic differentiation within species, little molecular parallelism between these species, and introgression from closely related species, concluding that standing genetic variation promotes parallel local adaptation.

The genetic basis of structural colour variation in mimetic Heliconius butterflies

Structural colours, produced by the reflection of light from ultrastructures, have evolved multiple times in butterflies. Unlike pigmentary colours and patterns, little is known about the genetic basis of these colours. Reflective structures on wing-scale ridges are responsible for iridescent structural colour in many butterflies, including the Müllerian mimics Heliconius erato and Heliconius melpomene. Here, we quantify aspects of scale ultrastructure variation and colour in crosses between iridescent and non-iridescent subspecies of both of these species and perform quantitative trait locus (QTL) mapping. We show that iridescent structural colour has a complex genetic basis in both species, with offspring from crosses having a wide variation in blue colour (both hue and brightness) and scale structure measurements. We detect two different genomic regions in each species that explain modest amounts of this variation, with a sex-linked QTL in H. erato but not H. melpomene. We also find differences between species in the relationships between structure and colour, overall suggesting that these species have followed different evolutionary trajectories in their evolution of structural colour. We then identify genes within the QTL intervals that are differentially expressed between subspecies and/or wing regions, revealing likely candidates for genes controlling structural colour formation.

This article is part of the theme issue ‘Genetic basis of adaptation and speciation: from loci to causative mutations’.

Tempo of degeneration across independently evolved non-recombining regions

Recombination is beneficial over the long term, allowing more effective selection. Despite long-term advantages of recombination, local recombination suppression can evolve and lead to genomic degeneration, in particular on sex chromosomes. Here, we investigated the tempo of degeneration in nonrecombining regions, that is, the function curve for the accumulation of deleterious mutations over time, leveraging on 22 independent events of recombination suppression identified on mating-type chromosomes of anther-smut fungi, including newly identified ones. Using previously available and newly generated high-quality genome assemblies of alternative mating types of 13 Microbotryum species, we estimated degeneration levels in terms of accumulation of nonoptimal codons and nonsynonymous substitutions in nonrecombining regions. We found a reduced frequency of optimal codons in the nonrecombining regions compared with autosomes, that was not due to less frequent GC-biased gene conversion or lower ancestral expression levels compared with recombining regions. The frequency of optimal codons rapidly decreased following recombination suppression and reached an asymptote after ca. 3 Ma. The strength of purifying selection remained virtually constant at d_N/d_S = 0.55, that is, at an intermediate level between purifying selection and neutral evolution. Accordingly, nonsynonymous differences between mating-type chromosomes increased linearly with stratum age, at a rate of 0.015 per My. We thus develop a method for disentangling effects of reduced selection efficacy from GC-biased gene conversion in the evolution of codon usage and we quantify the tempo of degeneration in nonrecombining regions, which is important for our knowledge on genomic evolution and on the maintenance of regions without recombination.

The impact of life stage and pigment source on the evolution of novel warning signal traits

Our understanding of how novel warning color traits evolve in natural populations is largely based on studies of reproductive stages and organisms with endogenously produced pigmentation. In these systems, genetic drift is often required for novel alleles to overcome strong purifying selection stemming from frequency‐dependent predation and positive assortative mating. Here, we integrate data from field surveys, predation experiments, population genomics, and phenotypic correlations to explain the origin and maintenance of geographic variation in a diet‐based larval pigmentation trait in the redheaded pine sawfly (Neodiprion lecontei), a pine‐feeding hymenopteran. Although our experiments confirm that N. lecontei larvae are indeed aposematic—and therefore likely to experience frequency‐dependent predation—our genomic data do not support a historical demographic scenario that would have facilitated the spread of an initially deleterious allele via drift. Additionally, significantly elevated differentiation at a known color locus suggests that geographic variation in larval color is currently maintained by selection. Together, these data suggest that the novel white morph likely spread via selection. However, white body color does not enhance aposematic displays, nor is it correlated with enhanced chemical defense or immune function. Instead, the derived white‐bodied morph is disproportionately abundant on a pine species with a reduced carotenoid content relative to other pine hosts, suggesting that bottom‐up selection via host plants may have driven divergence among populations. Overall, our results suggest that life stage and pigment source can have a substantial impact on the evolution of novel warning signals, highlighting the need to investigate diverse aposematic taxa to develop a comprehensive understanding of color variation in nature.

An accurate assignment test for extremely low-coverage whole-genome sequence data

Genomic assignment tests can provide important diagnostic biological characteristics, such as population of origin or ecotype. Yet, assignment tests often rely on moderate- to high-coverage sequence data that can be difficult to obtain for fields such as molecular ecology and ancient DNA. We have developed a novel approach that efficiently assigns biologically relevant information (i.e., population identity or structural variants such as inversions) in extremely low-coverage sequence data. First, we generate databases from existing reference data using a subset of diagnostic single nucleotide polymorphisms (SNPs) associated with a biological characteristic. Low-coverage alignment files are subsequently compared to these databases to ascertain allelic state, yielding a joint probability for each association. To assess the efficacy of this approach, we assigned haplotypes and population identity in Heliconius butterflies, Atlantic herring, and Atlantic cod using chromosomal inversion sites and whole-genome data. We scored both modern and ancient specimens, including the first whole-genome sequence data recovered from ancient Atlantic herring bones. The method accurately assigns biological characteristics, including population membership, using extremely low-coverage data (as low as 0.0001x) based on genome-wide SNPs. This approach will therefore increase the number of samples in evolutionary, ecological and archaeological research for which relevant biological information can be obtained.

Complex basis of hybrid female sterility and Haldane's rule in Heliconius butterflies: Z-linkage and epistasis

Hybrids between species are often sterile or inviable. Hybrid unfitness usually evolves first in the heterogametic sex-a pattern known as Haldane's rule. The genetics of Haldane's rule have been extensively studied in species where the male is the heterogametic (XX/XY) sex, but its basis in taxa where the female is heterogametic (ZW/ZZ), such as Lepidoptera and birds, is largely unknown. Here, we analyse a new case of female hybrid sterility between geographic subspecies of Heliconius pardalinus. The two subspecies mate freely in captivity, but female F1 hybrids in both directions of cross are sterile. Sterility is due to arrested development of oocytes after they become differentiated from nurse cells, but before yolk deposition. We backcrossed fertile male F1 hybrids to parental females and mapped quantitative trait loci (QTLs) for female sterility. We also identified genes differentially expressed in the ovary as a function of oocyte development. The Z chromosome has a major effect, similar to the 'large X effect' in Drosophila, with strong epistatic interactions between loci at either end of the Z chromosome, and between the Z chromosome and autosomal loci on chromosomes 8 and 20. By intersecting the list of genes within these QTLs with those differentially expressed in sterile and fertile hybrids, we identified three candidate genes with relevant phenotypes. This study is the first to characterize hybrid sterility using genome mapping in the Lepidoptera and shows that it is produced by multiple complex epistatic interactions often involving the sex chromosome, as predicted by the dominance theory of Haldane's rule.

Cortex cis-regulatory switches establish scale colour identity and pattern diversity in Heliconius

In Heliconius butterflies, wing colour pattern diversity and scale types are controlled by a few genes of large effect that regulate colour pattern switches between morphs and species across a large mimetic radiation. One of these genes, cortex, has been repeatedly associated with colour pattern evolution in butterflies. Here we carried out CRISPR knockouts in multiple Heliconius species and show that cortex is a major determinant of scale cell identity. Chromatin accessibility profiling and introgression scans identified cis-regulatory regions associated with discrete phenotypic switches. CRISPR perturbation of these regions in black hindwing genotypes recreated a yellow bar, revealing their spatially limited activity. In the H. melpomene/timareta lineage, the candidate CRE from yellow-barred phenotype morphs is interrupted by a transposable element, suggesting that cis-regulatory structural variation underlies these mimetic adaptations. Our work shows that cortex functionally controls scale colour fate and that its cis-regulatory regions control a phenotypic switch in a modular and pattern-specific fashion.

Many ways to make darker flies: Intra- and interspecific variation in Drosophila body pigmentation components

Body pigmentation is an evolutionarily diversified and ecologically relevant trait with substantial variation within and between species, and important roles in animal survival and reproduction. Insect pigmentation, in particular, provides some of the most compelling examples of adaptive evolution, including its ecological significance and genetic bases. Pigmentation includes multiple aspects of color and color pattern that may vary more or less independently, and can be under different selective pressures. We decompose Drosophila thorax and abdominal pigmentation, a valuable eco-evo-devo model, into distinct measurable traits related to color and color pattern. We investigate intra- and interspecific variation for those traits and assess its different sources. For each body part, we measured overall darkness, as well as four other pigmentation properties distinguishing between background color and color of the darker pattern elements that decorate each body part. By focusing on two standard D. melanogaster laboratory populations, we show that pigmentation components vary and covary in distinct manners depending on sex, genetic background, and temperature during development. Studying three natural populations of D. melanogaster along a latitudinal cline and five other Drosophila species, we then show that evolution of lighter or darker bodies can be achieved by changing distinct component traits. Our results paint a much more complex picture of body pigmentation variation than previous studies could uncover, including patterns of sexual dimorphism, thermal plasticity, and interspecific diversity. These findings underscore the value of detailed quantitative phenotyping and analysis of different sources of variation for a better understanding of phenotypic variation and diversification, and the ecological pressures and genetic mechanisms underlying them.

Extensive haplotypes are associated with population differentiation and environmental adaptability in Upland cotton (Gossypium hirsutum)

Three extensive eco-haplotypes associated with population differentiation and environmental adaptability in Upland cotton were identified, with A06_85658585, A08_43734499 and A06_113104285 considered the eco-loci for environmental adaptability. Population divergence is suggested to be the primary force driving the evolution of environmental adaptability in various species. Chromosome inversion increases reproductive isolation between subspecies and accelerates population divergence to adapt to new environments. Although modern cultivated Upland cotton (Gossypium hirsutum L.) has spread worldwide, the noticeable phenotypic differences still existed among cultivars grown in different areas. In recent years, the long-distance migration of cotton cultivation areas throughout China has demanded that breeders better understand the genetic basis of environmental adaptability in Upland cotton. Here, we integrated the genotypes of 419 diverse accessions, long-term environment-associated variables (EAVs) and environment-associated traits (EATs) to evaluate subgroup differentiation and identify adaptive loci in Upland cotton. Two highly divergent genomic regions were found on chromosomes A06 and A08, which likely caused by extensive chromosome inversions. The subgroups could be geographically classified based on distinct haplotypes in the divergent regions. A genome-wide association study (GWAS) also confirmed that loci located in these regions were significantly associated with environmental adaptability in Upland cotton. Our study first revealed the cause of population divergence in Upland cotton, as well as the consequences of variation in its environmental adaptability. These findings provide new insights into the genetic basis of environmental adaptability in Upland cotton, which could accelerate the development of molecular markers for adaptation to climate change in future cotton breeding.

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.

Drift and Directional Selection Are the Evolutionary Forces Driving Gene Expression Divergence in Eye and Brain Tissue of Heliconius Butterflies

Investigating gene expression evolution over micro- and macroevolutionary timescales will expand our understanding of the role of gene expression in adaptation and speciation. In this study, we characterized the evolutionary forces acting on gene expression levels in eye and brain tissue of five Heliconius butterflies with divergence times of ∼5-12 MYA. We developed and applied Brownian motion (BM) and Ornstein-Uhlenbeck (OU) models to identify genes whose expression levels are evolving through drift, stabilizing selection, or a lineage-specific shift. We found that 81% of the genes evolve under genetic drift. When testing for branch-specific shifts in gene expression, we detected 368 (16%) shift events. Genes showing a shift toward upregulation have significantly lower gene expression variance than those genes showing a shift leading toward downregulation. We hypothesize that directional selection is acting in shifts causing upregulation, since transcription is costly. We further uncovered through simulations that parameter estimation of OU models is biased when using small phylogenies and only becomes reliable with phylogenies having ≥ 50 taxa. Therefore, we developed a new statistical test based on BM to identify highly conserved genes (i.e., evolving under strong stabilizing selection), which comprised 3% of the orthoclusters. In conclusion, we found that drift is the dominant evolutionary force driving gene expression evolution in eye and brain tissue in Heliconius Nevertheless, the higher proportion of genes evolving under directional than under stabilizing selection might reflect species-specific selective pressures on vision and the brain that are necessary to fulfill species-specific requirements.

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.

OSF-Builder: A New Tool for Constructing and Representing Evolutionary Histories Involving Introgression

Introgression is an evolutionary process which provides an important source of innovation for evolution. Although various methods have been used to detect introgression, very few methods are currently available for constructing evolutionary histories involving introgression. In this article, we propose a new method for constructing such evolutionary histories whose starting point is a species forest (consisting of a collection of lineage trees, usually arising as a collection of clades or monophyletic groups in a species tree), and a gene tree for a specific allele of interest, or allele tree for short. Our method is based on representing introgression in terms of a certain "overlay" of the allele tree over the lineage trees, called an overlaid species forest (OSF). OSFs are similar to phylogenetic networks although a key difference is that they typically have multiple roots because each monophyletic group in the species tree has a different point of origin. Employing a new model for introgression, we derive an efficient algorithm for building OSFs called OSF-Builder that is guaranteed to return an optimal OSF in the sense that the number of potential introgression events is minimized. As well as using simulations to assess the performance of OSF-Builder, we illustrate its use on a butterfly data set in which introgression has been previously inferred. The OSF-Builder software is available for download from https://www.uea.ac.uk/computing/software/OSF-Builder.

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.

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.

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.

Ommochromes in invertebrates: biochemistry and cell biology

Ommochromes are widely occurring coloured molecules of invertebrates, arising from tryptophan catabolism through the so-called Tryptophan → Ommochrome pathway. They are mainly known to mediate compound eye vision, as well as reversible and irreversible colour patterning. Ommochromes might also be involved in cell homeostasis by detoxifying free tryptophan and buffering oxidative stress. These biological functions are directly linked to their unique chromophore, the phenoxazine/phenothiazine system. The most recent reviews on ommochrome biochemistry were published more than 30 years ago, since when new results on the enzymes of the ommochrome pathway, on ommochrome photochemistry as well as on their antiradical capacities have been obtained. Ommochromasomes are the organelles where ommochromes are synthesised and stored. Hence, they play an important role in mediating ommochrome functions. Ommochromasomes are part of the lysosome-related organelles (LROs) family, which includes other pigmented organelles such as vertebrate melanosomes. Ommochromasomes are unique because they are the only LRO for which a recycling process during reversible colour change has been described. Herein, we provide an update on ommochrome biochemistry, photoreactivity and antiradical capacities to explain their diversity and behaviour both in vivo and in vitro. We also highlight new biochemical techniques, such as quantum chemistry, metabolomics and crystallography, which could lead to major advances in their chemical and functional characterisation. We then focus on ommochromasome structure and formation by drawing parallels with the well-characterised melanosomes of vertebrates. The biochemical, genetic, cellular and microscopic tools that have been applied to melanosomes should provide important information on the ommochromasome life cycle. We propose LRO-based models for ommochromasome biogenesis and recycling that could be tested in the future. Using the context of insect compound eyes, we finally emphasise the importance of an integrated approach in understanding the biological functions of ommochromes.

Phenotypic variation in Heliconius erato crosses shows that iridescent structural colour is sex-linked and controlled by multiple genes

Bright, highly reflective iridescent colours can be seen across nature and are produced by the scattering of light from nanostructures. Heliconius butterflies have been widely studied for their diversity and mimicry of wing colour patterns. Despite iridescence evolving multiple times in this genus, little is known about the genetic basis of the colour and the development of the structures which produce it. Heliconius erato can be found across Central and South America, but only races found in western Ecuador and Colombia have developed blue iridescent colour. Here, we use crosses between iridescent and non-iridescent races of H. erato to study phenotypic variation in the resulting F₂ generation. Using measurements of blue colour from photographs, we find that iridescent structural colour is a quantitative trait controlled by multiple genes, with strong evidence for loci on the Z sex chromosome. Iridescence is not linked to the Mendelian colour pattern locus that also segregates in these crosses (controlled by the gene cortex). Small-angle X-ray scattering data show that spacing between longitudinal ridges on the scales, which affects the intensity of the blue reflectance, also varies quantitatively in F₂ crosses.

Evolution of flower color pattern through selection on regulatory small RNAs

Small RNAs (sRNAs) regulate genes in plants and animals. Here, we show that population-wide differences in color patterns in snapdragon flowers are caused by an inverted duplication that generates sRNAs. The complexity and size of the transcripts indicate that the duplication represents an intermediate on the pathway to microRNA evolution. The sRNAs repress a pigment biosynthesis gene, creating a yellow highlight at the site of pollinator entry. The inverted duplication exhibits steep clines in allele frequency in a natural hybrid zone, showing that the allele is under selection. Thus, regulatory interactions of evolutionarily recent sRNAs can be acted upon by selection and contribute to the evolution of phenotypic diversity.

Evolutionary Assembly of Communities in Butterfly Mimicry Rings

Species co-occurrence in ecological communities is thought to be influenced by multiple ecological and evolutionary processes, especially colonization and competition. However, effects of other interspecific interactions and evolutionary relationships are less explored. We examined evolutionary histories of community members and roles of mutualistic and parasitic interactions (Müllerian and Batesian mimicry, respectively) in the assembly of mimetic butterfly communities called mimicry rings in tropical forests of the Western Ghats, India. We found that Müllerian mimics were phylogenetically clustered, sharing aposematic signals due to common ancestry. On the other hand, Batesian mimics joined mimicry rings through convergent evolution and random phylogenetic assembly. Since the Western Ghats are a habitat island, we compared species diversity and composition in its mimicry rings with those of habitat mainland to test effects of biogeographic connectivity. The Western Ghats consisted of fewer mimicry rings and an overall smaller number of aposematic species and mimics compared to habitat mainland. The depauperate mimicry rings in the Western Ghats could have resulted from stochastic processes, reflecting their long temporal and spatial isolation and trickling colonization by the mimetic butterfly communities. These results highlight how evolutionary history, biogeographic isolation, and stochastic colonization influence the evolutionary assembly and diversity of ecological communities.

Developmental and evolutionary mechanisms shaping butterfly eyespots

Butterfly eyespots are visually compelling models to study the reciprocal interactions between evolutionary and developmental processes that shape phenotypic variation. They are evolutionarily diversified, ecologically relevant, and developmentally tractable, and have made key contributions to linking genotype, development, phenotype and fitness. Advances in the availability of analytical tools (e.g. gene editing and visualization techniques) and resources (e.g. genomic and transcriptomic data) are boosting the detailed dissection of the mechanisms underlying eyespot development and evolution. Here, we review current knowledge on the ecology, development, and evolution of butterfly eyespots, with focus on recent advances. We also highlight a number of unsolved mysteries in our understanding of the patterns and processes underlying the diversification of these structures.

Ecological Genetics: A Key Gene for Mimicry and Melanism

Mimicry and melanism in Lepidoptera provided the first convincing examples of natural selection in action. Genetic analysis has now shown that, surprisingly, mimicry in Heliconius butterflies and melanism in peppered moths are switched at precisely the same gene: cortex.