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

Lepidopteran scale cells derive from sensory organ precursors through a canonical lineage

The success of butterflies and moths is tightly linked to the origin of scales within the group. A long-standing hypothesis postulates that scales are homologous to the well-described mechanosensory bristles found in the fruit fly Drosophila melanogaster, as both derive from an epithelial precursor. Previous histological and candidate gene approaches identified parallels in genes involved in scale and bristle development. Here, we provide developmental and transcriptomic evidence that the differentiation of lepidopteran scales derives from the sensory organ precursor (SOP). Live imaging in lepidopteran pupae shows that SOP cells undergo two asymmetric divisions that first abrogate the neurogenic lineage, and then lead to a differentiated scale precursor and its associated socket cell. Single-nucleus RNA sequencing using early pupal wings revealed differential gene expression patterns that mirror SOP development, suggesting a shared developmental program. Additionally, we recovered a newly associated gene, the transcription factor pdm3, involved in the proper differentiation of butterfly wing scales. Altogether, these data open up avenues for understanding scale type specification and development, and illustrate how single-cell transcriptomics provide a powerful platform for understanding evolution of cell types.

The Genomics Revolution Drives a New Era in Entomology

Thanks to the fast development of sequencing techniques and bioinformatics tools, sequencing the genome of an insect species for specific research purposes has become an increasingly popular practice. Insect genomes not only provide sets of gene sequences but also represent a change in focus from reductionism to systemic biology in the field of entomology. Using insect genomes, researchers are able to identify and study the functions of all members of a gene family, pathway, or gene network associated with a trait of interest. Comparative genomics studies provide new insights into insect evolution, addressing long-lasting controversies in taxonomy. It is also now feasible to uncover the genetic basis of important traits by identifying variants using genome resequencing data of individual insects, followed by genome-wide association analysis. Here, we review the current progress in insect genome sequencing projects and the application of insect genomes in uncovering the phylogenetic relationships between insects and unraveling the mechanisms of important life-history traits. We also summarize the challenges in genome data sharing and possible solutions. Finally, we provide guidance for fully and deeply mining insect genome data.

Increased Evolutionary Rate in the Z chromosome of Sympatric and Allopatric Species of Morpho Butterflies

Divergent evolution of genomes among closely related species is shaped by both neutral processes and ecological forces, such as local adaptation and reinforcement. These factors can drive accelerated evolution of sex chromosomes relative to autosomes. Comparative genomic analyses between allopatric and sympatric species with overlapping or divergent ecological niches offer insights into reinforcement and ecological specialization on genome evolution. In the butterfly genus Morpho, several species coexist in sympatry, with specialization across forest strata and temporal niches. We analyzed the genomes of eight Morpho species, along with previously published genomes of three others, to compare chromosomal rearrangements and signs of positive selection in the Z chromosome vs. autosomes. We found extensive chromosomal rearrangements in Z chromosome, particularly in sympatric species with similar ecological niches, suggesting a role for inversions in restricting gene flow at a postzygotic level. Z-linked genes also exhibited significantly higher dN/dS ratios than autosomal genes across the genus, with pronounced differences in closely related species living in sympatry. Additionally, we examined the evolution of eight circadian clock genes, detecting positive selection in Period, located on the Z chromosome. Our findings suggest that the Z chromosome evolves more rapidly than autosomes, particularly among closely related species, raising questions about its role in prezygotic and postzygotic isolation mechanisms.

Developmental genetics of cuticular micro- and nano-structures in insects

Insect cuticle exhibits a wide array of micro- and nano-structures in terms of size, form, and function. However, the investigation of cellular mechanisms of morphogenesis has centered around a small number of structure types and organisms. The recent expansion of the taxa studied, and subsequent discoveries prompt us to revisit well-known models, like the one for bristle morphogenesis. In addition, common themes are emerging in the morphogenesis of cuticular structures, such as the polyploidy of precursor cells, the role of pigments and cuticular proteins in controlling chitin deposition in space and time, and the role of the apical extracellular matrix in defining the shape of the developing structure. Understanding how these structures are synthesized in biological systems holds promise for bioinspired design.

Cell membrane buckling governs early-stage ridge formation in butterfly wing scales

During the development of butterfly wing scales, ordered periodic cell membrane modulations occur at the upper surface of scale-forming cells, priming the formation of ridges. Ridges are critical for wing scale functionality, including structural color, wetting characteristics, and thermal performance. Here, we combine a morphoelastic model based on Föppl-von-Kármán plate theory with experimental observations to shed light on the biomechanical processes governing early-stage ridge formation in Painted Lady butterflies. By comparing the model predictions with time-resolved phase imaging data from live butterflies, we provide evidence that the onset of ridge formation is governed by a mechanical buckling transition induced by the interplay of membrane growth and confinement through association with regularly spaced actin bundles. Beyond ridge formation in Painted Lady scales, our theory offers a rationale for the absence of scale ridges in the lower lamina of many lepidopterans and for the alternating ridge pattern of other butterfly species.

Structural color in the bacterial domain: The ecogenomics of a 2-dimensional optical phenotype

Structural color is an optical phenomenon resulting from light interacting with nanostructured materials. Although structural color (SC) is widespread in the tree of life, the underlying genetics and genomics are not well understood. Here, we collected and sequenced a set of 87 structurally colored bacterial isolates and 30 related strains lacking SC. Optical analysis of colonies indicated that diverse bacteria from at least two different phyla (Bacteroidetes and Proteobacteria) can create two-dimensional packing of cells capable of producing SC. A pan-genome-wide association approach was used to identify genes associated with SC. The biosynthesis of uroporphyrin and pterins, as well as carbohydrate utilization and metabolism, was found to be involved. Using this information, we constructed a classifier to predict SC directly from bacterial genome sequences and validated it by cultivating and scoring 100 strains that were not part of the training set. We predicted that SCr is widely distributed within gram-negative bacteria. Analysis of over 13,000 assembled metagenomes suggested that SC is nearly absent from most habitats associated with multicellular organisms except macroalgae and is abundant in marine waters and surface/air interfaces. This work provides a large-scale ecogenomics view of SC in bacteria and identifies microbial pathways and evolutionary relationships that underlie this optical phenomenon.

The molecular basis of scale development highlighted by a single-cell atlas of Bicyclus anynana butterfly pupal forewings

Butterfly wings display a diversity of cell types, including large polyploid scale cells, yet the molecular basis of such diversity is poorly understood. To explore scale cell diversity at a transcriptomic level, we employ single-cell RNA sequencing of ∼5,200 large cells (>6 μm) from 22.5- to 25-h male pupal forewings of the butterfly Bicyclus anynana. Using unsupervised clustering, followed by in situ hybridization, immunofluorescence, and CRISPR-Cas9 editing of candidate genes, we annotate various cell types on the wing. We identify genes marking non-innervated scale cells, pheromone-producing glandular cells, and innervated sensory cell types. We show that senseless, a zinc-finger transcription factor, and HR38, a hormone receptor, determine the identity, size, and color of different scale cell types and are important regulators of scale cell differentiation. This dataset and the identification of various wing cell-type markers provide a foundation to compare and explore scale cell-type diversification across arthropod species.

The quantitative genetic basis of variation in sexual versus non-sexual butterfly wing colouration: autosomal, Z-linked, and maternal effects

Viability indicator traits are expected to be integrated extensively across the genome yet sex-limited to ensure that any benefits are sexually concordant. Understanding how such expectations are accommodated requires elucidating the quantitative genetic architecture of candidate traits in and across the sexes. Here we applied an animal modelling approach to partition the autosomal, allosomal, and direct maternal bases of variation in sexual versus non-sexual dorsal wing colouration in the butterfly Eurema hecabe. The sexual colour trait-coherently scattered ultraviolet that is under strong directional selection due to female choice-is brighter and more expansive in males, and overlays non-sexual pigmentary yellow markings that otherwise dominate both wing surfaces in each sex. Our modelling estimated high and sexually equivalent autosomal variances for ultraviolet reflectance (furnishing h2 ~ 0.58 overall and ~0.75 in males), accompanied by smaller but generally significant Z-linked and maternal components. By contrast, variation in non-sexual yellow was largely attributed to Z-linked sources. Intersexual genetic correlations based upon the major source of variation in each trait were high and not different from 1.0, implying regulation by a pool of genes common to each sex. An expansive autosomal basis for ultraviolet is consistent with its hypothesized role as a genome-wide viability indicator and ensures that both sons and daughters will inherit their father's attractiveness.

Adaptive introgression of a visual preference gene

Visual preferences are important drivers of mate choice and sexual selection, but little is known of how they evolve at the genetic level. In this study, we took advantage of the diversity of bright warning patterns displayed by Heliconius butterflies, which are also used during mate choice. Combining behavioral, population genomic, and expression analyses, we show that two Heliconius species have evolved the same preferences for red patterns by exchanging genetic material through hybridization. Neural expression of regucalcin1 correlates with visual preference across populations, and disruption of regucalcin1 with CRISPR-Cas9 impairs courtship toward conspecific females, providing a direct link between gene and behavior. Our results support a role for hybridization during behavioral evolution and show how visually guided behaviors contributing to adaptation and speciation are encoded within the genome.

The diversification of butterfly wing patterns: progress and prospects

Butterfly wings display rich phenotypic diversity and are associated with complex biological functions, thus serving as an important evolutionary system to address the genetic basis and evolution of phenotypic diversification. We review recent butterfly studies that revealed complex functions underlying diversified wing patterns and describe the genetic and environmental factors involved in wing pattern determinations. These factors lead to inter-specific divergence, genetic polymorphism, and phenotypic plasticity, which in many cases are decided by several key genes. We also summarize the research advances on gene co-option as an important origin of functional complexity and evolutionary novelty. These findings reveal a pattern of evolutionary innovation within a constrained developmental framework during butterfly wing morphogenesis, but further research is required to gain a systematic and comprehensive understanding.

Recording animal-view videos of the natural world using a novel camera system and software package

Plants, animals, and fungi display a rich tapestry of colors. Animals, in particular, use colors in dynamic displays performed in spatially complex environments. Although current approaches for studying colors are objective and repeatable, they miss the temporal variation of color signals entirely. Here, we introduce hardware and software that provide ecologists and filmmakers the ability to accurately record animal-perceived colors in motion. Specifically, our Python codes transform photos or videos into perceivable units (quantum catches) for animals of known photoreceptor sensitivity. The plans and codes necessary for end-users to capture animal-view videos are all open source and publicly available to encourage continual community development. The camera system and the associated software package will allow ecologists to investigate how animals use colors in dynamic behavioral displays, the ways natural illumination alters perceived colors, and other questions that remained unaddressed until now due to a lack of suitable tools. Finally, it provides scientists and filmmakers with a new, empirically grounded approach for depicting the perceptual worlds of nonhuman animals.

A meta-analysis of butterfly structural colors: their color range, distribution and biological production

Butterfly scales are among the richest natural sources of optical nanostructures, which produce structural color and iridescence. Several recurring nanostructure types have been described, such as ridge multilayers, gyroids and lower lamina thin films. While the optical mechanisms of these nanostructure classes are known, their phylogenetic distributions and functional ranges have not been described in detail. In this Review, we examine a century of research on the biological production of structural colors, including their evolution, development and genetic regulation. We have also created a database of more than 300 optical nanostructures in butterflies and conducted a meta-analysis of the color range, abundance and phylogenetic distribution of each nanostructure class. Butterfly structural colors are ubiquitous in short wavelengths but extremely rare in long wavelengths, especially red. In particular, blue wavelengths (around 450 nm) occur in more clades and are produced by more kinds of nanostructures than other hues. Nanostructure categories differ in prevalence, phylogenetic distribution, color range and brightness. For example, lamina thin films are the least bright; perforated lumen multilayers occur most often but are almost entirely restricted to the family Lycaenidae; and 3D photonic crystals, including gyroids, have the narrowest wavelength range (from about 450 to 550 nm). We discuss the implications of these patterns in terms of nanostructure evolution, physical constraint and relationships to pigmentary color. Finally, we highlight opportunities for future research, such as analyses of subadult and Hesperid structural colors and the identification of genes that directly build the nanostructures, with relevance for biomimetic engineering.

Hierarchical morphogenesis of swallowtail butterfly wing scale nanostructures

The study of color patterns in the animal integument is a fundamental question in biology, with many lepidopteran species being exemplary models in this endeavor due to their relative simplicity and elegance. While significant advances have been made in unraveling the cellular and molecular basis of lepidopteran pigmentary coloration, the morphogenesis of wing scale nanostructures involved in structural color production is not well understood. Contemporary research on this topic largely focuses on a few nymphalid model taxa (e.g., Bicyclus, Heliconius), despite an overwhelming diversity in the hierarchical nanostructural organization of lepidopteran wing scales. Here, we present a time-resolved, comparative developmental study of hierarchical scale nanostructures in Parides eurimedes and five other papilionid species. Our results uphold the putative conserved role of F-actin bundles in acting as spacers between developing ridges, as previously documented in several nymphalid species. Interestingly, while ridges are developing in P. eurimedes, plasma membrane manifests irregular mesh-like crossribs characteristic of Papilionidae, which delineate the accretion of cuticle into rows of planar disks in between ridges. Once the ridges have grown, disintegrating F-actin bundles appear to reorganize into a network that supports the invagination of plasma membrane underlying the disks, subsequently forming an extruded honeycomb lattice. Our results uncover a previously undocumented role for F-actin in the morphogenesis of complex wing scale nanostructures, likely specific to Papilionidae.

Genetics of yellow-orange color variation in a pair of sympatric sulphur butterflies

Continuous color polymorphisms can serve as a tractable model for the genetic and developmental architecture of traits. Here we investigated continuous color variation in Colias eurytheme and Colias philodice, two species of sulphur butterflies that hybridize in sympatry. Using quantitative trait locus (QTL) analysis and high-throughput color quantification, we found two interacting large-effect loci affecting orange-to-yellow chromaticity. Knockouts of red Malpighian tubules (red), likely involved in endosomal maturation, result in depigmented wing scales. Additionally, the transcription factor bric-a-brac can act as a modulator of orange pigmentation. We also describe the QTL architecture of other continuously varying traits, together supporting a large-X effect model where the genetic control of species-defining traits is enriched on sex chromosomes. This study sheds light on the range of possible genetic architectures that can underpin a continuously varying trait and illustrates the power of using automated measurement to score phenotypes that are not always conspicuous to the human eye.

Evidence for a single, ancient origin of a genus-wide alternative life history strategy

Understanding the evolutionary origins and factors maintaining alternative life history strategies (ALHS) within species is a major goal of evolutionary research. While alternative alleles causing discrete ALHS are expected to purge or fix over time, one-third of the ~90 species of Colias butterflies are polymorphic for a female-limited ALHS called Alba. Whether Alba arose once, evolved in parallel, or has been exchanged among taxa is currently unknown. Using comparative genome-wide association study (GWAS) and population genomic analyses, we placed the genetic basis of Alba in time-calibrated phylogenomic framework, revealing that Alba evolved once near the base of the genus and has been subsequently maintained via introgression and balancing selection. CRISPR-Cas9 mutagenesis was then used to verify a putative cis-regulatory region of Alba, which we identified using phylogenetic foot printing. We hypothesize that this cis-regulatory region acts as a modular enhancer for the induction of the Alba ALHS, which has likely facilitated its long evolutionary persistence.

Functional unit of supergene in female-limited Batesian mimicry of Papilio polytes

Supergenes are sets of genes and genetic elements that are inherited like a single gene and control complex adaptive traits, but their functional roles and units are poorly understood. In Papilio polytes, female-limited Batesian mimicry is thought to be regulated by a ∼130 kb inversion region (highly diversified region: HDR) containing 3 genes, UXT, U3X, and doublesex (dsx) which switches non-mimetic and mimetic types. To determine the functional unit, we here performed electroporation-mediated RNAi analyses (and further Crispr/Cas9 for UXT) of genes within and flanking the HDR in pupal hindwings. We first clarified that non-mimetic dsx-h had a function to form the non-mimetic trait in female and only dsx-H isoform 3 had an important function in the formation of mimetic traits. Next, we found that UXT was involved in making mimetic-type pale-yellow spots and adjacent gene sir2 in making red spots in hindwings, both of which refine more elaborate mimicry. Furthermore, downstream gene networks of dsx, U3X, and UXT screened by RNA sequencing showed that U3X upregulated dsx-H expression and repressed UXT expression. These findings demonstrate that a set of multiple genes, not only inside but also flanking HDR, can function as supergene members, which extends the definition of supergene unit than we considered before. Also, our results indicate that dsx functions as the switching gene and some other genes such as UXT and sir2 within the supergene unit work as the modifier gene.

Genome assembly of 3 Amazonian Morpho butterfly species reveals Z-chromosome rearrangements between closely related species living in sympatry

The genomic processes enabling speciation and species coexistence in sympatry are still largely unknown. Here we describe the whole-genome sequencing and assembly of 3 closely related species from the butterfly genus Morpho: Morpho achilles (Linnaeus, 1758), Morpho helenor (Cramer, 1776), and Morpho deidamia (Höbner, 1819). These large blue butterflies are emblematic species of the Amazonian rainforest. They live in sympatry in a wide range of their geographical distribution and display parallel diversification of dorsal wing color pattern, suggesting local mimicry. By sequencing, assembling, and annotating their genomes, we aim at uncovering prezygotic barriers preventing gene flow between these sympatric species. We found a genome size of  480 Mb for the 3 species and a chromosomal number ranging from 2n = 54 for M. deidamia to 2n = 56 for M. achilles and M. helenor. We also detected inversions on the sex chromosome Z that were differentially fixed between species, suggesting that chromosomal rearrangements may contribute to their reproductive isolation. The annotation of their genomes allowed us to recover in each species at least 12,000 protein-coding genes and to discover duplications of genes potentially involved in prezygotic isolation like genes controlling color discrimination (L-opsin). Altogether, the assembly and the annotation of these 3 new reference genomes open new research avenues into the genomic architecture of speciation and reinforcement in sympatry, establishing Morpho butterflies as a new eco-evolutionary model.

The scent chemistry of butterflies

Butterflies use structurally highly diverse volatile compounds for communication, in addition to visual signals. These compounds originate from plants or a formed de novo especially by male butterflies that possess specific scent organs.