Early origin and diverse phenotypic implementation of iridescent UV patterns for sexual signaling in pierid butterflies.

Iridescent ultraviolet (IUV) patterns on pierid butterfly wings are phenotypic adaptations commonly used as sexual signals, generated by scales with ultrastructural modifications. Pierid IUV patterns are sexually dichromatic, with reduced size in females, where conspicuous sexual signaling balances courtship against ecological predation. There have been no phylogenetic reconstructions of IUV within Pieridae and little morphological characterization of phenotypic diversity. Our genus-wide characterization of IUV revealed the uniform similarity of stacked lamellar ridges on the dorsal surface of cover scales. We tested a hypothesis of single versus multiple origins by reconstructing a phylogeny of 534 species (~43.2% described species), with all genera represented, and a trait matrix of 734 species (~59.4%) screened for IUV. A single, early dimorphic origin of IUV followed by several losses and gains received strong support, concluding that IUV patterns and structural coloration are old traits. Collectively, these results support the homology of IUV scales and patterns that diversified within several lineages, suggesting an interplay between female-mediated sexual selection and ecological predatory selection.

High level of novelty under the hood of convergent evolution.

Little is known about the extent to which species use homologous regulatory architectures to achieve phenotypic convergence. By characterizing chromatin accessibility and gene expression in developing wing tissues, we compared the regulatory architecture of convergence between a pair of mimetic butterfly species. Although a handful of color pattern genes are known to be involved in their convergence, our data suggest that different mutational paths underlie the integration of these genes into wing pattern development. This is supported by a large fraction of accessible chromatin being exclusive to each species, including the de novo lineage-specific evolution of a modular optix enhancer. These findings may be explained by a high level of developmental drift and evolutionary contingency that occurs during the independent evolution of mimicry.

Widespread Gene Expression Divergence in Butterfly Sensory Tissues Plays a Fundamental Role During Reproductive Isolation and Speciation.

Neotropical Heliconius butterflies are well known for their intricate behaviors and multiple instances of incipient speciation. Chemosensing plays a fundamental role in the life history of these groups of butterflies and in the establishment of reproductive isolation. However, chemical communication involves synergistic sensory and accessory functions, and it remains challenging to investigate the molecular mechanisms underlying behavioral differences. Here, we examine the gene expression profiles and genomic divergence of three sensory tissues (antennae, legs, and mouthparts) between sexes (females and males) and life stages (different adult stages) in two hybridizing butterflies, Heliconius melpomene and Heliconius cydno. By integrating comparative transcriptomic and population genomic approaches, we found evidence of widespread gene expression divergence, supporting a crucial role of sensory tissues in the establishment of species barriers. We also show that sensory diversification increases in a manner consistent with evolutionary divergence based on comparison with the more distantly related species Heliconius charithonia. The findings of our study strongly support the unique chemosensory function of antennae in all three species, the importance of the Z chromosome in interspecific divergence, and the nonnegligible role of nonchemosensory genes in the divergence of chemosensory tissues. Collectively, our results provide a genome-wide illustration of diversification in the chemosensory system under incomplete reproductive isolation, revealing strong molecular separation in the early stage of speciation. Here, we provide a unique perspective and relevant view of the genetic architecture (sensory and accessory functions) of chemosensing beyond the classic chemosensory gene families, leading to a better understanding of the magnitude and complexity of molecular changes in sensory tissues that contribute to the establishment of reproductive isolation and speciation.

A butterfly pan-genome reveals that a large amount of structural variation underlies the evolution of chromatin accessibility.

Despite insertions and deletions being the most common structural variants (SVs) found across genomes, not much is known about how much these SVs vary within populations and between closely related species, nor their significance in evolution. To address these questions, we characterized the evolution of indel SVs using genome assemblies of three closely related Heliconius butterfly species. Over the relatively short evolutionary timescales investigated, up to 18.0% of the genome was composed of indels between two haplotypes of an individual Heliconius charithonia butterfly and up to 62.7% included lineage-specific SVs between the genomes of the most distant species (11 Mya). Lineage-specific sequences were mostly characterized as transposable elements (TEs) inserted at random throughout the genome and their overall distribution was similarly affected by linked selection as single nucleotide substitutions. Using chromatin accessibility profiles (i.e., ATAC-seq) of head tissue in caterpillars to identify sequences with potential cis-regulatory function, we found that out of the 31,066 identified differences in chromatin accessibility between species, 30.4% were within lineage-specific SVs and 9.4% were characterized as TE insertions. These TE insertions were localized closer to gene transcription start sites than expected at random and were enriched for sites with significant resemblance to several transcription factor binding sites with known function in neuron development in Drosophila We also identified 24 TE insertions with head-specific chromatin accessibility. Our results show high rates of structural genome evolution that were previously overlooked in comparative genomic studies and suggest a high potential for structural variation to serve as raw material for adaptive evolution.

A genetic switch for male UV iridescence in an incipient species pair of sulphur butterflies.

Mating cues evolve rapidly and can contribute to species formation and maintenance. However, little is known about how sexual signals diverge and how this variation integrates with other barrier loci to shape the genomic landscape of reproductive isolation. Here, we elucidate the genetic basis of ultraviolet (UV) iridescence, a courtship signal that differentiates the males of Colias eurytheme butterflies from a sister species, allowing females to avoid costly heterospecific matings. Anthropogenic range expansion of the two incipient species established a large zone of secondary contact across the eastern United States with strong signatures of genomic admixtures spanning all autosomes. In contrast, Z chromosomes are highly differentiated between the two species, supporting a disproportionate role of sex chromosomes in speciation known as the large-X (or large-Z) effect. Within this chromosome-wide reproductive barrier, linkage mapping indicates that cis-regulatory variation of bric a brac (bab) underlies the male UV-iridescence polymorphism between the two species. Bab is expressed in all non-UV scales, and butterflies of either species or sex acquire widespread ectopic iridescence following its CRISPR knockout, demonstrating that Bab functions as a suppressor of UV-scale differentiation that potentiates mating cue divergence. These results highlight how a genetic switch can regulate a premating signal and integrate with other reproductive barriers during intermediate phases of speciation.

Novel Doublesex Duplication Associated with Sexually Dimorphic Development of Dogface Butterfly Wings.

Sexually dimorphic development is responsible for some of the most remarkable phenotypic variation found in nature. Alternative splicing of the transcription factor gene doublesex (dsx) is a highly conserved developmental switch controlling the expression of sex-specific pathways. Here, we leverage sex-specific differences in butterfly wing color pattern to characterize the genetic basis of sexually dimorphic development. We use RNA-seq, immunolocalization, and motif binding site analysis to test specific predictions about the role of dsx in the development of structurally based ultraviolet (UV) wing patterns in Zerene cesonia (Southern Dogface). Unexpectedly, we discover a novel duplication of dsx that shows a sex-specific burst of expression associated with the sexually dimorphic UV coloration. The derived copy consists of a single exon that encodes a DNA binding but no protein-binding domain and has experienced rapid amino-acid divergence. We propose the novel dsx paralog may suppress UV scale differentiation in females, which is supported by an excess of Dsx-binding sites at cytoskeletal and chitin-related genes with sex-biased expression. These findings illustrate the molecular flexibility of the dsx gene in mediating the differentiation of secondary sexual characteristics.

Selection and isolation define a heterogeneous divergence landscape between hybridizing Heliconius butterflies.

Hybridizing species provide a powerful system to identify the processes that shape genomic variation and maintain species boundaries. However, complex histories of isolation, gene flow, and selection often generate heterogeneous genomic landscapes of divergence that complicate reconstruction of the speciation history. Here, we explore patterns of divergence to reconstruct recent speciation in the erato clade of Heliconius butterflies. We focus on the genomic landscape of divergence across three contact zones of the species H. erato and H. himera. We show that these hybridizing species have an intermediate level of divergence in the erato clade, which fits with their incomplete levels of reproductive isolation. Using demographic modeling and the relationship between admixture and divergence with recombination rate variation, we reconstruct histories of gene flow, selection, and demographic change that explain the observed patterns of genomic divergence. We find that periods of isolation and selection within populations, followed by secondary contact with asymmetrical gene flow are key factors in shaping the heterogeneous genomic landscapes. Collectively, these results highlight the effectiveness of demographic modeling and recombination rate estimates to disentangling the distinct contributions of gene flow and selection to patterns of genomic divergence.

Balanced polymorphisms and their divergence in a Heliconius butterfly.

The evolution of mimicry in similarly defended prey is well described by the Müllerian mimicry theory, which predicts the convergence of warning patterns in order to gain the most protection from predators. However, despite this prediction, we can find great diversity of color patterns among Müllerian mimics such as Heliconius butterflies in the neotropics. Furthermore, some species have evolved the ability to maintain multiple distinct warning patterns in single populations, a phenomenon known as polymorphic mimicry. The adaptive benefit of these polymorphisms is questionable since variation from the most common warning patterns is expected to be disadvantageous as novel signals are punished by predators naive to them. In this study, we use artificial butterfly models throughout Central and South America to characterize the selective pressures maintaining polymorphic mimicry in Heliconius doris. Our results highlight the complexity of positive frequency-dependent selection, the principal selective pressure driving convergence among Müllerian mimics, and its impacts on interspecific variation of mimetic warning coloration. We further show how this selection regime can both limit and facilitate the diversification of mimetic traits.

Perfect mimicry between Heliconius butterflies is constrained by genetics and development.

Müllerian mimicry strongly exemplifies the power of natural selection. However, the exact measure of such adaptive phenotypic convergence and the possible causes of its imperfection often remain unidentified. Here, we first quantify wing colour pattern differences in the forewing region of 14 co-mimetic colour pattern morphs of the butterfly species Heliconius erato and Heliconius melpomene and measure the extent to which mimicking colour pattern morphs are not perfectly identical. Next, using gene-editing CRISPR/Cas9 KO experiments of the gene WntA, which has been mapped to colour pattern diversity in these butterflies, we explore the exact areas of the wings in which WntA affects colour pattern formation differently in H. erato and H. melpomene. We find that, while the relative size of the forewing pattern is generally nearly identical between co-mimics, the CRISPR/Cas9 KO results highlight divergent boundaries in the wing that prevent the co-mimics from achieving perfect mimicry. We suggest that this mismatch may be explained by divergence in the gene regulatory network that defines wing colour patterning in both species, thus constraining morphological evolution even between closely related species.

Genome Assembly of the Dogface Butterfly Zerene cesonia.

Comparisons of high-quality, reference butterfly, and moth genomes have been instrumental to advancing our understanding of how hybridization, and natural selection drive genomic change during the origin of new species and novel traits. Here, we present a genome assembly of the Southern Dogface butterfly, Zerene cesonia (Pieridae) whose brilliant wing colorations have been implicated in developmental plasticity, hybridization, sexual selection, and speciation. We assembled 266,407,278 bp of the Z. cesonia genome, which accounts for 98.3% of the estimated 271 Mb genome size. Using a hybrid approach involving Chicago libraries with Hi-Rise assembly and a diploid Meraculous assembly, the final haploid genome was assembled. In the final assembly, nearly all autosomes and the Z chromosome were assembled into single scaffolds. The largest 29 scaffolds accounted for 91.4% of the genome assembly, with the remaining ∼8% distributed among another 247 scaffolds and overall N50 of 9.2 Mb. Tissue-specific RNA-seq informed annotations identified 16,442 protein-coding genes, which included 93.2% of the arthropod Benchmarking Universal Single-Copy Orthologs (BUSCO). The Z. cesonia genome assembly had ∼9% identified as repetitive elements, with a transposable element landscape rich in helitrons. Similar to other Lepidoptera genomes, Z. cesonia showed a high conservation of chromosomal synteny. The Z. cesonia assembly provides a high-quality reference for studies of chromosomal arrangements in the Pierid family, as well as for population, phylo, and functional genomic studies of adaptation and speciation.

Genomic architecture and introgression shape a butterfly radiation.

We used 20 de novo genome assemblies to probe the speciation history and architecture of gene flow in rapidly radiating Heliconius butterflies. Our tests to distinguish incomplete lineage sorting from introgression indicate that gene flow has obscured several ancient phylogenetic relationships in this group over large swathes of the genome. Introgressed loci are underrepresented in low-recombination and gene-rich regions, consistent with the purging of foreign alleles more tightly linked to incompatibility loci. Here, we identify a hitherto unknown inversion that traps a color pattern switch locus. We infer that this inversion was transferred between lineages by introgression and is convergent with a similar rearrangement in another part of the genus. These multiple de novo genome sequences enable improved understanding of the importance of introgression and selective processes in adaptive radiation.

Interplay between Developmental Flexibility and Determinism in the Evolution of Mimetic Heliconius Wing Patterns.

To what extent can we predict how evolution occurs? Do genetic architectures and developmental processes canalize the evolution of similar outcomes in a predictable manner? Or do historical contingencies impose alternative pathways to answer the same challenge? Examples of Müllerian mimicry between distantly related butterfly species provide natural replicates of evolution, allowing us to test whether identical wing patterns followed parallel or novel trajectories. Here, we explore the role that the signaling ligand WntA plays in generating mimetic wing patterns in Heliconius butterflies, a group with extraordinary mimicry-related wing pattern diversity. The radiation is relatively young, and numerous cases of wing pattern mimicry have evolved within the last 2.5-4.5 Ma. WntA is an important target of natural selection and is one of four major effect loci that underlie much of the pattern variation in the group. We used CRISPR/Cas9 targeted mutagenesis to generate WntA-deficient wings in 12 species and a further 10 intraspecific variants, including three co-mimetic pairs. In all tested butterflies, WntA knockouts affect pattern broadly and cause a shift among every possible scale cell type. Interestingly, the co-mimics lacking WntA were very different, suggesting that the gene networks that pattern a wing have diverged considerably among different lineages. Thus, although natural selection channeled phenotypic convergence, divergent developmental contexts between the two major Heliconius lineages opened different developmental routes to evolve resemblance. Consequently, even under very deterministic evolutionary scenarios, our results underscore a surprising unpredictability in the developmental paths underlying convergence in a recent radiation.

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.

Comparative Transcriptomics Provides Insights into Reticulate and Adaptive Evolution of a Butterfly Radiation.

Butterfly eyes are complex organs that are composed of a diversity of proteins and they play a central role in visual signaling and ultimately, speciation, and adaptation. Here, we utilized the whole eye transcriptome to obtain a more holistic view of the evolution of the butterfly eye while accounting for speciation events that co-occur with ancient hybridization. We sequenced and assembled transcriptomes from adult female eyes of eight species representing all major clades of the Heliconius genus and an additional outgroup species, Dryas iulia. We identified 4,042 orthologous genes shared across all transcriptome data sets and constructed a transcriptome-wide phylogeny, which revealed topological discordance with the mitochondrial phylogenetic tree in the Heliconius pupal mating clade. We then estimated introgression among lineages using additional genome data and found evidence for ancient hybridization leading to the common ancestor of Heliconius hortense and Heliconius clysonymus. We estimated the Ka/Ks ratio for each orthologous cluster and performed further tests to demonstrate genes showing evidence of adaptive protein evolution. Furthermore, we characterized patterns of expression for a subset of these positively selected orthologs using qRT-PCR. Taken together, we identified candidate eye genes that show signatures of adaptive molecular evolution and provide evidence of their expression divergence between species, tissues, and sexes. Our results demonstrate: 1) greater evolutionary changes in younger Heliconius lineages, that is, more positively selected genes in the cydno-melpomene-hecale group as opposed to the sara-hortense-erato group, and 2) suggest an ancient hybridization leading to speciation among Heliconius pupal-mating species.

patternize: An R package for quantifying colour pattern variation.

The use of image data to quantify, study and compare variation in the colors and patterns of organisms requires the alignment of images to establish homology, followed by color-based segmentation of images. Here we describe an R package for image alignment and segmentation that has applications to quantify color patterns in a wide range of organisms. patternize is an R package that quantifies variation in color patterns obtained from image data. patternize first defines homology between pattern positions across specimens either through manually placed homologous landmarks or automated image registration. Pattern identification is performed by categorizing the distribution of colors using an RGB threshold, k-means clustering or watershed transformation.We demonstrate that patternize can be used for quantification of the color patterns in a variety of organisms by analyzing image data for butterflies, guppies, spiders and salamanders. Image data can be compared between sets of specimens, visualized as heatmaps and analyzed using principal component analysis (PCA). patternize has potential applications for fine scale quantification of color pattern phenotypes in population comparisons, genetic association studies and investigating the basis of color pattern variation across a wide range of organisms.

Patterns of Z chromosome divergence among Heliconius species highlight the importance of historical demography.

Sex chromosomes are disproportionately involved in reproductive isolation and adaptation. In support of such a "large-X" effect, genome scans between recently diverged populations and species pairs often identify distinct patterns of divergence on the sex chromosome compared to autosomes. When measures of divergence between populations are higher on the sex chromosome compared to autosomes, such patterns could be interpreted as evidence for faster divergence on the sex chromosome, that is "faster-X", barriers to gene flow on the sex chromosome. However, demographic changes can strongly skew divergence estimates and are not always taken into consideration. We used 224 whole-genome sequences representing 36 populations from two Heliconius butterfly clades (H. erato and H. melpomene) to explore patterns of Z chromosome divergence. We show that increased divergence compared to equilibrium expectations can in many cases be explained by demographic change. Among Heliconius erato populations, for instance, population size increase in the ancestral population can explain increased absolute divergence measures on the Z chromosome compared to the autosomes, as a result of increased ancestral Z chromosome genetic diversity. Nonetheless, we do identify increased divergence on the Z chromosome relative to the autosomes in parapatric or sympatric species comparisons that imply postzygotic reproductive barriers. Using simulations, we show that this is consistent with reduced gene flow on the Z chromosome, perhaps due to greater accumulation of incompatibilities. Our work demonstrates the importance of taking demography into account to interpret patterns of divergence on the Z chromosome, but nonetheless provides evidence to support the Z chromosome as a strong barrier to gene flow in incipient Heliconius butterfly species.

Complex modular architecture around a simple toolkit of wing pattern genes.

Identifying the genomic changes that control morphological variation and understanding how they generate diversity is a major goal of evolutionary biology. In Heliconius butterflies, a small number of genes control the development of diverse wing color patterns. Here, we used full genome sequencing of individuals across the Heliconius erato radiation and closely related species to characterize genomic variation associated with wing pattern diversity. We show that variation around color pattern genes is highly modular, with narrow genomic intervals associated with specific differences in color and pattern. This modular architecture explains the diversity of color patterns and provides a flexible mechanism for rapid morphological diversification.

Demographic Variation of Wolbachia Infection in the Endangered Mitchell's Satyr Butterfly.

The Mitchell's satyr, Neonympha mitchellii, is an endangered species that is limited to highly isolated habitats in the northern and southern United States. Conservation strategies for isolated endangered species often implement captive breeding and translocation programs for repopulation. However, these programs risk increasing the spread of harmful pathogens, such as the bacterial endosymbiont Wolbachia. Wolbachia can manipulate the host's reproduction leading to incompatibilities between infected and uninfected hosts. This study uses molecular methods to screen for Wolbachia presence across the distribution of the Mitchell's satyr and its subspecies, St. Francis satyr, which are both federally listed as endangered and are considered two of the rarest butterflies in North America. The screens confirmed the presence of Wolbachia in the northern and newly discovered southern populations of the Mitchell's satyr, but not in the St. Francis satyr population. These results combined with previous reports of Wolbachia in N. mitchellii, highlight that Wolbachia infection varies both geographically and temporally in satyr populations. The temporal variance shows the importance of continued monitoring of Wolbachia infection during conservation programs. To reduce the risk of reproductive incompatibilities, it is advised that all individuals collected for conservation purposes be screened for Wolbachia and recommended to avoid the use of infected individuals for captive breeding and translocation programs.