Cis-regulatory modes of Ultrabithorax inactivation in butterfly forewings.

Hox gene clusters encode transcription factors that drive regional specialization during animal development: for example the Hox factor Ubx is expressed in the insect metathoracic (T3) wing appendages and differentiates them from T2 mesothoracic identities. Hox transcriptional regulation requires silencing activities that prevent spurious activation and regulatory crosstalks in the wrong tissues, but this has seldom been studied in insects other than Drosophila, which shows a derived Hox dislocation into two genomic clusters that disjoined Antennapedia (Antp) and Ultrabithorax (Ubx). Here, we investigated how Ubx is restricted to the hindwing in butterflies, amidst a contiguous Hox cluster. By analysing Hi-C and ATAC-seq data in the butterfly Junonia coenia, we show that a Topologically Associated Domain (TAD) maintains a hindwing-enriched profile of chromatin opening around Ubx. This TAD is bordered by a Boundary Element (BE) that separates it from a region of joined wing activity around the Antp locus. CRISPR mutational perturbation of this BE releases ectopic Ubx expression in forewings, inducing homeotic clones with hindwing identities. Further mutational interrogation of two non-coding RNA encoding regions and one putative cis-regulatory module within the Ubx TAD cause rare homeotic transformations in both directions, indicating the presence of both activating and repressing chromatin features. We also describe a series of spontaneous forewing homeotic phenotypes obtained in Heliconius butterflies, and discuss their possible mutational basis. By leveraging the extensive wing specialization found in butterflies, our initial exploration of Ubx regulation demonstrates the existence of silencing and insulating sequences that prevent its spurious expression in forewings.

Evolutionary dynamics of genome size and content during the adaptive radiation of Heliconiini butterflies.

Heliconius butterflies, a speciose genus of Müllerian mimics, represent a classic example of an adaptive radiation that includes a range of derived dietary, life history, physiological and neural traits. However, key lineages within the genus, and across the broader Heliconiini tribe, lack genomic resources, limiting our understanding of how adaptive and neutral processes shaped genome evolution during their radiation. Here, we generate highly contiguous genome assemblies for nine Heliconiini, 29 additional reference-assembled genomes, and improve 10 existing assemblies. Altogether, we provide a dataset of annotated genomes for a total of 63 species, including 58 species within the Heliconiini tribe. We use this extensive dataset to generate a robust and dated heliconiine phylogeny, describe major patterns of introgression, explore the evolution of genome architecture, and the genomic basis of key innovations in this enigmatic group, including an assessment of the evolution of putative regulatory regions at the Heliconius stem. Our work illustrates how the increased resolution provided by such dense genomic sampling improves our power to generate and test gene-phenotype hypotheses, and precisely characterize how genomes evolve.

Frizzled2 receives WntA signaling during butterfly wing pattern formation.

Butterfly color patterns provide visible and biodiverse phenotypic readouts of the patterning processes. Although the secreted ligand WntA has been shown to instruct the color pattern formation in butterflies, its mode of reception remains elusive. Butterfly genomes encode four homologs of the Frizzled-family of Wnt receptors. Here, we show that CRISPR mosaic knockouts of frizzled2 (fz2) phenocopy the color pattern effects of WntA loss of function in multiple nymphalids. Whereas WntA mosaic clones result in intermediate patterns of reduced size, fz2 clones are cell-autonomous, consistent with a morphogen function. Shifts in expression of WntA and fz2 in WntA crispant pupae show that they are under positive and negative feedback, respectively. Fz1 is required for Wnt-independent planar cell polarity in the wing epithelium. Fz3 and Fz4 show phenotypes consistent with Wnt competitive-antagonist functions in vein formation (Fz3 and Fz4), wing margin specification (Fz3), and color patterning in the Discalis and Marginal Band Systems (Fz4). Overall, these data show that the WntA/Frizzled2 morphogen-receptor pair forms a signaling axis that instructs butterfly color patterning and shed light on the functional diversity of insect Frizzled receptors.

Deep cis-regulatory homology of the butterfly wing pattern ground plan.

Butterfly wing patterns derive from a deeply conserved developmental ground plan yet are diverse and evolve rapidly. It is poorly understood how gene regulatory architectures can accommodate both deep homology and adaptive change. To address this, we characterized the cis-regulatory evolution of the color pattern gene WntA in nymphalid butterflies. Comparative assay for transposase-accessible chromatin using sequencing (ATAC-seq) and in vivo deletions spanning 46 cis-regulatory elements across five species revealed deep homology of ground plan-determining sequences, except in monarch butterflies. Furthermore, noncoding deletions displayed both positive and negative regulatory effects that were often broad in nature. Our results provide little support for models predicting rapid enhancer turnover and suggest that deeply ancestral, multifunctional noncoding elements can underlie rapidly evolving trait systems.

Erratum: Mapping and CRISPR homology-directed repair of a recessive white eye mutation in Plodia moths.

[This corrects the article DOI: 10.1016/j.isci.2022.103885.].

Mapping and CRISPR homology-directed repair of a recessive white eye mutation in Plodia moths.

The pantry moth Plodia interpunctella is a worldwide pest of stored food products and a promising laboratory model system for lepidopteran functional genomics. Here we describe efficient methods for precise genome editing in this insect. A spontaneous recessive white-eyed phenotype maps to a frameshift deletion (c.737delC) in the white gene. CRISPR NHEJ mutagenesis of white replicates this phenotype with high rates of somatic biallelic knockout. G0 individuals with mutant clones on both eyes produced 100% mutant progeny, making white an ideal marker for co-conversion when targeting other genes. CRISPR HDR experiments corrected c.737delC and reverted white eyes to a pigmented state in 37% of G0 mosaic adults. These repaired alleles showed practical rates of germline transmission in backcrosses, demonstrating the potential of the technique for precise genome editing. Plodia offers a promising avenue for research in this taxon because of its lab-ready features, egg injectability, and editability.

Efficient hyperactive piggyBac transgenesis in Plodia pantry moths.

While piggyBac transposon-based transgenesis is widely used in various emerging model organisms, its relatively low transposition rate in butterflies and moths has hindered its use for routine genetic transformation in Lepidoptera. Here, we tested the suitability of a codon-optimized hyperactive piggyBac transposase (hyPBase) in mRNA form to deliver and integrate transgenic cassettes into the genome of the pantry moth Plodia interpunctella. Co-injection of hyPBase mRNA with donor plasmids successfully integrated 1.5-4.4 kb expression cassettes driving the fluorescent markers EGFP, DsRed, or EYFP in eyes and glia with the 3xP3 promoter. Somatic integration and expression of the transgene in the G0 injected generation was detectable from 72-h embryos and onward in larvae, pupae and adults carrying a recessive white-eyed mutation. Overall, 2.5% of injected eggs survived into transgene-bearing adults with mosaic fluorescence. Subsequent outcrossing of fluorescent G0 founders transmitted single-insertion copies of 3xP3::EGFP and 3xP3::EYFP and generated stable isogenic lines. Random in-crossing of a small cohort of G0 founders expressing 3xP3::DsRed yielded a stable transgenic line segregating for more than one transgene insertion site. We discuss how hyPBase can be used to generate stable transgenic resources in Plodia and other moths.

Molecular Phylogeny and Evolution of Amazon Parrots in the Greater Antilles.

Amazon parrots (Amazona spp.) colonized the islands of the Greater Antilles from the Central American mainland, but there has not been a consensus as to how and when this happened. Today, most of the five remaining island species are listed as endangered, threatened, or vulnerable as a consequence of human activity. We sequenced and annotated full mitochondrial genomes of all the extant Amazon parrot species from the Greater Antillean (A. leucocephala (Cuba), A. agilis, A. collaria (both from Jamaica), A. ventralis (Hispaniola), and A. vittata (Puerto Rico)), A. albifrons from mainland Central America, and A. rhodocorytha from the Atlantic Forest in Brazil. The assembled and annotated mitogenome maps provide information on sequence organization, variation, population diversity, and evolutionary history for the Caribbean species including the critically endangered A. vittata. Despite the larger number of available samples from the Puerto Rican Parrot Recovery Program, the sequence diversity of the A. vittata population in Puerto Rico was the lowest among all parrot species analyzed. Our data support the stepping-stone dispersal and speciation hypothesis that has started approximately 3.47 MYA when the ancestral population arrived from mainland Central America and led to diversification across the Greater Antilles, ultimately reaching the island of Puerto Rico 0.67 MYA. The results are presented and discussed in light of the geological history of the Caribbean and in the context of recent parrot evolution, island biogeography, and conservation. This analysis contributes to understating evolutionary history and empowers subsequent assessments of sequence variation and helps design future conservation efforts in the Caribbean.

Genomic architecture of a genetically assimilated seasonal color pattern.

Developmental plasticity allows genomes to encode multiple distinct phenotypes that can be differentially manifested in response to environmental cues. Alternative plastic phenotypes can be selected through a process called genetic assimilation, although the mechanisms are still poorly understood. We assimilated a seasonal wing color phenotype in a naturally plastic population of butterflies (Junonia coenia) and characterized three responsible genes. Endocrine assays and chromatin accessibility and conformation analyses showed that the transition of wing coloration from an environmentally determined trait to a predominantly genetic trait occurred through selection for regulatory alleles of downstream wing-patterning genes. This mode of genetic evolution is likely favored by selection because it allows tissue- and trait-specific tuning of reaction norms without affecting core cue detection or transduction mechanisms.

Multiple roles forlaccase2 in butterfly wing pigmentation, scale development, and cuticle tanning.

Lepidopteran wing scales play important roles in a number of functions including color patterning and thermoregulation. Despite the importance of wing scales, however, we still have a limited understanding of the genetic mechanisms that underlie scale patterning, development, and coloration. Here, we explore the function of the phenoloxidase-encoding gene laccase2 in wing and scale development in the nymphalid butterfly Vanessa cardui. Somatic deletion mosaics of laccase2 generated by CRISPR/Cas9 genome editing presented several distinct mutant phenotypes. Consistent with the work in other nonlepidopteran insect groups, we observed reductions in melanin pigmentation and defects in cuticle formation. We were also surprised, however, to see distinct effects on scale development including complete loss of wing scales. This study highlights laccase2 as a gene that plays multiple roles in wing and scale development and provides new insight into the evolution of lepidopteran wing coloration.

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.

Single master regulatory gene coordinates the evolution and development of butterfly color and iridescence.

The optix gene has been implicated in butterfly wing pattern adaptation by genetic association, mapping, and expression studies. The actual developmental function of this gene has remained unclear, however. Here we used CRISPR/Cas9 genome editing to show that optix plays a fundamental role in nymphalid butterfly wing pattern development, where it is required for determination of all chromatic coloration. optix knockouts in four species show complete replacement of color pigments with melanins, with corresponding changes in pigment-related gene expression, resulting in black and gray butterflies. We also show that optix simultaneously acts as a switch gene for blue structural iridescence in some butterflies, demonstrating simple regulatory coordination of structural and pigmentary coloration. Remarkably, these optix knockouts phenocopy the recurring "black and blue" wing pattern archetype that has arisen on many independent occasions in butterflies. Here we demonstrate a simple genetic basis for structural coloration, and show that optix plays a deeply conserved role in butterfly wing pattern development.

Macroevolutionary shifts of WntA function potentiate butterfly wing-pattern diversity.

Butterfly wing patterns provide a rich comparative framework to study how morphological complexity develops and evolves. Here we used CRISPR/Cas9 somatic mutagenesis to test a patterning role for WntA, a signaling ligand gene previously identified as a hotspot of shape-tuning alleles involved in wing mimicry. We show that WntA loss-of-function causes multiple modifications of pattern elements in seven nymphalid butterfly species. In three butterflies with a conserved wing-pattern arrangement, WntA is necessary for the induction of stripe-like patterns known as symmetry systems and acquired a novel eyespot activator role specific to Vanessa forewings. In two Heliconius species, WntA specifies the boundaries between melanic fields and the light-color patterns that they contour. In the passionvine butterfly Agraulis, WntA removal shows opposite effects on adjacent pattern elements, revealing a dual role across the wing field. Finally, WntA acquired a divergent role in the patterning of interveinous patterns in the monarch, a basal nymphalid butterfly that lacks stripe-like symmetry systems. These results identify WntA as an instructive signal for the prepatterning of a biological system of exuberant diversity and illustrate how shifts in the deployment and effects of a single developmental gene underlie morphological change.

ChIP-Seq-Annotated Heliconius erato Genome Highlights Patterns of cis-Regulatory Evolution in Lepidoptera.

Uncovering phylogenetic patterns of cis-regulatory evolution remains a fundamental goal for evolutionary and developmental biology. Here, we characterize the evolution of regulatory loci in butterflies and moths using chromatin immunoprecipitation sequencing (ChIP-seq) annotation of regulatory elements across three stages of head development. In the process we provide a high-quality, functionally annotated genome assembly for the butterfly, Heliconius erato. Comparing cis-regulatory element conservation across six lepidopteran genomes, we find that regulatory sequences evolve at a pace similar to that of protein-coding regions. We also observe that elements active at multiple developmental stages are markedly more conserved than elements with stage-specific activity. Surprisingly, we also find that stage-specific proximal and distal regulatory elements evolve at nearly identical rates. Our study provides a benchmark for genome-wide patterns of regulatory element evolution in insects, and it shows that developmental timing of activity strongly predicts patterns of regulatory sequence evolution.