Genetics of flight in spongy moths (Lymantria dispar ssp.): functionally integrated profiling of a complex invasive trait

BACKGROUND

Flight can drastically enhance dispersal capacity and is a key trait defining the potential of exotic insect species to spread and invade new habitats. The phytophagous European spongy moths (ESM, Lymantria dispar dispar) and Asian spongy moths (ASM; a multi-species group represented here by L. d. asiatica and L. d. japonica), are globally invasive species that vary in adult female flight capability-female ASM are typically flight capable, whereas female ESM are typically flightless. Genetic markers of flight capability would supply a powerful tool for flight profiling of these species at any intercepted life stage. To assess the functional complexity of spongy moth flight and to identify potential markers of flight capability, we used multiple genetic approaches aimed at capturing complementary signals of putative flight-relevant genetic divergence between ESM and ASM: reduced representation genome-wide association studies, whole genome sequence comparisons, and developmental transcriptomics. We then judged the candidacy of flight-associated genes through functional analyses aimed at addressing the proximate demands of flight and salient features of the ecological context of spongy moth flight evolution.

RESULTS

Candidate gene sets were typically non-overlapping across different genetic approaches, with only nine gene annotations shared between any pair of approaches. We detected an array of flight-relevant functional themes across gene sets that collectively suggest divergence in flight capability between European and Asian spongy moth lineages has coincided with evolutionary differentiation in multiple aspects of flight development, execution, and surrounding life history. Overall, our results indicate that spongy moth flight evolution has shaped or been influenced by a large and functionally broad network of traits.

CONCLUSIONS

Our study identified a suite of flight-associated genes in spongy moths suited to exploration of the genetic architecture and evolution of flight, or validation for flight profiling purposes. This work illustrates how complementary genetic approaches combined with phenotypically targeted functional analyses can help to characterize genetically complex traits.

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.

Spatial and temporal regulation of Wnt signaling pathway members in the development of butterfly wing patterns

Wnt signaling members are involved in the differentiation of cells associated with eyespot and band color patterns on the wings of butterflies, but the identity and spatio-temporal regulation of specific Wnt pathway members remains unclear. Here, we explore the localization and function of Armadillo/β-catenin dependent (canonical) and Armadillo/β-catenin independent (noncanonical) Wnt signaling in eyespot and band development in Bicyclus anynana by localizing Armadillo (Arm), the expression of all eight Wnt ligand and four frizzled receptor transcripts present in the genome of this species and testing the function of some of the ligands and receptors using CRISPR-Cas9. We show that distinct Wnt signaling pathways are essential for eyespot and band patterning in butterflies and are likely interacting to control their active domains.

Ultrabithorax modifies a regulatory network of genes essential for butterfly eyespot development in a wing sector-specific manner

Nymphalid butterfly species often have a different number of eyespots in forewings and hindwings, but how the hindwing identity gene Ultrabithorax (Ubx) drives this asymmetry is not fully understood. We examined a three-gene regulatory network for eyespot development in the hindwings of Bicyclus anynana butterflies and compared it with the same network previously described for forewings. We also examined how Ubx interacts with each of these three eyespot-essential genes. We found similar genetic interactions between the three genes in fore- and hindwings, but we discovered three regulatory differences: Antennapedia (Antp) merely enhances spalt (sal) expression in the eyespot foci in hindwings, but is not essential for sal activation, as in forewings; Ubx upregulates Antp in all hindwing eyespot foci but represses Antp outside these wing regions; and Ubx regulates sal in a wing sector-specific manner, i.e. it activates sal expression only in the sectors that have hindwing-specific eyespots. We propose a model for how the regulatory connections between these four genes evolved to produce wing- and sector-specific variation in eyespot number.

Laser Microdissection-Mediated Isolation of Butterfly Wing Tissue for Spatial Transcriptomics

The assignment of specific patterns of gene expression to specific cells in a complex tissue facilitates the connection between genotype and phenotype. Single-cell sequencing of whole tissues produces single-cell transcript resolution but lacks the spatial information of the derivation of each cell, whereas techniques such as multiplex FISH localize transcripts to specific cells in a tissue but require a priori information of the target transcripts to examine. Laser dissection of tissues followed by transcriptome analysis is an efficient and cost-effective technique that provides both unbiased gene expression discovery together with spatial information. Here, we detail a laser dissection protocol for total RNA extraction from butterfly larval and pupal wing tissues, without the need of paraffin embedding or the use of a microtome, that could be useful to researchers interested in the transcriptome of specific areas of the wing during development. This protocol can bypass difficulties in extracting high quality RNA from thick fixed tissues for sequencing applications.

Distal-less and spalt are distal organisers of pierid wing patterns

Two genes, Distal-less (Dll) and spalt (sal), are known to be involved in establishing nymphalid butterfly wing patterns. They function in several ways: in the differentiation of the eyespot’s central signalling cells, or foci; in the differentiation of the surrounding black disc; in overall scale melanisation (Dll); and in elaborating marginal patterns, such as parafocal elements. However, little is known about the functions of these genes in the development of wing patterns in other butterfly families. Here, we study the expression and function of Dll and sal in the development of spots and other melanic wing patterns of the Indian cabbage white, Pieris canidia, a pierid butterfly. In P. canidia, both Dll and Sal proteins are expressed in the scale-building cells at the wing tips, in chevron patterns along the pupal wing margins, and in areas of future scale melanisation. Additionally, Sal alone is expressed in the future black spots. CRISPR knockouts of Dll and sal showed that each gene is required for the development of melanic wing pattern elements, and repressing pteridine granule formation, in the areas where they are expressed. We conclude that both genes likely play ancestral roles in organising distal butterfly wing patterns, across pierid and nymphalid butterflies, but are unlikely to be differentiating signalling centres in pierids black spots. The genetic and developmental mechanisms that set up the location of spots and eyespots are likely distinct in each lineage.

The development of shape. Modular control of growth in the lepidopteran forewing

The mechanisms by which tissues and organs achieve their final size and shape during development are largely unknown. Although we have learned much about the mechanisms that control growth, little is known about how those play out to achieve a structure's specific final size and shape. The wings of insects are attractive systems for the study of the control of morphogenesis, because they are perfectly flat and two-dimensional, composed of two closely appressed cellular monolayers in which morphogenetic processes can be easily visualized. The wings of Lepidoptera arise from imaginal disks whose structure is always perfectly congruent with that of the adult wing, so that it is possible to fate-map corresponding positions on the larval disk to those of the adult wing. Here we show that the forewing imaginal disks of Junonia coenia are subdivided into four domains, with characteristic patterns of expression of known patterning genes Spalt (Sal), Engrailed (En), and Cubitus interruptus (Ci). We show that DNA and protein synthesis, as well as mitoses, are spatially patterned in a domain-specific way. Knockdown of Sal and En using produced domain-specific reductions in the shape of the forewing. Knockdown of signaling pathways involved in the regulation of growth likewise altered the shape of the forewing in a domain-specific way. Our results reveal a multi-level regulation of forewing shape involving hormones and growth-regulating genes.

Butterfly eyespots evolved via cooption of an ancestral gene-regulatory network that also patterns antennae, legs, and wings

Where do butterfly eyespots come from? One of the long-standing questions in the field of evolution concerns addressing where novel complex traits come from. Here we show that butterfly eyespots, a novel complex trait, likely originated from the redeployment of a preexisting gene-regulatory network regulating antennae, legs, and wings, to novel locations on the wing.

Butterfly eyespots are beautiful novel traits with an unknown developmental origin. Here we show that eyespots likely originated via cooption of parts of an ancestral appendage gene-regulatory network (GRN) to novel locations on the wing. Using comparative transcriptome analysis, we show that eyespots cluster most closely with antennae, relative to multiple other tissues. Furthermore, three genes essential for eyespot development, Distal-less (Dll), spalt (sal), and Antennapedia (Antp), share similar regulatory connections as those observed in the antennal GRN. CRISPR knockout of cis-regulatory elements (CREs) for Dll and sal led to the loss of eyespots, antennae, legs, and also wings, demonstrating that these CREs are highly pleiotropic. We conclude that eyespots likely reused an ancient GRN for their development, a network also previously implicated in the development of antennae, legs, and wings.

Sob gene is critical to wing development in Bombyx mori and Tribolium castaneum

The development of insect appendages requires the expression of multiple genes in a strict spatial and temporal order. The odd-skipped family genes are vital transcriptional factors involved in embryonic development. The development and morphogenesis of the insect wing requires multiple transcription factors to regulate the expression of wing patterning genes at the transcriptional level. However, the function of odd-related genes in insect wing morphogenesis and development during postembryonic stages is unclear. We focused on the roles of the sister of odd and bowl (sob) gene, a member of odd-skipped family genes, during the wing morphopoiesis in Bombyx mori using the clustered regularly interspaced palindromic repeats (CRISPR)/CRISPR-associated protein 9 system and in Tribolium castaneum by RNA interference. The results showed that the wings were significantly smaller and degenerated, and wing veins were indistinct in the sob gene loss-of-function group in both B. mori and T. castaneum. Quantitative real-time polymerase chain reaction revealed that the Tcsob gene regulated the expression of wing development genes, such as the cht 7 and the vg gene. The findings suggest the importance of sob gene in insect wing morphology formation during postembryonic stages.