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

Comprehensive Identification and Characterization of HML-9 Group in Chimpanzee Genome

Endogenous retroviruses (ERVs) are related to long terminal repeat (LTR) retrotransposons, comprising gene sequences of exogenous retroviruses integrated into the host genome and inherited according to Mendelian law. They are considered to have contributed greatly to the evolution of host genome structure and function. We previously characterized HERV-K HML-9 in the human genome. However, the biological function of this type of element in the genome of the chimpanzee, which is the closest living relative of humans, largely remains elusive. Therefore, the current study aims to characterize HML-9 in the chimpanzee genome and to compare the results with those in the human genome. Firstly, we report the distribution and genetic structural characterization of the 26 proviral elements and 38 solo LTR elements of HML-9 in the chimpanzee genome. The results showed that the distribution of these elements displayed a non-random integration pattern, and only six elements maintained a relatively complete structure. Then, we analyze their phylogeny and reveal that the identified elements all cluster together with HML-9 references and with those identified in the human genome. The HML-9 integration time was estimated based on the 2-LTR approach, and the results showed that HML-9 elements were integrated into the chimpanzee genome between 14 and 36 million years ago and into the human genome between 18 and 49 mya. In addition, conserved motifs, cis-regulatory regions, and enriched PBS sequence features in the chimpanzee genome were predicted based on bioinformatics. The results show that pathways significantly enriched for ERV LTR-regulated genes found in the chimpanzee genome are closely associated with disease development, including neurological and neurodevelopmental psychiatric disorders. In summary, the identification, characterization, and genomics of HML-9 presented here not only contribute to our understanding of the role of ERVs in primate evolution but also to our understanding of their biofunctional significance.

Parthenogenetic Stick Insects Exhibit Signatures of Preservation in the Molecular Architecture of Male Reproduction

After the loss of a trait, theory predicts that the molecular machinery underlying its phenotypic expression should decay. Yet, empirical evidence is contrasting. Here, we test the hypotheses that (i) the molecular ground plan of a lost trait could persist due to pleiotropic effects on other traits and (ii) that gene co-expression network architecture could constrain individual gene expression. Our testing ground has been the Bacillus stick insect species complex, which contains close relatives that are either bisexual or parthenogenetic. After the identification of genes expressed in male reproductive tissues in a bisexual species, we investigated their gene co-expression network structure in two parthenogenetic species. We found that gene co-expression within the male gonads was partially preserved in parthenogens. Furthermore, parthenogens did not show relaxed selection on genes upregulated in male gonads in the bisexual species. As these genes were mostly expressed in female gonads, this preservation could be driven by pleiotropic interactions and an ongoing role in female reproduction. Connectivity within the network also played a key role, with highly connected-and more pleiotropic-genes within male gonad also having a gonad-biased expression in parthenogens. Our findings provide novel insight into the mechanisms which could underlie the production of rare males in parthenogenetic lineages; more generally, they provide an example of the cryptic persistence of a lost trait molecular architecture, driven by gene pleiotropy on other traits and within their co-expression network.

A new dawn for the study of cell type evolution

Animal evolution is influenced by the emergence of new cell types, yet our understanding of this process remains elusive. This prompts the need for a broader exploration across diverse research organisms, facilitated by recent breakthroughs, such as gene editing tools and single-cell genomics. Essential to our understanding of cell type evolution is the accurate identification of homologous cells. We delve into the significance of considering developmental ontogeny and potential pitfalls when drawing conclusions about cell type homology. Additionally, we highlight recent discoveries in the study of cell type evolution through the application of single-cell transcriptomics and pinpoint areas ripe for further exploration.

Protein Delivery to Insect Epithelial Cells In Vivo: Potential Application to Functional Molecular Analysis of Proteins in Butterfly Wing Development

Protein delivery to cells in vivo has great potential for the functional analysis of proteins in nonmodel organisms. In this study, using the butterfly wing system, we investigated a method of protein delivery to insect epithelial cells that allows for easy access, treatment, and observation in real time in vivo. Topical and systemic applications (called the sandwich and injection methods, respectively) were tested. In both methods, green/orange fluorescent proteins (GFP/OFP) were naturally incorporated into intracellular vesicles and occasionally into the cytosol from the apical surface without any delivery reagent. However, the antibodies were not delivered by the sandwich method at all, and were delivered only into vesicles by the injection method. A membrane-lytic peptide, L17E, appeared to slightly improve the delivery of GFP/OFP and antibodies. A novel peptide reagent, ProteoCarry, successfully promoted the delivery of both GFP/OFP and antibodies into the cytosol via both the sandwich and injection methods. These protein delivery results will provide opportunities for the functional molecular analysis of proteins in butterfly wing development, and may offer a new way to deliver proteins into target cells in vivo in nonmodel organisms.

Chromosome-Level Genome Assembly of Papilio elwesi Leech, 1889 (Lepidoptera: Papilionidae)

A rarely seen butterfly species, the large swallowtail butterfly Papilio elwesi Leech, 1889 (Lepidoptera: Papilionidae), endemic to the Chinese mainland, has been declared a state-protected animal in China since 2000, but its genome is not yet available. To obtain high-quality genome assembly and annotation, we sequenced the genome and transcriptome of P. elwesi using the PacBio and PromethION platforms, respectively. The final assembled genome was 358.51 Mb, of which 97.59% was anchored to chromosomes (30 autosomes and 1 Z sex chromosome), with a contig/scaffold N50 length of 6.79/12.32 Mb and 99.0% (n = 1367) BUSCO completeness. The genome annotation pointed to 36.82% (131.99 Mb) repetitive elements and 1296 non-coding RNAs in the genome, along with 13,681 protein-coding genes that cover 98.6% (1348) of the BUSCO genes. Among the 11,499 identified gene families, 104 underwent significantly rapid expansions or contractions, and these rapidly expanding families play roles in detoxification and metabolism. Additionally, strong synteny exists between the chromosomes of P. elwesi and P. machaon. The chromosome-level genome of P. elwesi could serve as an important genomic resource for furthering our understanding of butterfly evolution and for more in-depth genomic analyses.

High level of novelty under the hood of convergent evolution

Little is known about the extent to which species use homologous regulatory architectures to achieve phenotypic convergence. By characterizing chromatin accessibility and gene expression in developing wing tissues, we compared the regulatory architecture of convergence between a pair of mimetic butterfly species. Although a handful of color pattern genes are known to be involved in their convergence, our data suggest that different mutational paths underlie the integration of these genes into wing pattern development. This is supported by a large fraction of accessible chromatin being exclusive to each species, including the de novo lineage-specific evolution of a modular optix enhancer. These findings may be explained by a high level of developmental drift and evolutionary contingency that occurs during the independent evolution of mimicry.

The Genetic Mechanisms Underlying the Concerted Expression of the yellow and tan Genes in Complex Patterns on the Abdomen and Wings of Drosophila guttifera

How complex morphological patterns form is an intriguing question in developmental biology. However, the mechanisms that generate complex patterns remain largely unknown. Here, we sought to identify the genetic mechanisms that regulate the tan (t) gene in a multi-spotted pigmentation pattern on the abdomen and wings of Drosophila guttifera. Previously, we showed that yellow (y) gene expression completely prefigures the abdominal and wing pigment patterns of this species. In the current study, we demonstrate that the t gene is co-expressed with the y gene in nearly identical patterns, both transcripts foreshadowing the adult abdominal and wing melanin spot patterns. We identified cis-regulatory modules (CRMs) of t, one of which drives reporter expression in six longitudinal rows of spots on the developing pupal abdomen, while the second CRM activates the reporter gene in a spotted wing pattern. Comparing the abdominal spot CRMs of y and t, we found a similar composition of putative transcription factor binding sites that are thought to regulate the complex expression patterns of both terminal pigmentation genes y and t. In contrast, the y and t wing spots appear to be regulated by distinct upstream factors. Our results suggest that the D. guttifera abdominal and wing melanin spot patterns have been established through the co-regulation of y and t, shedding light on how complex morphological traits may be regulated through the parallel coordination of downstream target genes.

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.

Butterfly Wing Color Pattern Modification Inducers May Act on Chitin in the Apical Extracellular Site: Implications in Morphogenic Signals for Color Pattern Determination

Butterfly wing color patterns are modified by various treatments, such as temperature shock, injection of chemical inducers, and covering materials on pupal wing tissue. Their mechanisms of action have been enigmatic. Here, we investigated the mechanisms of color pattern modifications usingthe blue pansy butterfly Junoniaorithya. We hypothesized that these modification-inducing treatments act on the pupal cuticle or extracellular matrix (ECM). Mechanical load tests revealed that pupae treated with cold shock or chemical inducers were significantly less rigid, suggesting that these treatments made cuticle formation less efficient. A known chitin inhibitor, FB28 (fluorescent brightener 28), was discovered to efficiently induce modifications. Taking advantage of its fluorescent character, fluorescent signals from FB28 were observed in live pupae in vivo from the apical extracellular side and were concentrated at the pupal cuticle focal spots immediately above the eyespot organizing centers. It was shown that chemical modification inducers and covering materials worked additively. Taken together, various modification-inducing treatments likely act extracellularly on chitin or other polysaccharides to inhibit pupal cuticle formation or ECM function, which probably causes retardation of morphogenic signals. It is likely that an interactive ECM is required for morphogenic signals for color pattern determination to travel long distances.

Melanin Synthesis Pathway Interruption: CRISPR/Cas9-mediated Knockout of dopa decarboxylase (DDC) in Harmonia axyridis (Coleoptera: Coccinellidae)

CRISPR/Cas9 technology is a very powerful genome editing tool and has been used in many insect species for functional genomics studies through targeted gene mutagenesis. Here, we successfully established CRISPR/Cas9 research platform in Asian multi-colored ladybird beetle, Harmonia axyridis, an important natural enemy in biological control. In this study, one pivotal gene dopa decarboxylase (DDC) in melanin synthesis was targeted by CRISPR/Cas9 to generate mutants in H. axyridis by CRISPR/Cas9 technology. Our results showed that injection of single guide RNA of the DDC and Cas9 protein into preblastoderm eggs induced one insertion and four deletion (indels) mutant H. axyridis. Mutations of HaDDC gene generated 25% mutant rate with melanin missing phenotype in larva, pupa,l and adult stage. The predation ability of the fourth instar larvae has no significant difference between wild (control) and mutant H. axyridis (G0), while these mutant fourth instar larvae had longer developmental period than that of the wild type. Consequently, the total predation of the fourth instar larvae was significantly increased in H. axyridis mutants comparing with the wild type. These results indicated that the success of CRISPR/Cas9 gene editing in H. axyridis. The gene editing platform in H. axyridis would facilitate the gene function research and promote special strain of predatory ladybird beetle generation.

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.

Microinjection-based CRISPR/Cas9 mutagenesis in the decapoda crustaceans, Neocaridina heteropoda and Eriocheir sinensis

CRISPR/Cas9 technology has been applied to many arthropods. However, application of this technology to crustaceans remains limited due to the unique characteristics of embryos. Our group has developed a microinjection system to introduce the CRISPR/Cas9 system into Neocaridina heteropoda embryos (one-cell stage). Using the developed method, we mutated the target gene Nh-scarlet (N. heteropoda scarlet), which functions in eye development and pigmentation. The results showed that both eye color and shape were altered in individuals in which Nh-scarlet was knocked out. Furthermore, this system was also successfully applied to another decapod crustacean, Eriocheir sinensis. DNA sequencing revealed that the zoeae with red eyes had an edited version of Es-scarlet. This study provides a stable microinjection method for freshwater crustaceans, and will contribute to functional genomics studies in various decapods.

The genome sequence of the painted lady, Vanessa cardui Linnaeus 1758

We present a genome assembly from an individual female Vanessa cardui (the painted lady; Arthropoda; Insecta; Lepidoptera; Nymphalidae). The genome sequence is 425 megabases in span. The majority of the assembly is scaffolded into 32 chromosomal pseudomolecules, with the W and Z sex chromosome assembled. Gene annotation of this assembly on Ensembl has identified 12,821 protein coding genes.

Transgenesis in the acoel worm Hofstenia miamia

The acoel worm Hofstenia miamia, which can replace tissue lost to injury via differentiation of a population of stem cells, has emerged as a new research organism for studying regeneration. To enhance the depth of mechanistic studies in this system, we devised a protocol for microinjection into embryonic cells that resulted in stable transgene integration into the genome and generated animals with tissue-specific fluorescent transgene expression in epidermis, gut, and muscle. We demonstrate that transgenic Hofstenia are amenable to the isolation of specific cell types, investigations of regeneration, tracking of photoconverted molecules, and live imaging. Further, our stable transgenic lines revealed insights into the biology of Hofstenia, including a high-resolution three-dimensional view of cell morphology and the organization of muscle as a cellular scaffold for other tissues. Our work positions Hofstenia as a powerful system with multiple toolkits for mechanistic investigations of development, whole-body regeneration, and stem cell biology.

Genetic knockout and general odorant-binding/chemosensory protein interactions: Revealing the function and importance of GOBP2 in the yellow peach moth's olfactory system

The GOBP2 protein has a unique function in the yellow peach moth (Conogethes punctiferalis (Guenée)). Several general odorant-binding proteins (GOBPs) have been identified in various lepidopteran species, but the functional difference between GOBP1 and GOBP2 in recognition of host plant odorants is still unknown. The functions of GOBP1 and GOBP2 in the yellow peach moth were evaluated in this study by using the CRISPR-Cas9 system. The results revealed the importance of GOBP2 in the olfaction mechanism in the yellow peach moth. The perception of the GOBP1-knockout larvae toward feeding decreased but did not reach a significant level while knocking out the GOBP2 and GOBP1/2 genes resulted in huge differences. On the other hand, electroantennograms (EAGs) and wind tunnel tests showed that the sensitivity of GOBP2 knockout adults to odorants decreased more than that of GOBP1 knockout individuals. The results of STRING database text mining grabbed our attention in protein-protein interaction studies. In this research, we first proved the existence of physical interactions between GOBPs and chemosensory proteins (CSPs) through the surface plasmon resonance (SPR) and isothermal titration calorimetry (ITC) methods. Interestingly, GOBP1 and GOBP2 could not interact with each other, but they could interact with CSPs. The interaction results indicated that GOBP2 could physically interact with CSP15, CSP5, and OBP17, whereas GOBP1 could bind only with CSP5 and CSP10, and its association constant (k_a) was also more substantial than that of GOBP1. These results strongly suggest the importance of the function of GOBP2 in the perception of host plant odorants by the yellow peach moth.

Non-molting dwarf (nm-d) as a mutant of Bombyx mori with a defect in purine synthesis

There are several known non-molting mutations of the silkworm, Bombyx mori, including non-molting dwarf (nm-d). Larvae with this mutation hatch normally and start eating leaves, but die before the completion of the first ecdysis. Genetic analysis of the nm-d mutation would contribute to the isolation of essential genes for the larval development of lepidopteran insects. To identify the causative gene of the nm-d locus, we conducted RNA-seq based rough mapping. Using two sets of RNA-seq data, one from a pooled sample of normal larvae, and one from a pooled sample of nm-d larvae, the nm-d locus was narrowed to a 500 kb region. Among the genes located in this region, a nm-d-specific exon loss was identified in the Bombyx homolog of the ATIC (5-aminoimidazole-4-carboxamide ribonucleotide transformylase/Inosine 5'-monophosphate cyclohydrolase) (BmATIC) gene, which catalyzes the final two steps of the de novo purine biosynthetic pathway in mammals. PCR and subsequent sequencing analysis revealed that a region containing exon 9 of the BmATIC gene is deleted in the nm-d larvae. A knockout allele of the BmATIC gene (BmATIC^KO), that was generated using the CRISPR/Cas9 system, revealed that first instar knockout larvae died while exhibiting the dark brown larval body that is a typical feature of mutants that lack uric acid in the integument. Lethal larvae resulted from crosses between +/BmATIC^KO moths. The uric acid content in the whole-body of the first instar was drastically reduced in the nm-d larvae compared to normal larvae. These results indicated that the BmATIC gene is responsible for the nm-d phenotype, and that nm-d larvae have a defect in purine biosynthesis, including uric acid. We also discuss the possibility that the BmATIC mRNA is maternally transmitted to eggs. Our results indicated that RNA-seq based mapping using pooled samples is a practical method for the identification of the causative genes of lethal mutations.

Comparative transcriptome analysis at the onset of speciation in a mimetic butterfly-The Ithomiini Melinaea marsaeus

Ecological speciation entails divergent selection on specific traits and ultimately on the developmental pathways responsible for these traits. Selection can act on gene sequences but also on regulatory regions responsible for gene expression. Mimetic butterflies are a relevant system for speciation studies because wing colour pattern (WCP) often diverges between closely related taxa and is thought to drive speciation through assortative mating and increased predation on hybrids. Here, we generate the first transcriptomic resources for a mimetic butterfly of the tribe Ithomiini, Melinaea marsaeus, to examine patterns of differential expression between two subspecies and between tissues that express traits that likely drive reproductive isolation; WCP and chemosensory genes. We sequenced whole transcriptomes of three life stages to cover a large catalogue of transcripts, and we investigated differential expression between subspecies in pupal wing discs and antennae. Eighteen known WCP genes were expressed in wing discs and 115 chemosensory genes were expressed in antennae, with a remarkable diversity of chemosensory protein genes. Many transcripts were differentially expressed between subspecies, including two WCP genes and one odorant receptor. Our results suggest that in M. marsaeus the same genes as in other mimetic butterflies are involved in traits causing reproductive isolation, and point at possible candidates for the differences in those traits between subspecies. Differential expression analyses of other developmental stages and body organs and functional studies are needed to confirm and expand these results. Our work provides key resources for comparative genomics in mimetic butterflies, and more generally in Lepidoptera.

Hox genes are essential for the development of eyespots in Bicyclus anynana butterflies

The eyespot patterns found on the wings of nymphalid butterflies are novel traits that originated first in hindwings and subsequently in forewings, suggesting that eyespot development might be dependent on Hox genes. Hindwings differ from forewings in the expression of Ultrabithorax (Ubx), but the function of this Hox gene in eyespot development as well as that of another Hox gene Antennapedia (Antp), expressed specifically in eyespots centers on both wings, are still unclear. We used CRISPR-Cas9 to target both genes in Bicyclus anynana butterflies. We show that Antp is essential for eyespot development on the forewings and for the differentiation of white centers and larger eyespots on hindwings, whereas Ubx is essential not only for the development of at least some hindwing eyespots but also for repressing the size of other eyespots. Additionally, Antp is essential for the development of silver scales in male wings. In summary, Antp and Ubx, in addition to their conserved roles in modifying serially homologous segments along the anterior-posterior axis of insects, have acquired a novel role in promoting the development of a new set of serial homologs, the eyespot patterns, in both forewings (Antp) and hindwings (Antp and Ubx) of B. anynana butterflies. We propose that the peculiar pattern of eyespot origins on hindwings first, followed by forewings, could be due to an initial co-option of Ubx into eyespot development followed by a later, partially redundant, co-option of Antp into the same network.

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.

Developmental Transcriptomics Reveals a Gene Network Driving Mimetic Color Variation in a Bumble Bee

A major goal of evolutionary genetics and evo-devo is to understand how changes in genotype manifest as changes in phenotype. Bumble bees display remarkable color pattern diversity while converging onto numerous regional Müllerian mimicry patterns, thus enabling exploration of the genetic mechanisms underlying convergent phenotypic evolution. In western North America, multiple bumble bee species converge onto local mimicry patterns through parallel shifts of midabdominal segments from red to black. It was previously demonstrated that a Hox gene, Abd-B, is the key regulator of the phenotypic switch in one of these species, Bombus melanopygus, however, the mechanism by which Abd-B regulates color differentiation remains unclear. Using tissue/stage-specific transcriptomic analysis followed by qRT–PCR validation, this study reveals a suite of genes potentially involved downstream of Abd-B during color pattern differentiation. The data support differential genes expression of not only the first switch gene Abd-B, but also an intermediate developmental gene nubbin, and a whole suite of downstream melanin and redox genes that together reinforce the observed eumelanin (black)-pheomelanin (red) ratios. These include potential genes involved in the production of insect pheomelanins, a pigment until recently not thought to occur in insects and thus lacking known regulatory enzymes. The results enhance understanding of pigmentation gene networks involved in bumble bee color pattern development and diversification, while providing insights into how upstream regulators such as Hox genes interact with downstream morphogenic players to facilitate this adaptive phenotypic radiation.

The Tiger Rattlesnake genome reveals a complex genotype underlying a simple venom phenotype

Variation in gene regulation is ubiquitous, yet identifying the mechanisms producing such variation, especially for complex traits, is challenging. Snake venoms provide a model system for studying the phenotypic impacts of regulatory variation in complex traits because of their genetic tractability. Here, we sequence the genome of the Tiger Rattlesnake, which possesses the simplest and most toxic venom of any rattlesnake species, to determine whether the simple venom phenotype is the result of a simple genotype through gene loss or a complex genotype mediated through regulatory mechanisms. We generate the most contiguous snake-genome assembly to date and use this genome to show that gene loss, chromatin accessibility, and methylation levels all contribute to the production of the simplest, most toxic rattlesnake venom. We provide the most complete characterization of the venom gene-regulatory network to date and identify key mechanisms mediating phenotypic variation across a polygenic regulatory network.

The Fractal Geometry of the Nymphalid Groundplan: Self-Similar Configuration of Color Pattern Symmetry Systems in Butterfly Wings

Simple Summary

Highly diverse color patterns of butterfly wings can be explained as modifications of an archetypical color pattern of nymphalid butterflies called the nymphalid groundplan. The nymphalid groundplan contains three major symmetry systems and a discal symmetry system, but their relationships have been elusive. Here, the morphological and spatial relationships among these symmetry systems were studied based on cross-species color-pattern comparisons of the hindwings in nymphalid butterflies. It was shown that all symmetry systems can be expressed as various structures, suggesting the equivalence (homology) of these systems in developmental potential. In some cases, the discal symmetry system is circularly surrounded by the central symmetry system, which may then be surrounded by the border and basal symmetry systems, indicating a unified supersymmetry system covering the entire wing. These results suggest that butterfly color patterns are hierarchically constructed; one system is nested within another system, which is a self-similar relationship that achieves the fractal geometry. This self-similarity is likely mediated by the serial induction of organizers during development, and a possible mechanism is proposed for symmetry breaking of the system morphology, which contributes to the diversity of butterfly wing color patterns.

Abstract

The nymphalid groundplan is an archetypical color pattern of nymphalid butterflies involving three major symmetry systems and a discal symmetry system, which share the basic morphogenesis unit. Here, the morphological and spatial relationships among these symmetry systems were studied based on cross-species comparisons of nymphalid hindwings. Based on findings in Neope and Symbrenthia, all three major symmetry systems can be expressed as bands, spots, or eyespot-like structures, suggesting equivalence (homology) of these systems in developmental potential. The discal symmetry system can also be expressed as various structures. The discal symmetry system is circularly surrounded by the central symmetry system, which may then be surrounded by the border and basal symmetry systems, based mainly on findings in Agrias, indicating a unified supersymmetry system covering the entire wing. The border symmetry system can occupy the central part of the wing when the central symmetry system is compromised, as seen in Callicore. These results suggest that butterfly color patterns are hierarchically constructed in a self-similar fashion, as the fractal geometry of the nymphalid groundplan. This self-similarity is likely mediated by the serial induction of organizers, and symmetry breaking of the system morphology may be generated by the collision of opposing signals during development.

The modular expression patterns of three pigmentation genes prefigure unique abdominal morphologies seen among three Drosophila species

To understand how novel animal body colorations emerged, one needs to ask how the development of color patterns differs among closely related species. Here we examine three species of fruit flies - Drosophila guttifera (D. guttifera), D. palustris, and D. subpalustris - displaying a varying number of abdominal spot rows. Through in situ hybridization experiments, we examine the mRNA expression patterns for the pigmentation genes Dopa decarboxylase (Ddc), tan (t), and yellow (y) during pupal development. Our results show that Ddc, t, and y are co-expressed in modular, identical patterns, each foreshadowing the adult abdominal spots in D. guttifera, D. palustris, and D. subpalustris. We suggest that differences in the expression patterns of these three genes partially underlie the morphological diversity of the quinaria species group.

A defect in purine nucleotide metabolism in the silkworm, Bombyx mori, causes a translucent larval integument and male infertility

p-oily (op) is a novel mutant of Bombyx mori exhibiting translucent larval integument and male infertility. Elucidation of the causative gene of the op mutant will help understand the genetic mechanism underlying larval integument coloration and male fertility. Using polymorphisms between B. mori and B. mandarina, the op locus was narrowed down to a 375-kb region. Using RNA-seq analysis, we found that op mutants have a frameshift mutation in the KWMTBOMO13770 gene located in the 375-kb region. A database search indicated that this gene is the human cytosolic 5'-nucleotidase II gene (cN-II) homolog in Bombyx, which mediates the conversion of inosine monophosphate (IMP) to inosine, a precursor of uric acid. CRISPR/Cas9-mediated knockout mutants of the Bm-cN-II gene showed translucent integuments, and there appeared translucent larvae in the crosses between knockout moths and +/op moths. Moreover, the translucent phenotype of, and decreased uric acid content in the larval integument caused by the mutations in the Bm-cN-II gene were rescued by oral administration of inosine. These results indicated that the Bm-cN-II gene is responsible for the op phenotype and that the molecular function of the Bm-cN-II gene is the conversion of IMP to inosine. We also discuss the genetic relationship between the Bm-cN-II gene and male fertility.

A release from developmental bias accelerates morphological diversification in butterfly eyespots

Development can bias the independent evolution of traits sharing ontogenetic pathways, making certain evolutionary changes less likely. The eyespots commonly found on butterfly wings each have concentric rings of differing colors, and these serially repeated pattern elements have been a focus for evo-devo research. In the butterfly family Nymphalidae, eyespots have been shown to function in startling or deflecting predators and to be involved in sexual selection. Previous work on a model species of Mycalesina butterfly, Bicyclus anynana, has provided insights into the developmental control of the size and color composition of individual eyespots. Experimental evolution has also shown that the relative size of a pair of eyespots on the same wing surface is highly flexible, whereas they are resistant to diverging in color composition, presumably due to the underlying shared developmental process. This fixed color composition has been considered as a prime example of developmental bias with significant consequences for wing pattern evolution. Here, we test this proposal by surveying eyespots across the whole subtribe of Mycalesina butterflies and demonstrate that developmental bias shapes evolutionary diversification except in the genus Heteropsis which has gained independent control of eyespot color composition. Experimental manipulations of pupal wings reveal that the bias has been released through a novel regional response of the wing tissue to a conserved patterning signal. Our study demonstrates that development can bias the evolutionary independence of traits, but it also shows how bias can be released through developmental innovations, thus, allowing rapid morphological change, facilitating evolutionary diversification.

Speciation in North American Junonia from a genomic perspective

Delineating species boundaries in phylogenetic groups undergoing recent radiation is a daunting challenge akin to discretizing continuity. Here, we propose a general approach exemplified by American butterflies from the genus Junonia Hübner notorious for the variety of similar phenotypes, ease of hybridization, and the lack of consensus about their classification. We obtain whole-genome shotgun sequences of about 200 specimens. We reason that discreteness emerges from continuity by means of a small number of key players, and search for the proteins that diverged markedly between sympatric populations of different species, while keeping low polymorphism within these species. Being 0.25% of the total number, these three dozen 'speciation' proteins indeed partition pairs of Junonia populations into two clusters with a prominent break in between, while all proteins taken together fail to reveal this discontinuity. Populations with larger divergence from each other, comparable to that between two sympatric species, form the first cluster and correspond to different species. The other cluster is characterized by smaller divergence, similar to that between allopatric populations of the same species and comprise conspecific pairs. Using this method, we conclude that J. genoveva (Cramer), J. litoralis Brévignon, J. evarete (Cramer), and J. divaricata C. & R. Felder are restricted to South America. We find that six species of Junonia are present in the United States, one of which is new: Junonia stemosa Grishin, sp.n. (i), found in south Texas and phenotypically closest to J. nigrosuffusa W. Barnes & McDunnough (ii) in its dark appearance. In the pale nudum of the antennal club, these two species resemble J. zonalis C. & R. Felder (iii) from Florida and the Caribbean Islands. The pair of sister species, J. grisea Austin & J. Emmel (iv) and J. coenia Hübner (v), represent the classic west/east U.S.A. split. The mangrove feeder (as caterpillar), dark nudum J. neildi Brévignon (vi) enters south Texas as a new subspecies Junonia neildi varia Grishin ssp.n. characterized by more extensive hybridization with and introgression from J. coenia, and, as a consequence, more variable wing patterns compared with the nominal J. n. neildi in Florida. Furthermore, a new mangrove-feeding species from the Pacific Coast of Mexico is described as Junonia pacoma Grishin sp.n. Finally, genomic analysis suggests that J. nigrosuffusa may be a hybrid species formed by the ancestors of J. grisea and J. stemosa sp.n.

Highly Efficient Knockout of a Squid Pigmentation Gene

Seminal studies using squid as a model led to breakthroughs in neurobiology. The squid giant axon and synapse, for example, laid the foundation for our current understanding of the action potential [1], ionic gradients across cells [2], voltage-dependent ion channels [3], molecular motors [4-7], and synaptic transmission [8-11]. Despite their anatomical advantages, the use of squid as a model receded over the past several decades as investigators turned to genetically tractable systems. Recently, however, two key advances have made it possible to develop techniques for the genetic manipulation of squid. The first is the CRISPR-Cas9 system for targeted gene disruption, a largely species-agnostic method [12, 13]. The second is the sequencing of genomes for several cephalopod species [14-16]. If made genetically tractable, squid and other cephalopods offer a wealth of biological novelties that could spur discovery. Within invertebrates, not only do they possess by far the largest brains, they also express the most sophisticated behaviors [17]. In this paper, we demonstrate efficient gene knockout in the squid Doryteuthis pealeii using CRISPR-Cas9. Ommochromes, the pigments found in squid retinas and chromatophores, are derivatives of tryptophan, and the first committed step in their synthesis is normally catalyzed by Tryptophan 2,3 Dioxygenase (TDO [18-20]). Knocking out TDO in squid embryos efficiently eliminated pigmentation. By precisely timing CRISPR-Cas9 delivery during early development, the degree of pigmentation could be finely controlled. Genotyping revealed knockout efficiencies routinely greater than 90%. This study represents a critical advancement toward making squid genetically tractable.

Divergent sensory investment mirrors potential speciation via niche partitioning across Drosophila

The examination of phylogenetic and phenotypic characteristics of the nervous system, such as behavior and neuroanatomy, can be utilized as a means to assess speciation. Recent studies have proposed a fundamental tradeoff between two sensory organs, the eye and the antenna. However, the identification of ecological mechanisms for this observed tradeoff have not been firmly established. Our current study examines several monophyletic species within the obscura group, and asserts that despite their close relatedness and overlapping ecology, they deviate strongly in both visual and olfactory investment. We contend that both courtship and microhabitat preferences support the observed inverse variation in these sensory traits. Here, this variation in visual and olfactory investment seems to provide relaxed competition, a process by which similar species can use a shared environment differently and in ways that help them coexist. Moreover, that behavioral separation according to light gradients occurs first, and subsequently, courtship deviations arise.

NPC1b as a novel target in controlling the cotton bollworm, Helicoverpa armigera

BACKGROUND

Insects cannot synthesize sterols and must acquire them from food. The mechanisms underlying how insects uptake dietary sterols are largely unknown except that NPC1b, an integral membrane protein, has been shown to be responsible for dietary cholesterol uptake in Drosophila melanogaster. However, whether NPC1b orthologs in other insect species, particularly the economically important pests, function similarly remains to be determined.

RESULTS

In this study, we characterized the function of NPC1b in Helicoverpa armigera, a global pest that causes severe yield losses to many important crops. Limiting dietary cholesterol uptake to insects significantly inhibited food ingestion and weight gain. Compared to the wild-type H. armigera, the CRISPR/Cas9-edited NPC1b mutant larvae were incapable of getting adequate cholesterol and died in their early life stage. Gene expression profile and in situ hybridization analyses indicated that NPC1b was mainly expressed in the midgut where dietary cholesterol was absorbed. Expression of NPC1b was also correlated with the feeding life stages and was especially upregulated during early larval instars. Protein-ligand docking and sequence similarity analyses further demonstrated that NPC1b proteins of lepidopteran insects shared a relatively conserved cholesterol binding region, NPC1b_NTD, which, however, was highly divergent from bees-derived sequences.

CONCLUSION

NPC1b was crucial for dietary cholesterol uptake and growth of H. armigera, and therefore could serve as an insecticide target for the development of a novel pest-management approach to control this economically significant insect pest with little off-target effect on bees and sterol-autotrophic animals. © 2020 Society of Chemical Industry.

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.

Transcriptome analysis reveals wingless regulates neural development and signaling genes in the region of wing pigmentation of a polka-dotted fruit fly

How evolutionary novelties have arisen is one of the central questions in evolutionary biology. Preexisting gene regulatory networks or signaling pathways have been shown to be co-opted for building novel traits in several organisms. However, the structure of entire gene regulatory networks and evolutionary events of gene co-option for emergence of a novel trait are poorly understood. In this study, to explore the genetic and molecular bases of the novel wing pigmentation pattern of a polka-dotted fruit fly (Drosophila guttifera), we performed de novo genome sequencing and transcriptome analyses. As a result, we comprehensively identified the genes associated with the pigmentation pattern. Furthermore, we revealed that 151 of these associated genes were positively or negatively regulated by wingless, a master regulator of wing pigmentation. Genes for neural development, Wnt signaling, Dpp signaling, and effectors (such as enzymes) for melanin pigmentation were included among these 151 genes. None of the known regulatory genes that regulate pigmentation pattern formation in other fruit fly species were included. Our results suggest that the novel pigmentation pattern of a polka-dotted fruit fly might have emerged through multistep co-options of multiple gene regulatory networks, signaling pathways, and effector genes, rather than recruitment of one large gene circuit.

Evolution of wing pigmentation in Drosophila: Diversity, physiological regulation, and cis-regulatory evolution

Fruit flies (Drosophila and its close relatives, or “drosophilids”) are a group that includes an important model organism, Drosophila melanogaster, and also very diverse species distributed worldwide. Many of these species have black or brown pigmentation patterns on their wings, and have been used as material for evo‐devo research. Pigmentation patterns are thought to have evolved rapidly compared with body plans or body shapes; hence they are advantageous model systems for studying evolutionary gains of traits and parallel evolution. Various groups of drosophilids, including genus Idiomyia (Hawaiian Drosophila), have a variety of pigmentations, ranging from simple black pigmentations around crossveins to a single antero‐distal spot and a more complex mottled pattern. Pigmentation patterns are sometimes obviously used for sexual displays; however, in some cases they may have other functions. The process of wing formation in Drosophila, the general mechanism of pigmentation formation, and the transport of substances necessary for pigmentation, including melanin precursors, through wing veins are summarized here. Lastly, the evolution of the expression of genes regulating pigmentation patterns, the role of cis‐regulatory regions, and the conditions required for the evolutionary emergence of pigmentation patterns are discussed. Future prospects for research on the evolution of wing pigmentation pattern formation in drosophilids are presented, particularly from the point of view of how they compare with other studies of the evolution of new traits.

Multiple Loci Control Eyespot Number Variation on the Hindwings of Bicyclus anynana Butterflies

The underlying genetic changes that regulate the appearance and disappearance of repeated traits, or serial homologs, remain poorly understood. One hypothesis is that variation in genomic regions flanking master regulatory genes, also known as input-output genes, controls variation in trait number, making the locus of evolution almost predictable. Another hypothesis implicates genetic variation in up- or downstream loci of master control genes. Here, we use the butterfly Bicyclus anynana, a species that exhibits natural variation in eyespot number on the dorsal hindwing, to test these two hypotheses. We first estimated the heritability of dorsal hindwing eyespot number by breeding multiple butterfly families differing in eyespot number and regressing eyespot numbers of offspring on midparent values. We then estimated the number and identity of independent genetic loci contributing to eyespot number variation by performing a genome-wide association study with restriction site-associated DNA sequencing from multiple individuals varying in number of eyespots sampled across a freely breeding laboratory population. We found that dorsal hindwing eyespot number has a moderately high heritability of ∼0.50 and is characterized by a polygenic architecture. Previously identified genomic regions involved in eyespot development, and novel ones, display high association with dorsal hindwing eyespot number, suggesting that homolog number variation is likely determined by regulatory changes at multiple loci that build the trait, and not by variation at single master regulators or input-output genes.

Butterfly eyespot color pattern formation requires physical contact of the pupal wing epithelium with extracellular materials for morphogenic signal propagation

BACKGROUND

Eyespot color pattern formation on butterfly wings is sensitive to physical damage and physical distortion as well as physical contact with materials on the surface of wing epithelial tissue at the pupal stage. Contact-mediated eyespot color pattern changes may imply a developmental role of the extracellular matrix in morphogenic signal propagation. Here, we examined eyespot responses to various contact materials, focusing on the hindwing posterior eyespots of the blue pansy butterfly, Junonia orithya.

RESULTS

Contact with various materials, including both nonbiological and biological materials, induced eyespot enlargement, reduction, or no change in eyespot size, and each material was characterized by a unique response profile. For example, silicone glassine paper almost always induced a considerable reduction, while glass plates most frequently induced enlargement, and plastic plates generally produced no change. The biological materials tested here (fibronectin, polylysine, collagen type I, and gelatin) resulted in various responses, but polylysine induced more cases of enlargement, similar to glass plates. The response profile of the materials was not readily predictable from the chemical composition of the materials but was significantly correlated with the water contact angle (water repellency) of the material surface, suggesting that the surface physical chemistry of materials is a determinant of eyespot size. When the proximal side of a prospective eyespot was covered with a size-reducing material (silicone glassine paper) and the distal side and the organizer were covered with a material that rarely induced size reduction (plastic film), the proximal side of the eyespot was reduced in size in comparison with the distal side, suggesting that signal propagation but not organizer activity was inhibited by silicone glassine paper.

CONCLUSIONS

These results suggest that physical contact with an appropriate hydrophobic surface is required for morphogenic signals from organizers to propagate normally. The binding of the apical surface of the epithelium with an opposing surface may provide mechanical support for signal propagation. In addition to conventional molecular morphogens, there is a possibility that mechanical distortion of the epithelium that is propagated mechanically serves as a nonmolecular morphogen to induce subsequent molecular changes, in accordance with the distortion hypothesis for butterfly wing color pattern formation.

Bringing immersive science to undergraduate laboratory courses using CRISPR gene knockouts in frogs and butterflies

The use of CRISPR/Cas9 for gene editing offers new opportunities for biology students to perform genuine research exploring the gene-to-phenotype relationship. It is important to introduce the next generation of scientists, health practitioners and other members of society to the technical and ethical aspects of gene editing. Here, we share our experience leading hands-on undergraduate laboratory classes, where students formulate hypotheses regarding the roles of candidate genes involved in development, perform loss-of-function experiments using programmable nucleases and analyze the phenotypic effects of mosaic mutant animals. This is enabled by the use of the amphibian Xenopus laevis and the butterfly Vanessa cardui, two organisms that reliably yield hundreds of large and freshly fertilized eggs in a scalable manner. Frogs and butterflies also present opportunities to teach key biological concepts about gene regulation and development. To complement these practical aspects, we describe learning activities aimed at equipping students with a broad understanding of genome editing techniques, their application in fundamental and translational research, and the bioethical challenges they raise. Overall, our work supports the introduction of CRISPR technology into undergraduate classrooms and, when coupled with classroom undergraduate research experiences, enables hypothesis-driven research by undergraduates.

How to study enhancers in non-traditional insect models

Transcriptional enhancers are central to the function and evolution of genes and gene regulation. At the organismal level, enhancers play a crucial role in coordinating tissue- and context-dependent gene expression. At the population level, changes in enhancers are thought to be a major driving force that facilitates evolution of diverse traits. An amazing array of diverse traits seen in insect morphology, physiology and behavior has been the subject of research for centuries. Although enhancer studies in insects outside of Drosophila have been limited, recent advances in functional genomic approaches have begun to make such studies possible in an increasing selection of insect species. Here, instead of comprehensively reviewing currently available technologies for enhancer studies in established model organisms such as Drosophila, we focus on a subset of computational and experimental approaches that are likely applicable to non-Drosophila insects, and discuss the pros and cons of each approach. We discuss the importance of validating enhancer function and evaluate several possible validation methods, such as reporter assays and genome editing. Key points and potential pitfalls when establishing a reporter assay system in non-traditional insect models are also discussed. We close with a discussion of how to advance enhancer studies in insects, both by improving computational approaches and by expanding the genetic toolbox in various insects. Through these discussions, this Review provides a conceptual framework for studying the function and evolution of enhancers in non-traditional insect models.

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.

Morphological Murals: The Scaling and Allometry of Butterfly Wing Patterns

The color patterns of butterflies moths are exceptionally diverse, but are very stable within a species, so that most species can be identified on the basis of their color pattern alone. The color pattern is established in the wing imaginal disc during a prolonged period of growth and differentiation, beginning during the last larval instar and ending during the first few days of the pupal stage. During this period, a variety of diffusion and reaction-diffusion signaling mechanisms determine the positions and sizes of the various elements that make up the overall color pattern. The patterning occurs while the wing is growing from a small imaginal disc to a very large pupal wing. One would therefore expect that some or all aspects of the color pattern would be sensitive to the size of the developmental field on which pattern formation takes place. To study this possibility, we analyzed the color patterns of Junonia coenia from animals whose growth patterns were altered by periodic starvation during larval growth, which produced individuals with a large range of variation in body size and wing size. Analyses of the color patterns showed that the positions and size of the pattern elements scaled perfectly isometrically with wing size. This is a puzzling finding and suggests the operation of a homeostatic or robustness mechanism that stabilizes pattern in spite of variation in the growth rate and final size of the wing.

Chromosomal-level reference genome of Chinese peacock butterfly (Papilio bianor) based on third-generation DNA sequencing and Hi-C analysis

BACKGROUND

Papilio bianor Cramer, 1777 (commonly known as the Chinese peacock butterfly) (Insecta, Lepidoptera, Papilionidae) is a widely distributed swallowtail butterfly with a wide number of geographic populations ranging from the southeast of Russia to China, Japan, India, Vietnam, Myanmar, and Thailand. Its wing color consists of both pigmentary colored scales (black, reddish) and structural colored scales (iridescent blue or green dust). A high-quality reference genome of P. bianor is an important foundation for investigating iridescent color evolution, phylogeography, and the evolution of swallowtail butterflies.

FINDINGS

We obtained a chromosome-level de novo genome assembly of the highly heterozygous P. bianor using long Pacific Biosciences sequencing reads and high-throughput chromosome conformation capture technology. The final assembly is 421.52 Mb on 30 chromosomes (29 autosomes and 1 Z sex chromosome) with 13.12 Mb scaffold N50. In total, 15,375 protein-coding genes and 233.09 Mb of repetitive sequences were identified. Phylogenetic analyses indicated that P. bianor separated from a common ancestor of swallowtails ∼23.69-36.04 million years ago. Demographic history suggested that the population expansion of this species from the last interglacial period to the last glacial maximum possibly resulted from its decreased natural enemies and its adaptation to climate change during the glacial period.

CONCLUSIONS

We present a high-quality chromosome-level reference genome of P. bianor using long-read single-molecule sequencing and Hi-C-based chromatin interaction maps. Our results lay the foundation for exploring the genetic basis of special biological features of P. bianor and also provide a useful data source for comparative genomics and phylogenomics among butterflies and moths.

Modular cis-regulatory logic of yellow gene expression in silkmoth larvae

Colour patterns in butterflies and moths are crucial traits for adaptation. Previous investigations have highlighted genes responsible for pigmentation (ie yellow and ebony). However, the mechanisms by which these genes are regulated in lepidopteran insects remain poorly understood. To elucidate this, molecular studies involving dipterans have largely analysed the cis-regulatory regions of pigmentation genes and have revealed cis-regulatory modularity. Here, we used well-developed transgenic techniques in Bombyx mori and demonstrated that cis-regulatory modularity controls tissue-specific expression of the yellow gene. We first identified which body parts are regulated by the yellow gene via black pigmentation. We then isolated three discrete regulatory elements driving tissue-specific gene expression in three regions of B. mori larvae. Finally, we found that there is no apparent sequence conservation of cis-regulatory regions between B. mori and Drosophila melanogaster, and no expression driven by the regulatory regions of one species when introduced into the other species. Therefore, the trans-regulatory landscapes of the yellow gene differ significantly between the two taxa. The results of this study confirm that lepidopteran species use cis-regulatory modules to control gene expression related to pigmentation, and represent a powerful cadre of transgenic tools for studying evolutionary developmental mechanisms.

Conservation and flexibility in the gene regulatory landscape of heliconiine butterfly wings

BACKGROUND

Many traits evolve by cis-regulatory modification, by which changes to noncoding sequences affect the binding affinity for available transcription factors and thus modify the expression profile of genes. Multiple examples of cis-regulatory evolution have been described at pattern switch genes responsible for butterfly wing pattern polymorphism, including in the diverse neotropical genus Heliconius, but the identities of the factors that can regulate these switch genes have not been identified.

RESULTS

We investigated the spatial transcriptomic landscape across the wings of three closely related butterfly species, two of which have a convergently evolved co-mimetic pattern and the other having a divergent pattern. We identified candidate factors for regulating the expression of wing patterning genes, including transcription factors with a conserved expression profile in all three species, and others, including both transcription factors and Wnt pathway genes, with markedly different profiles in each of the three species. We verified the conserved expression profile of the transcription factor homothorax by immunofluorescence and showed that its expression profile strongly correlates with that of the selector gene optix in butterflies with the Amazonian forewing pattern element 'dennis.'

CONCLUSION

Here we show that, in addition to factors with conserved expression profiles like homothorax, there are also a variety of transcription factors and signaling pathway components that appear to vary in their expression profiles between closely related butterfly species, highlighting the importance of genome-wide regulatory evolution between species.

Activation of butterfly eyespots by Distal-less is consistent with a reaction-diffusion process

Eyespots on the wings of nymphalid butterflies represent colorful examples of pattern formation, yet the developmental origins and mechanisms underlying eyespot center differentiation are still poorly understood. Using CRISPR-Cas9 we re-examine the function of Distal-less (Dll) as an activator or repressor of eyespots, a topic that remains controversial. We show that the phenotypic outcome of CRISPR mutations depends upon which specific exon is targeted. In Bicyclus anynana, exon 2 mutations are associated with both missing and ectopic eyespots, and also exon skipping. Exon 3 mutations, which do not lead to exon skipping, produce only null phenotypes, including missing eyespots, lighter wing coloration and loss of scales. Reaction-diffusion modeling of Dll function, using Wnt and Dpp as candidate morphogens, accurately replicates these complex crispant phenotypes. These results provide new insight into the function of Dll as a potential activator of eyespot development, scale growth and melanization, and suggest that the tuning of Dll expression levels can generate a diversity of eyespot phenotypes, including their appearance on the wing.This article has an associated 'The people behind the papers' interview.

Molecular Interactions Between Plants and Insect Herbivores

Diverse molecular processes regulate the interactions between plants and insect herbivores. Here, we review genes and proteins that are involved in plant-herbivore interactions and discuss how their discovery has structured the current standard model of plant-herbivore interactions. Plants perceive damage-associated and, possibly, herbivore-associated molecular patterns via receptors that activate early signaling components such as Ca²⁺, reactive oxygen species, and MAP kinases. Specific defense reprogramming proceeds via signaling networks that include phytohormones, secondary metabolites, and transcription factors. Local and systemic regulation of toxins, defense proteins, physical barriers, and tolerance traits protect plants against herbivores. Herbivores counteract plant defenses through biochemical defense deactivation, effector-mediated suppression of defense signaling, and chemically controlled behavioral changes. The molecular basis of plant-herbivore interactions is now well established for model systems. Expanding molecular approaches to unexplored dimensions of plant-insect interactions should be a future priority.

A genetic switch for worker nutrition-mediated traits in honeybees

Highly social insects are characterized by caste dimorphism, with distinct size differences of reproductive organs between fertile queens and the more or less sterile workers. An abundance of nutrition or instruction via diet-specific compounds has been proposed as explanations for the nutrition-driven queen and worker polyphenism. Here, we further explored these models in the honeybee (Apis mellifera) using worker nutrition rearing and a novel mutational screening approach using the clustered regularly interspaced short palindromic repeats/CRISPR-associated protein 9 (CRISPR/Cas9) method. The worker nutrition-driven size reduction of reproductive organs was restricted to the female sex, suggesting input from the sex determination pathway. Genetic screens on the sex determination genes in genetic females for size polyphenism revealed that doublesex (dsx) mutants display size-reduced reproductive organs irrespective of the sexual morphology of the organ tissue. In contrast, feminizer (fem) mutants lost the response to worker nutrition-driven size control. The first morphological worker mutants in honeybees demonstrate that the response to nutrition relies on a genetic program that is switched “ON” by the fem gene. Thus, the genetic instruction provided by the fem gene provides an entry point to genetically dissect the underlying processes that implement the size polyphenism.

In honeybees, nutrition drives dimorphic size development of reproductive organs in fertile queens and sterile workers. A study using the first induced morphological mutants in honeybees demonstrates that this developmental plasticity requires a genetic program that is switched on by the “feminizer” gene.

In honeybees, nutrition drives dimorphic size development of reproductive organs in fertile queens and sterile workers. The first induced morphological mutants in honeybees demonstrate that this developmental plasticity requires a genetic program that is switched “ON” by the feminizer (fem) gene.

Identification of a novel function of the silkworm integument in nitrogen metabolism: Uric acid is synthesized within the epidermal cells in B. mori

During nitrogen metabolism, animals convert toxic ammonia to less toxic forms. Uric acid (UA) is an end product of this process in terrestrial insects. In lepidopteran larvae, a large amount of UA is stored in the integument via a phenomenon known as storage excretion. Physiologically, integumental UA plays crucial roles as a barrier against sunlight and as a white pigment for larval pigmentation patterns. Conventionally, UA is thought to be synthesized in the fat body, the insect equivalent of the liver of vertebrates, and to be transported to the epidermis via the hemolymph. Here, we reconsidered the conventional theory by a mosaic analysis targeting genes governing UA synthesis, using CRISPR/Cas9 mutagenesis and a traditional genetic method in Bombyx mori. Notably, we observed mosaic larvae in which the integument comprised both UA-containing white and UA-lacking translucent areas, indicating that UA synthesis in the epidermis is indispensable to the accumulation of a large amount of highly insoluble UA in the epidermis. Our results thus provide a genetic basis for storage excretion wherein lepidopteran insects use nitrogenous waste to adapt to their environment.

Long-Range Effects of Wing Physical Damage and Distortion on Eyespot Color Patterns in the Hindwing of the Blue Pansy Butterfly Junonia orithya

Butterfly eyespot color patterns have been studied using several different approaches, including applications of physical damage to the forewing. Here, damage and distortion experiments were performed, focusing on the hindwing eyespots of the blue pansy butterfly Junonia orithya. Physical puncture damage with a needle at the center of the eyespot reduced the eyespot size. Damage at the eyespot outer rings not only deformed the entire eyespot, but also diminished the eyespot core disk size, despite the distance from the damage site to the core disk. When damage was inflicted near the eyespot, the eyespot was drawn toward the damage site. The induction of an ectopic eyespot-like structure and its fusion with the innate eyespots were observed when damage was inflicted in the background area. When a small stainless ball was placed in close proximity to the eyespot using the forewing-lift method, the eyespot deformed toward the ball. Taken together, physical damage and distortion elicited long-range inhibitory, drawing (attracting), and inducing effects, suggesting that the innate and induced morphogenic signals travel long distances and interact with each other. These results are consistent with the distortion hypothesis, positing that physical distortions of wing tissue contribute to color pattern determination in butterfly wings.

CRISPR/Cas9 deletions in a conserved exon of Distal-less generates gains and losses in a recently acquired morphological novelty in flies

Distal-less has been repeatedly co-opted for the development of many novel traits. Here, we document its curious role in the development of a novel abdominal appendage (“sternite brushes”) in sepsid flies. CRISPR/Cas9 deletions in the homeodomain result in losses of sternite brushes, demonstrating that Distal-less is necessary for their development. However, deletions in the upstream coding exon (Exon 2) produce losses or gains of brushes. A dissection of Exon 2 reveals that the likely mechanism for gains involves a deletion in an exon-splicing enhancer site that leads to exon skipping. Such contradictory phenotypes are also observed in butterflies, suggesting that mutations in the conserved upstream regions have the potential to generate phenotypic variability in insects that diverged 300 million years ago. Our results demonstrate the importance of Distal-less for the development of a novel abdominal appendage in insects and highlight how site-specific mutations in the same exon can produce contradictory phenotypes.

•

Distal-less is necessary for the development of a novel abdominal appendage

•

CRISPR/Cas9 editing produced both losses and gains of novel abdominal appendages

•

Gains of appendages result from mutations in exonic splicing enhancer (ESEs) sites

•

ESE mutations likely led to exon skipping and an altered Distal-less protein

Developmental dynamics of butterfly wings: real-time in vivo whole-wing imaging of twelve butterfly species

Colour pattern development of butterfly wings has been studied from several different approaches. However, developmental changes in the pupal wing tissues have rarely been documented visually. In this study, we recorded real-time developmental changes of the pupal whole wings of 9 nymphalid, 2 lycaenid, and 1 pierid species in vivo, from immediately after pupation to eclosion, using the forewing-lift method. The developmental period was roughly divided into four sequential stages. At the very early stage, the wing tissue was transparent, but at the second stage, it became semi-transparent and showed dynamic peripheral adjustment and slow low-frequency contractions. At this stage, the wing peripheral portion diminished in size, but simultaneously, the ventral epithelium expanded in size. Likely because of scale growth, the wing tissue became deeply whitish at the second and third stages, followed by pigment deposition and structural colour expression at the fourth stage. Some red or yellow (light-colour) areas that emerged early were “overpainted” by expanding black areas, suggesting the coexistence of two morphogenic signals in some scale cells. The discal spot emerged first in some nymphalid species, as though it organised the entire development of colour patterns. These results indicated the dynamic wing developmental processes common in butterflies.