The ivory lncRNA regulates seasonal color patterns in buckeye butterflies.

Long noncoding RNAs (lncRNAs) are transcribed elements increasingly recognized for their roles in regulating gene expression. Thus far, however, we have little understanding of how lncRNAs contribute to evolution and adaptation. Here, we show that a conserved lncRNA, ivory, is an important color patterning gene in the buckeye butterfly Junonia coenia. ivory overlaps with cortex, a locus linked to multiple cases of crypsis and mimicry in Lepidoptera. Along with a companion paper by Livraghi et al., we argue that ivory, not cortex, is the color pattern gene of interest at this locus. In J. coenia, a cluster of cis-regulatory elements (CREs) in the first intron of ivory are genetically associated with natural variation in seasonal color pattern plasticity, and targeted deletions of these CREs phenocopy seasonal phenotypes. Deletions of different ivory CREs produce other distinct phenotypes as well, including loss of melanic eyespot rings, and positive and negative changes in overall wing pigmentation. We show that the color pattern transcription factors Spineless, Bric-a-brac, and Ftz-f1 bind to the ivory promoter during wing pattern development, suggesting that they directly regulate ivory. This case study demonstrates how cis-regulation of a single noncoding RNA can exert diverse and nuanced effects on the evolution and development of color patterns, including modulating seasonally plastic color patterns.

A long noncoding RNA at the cortex locus controls adaptive coloration in butterflies.

Evolutionary variation in the wing pigmentation of butterflies and moths offers striking examples of adaptation by crypsis and mimicry. The cortex locus has been independently mapped as the locus controlling color polymorphisms in 15 lepidopteran species, suggesting that it acts as a genomic hotspot for the diversification of wing patterns, but functional validation through protein-coding knockouts has proven difficult to obtain. Our study unveils the role of a long noncoding RNA (lncRNA) which we name ivory, transcribed from the cortex locus, in modulating color patterning in butterflies. Strikingly, ivory expression prefigures most melanic patterns during pupal development, suggesting an early developmental role in specifying scale identity. To test this, we generated CRISPR mosaic knock-outs in five nymphalid butterfly species and show that ivory mutagenesis yields transformations of dark pigmented scales into white or light-colored scales. Genotyping of Vanessa cardui germline mutants associates these phenotypes to small on-target deletions at the conserved first exon of ivory. In contrast, cortex germline mutant butterflies with confirmed null alleles lack any wing phenotype and exclude a color patterning role for this adjacent gene. Overall, these results show that a lncRNA gene acts as a master switch of color pattern specification and played key roles in the adaptive diversification of wing patterns in butterflies.

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.