No abstract

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.

22The wing patterns of butterflies are an excellent system with which to study phenotypic evolution. The 23 incredibly diverse patterns are generated from an array of pigmented scales on a largely two-24 dimensional surface, resulting in a visibly tractable system for studying the evolution of pigmentation. 25In Heliconius butterflies, much of this diversity is controlled by a few genes of large effect that regulate 26 pattern switches between races and species across a large mimetic radiation. One of these genescortex 27 -has been repeatedly mapped in association with colour pattern evolution in both Heliconius and other 28Lepidoptera, but we lack functional data supporting its role in modulating wing patterns. Here we 29 carried out CRISPR knock-outs in multiple Heliconius species and show that cortex is a major 30 determinant of scale cell identity. Mutant wing clones lacking cortex showed shifts in colour identity, 31 with melanic and red scales acquiring a yellow or white state. These homeotic transformations include 32 changes in both pigmentation and scale ultrastructure, suggesting that cortex acts during early stages of 33 scale cell fate specification rather than during the deployment of effector genes. In addition, mutant 34 clones were observed across the entire wing surface, contrasting with other known Heliconius mimicry 35 loci that act in specific patterns. Cortex is known as a cell-cycle regulator that modulates mitotic entry 36in Drosophila, and we found the Cortex protein to accumulate in the nuclei of the polyploid scale 37 building cells of the butterfly wing epithelium, speculatively suggesting a connection between scale cell 38 endocycling and colour identity. In summary, and while its molecular mode of action remains 39 mysterious, we conclude that cortex played key roles in the diversification of lepidopteran wing patterns 40 in part due to its switch-like effects in scale identity across the entire wing surface. 41 42 3

Chromosomes are a central unit of genome organisation. One tenth of all described species on Earth are Lepidoptera, butterflies and moths, and these generally possess 31 holocentric chromosomes. However, a subset of lepidopteran species display dramatic variation in chromosome counts. By analysing 210 chromosomally-complete lepidopteran genomes, the largest analysis of eukaryotic chromosomal-level reference genomes to date, we show that the diverse karyotypes of extant species are derived from 32 ancestral linkage groups, which we term Merian elements. Merian elements have remained largely intact across 250 million years of evolution and diversification. Against this stable background, we identify eight independent lineages that have evaded constraint and undergone extensive reorganisation - either by numerous fissions or a combination of fusion and fission events. Outside these lineages, fusions are rare and fissions are rarer still. Fusions tend to involve small, repeat-rich Merian elements and/or the Z chromosome. Together, our results reveal the constraints on genome architecture in Lepidoptera and enable a deeper understanding of the importance of chromosomal rearrangements in shaping the evolution of eukaryotic genomes.

We present a genome assembly from an individual female Vanessa cardui (the painted lady; Arthropoda; Insecta; Lepidoptera; Nymphalidae). The genome sequence is 425 megabases in span. The majority of the assembly is scaffolded into 32 chromosomal pseudomolecules, with the W and Z sex chromosome assembled. Gene annotation of this assembly on Ensembl has identified 12,821 protein coding genes.

Butterflies and moths (Lepidoptera) are one of the most ecologically diverse and speciose insect orders, with more than 157,000 described species. However, the abundance and diversity of Lepidoptera are declining worldwide at an alarming rate. As few Lepidoptera are explicitly recognised as at risk globally, the need for conservation is neither mandated nor well-evidenced. Large-scale biodiversity genomics projects that take advantage of the latest developments in long-read sequencing technologies offer a valuable source of information. We here present a comprehensive, reference-free, whole-genome, multiple sequence alignment of 88 species of Lepidoptera. We show that the accuracy and quality of the alignment is influenced by the contiguity of the reference genomes analysed. We explored genomic signatures that might indicate conservation concern in these species. In our dataset, which is largely from Britain, many species, in particular moths, display low heterozygosity and a high level of inbreeding, reflected in medium (0.1 - 1 Mb) and long (> 1 Mb) runs of homozygosity. Many species with low inbreeding display a higher masked load, estimated from the sum of rejected substitution scores at heterozygous sites. Our study shows that the analysis of a single diploid genome in a comparative phylogenetic context can provide relevant genetic information to prioritise species for future conservation investigation, particularly for those with an unknown conservation status.

The Xerces Blue (Glaucopsyche xerces) is considered to be the first butterfly to become extinct at global scale in historical times. It was notable for its chalky lavender wings with conspicuous white spots on the ventral wings. The last individuals were collected in their restricted habitat, in the dunes near the Presidio military base in San Francisco, in 1941. We sequenced the genomes of four 80 to 100-year-old Xerces Blue, and seven historical and one modern specimens of its closest relative, the Silvery Blue (G. lygdamus). We compared these to a novel annotated genome of the Green-Underside Blue (G. alexis). Phylogenetic relationships inferred from complete mitochondrial genomes indicate that Xerces Blue was a distinct species that diverged from the Silvery Blue lineage at least 850,000 years ago. Using nuclear genomes, both species experienced population growth during the Eemian interglacial period, but the Xerces Blue decreased to a very low effective population size subsequently, a trend opposite to that observed in the Silvery Blue. Runs of homozygosity and deleterious load in the Xerces Blue were significantly greater than in the Silvery Blue, suggesting a higher incidence of inbreeding. These signals of population decline observed in Xerces Blue could be used to identify and monitor other insects threatened by human activities, whose extinction patterns are still not well known.

The Xerces Blue (Glaucopsyche xerces) is considered to be the first butterfly to become extinct at global scale in historical times. It was notable for its chalky lavender wings with conspicuous white spots on the ventral wings. The last individuals were collected in their restricted habitat, in the dunes near the Presidio military base in San Francisco, in 1941. To explore the demographic history of this iconic butterfly and to better understand why it went extinct, we sequenced at medium coverage the genomes of four 80 to 100-year-old Xerces Blue specimens and seven historic specimens of its closest relative, the Silvery Blue (G. lygdamus). We compared these to a novel annotated genome of the Green-Underside Blue (G. alexis). Phylogenetic relationships inferred from complete mitochondrial genomes indicate that Xerces Blue was a distinct species that diverged from the Silvery Blue lineage at least 850,000 years ago. Using nuclear genomes, we show that both species experienced population growth during the MIS 7 interglacial period, but the Xerces Blue decreased to a very low effective population size subsequently, a trend opposite to that observed in the Silvery Blue. Runs of homozygosity in the Xerces Blue were significantly greater than in the Silvery Blue, suggesting a higher incidence of inbreeding. In addition, the Xerces Blue carried a higher proportion of derived, putatively deleterious amino acid-changing alleles than the Silvery Blue. These results demonstrate that the Xerces Blue experienced more than 100 thousand years of population decline, prior to its human-induced final extinction.