Patterns of selection across gene regulatory networks.

Gene regulatory networks (GRNs) are the core engine of organismal development. If we would like to understand the origin and diversification of phenotypes, it is necessary to consider the structure of GRNs in order to reconstruct the links between genetic mutations and phenotypic change. Much of the progress in evolutionary developmental biology, however, has occurred without a nuanced consideration of the evolution of functional relationships between genes, especially in the context of their broader network interactions. Characterizing and comparing GRNs across traits and species in a more detailed way will allow us to determine how network position influences what genes drive adaptive evolution. In this perspective paper, we consider the architecture of developmental GRNs and how positive selection strength may vary across a GRN. We then propose several testable models for these patterns of selection and experimental approaches to test these models.

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.

High-Quality Genome Assembly and Comprehensive Transcriptome of the Painted Lady Butterfly Vanessa cardui.

The painted lady butterfly, Vanessa cardui, has the longest migration routes, the widest hostplant diversity, and one of the most complex wing patterns of any insect. Due to minimal culturing requirements, easily characterized wing pattern elements, and technical feasibility of CRISPR/Cas9 genome editing, V. cardui is emerging as a functional genomics model for diverse research programs. Here, we report a high-quality, annotated genome assembly of the V. cardui genome, generated using 84× coverage of PacBio long-read data, which we assembled into 205 contigs with a total length of 425.4 Mb (N50 = 10.3 Mb). The genome was very complete (single-copy complete Benchmarking Universal Single-Copy Orthologs [BUSCO] 97%), with contigs assembled into presumptive chromosomes using synteny analyses. Our annotation used embryonic, larval, and pupal transcriptomes, and 20 transcriptomes across five different wing developmental stages. Gene annotations showed a high level of accuracy and completeness, with 14,437 predicted protein-coding genes. This annotated genome assembly constitutes an important resource for diverse functional genomic studies ranging from the developmental genetic basis of butterfly color pattern, to coevolution with diverse hostplants.

The Dryas iulia Genome Supports Multiple Gains of a W Chromosome from a B Chromosome in Butterflies.

In butterflies and moths, which exhibit highly variable sex determination mechanisms, the homogametic Z chromosome is deeply conserved and is featured in many genome assemblies. The evolution and origin of the female W sex chromosome, however, remains mostly unknown. Previous studies have proposed that a ZZ/Z0 sex determination system is ancestral to Lepidoptera, and that W chromosomes may originate from sex-linked B chromosomes. Here, we sequence and assemble the female Dryas iulia genome into 32 highly contiguous ordered and oriented chromosomes, including the Z and W sex chromosomes. We then use sex-specific Hi-C, ATAC-seq, PRO-seq, and whole-genome DNA sequence data sets to test if features of the D. iulia W chromosome are consistent with a hypothesized B chromosome origin. We show that the putative W chromosome displays female-associated DNA sequence, gene expression, and chromatin accessibility to confirm the sex-linked function of the W sequence. In contrast with expectations from studies of homologous sex chromosomes, highly repetitive DNA content on the W chromosome, the sole presence of domesticated repetitive elements in functional DNA, and lack of sequence homology with the Z chromosome or autosomes is most consistent with a B chromosome origin for the W, although it remains challenging to rule out extensive sequence divergence. Synteny analysis of the D. iulia W chromosome with other female lepidopteran genome assemblies shows no homology between W chromosomes and suggests multiple, independent origins of the W chromosome from a B chromosome likely occurred in butterflies.

Chromosome Fusion Affects Genetic Diversity and Evolutionary Turnover of Functional Loci but Consistently Depends on Chromosome Size.

Major changes in chromosome number and structure are linked to a series of evolutionary phenomena, including intrinsic barriers to gene flow or suppression of recombination due to chromosomal rearrangements. However, chromosome rearrangements can also affect the fundamental dynamics of molecular evolution within populations by changing relationships between linked loci and altering rates of recombination. Here, we build chromosome-level assembly Eueides isabella and, together with a recent chromosome-level assembly of Dryas iulia, examine the evolutionary consequences of multiple chromosome fusions in Heliconius butterflies. These assemblies pinpoint fusion points on 10 of the 20 autosomal chromosomes and reveal striking differences in the characteristics of fused and unfused chromosomes. The ten smallest autosomes in D. iulia and E. isabella, which have each fused to a longer chromosome in Heliconius, have higher repeat and GC content, and longer introns than predicted by their chromosome length. When fused, these characteristics change to become more in line with chromosome length. The fusions also led to reduced diversity, which likely reflects increased background selection and selection against introgression between diverging populations, following a reduction in per-base recombination rate. We further show that chromosome size and fusion impact turnover rates of functional loci at a macroevolutionary scale. Together these results provide further evidence that chromosome fusion in Heliconius likely had dramatic effects on population level processes shaping rates of neutral and adaptive divergence. These effects may have impacted patterns of diversification in Heliconius, a classic example of an adaptive radiation.

Seasonal plasticity: how do butterfly wing pattern traits evolve environmental responsiveness?

Phenotypic plasticity in response to environmental cues is common in butterflies, and is a major driver of butterfly wing pattern diversity. The endocrine signal ecdysone has been revealed as a major modulator of plasticity in butterflies. External cues such as day length or temperature are translated internally into variation in ecdysone titers, which in turn lead to alternate phenotypes such as seasonal wing patterns. Here we review the evidence showing that ecdysone-mediated plasticity of different wing pattern features such as wing color and eyespot size can evolve independently. Recent studies show that ecdysone regulates gene expression in Drosophila melanogaster via a chromatin remodeling mechanism. We thus propose that environmentally responsive ecdysone titers in butterflies may also function via chromatin regulation to promote different seasonal phenotypes. We present a model of ecdysone response evolution that integrates both gene regulatory architecture and organismal development, and propose a set of testable mechanistic hypotheses for how plastic response profiles of specific genes can evolve.

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.

Many functionally connected loci foster adaptive diversification along a neotropical hybrid zone.

Characterizing the genetic complexity of adaptation and trait evolution is a major emphasis of evolutionary biology and genetics. Incongruent findings from genetic studies have resulted in conceptual models ranging from a few large-effect loci to massively polygenic architectures. Here, we combine chromatin immunoprecipitation sequencing, Hi-C, RNA sequencing, and 40 whole-genome sequences from Heliconius butterflies to show that red color pattern diversification occurred via many genomic loci. We find that the red wing pattern master regulatory transcription factor Optix binds dozens of loci also under selection, which frequently form three-dimensional adaptive hubs with selection acting on multiple physically interacting genes. Many Optix-bound genes under selection are tied to pigmentation and wing development, and these loci collectively maintain separation between adaptive red color pattern phenotypes in natural populations. We propose a model of trait evolution where functional connections between loci may resolve much of the disparity between large-effect and polygenic evolutionary models.

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.

Transcription factors underlying wing margin color patterns and pupal cuticle markings in butterflies.

BACKGROUND: The diversity of butterfly color patterns can be attributed to a relatively small number of pattern elements that are homologous across Lepidoptera. Although genes involved in patterning some of these elements have been identified, the development of several major elements remains poorly understood. To identify genes underlying wing pupal cuticle markings and wing margin color patterns, we examined expression of the candidate transcription factors Engrailed/Invected (En/Inv), Distal-less (Dll), Cubitus interruptus (Ci), and Spalt in two nymphalids: Junonia coenia and Bicyclus anynana. RESULTS: We found that En/Inv, Dll, and Ci mark domains on the J. coenia last-instar forewing disc that closely correspond to the position and shape of pupal cuticle markings. We also found that Spalt demarcates wing margin color patterns in both J. coenia and B. anynana, and that CRISPR/Cas9 deletions in the spalt gene result in reduction and loss of wing margin color patterns in J. coenia. These data demonstrate a role for spalt in promoting wing margin color patterning, in addition to its previously described role in eyespot patterning. CONCLUSION: Our observations support the model that a core set of regulatory genes are redeployed multiple times, and in multiple roles, during butterfly wing pattern development. Of these genes, spalt is of special interest as it plays a dual role in both eyespot and margin color pattern development.

Dichotomy of Dosage Compensation along the Neo Z Chromosome of the Monarch Butterfly.

Mechanisms of sex chromosome dosage compensation (SCDC) differ strikingly among animals. In Drosophila flies, chromosome-wide transcription is doubled from the single X chromosome in hemizygous (XY) males, whereas in Caenorhabditis nematodes, expression is halved for both X copies in homozygous (XX) females [1, 2]. Unlike other female-heterogametic (WZ female and ZZ male) animals, moths and butterflies exhibit sex chromosome dosage compensation patterns typically seen only in male-heterogametic species [3]. The monarch butterfly carries a newly derived Z chromosome segment that arose from an autosomal fusion with the ancestral Z [4]. Using a highly contiguous genome assembly, we show that gene expression is balanced between sexes along the entire Z chromosome but with distinct modes of compensation on the two segments. On the ancestral Z segment, depletion of H4K16ac corresponds to nearly halving of biallelic transcription in males, a pattern convergent to nematodes. Conversely, the newly derived Z segment shows a Drosophila-like mode of compensation, with enriched H4K16ac levels corresponding to doubled monoallelic transcription in females. Our work reveals that, contrary to the expectation of co-opting regulatory mechanisms readily in place, the evolution of plural modes of dosage compensation is also possible along a single sex chromosome within a species.

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.

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.

Contrasting Roles of Transcription Factors Spineless and EcR in the Highly Dynamic Chromatin Landscape of Butterfly Wing Metamorphosis.

Development requires highly coordinated changes in chromatin accessibility in order for proper gene regulation to occur. Here, we identify factors associated with major, discrete changes in chromatin accessibility during butterfly wing metamorphosis. By combining mRNA sequencing (mRNA-seq), assay for transposase-accessible chromatin using sequencing (ATAC-seq), and machine learning analysis of motifs, we show that distinct sets of transcription factors are predictive of chromatin opening at different developmental stages. Our data suggest an important role for nuclear hormone receptors early in metamorphosis, whereas PAS-domain transcription factors are strongly associated with later chromatin opening. Chromatin immunoprecipitation sequencing (ChIP-seq) validation of select candidate factors showed spineless binding to be a major predictor of opening chromatin. Surprisingly, binding of ecdysone receptor (EcR), a candidate accessibility factor in Drosophila, was not predictive of opening but instead marked persistent sites. This work characterizes the chromatin dynamics of insect wing metamorphosis, identifies candidate chromatin remodeling factors in insects, and presents a genome assembly of the model butterfly Junonia coenia.

Genome-Wide Regulatory Adaptation Shapes Population-Level Genomic Landscapes in Heliconius.

Cis-regulatory evolution is an important engine of organismal diversification. Although recent studies have looked at genomic patterns of regulatory evolution between species, we still have a poor understanding of the magnitude and nature of regulatory variation within species. Here, we examine the evolution of regulatory element activity over wing development in three Heliconius erato butterfly populations to determine how regulatory variation is associated with population structure. We show that intraspecific divergence in chromatin accessibility and regulatory activity is abundant, and that regulatory variants are spatially clustered in the genome. Regions with strong population structure are highly enriched for regulatory variants, and enrichment patterns are associated with developmental stage and gene expression. We also found that variable regulatory elements are particularly enriched in species-specific genomic regions and long interspersed nuclear elements. Our findings suggest that genome-wide selection on chromatin accessibility and regulatory activity is an important force driving patterns of genomic divergence within Heliconius species. This work also provides a resource for the study of gene regulatory evolution in H. erato and other heliconiine butterflies.

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.

Genetic Basis of Melanin Pigmentation in Butterfly Wings.

Despite the variety, prominence, and adaptive significance of butterfly wing patterns, surprisingly little is known about the genetic basis of wing color diversity. Even though there is intense interest in wing pattern evolution and development, the technical challenge of genetically manipulating butterflies has slowed efforts to functionally characterize color pattern development genes. To identify candidate wing pigmentation genes, we used RNA sequencing to characterize transcription across multiple stages of butterfly wing development, and between different color pattern elements, in the painted lady butterfly Vanessa cardui This allowed us to pinpoint genes specifically associated with red and black pigment patterns. To test the functions of a subset of genes associated with presumptive melanin pigmentation, we used clustered regularly interspaced short palindromic repeats (CRISPR)/Cas9 genome editing in four different butterfly genera. pale, Ddc, and yellow knockouts displayed reduction of melanin pigmentation, consistent with previous findings in other insects. Interestingly, however, yellow-d, ebony, and black knockouts revealed that these genes have localized effects on tuning the color of red, brown, and ochre pattern elements. These results point to previously undescribed mechanisms for modulating the color of specific wing pattern elements in butterflies, and provide an expanded portrait of the insect melanin pathway.

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.

Genome editing in butterflies reveals that spalt promotes and Distal-less represses eyespot colour patterns.

Butterfly eyespot colour patterns are a key example of how a novel trait can appear in association with the co-option of developmental patterning genes. Little is known, however, about how, or even whether, co-opted genes function in eyespot development. Here we use CRISPR/Cas9 genome editing to determine the roles of two co-opted transcription factors that are expressed during early eyespot determination. We found that deletions in a single gene, spalt, are sufficient to reduce or completely delete eyespot colour patterns, thus demonstrating a positive regulatory role for this gene in eyespot determination. Conversely, and contrary to previous predictions, deletions in Distal-less (Dll) result in an increase in the size and number of eyespots, illustrating a repressive role for this gene in eyespot development. Altogether our results show that the presence, absence and shape of butterfly eyespots can be controlled by the activity of two co-opted transcription factors.