Periodic Wnt1 expression in response to ecdysteroid generates twin-spot markings on caterpillars

Musashi orchestrates melanism in Laodelphax striatellus

In insects, melanism, a fundamental pigmentation process, is of significant importance in evolutionary biology due to its complex genetic foundation. We investigated the role of the RNA-binding gene Musashi (msi) in melanism in Laodelphax striatellus, a Hemiptera species. We identified a single L. striatellus msi homolog, Lsmsi, encoding a 357 amino acid protein with 2 RNA recognition motifs. RNA interference-mediated knockdown of LsMsi resulted in complete body melanism and increased cuticular permeability. Additionally, we found the involvement of G protein-coupled receptor A42 and tyrosine hydroxylase (Th) in L. striatellus melanism. Knockdown of LsTh lightened the epidermis, showing dehydration signs, while LsA42 knockdown enhanced LsTh expression, leading to melanism. Surprisingly, Lsmsi knockdown decreased both LsA42 and LsTh expression, which was expected to cause whitening but resulted in melanism. Further, we found that Lsmsi influenced downstream genes like phenoloxidase homolog LsPo and dopa decarboxylase (Ddc) homolog LsDdc in the tyrosine-mediated melanism pathway. Extending to Nilaparvata lugens and Sogatella furcifera, we demonstrated the conserved role of msi in melanism among Delphacidae. Given MSI proteins' roles in cancer and tumors in vertebrates, our study is the first to link msi in insects to Delphacidae body color melanization via the tyrosine-mediated pathway, offering fresh perspectives on the genetic basis of insect melanism and msi functions.

The BTB-ZF gene Bm-mamo regulates pigmentation in silkworm caterpillars

The color pattern of insects is one of the most diverse adaptive evolutionary phenotypes. However, the molecular regulation of this color pattern is not fully understood. In this study, we found that the transcription factor Bm-mamo is responsible for black dilute (bd) allele mutations in the silkworm. Bm-mamo belongs to the BTB zinc finger family and is orthologous to mamo in Drosophila melanogaster. This gene has a conserved function in gamete production in Drosophila and silkworms and has evolved a pleiotropic function in the regulation of color patterns in caterpillars. Using RNAi and clustered regularly interspaced short palindromic repeats (CRISPR) technology, we showed that Bm-mamo is a repressor of dark melanin patterns in the larval epidermis. Using in vitro binding assays and gene expression profiling in wild-type and mutant larvae, we also showed that Bm-mamo likely regulates the expression of related pigment synthesis and cuticular protein genes in a coordinated manner to mediate its role in color pattern formation. This mechanism is consistent with the dual role of this transcription factor in regulating both the structure and shape of the cuticle and the pigments that are embedded within it. This study provides new insight into the regulation of color patterns as well as into the construction of more complex epidermal features in some insects.

High-throughput and genome-scale targeted mutagenesis using CRISPR in a nonmodel multicellular organism, Bombyx mori

Large-scale genetic mutant libraries are powerful approaches to interrogating genotype-phenotype correlations and identifying genes responsible for certain environmental stimuli, both of which are the central goal of life science study. We produced the first large-scale CRISPR-Cas9-induced library in a nonmodel multicellular organism, Bombyx mori We developed a piggyBac-delivered binary genome editing strategy, which can simultaneously meet the requirements of mixed microinjection, efficient multipurpose genetic operation, and preservation of growth-defect lines. We constructed a single-guide RNA (sgRNA) plasmid library containing 92,917 sgRNAs targeting promoters and exons of 14,645 protein-coding genes, established 1726 transgenic sgRNA lines following microinjection of 66,650 embryos, and generated 300 mutant lines with diverse phenotypic changes. Phenomic characterization of mutant lines identified a large set of genes responsible for visual phenotypic or economically valuable trait changes. Next, we performed pooled context-specific positive screens for tolerance to environmental pollutant cadmium exposure, and identified KWMTBOMO12902 as a strong candidate gene for breeding applications in sericulture industry. Collectively, our results provide a novel and versatile approach for functional B. mori genomics, as well as a powerful resource for identifying the potential of key candidate genes for improving various economic traits. This study also shows the effectiveness, practicality, and convenience of large-scale mutant libraries in other nonmodel organisms.

SilkMeta: a comprehensive platform for sharing and exploiting pan-genomic and multi-omic silkworm data

The silkworm Bombyx mori is a domesticated insect that serves as an animal model for research and agriculture. The silkworm super-pan-genome dataset, which we published last year, is a unique resource for the study of global genomic diversity and phenotype-genotype association. Here we present SilkMeta (http://silkmeta.org.cn), a comprehensive database covering the available silkworm pan-genome and multi-omics data. The database contains 1082 short-read genomes, 546 long-read assembled genomes, 1168 transcriptomes, 294 phenotype characterizations (phenome), tens of millions of variations (variome), 7253 long non-coding RNAs (lncRNAs), 18 717 full length transcripts and a set of population statistics. We have compiled publications on functional genomics research and genetic stock deciphering (mutant map). A range of bioinformatics tools is also provided for data visualization and retrieval. The large batch of omics data and tools were integrated in twelve functional modules that provide useful strategies and data for comparative and functional genomics research. The interactive bioinformatics platform SilkMeta will benefit not only the silkworm but also the insect biology communities.

Disruption of a BTB-ZF transcription factor causes female sterility and melanization in the larval body of the silkworm, Bombyx mori

The dilute black (bd) of the silkworm Bombyx mori is a recessive mutant that produces a grayish-black color in the larval integument, instead of the characteristic white color found in wild-type larvae. In addition, eggs produced by bd females are sterile due to a deficiency in the micropylar apparatus. We identified candidate genes responsible for the bd phenotype using publicly available RNA-seq data. One of these candidate genes was homologous to the maternal gene required for meiosis (mamo) of Drosophila melanogaster, which encodes a broad-complex, tramtrack, and bric-à-brac-zinc finger (BTB-ZF) transcription factor essential for female fertility. In three independent bd strains, the expression of the B. mori mamo (Bmmamo) was downregulated in the larval integument. Using a CRISPR/Cas9-mediated knockout strategy, we found that Bmmamo knockout mutants exhibit a grayish-black color in the larval integument and female infertility. Moreover, larvae obtained from the complementation cross between bd/+ mutants and heterozygous knockouts for the Bmmamo also exhibited a grayish-black color, indicating that Bmmamo is responsible for the bd phenotype. Gene expression analysis using Bmmamo knockout mutants suggested that the BmMamo protein suppresses the expression of melanin synthesis genes. Previous comparative genome analysis revealed that the Bmmamo was selected during silkworm domestication, and we found that Bmmamo expression in the larval integument is higher in B. mori than in the wild silkworm B. mandarina, suggesting that the Bmmamo is involved in domestication-associated pigmentation changes of the silkworm.

A high-quality, long-read genome assembly of the whitelined sphinx moth (Lepidoptera: Sphingidae: Hyles lineata) shows highly conserved melanin synthesis pathway genes

The sphinx moth genus Hyles comprises 29 described species inhabiting all continents except Antarctica. The genus diverged relatively recently (40-25 MYA), arising in the Americas and rapidly establishing a cosmopolitan distribution. The whitelined sphinx moth, Hyles lineata, represents the oldest extant lineage of this group and is one of the most widespread and abundant sphinx moths in North America. Hyles lineata exhibits the large body size and adept flight control characteristic of the sphinx moth family (Sphingidae), but it is unique in displaying extreme larval color variation and broad host plant use. These traits, in combination with its broad distribution and high relative abundance within its range, have made H. lineata a model organism for studying phenotypic plasticity, plant-herbivore interactions, physiological ecology, and flight control. Despite being one of the most well-studied sphinx moths, little data exist on genetic variation or regulation of gene expression. Here, we report a high-quality genome showing high contiguity (N50 of 14.2 Mb) and completeness (98.2% of Lepidoptera BUSCO genes), an important first characterization to facilitate such studies. We also annotate the core melanin synthesis pathway genes and confirm that they have high sequence conservation with other moths and are most similar to those of another, well-characterized sphinx moth, the tobacco hornworm (Manduca sexta).

Abnormal overexpression of SoxD enhances melanin synthesis in the Ursa mutant of Bombyx mori

The pigment and structural color of insects play crucial roles in body protection, ecological adaptation, and signal communication. Epidermal melanization is a common and main coloring pattern, which results in broad phenotypic diversity. Melanin is one of the compounds contributing to dark brown-black pigmentation, which is synthesized from dopamine and tyrosine by the melanin metabolism pathway. The Ursa mutant of the silkworm Bombyx mori is a body-color mutant characterized by excessive melanin pigmentation in the larval epidermis. However, the exact gene responsible for this phenotype remains unclear. Here, we performed positional cloning of the gene responsible for Ursa, which was mapped to an 83-kb region on chromosome 14. The genomic region contains a protein-coding gene encoding a transcription factor, which was designated BmSoxD. The mutation site was determined by analysis of nucleotide sequences of the genomic region corresponding to BmSoxD, which identified a 449-bp transposable sequence similar to that of the B. mori transposon Helitron inserted into the sixth intron. BmSoxD was dramatically overexpressed in the epidermis of Ursa at the end of the molting stage compared with that of wild-type B. mori. Overexpression of BmSoxD led to upregulation of genes involved in the melanin metabolism pathway, whereas knocking down BmSoxD via small interfering RNAs blocked melanin pigment production in the larval epidermis. These data indicate that the mutation in BmSoxD is responsible for the Ursa mutant phenotype. We propose that the transposable sequence insertion causes abnormal overexpression of BmSoxD at the molting stage in the Ursa mutant, resulting in excessive melanin synthesis and its accumulation in epidermal cells.

High-resolution silkworm pan-genome provides genetic insights into artificial selection and ecological adaptation

The silkworm Bombyx mori is an important economic insect for producing silk, the “queen of fabrics”. The currently available genomes limit the understanding of its genetic diversity and the discovery of valuable alleles for breeding. Here, we deeply re-sequence 1,078 silkworms and assemble long-read genomes for 545 representatives. We construct a high-resolution pan-genome dataset representing almost the entire genomic content in the silkworm. We find that the silkworm population harbors a high density of genomic variants and identify 7308 new genes, 4260 (22%) core genes, and 3,432,266 non-redundant structure variations (SVs). We reveal hundreds of genes and SVs that may contribute to the artificial selection (domestication and breeding) of silkworm. Further, we focus on four genes responsible, respectively, for two economic (silk yield and silk fineness) and two ecologically adaptive traits (egg diapause and aposematic coloration). Taken together, our population-scale genomic resources will promote functional genomics studies and breeding improvement for silkworm.

The Landscapes of Full-Length Transcripts and Splice Isoforms as Well as Transposons Exonization in the Lepidopteran Model System, Bombyx mori

The domesticated silkworm, Bombyx mori, is an important model system for the order Lepidoptera. Currently, based on third-generation sequencing, the chromosome-level genome of Bombyx mori has been released. However, its transcripts were mainly assembled by using short reads of second-generation sequencing and expressed sequence tags which cannot explain the transcript profile accurately. Here, we used PacBio Iso-Seq technology to investigate the transcripts from 45 developmental stages of Bombyx mori. We obtained 25,970 non-redundant high-quality consensus isoforms capturing ∼60% of previous reported RNAs, 15,431 (∼47%) novel transcripts, and identified 7,253 long non-coding RNA (lncRNA) with a large proportion of novel lncRNA (∼56%). In addition, we found that transposable elements (TEs) exonization account for 11,671 (∼45%) transcripts including 5,980 protein-coding transcripts (∼32%) and 5,691 lncRNAs (∼79%). Overall, our results expand the silkworm transcripts and have general implications to understand the interaction between TEs and their host genes. These transcripts resource will promote functional studies of genes and lncRNAs as well as TEs in the silkworm.

The evolution and genetics of lepidopteran egg and caterpillar coloration

Insect colors and color patterns have fascinated biologists for centuries. While extensive research has focused on the adult colors of Drosophila and butterflies, our understanding of how colors are generated and diversified in embryonic and larval stages remains limited, especially, the genetics behind the protective coloration of the immobile embryonic and larval stages. Lepidoptera, one of the most widespread and species-rich insect orders, are extremely helpful uncovering those mechanisms due to their remarkable diverse colors in eggs and caterpillars within or among species, and these colors usually are variable in different developmental stages or in response to different environments. Here we review the recent progress on coloration of lepidopteran eggs and caterpillars, focusing on the genetic basis, developmental mechanisms, ecology, and evolution underlying the remarkable color diversity.

The Number of Larval Molts Is Controlled by Hox in Caterpillars

Animals with exoskeletons molt for further growth. In insects, the number of larval (or nymphal) molts varies inter- and intra-specifically, and it is widely accepted that the variation in the number of larval molts is an adaptive response to diverse environmental conditions.^1-5 However, the molecular mechanism that underlies the variety and plasticity in the number of larval molts is largely unknown. In the silkworm, Bombyx mori, there are strains that molt three, four, or five times, and these numbers are determined by allelic variation at a single autosomal locus, Moltinism (M).^6-9 Here, we demonstrate that the Hox gene Sex combs reduced (Scr) is responsible for the phenotypes of the M locus. Scr is selectively expressed in the larval prothoracic gland (PG), an endocrine organ that produces molting hormones.²Scr represses the biosynthesis of molting hormones in the PG, thereby regulating the incremental increase in body size during each larval instar. Our experiments consistently suggest that the differential expression levels of Scr among the three M alleles result in different growth ratios that ultimately lead to the different number of larval molts. Although the role of Hox genes in conferring segmental identity along the body axis and in molding segment-specific structure later in development has been well established,^10-13 the present study identifies an unexpected role of Hox gene in hormone biosynthesis. This new role means that, in addition to shaping segment-specific morphology, Hox genes also drive the evolution of life history traits by regulating animal physiology.

Functional conservation and diversification of yellow-y in lepidopteran insects

The diverse colors and patterns found in Lepidoptera are important for success of these species. Similar to the wings of adult butterflies, lepidopteran larvae exhibit diverse color variations to adapt to their habitats. Compared with butterfly wings, however, less attention has been paid to larval body colorations and patterns. In the present study, we focus on the yellow-y gene, which participates in the melanin synthesis pathway. We conducted CRISPR/Cas9-mediated targeted mutagenesis of yellow-y in the tobacco cutworm Spodoptera litura. We analyzed the role of S. litura yellow-y in pigmentation by morphological observation and discovered that yellow-y is necessary for normal black pigmentation in S. litura. We also showed species- and tissue-specific requirements of yellow-y in pigmentation in comparison with those of Bombyx mori yellow-y mutants. Furthermore, we found that almost none of the yellow-y mutant embryos hatched unaided. We provide evidence that S. litura yellow-y has a novel important function in egg hatching, in addition to pigmentation. The present study will enable a greater understanding of the functions and diversification of the yellow-y gene in insects.

Electroporation-mediated nucleic acid delivery during non-embryonic stages for gene-function analysis in Anopheles sinensis

The delivery of exogenous nucleic acids to eggs or non-embryonic individuals by microinjection is a vital reverse genetics technique used to determine gene function in mosquitoes. However, DNA delivery to eggs is complex and time-consuming, and conventional, non-embryonic-injection techniques may result in unobvious phenotypes caused by insufficient absorption of nucleic acid fragments by cells at target body parts or tissues. In this study, we developed a set of electroporation-mediated non-embryonic microinjections for the delivery of exogenous nucleic acids in Anopheles sinensis. Gene silencing using this method led to down-regulation of target gene expression (AsCPR128) by 77% in targeted body parts, compared with only 10% in non-targeted body parts, thus increasing the defect-phenotype rate in the target area by 5.3-fold, compared with non-shock injected controls. Electroporation-mediated somatic transgenesis resulted in stable phenotypic characteristics of the reporter gene at the shocked body parts during the pupal-adult stages in about 69% of individuals. Furthermore, injecting plasmid DNA near the ovaries of female mosquitoes after a blood meal followed by electric shock produced three germline G₁ transgenic lines, with a transformation rate of about 11.1% (calculated from ovulatory G₀ females). Among the positive G₁ lines, 42%, 40%, and 31% of individuals emitted red fluorescence in the larval stage. When the red fluorescent larvae developed into adults, green fluorescence was emitted from the ovaries of the females upon feeding. These results suggest that electroporation-mediated non-embryonic microinjection can be an efficient, rapid, and simple technique for analyzing gene function in non-model mosquitoes or other small insects.

Imaginal disc growth factor maintains cuticle structure and controls melanization in the spot pattern formation of Bombyx mori

The complex stripes and patterns of insects play key roles in behavior and ecology. However, the fine-scale regulation mechanisms underlying pigment formation and morphological divergence remain largely unelucidated. Here we demonstrated that imaginal disc growth factor (IDGF) maintains cuticle structure and controls melanization in spot pattern formation of Bombyx mori. Moreover, our knockout experiments showed that IDGF is suggested to impact the expression levels of the ecdysone inducible transcription factor E75A and pleiotropic factors apt-like and Toll8/spz3, to further control the melanin metabolism. Furthermore, the untargeted metabolomics analyses revealed that BmIDGF significantly affected critical metabolites involved in phenylalanine, beta-alanine, purine, and tyrosine metabolism pathways. Our findings highlighted not only the universal function of IDGF to the maintenance of normal cuticle structure but also an underexplored space in the gene function affecting melanin formation. Therefore, this study furthers our understanding of insect pigment metabolism and melanin pattern polymorphisms.

The diverse stripe patterns of animals are usually used for warning or camouflage. However, the actual mechanisms underlying diverse stripe pattern formation remains largely unknown. This study provides direct evidence that imaginal disc growth factor (IDGF) maintains cuticle structure and controls melanization in the spot pattern formation. Our exhaustive knockout experiments reveal that BmIDGF is involved in the melanin pigmentation of Bombyx mori. We demonstrate that IDGF impacts the expression levels of the 20E-inducible transcription factor E75A and pleiotropic factors apt-like and Toll8/spz3, to further affect the melanin metabolism. Furthermore, the metabolome of BmIDGF gene deletion connects metabolism to gene function. Thus, this study shed light on not only the unique function of IDGF but also the molecular mechanism of spot pattern formation.

Notch and Delta Control the Switch and Formation of Camouflage Patterns in Caterpillars

In most Papilio species, a younger larva mimics bird droppings but changes its pattern to match host plant colors in its final instar. This change is determined by juvenile hormone (JH) during the JH-sensitive period (JHSP) early in the fourth instar. Recently, we found that homeobox genes control the pre-pattern formation specifically during JHSP, but the molecular mechanisms underlying final patterning and pigmentation at molt are unknown. By knockdown of Delta and Notch in Papilio xuthus larvae, we here showed that these genes define the edge and pigmentation area in final patterns, during and even after JHSP, suggesting that they bridge the JHSP and molt. Knockdown of Delta in Papilio machaon led to similar phenotypic changes, and knockdown of Notch caused pigmentation loss in twin spots of the silkworm Multilunar (L) mutant. Our findings suggest the importance of the Notch signaling pathway in caterpillars' adaptive evolution of color pattern formation.

•

Notch and its ligand Delta regulate camouflage patterns of caterpillars

•

They define edge and pigmentation area in Papilio xuthus final larval patterns

•

They are suggested to bridge the juvenile hormone response period and final molt

•

Notch signaling pathway is important for caterpillars' color pattern evolution

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.

CRISPR/Cas9-based mutation reveals Argonaute 1 is essential for pigmentation in Ostrinia furnacalis

Ostrinia furnacalis (Lepidoptera: Pyralidae) is one of the most destructive agricultural pests in Asia. Traditional pest-management methods include sex pheromone capture, transgenic crops that produce Bacillus thuringiensis toxin, and pesticides. Although these strategies control pest populations effectively, they also cause negative side effects, including dramatically increased pesticide resistance, severe pollution, and hazards for human health. Recently developed genome editing tools provide new prospects for pest management and have been successfully used in several species. However, few examples have been reported in the agricultural pest O. furnacalis due to a lack in genomic information. In this report, we identified only one transcript of O. furnacalis Argonaute 1 (OfAgo1) gene from the genome and cloned the open reading frame. OfAgo1 presented the maximum expression at the embryo stage or in the fat body during the larval stages. To understand its function, an OfAgo1 mutant was constructed using the Clustered Regularly Interspaced Short Palindromic Repeat/RNA-guided Cas9 nuclease (CRISPR/Cas9). Mutagenesis of OfAgo1 disrupted cuticle pigmentation by down-regulating micro RNAs and pigmentation-related genes. This is the first report for the cloning and functional analysis of OfAgo1, revealing a role of OfAgo1 in cuticle pigmentation. The current report also established a CRISPR/Cas9 system in O. furnacalis, providing a new insight for pest management.

The silkworm (Bombyx mori) neuropeptide orcokinin is involved in the regulation of pigmentation

The natural colorful cuticles of insects play important roles in many physiological processes. Pigmentation is a physiological process with a complex regulatory network whose regulatory mechanism remains unclear. Bombyx mori pigmentation mutants are ideal materials for research on pigmentation mechanisms. The purple quail-like (q-l^p) and brown quail-like (q-l^b) mutants originated from plain silkworm breeds 932VR and 0223JH respectively exhibit similar cuticle pigmentation to that of the quail mutant. The q-l^p mutant also presents a developmental abnormality. In this study, genes controlling q-l^p and q-l^b mutants were located on chromosome 8 by positional cloning. Then the neuropeptide gene orcokinin (OK) was identified to be the major gene responsible for two quail-like mutants. The B. mori orcokinin gene (BommoOK) produces two transcripts, BommoOKA and BommoOKB, by alternative splicing. The CRISPR/Cas9 system and orcokinin peptides injection were used for further functional verification. We show a novel function of BommoOKA in inhibiting pigmentation, and one mature peptide of orcokinin A, OKA_type2, is the key factor in pigmentation inhibition. These results provide a reference for studying the function of orcokinin and are of theoretical importance for studying the regulatory mechanism of pigmentation.

Comparative analysis of integument transcriptomes identifies genes that participate in marking pattern formation in three allelic mutants of silkworm, Bombyx mori

The diversity markings and pigment patterns in insects are outcomes of adaptive evolution. The elucidation of the molecular mechanism underlying variations in pigment patterns may improve our understanding of the origin and evolution of these spectacular diverse phenotypes. Melanin, ommochrome, and pteridine are the three main types of insect pigments, and the genes that directly participate in pigment biosynthesis have been extensively studied. However, available information on gene interactions and the whole pigment regulatory network is limited. In this study, we performed integument transcriptome sequencing to analyze three larval marking allelic mutants, namely, multi lunar (L), L^C, and L^Ca, which have similar twin-spot markings on the dorsal side of multiple segments. Further analysis identified 336 differentially expressed genes (DEGs) between L and Dazao (wild type which exhibits normal markings), 68 DEGs between L^C/+ and +^LC/+^LC, and 188 DEGs between L^Ca/+ and +^LCa/+^LCa. Gene Ontology (GO) analysis indicated a significant DEG enrichment of the functional terms catalytic activity, binding, metabolic process, and cellular process. Furthermore, three mutants share six common enriched KEGG pathways. We finally identified eight common DEGs among three pairwise comparisons, including Krueppel-like factor, TATA-binding protein, protein patched, UDP-glycosyltransferase, an unknown secreted protein, and three cuticular proteins. Microarray-based gene expression analysis revealed that the eight genes are upregulated during molting, which coincides with marking formation, and are significantly differentially expressed between marking and non-marking regions. The results suggest that the eight common genes are involved in the construction of the multiple twin-spot marking patterns in the three mutants.

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.

The mimetic wing pattern of Papilio polytes butterflies is regulated by a doublesex-orchestrated gene network

The swallowtail butterfly Papilio polytes is sexually dimorphic and exhibits female-limited Batesian mimicry. This species also has two female forms, a non-mimetic form with male-like wing patterns, and a mimetic form resembling an unpalatable model, Pachliopta aristolochiae. The mimicry locus H constitutes a dimorphic Mendelian 'supergene', including a transcription factor gene doublesex (dsx). However, how the mimetic-type dsx (dsx-H) orchestrates the downstream gene network and causes the mimetic traits remains unclear. Here we performed RNA-seq-based gene screening and found that Wnt1 and Wnt6 are up-regulated by dsx-H during the early pupal stage and are involved in the red/white pigmentation and patterning of mimetic female wings. In contrast, a homeobox gene abdominal-A is repressed by dsx-H and involved in the non-mimetic colouration pattern. These findings suggest that dual regulation by dsx-H, induction of mimetic gene networks and repression of non-mimetic gene networks, is essential for the switch from non-mimetic to mimetic pattern in mimetic female wings.

Prepatterning of Papilio xuthus caterpillar camouflage is controlled by three homeobox genes: clawless, abdominal-A, and Abdominal-B

Color patterns often function as camouflage to protect insects from predators. In most swallowtail butterflies, younger larvae mimic bird droppings but change their pattern to mimic their host plants during their final molt. This pattern change is determined during the early fourth instar by juvenile hormone (JH-sensitive period), but it remains unclear how the prepatterning process is controlled. Using Papilio xuthus larvae, we performed transcriptome comparisons to identify three camouflage pattern-associated homeobox genes [clawless, abdominal-A, and Abdominal-B (Abd-B)] that are up-regulated during the JH-sensitive period in a region-specific manner. Electroporation-mediated knockdown of each gene at the third instar caused loss or change of original fifth instar patterns, but not the fourth instar mimetic pattern, and knockdown of Abd-B after the JH-sensitive period had no effect on fifth instar patterns. These results indicate the role of these genes during the JH-sensitive period and in the control of the prepatterning gene network.

Genome-Wide Identification and Expression Profiling of Wnt Family Genes in the Silkworm, Bombyx mori

Wnt is a family of conserved glycoproteins that participate in a variety of important biological processes including embryo development, cell proliferation and differentiation, and tissue regeneration. The Wnt family is a metazoan novelty found in all animal phyla. Studies have revealed that the number of Wnt genes varies among species, presumably due to reproduction and loss of genes during evolution. However, a comprehensive inventory of Wnt genes in Lepidoptera is lacking. In this study, we identified the repertoire of Wnt genes in the silkworm and seven other species of Lepidoptera and obtained eight Wnt genes (Wnt1, Wnt5–Wnt7, Wnt9–Wnt11, and WntA) in each species. Four of these Wnt genes are clustered in two orientations (5′-Wnt9-Wnt1-Wnt6-Wnt10-3′ and 5′-Wnt10-Wnt6-Wnt1-Wnt9-3′) in both moths and butterflies. Transcript analysis of Wnt in silkworm embryonic stages showed that each BmWnt gene had a unique expression pattern during embryological development. Analysis of a larval stage revealed differential expression of Wnt family members in diverse tissues. Our study provides an overview of the Wnt family in Lepidoptera and will inspire further functional study of the Wnt genes in the silkworm.

Proteomic Analysis of Larval Integument in a Dominant Obese Translucent (Obs) Silkworm Mutant

The dominant obese translucent (Obs) mutant of the silkworm (Bombyx mori) results in a short and stout larval body, translucent phenotype, and abnormal pigmentation in the integument. The Obs mutant also displays deficiency in ecdysis and metamorphosis. In the present study, to gain an understanding of multiple Obs phenotypes, we investigated the phenotypes of Obs and performed a comparative analysis of the larval integument proteomes of Obs and normal silkworms. The phenotypic analysis revealed that the Obs larvae were indeed short and fat, and that chitin and uric acid content were lower but melanin content was higher in the Obs mutant. Proteomic analysis revealed that 244 proteins were significantly differentially expressed between Obs and normal silkworms, some of which were involved in uric acid metabolism and melanin pigmentation. Twenty-six proteins were annotated as cuticular proteins, including RR motif-rich cuticular proteins (CPR), glycine-rich cuticular protein (CPG), hypothetical cuticular protein (CPH), cuticular protein analogous to peritrophins (CPAPs), and the chitin_bind_3 motif proteins, and accounted for over 84% of the abundance of the total significantly differentially expressed proteins. Moreover, 22 of the 26 cuticular proteins were downregulated in the Obs mutant. Comparative proteomic analysis suggested that the multiple phenotypes of the Obs mutant might be related to changes in the expression of proteins that participate in cuticular formation, uric acid metabolism, and melanin pigmentation. These results could lay a basis for further identification of the gene responsible for the Obs mutant. The data have been deposited to ProteomeXchange with identifier PXD010998.

Macroevolutionary shifts of WntA function potentiate butterfly wing-pattern diversity

Butterfly wing patterns provide a rich comparative framework to study how morphological complexity develops and evolves. Here we used CRISPR/Cas9 somatic mutagenesis to test a patterning role for WntA, a signaling ligand gene previously identified as a hotspot of shape-tuning alleles involved in wing mimicry. We show that WntA loss-of-function causes multiple modifications of pattern elements in seven nymphalid butterfly species. In three butterflies with a conserved wing-pattern arrangement, WntA is necessary for the induction of stripe-like patterns known as symmetry systems and acquired a novel eyespot activator role specific to Vanessa forewings. In two Heliconius species, WntA specifies the boundaries between melanic fields and the light-color patterns that they contour. In the passionvine butterfly Agraulis, WntA removal shows opposite effects on adjacent pattern elements, revealing a dual role across the wing field. Finally, WntA acquired a divergent role in the patterning of interveinous patterns in the monarch, a basal nymphalid butterfly that lacks stripe-like symmetry systems. These results identify WntA as an instructive signal for the prepatterning of a biological system of exuberant diversity and illustrate how shifts in the deployment and effects of a single developmental gene underlie morphological change.

Body Shape and Coloration of Silkworm Larvae Are Influenced by a Novel Cuticular Protein

The genetic basis of body shape and coloration patterns on caterpillars is often assumed to be regulated separately, but it is possible that common molecules affect both types of trait simultaneously. Here we examine the genetic basis of a spontaneous cuticle defect in silkworm, where larvae exhibit a bamboo-like body shape and decreased pigmentation. We performed linkage mapping and mutation screening to determine the gene product that affects body shape and coloration simultaneously. In these mutant larvae we identified a null mutation in BmorCPH24, a gene encoding a cuticular protein with low complexity sequence. Spatiotemporal expression analyses showed that BmorCPH24 is expressed in the larval epidermis postecdysis. RNAi-mediated knockdown and CRISPR/Cas9-mediated knockout of BmorCPH24 produced the abnormal body shape and the inhibited pigment typical of the mutant phenotype. In addition, our results showed that BmorCPH24 may be involved in the synthesis of endocuticle and its disruption-induced apoptosis of epidermal cells that accompanied the reduced expression of R&R-type larval cuticle protein genes and pigmentation gene Wnt1 Strikingly, BmorCPH24, a fast-evolving gene, has evolved a new function responsible for the assembly of silkworm larval cuticle and has evolved to be an indispensable factor maintaining the larval body shape and its coloration pattern. This is the first study to identify a molecule whose pleiotropic function affects the development of body shape and color patterns in insect larvae.

Progress and Prospects of CRISPR/Cas Systems in Insects and Other Arthropods

Clustered regularly interspaced short palindromic repeats (CRISPR) and the CRISPR-associated gene Cas9 represent an invaluable system for the precise editing of genes in diverse species. The CRISPR/Cas9 system is an adaptive mechanism that enables bacteria and archaeal species to resist invading viruses and phages or plasmids. Compared with zinc finger nucleases and transcription activator-like effector nucleases, the CRISPR/Cas9 system has the advantage of requiring less time and effort. This efficient technology has been used in many species, including diverse arthropods that are relevant to agriculture, forestry, fisheries, and public health; however, there is no review that systematically summarizes its successful application in the editing of both insect and non-insect arthropod genomes. Thus, this paper seeks to provide a comprehensive and impartial overview of the progress of the CRISPR/Cas9 system in different arthropods, reviewing not only fundamental studies related to gene function exploration and experimental optimization but also applied studies in areas such as insect modification and pest control. In addition, we also describe the latest research advances regarding two novel CRISPR/Cas systems (CRISPR/Cpf1 and CRISPR/C2c2) and discuss their future prospects for becoming crucial technologies in arthropods.

Toll ligand Spätzle3 controls melanization in the stripe pattern formation in caterpillars

A stripe pattern is an aposematic or camouflage coloration often observed among various caterpillars. However, how this ecologically important pattern is formed is largely unknown. The silkworm dominant mutant Zebra (Ze) has a black stripe in the anterior margin of each dorsal segment. Here, fine linkage mapping of 3,135 larvae revealed a 63-kbp region responsible for the Ze locus, which contained three candidate genes, including the Toll ligand gene spätzle3 (spz-3). Both electroporation-mediated ectopic expression and RNAi analyses showed that, among candidate genes, only processed spz-3 induced melanin pigmentation and that Toll-8 was the candidate receptor gene of spz-3 This Toll ligand/receptor set is also involved in melanization of other mutant Striped (p^S ), which has broader stripes. Additional knockdown of 5 other spz family and 10 Toll-related genes caused no drastic change in the pigmentation of either mutant, suggesting that only spz-3/Toll-8 is mainly involved in the melanization process rather than pattern formation. The downstream pigmentation gene yellow was specifically up-regulated in the striped region of the Ze mutant, but spz-3 showed no such region-specific expression. Toll signaling pathways are known to be involved in innate immunity, dorsoventral axis formation, and neurotrophic functions. This study provides direct evidence that a Toll signaling pathway is co-opted to control the melanization process and adaptive striped pattern formation in caterpillars.

Wingless is a positive regulator of eyespot color patterns in Bicyclus anynana butterflies

Eyespot patterns of nymphalid butterflies are an example of a novel trait yet, the developmental origin of eyespots is still not well understood. Several genes have been associated with eyespot development but few have been tested for function. One of these genes is the signaling ligand, wingless, which is expressed in the eyespot centers during early pupation and may function in eyespot signaling and color ring differentiation. Here we tested the function of wingless in wing and eyespot development by down-regulating it in transgenic Bicyclus anynana butterflies via RNAi driven by an inducible heat-shock promoter. Heat-shocks applied during larval and early pupal development led to significant decreases in wingless mRNA levels and to decreases in eyespot size and wing size in adult butterflies. We conclude that wingless is a positive regulator of eyespot and wing development in B. anynana butterflies.

Functional analysis of genes involved in color pattern formation in Lepidoptera

In addition to the genome editing technology, novel functional analyses using electroporation are powerful tools to reveal the gene function in the color pattern formation. Using these methods, several genes involved in various larval color pattern formation are clarified in the silkworm Bombyx mori and some Papilio species. Furthermore, the coloration pattern mechanism underlying the longtime mystery of female-limited Batesian mimicry of Papilio polytes has been recently revealed. This review presents the recent progress on the molecular mechanisms and evolutionary process of coloration patterns contributing to various mimicry in Lepidoptera, especially focusing on the gene function in the silkworm and Papilio species.

Identification of the silkworm quail gene reveals a crucial role of a receptor guanylyl cyclase in larval pigmentation

Diverse color patterns on the integument of lepidopteran larvae play important roles in their survival through camouflage, mimicry, sexual signaling, and aposematism. In the silkworm Bombyx mori, many color pattern variations have been preserved in inbred strains making them a good model for elucidating the molecular mechanisms that underlie color pattern formation. In this study, we focused on the silkworm quail (q) mutant, which exhibits abnormalities in multiple pigment biosynthesis pathways. Positional cloning of the q gene revealed that disruption of a guanylyl cyclase gene, BmGC-I, is responsible for its abnormal pigmentation. In q mutants, we identified a 16-bp deletion in the BmGC-I transcript, resulting in the production of a premature stop codon. Knockout of the BmGC-I gene resulted in the q-like abnormal pigmentation, thereby demonstrating that the BmGC-I gene is involved in the pigment biosynthesis pathway in the integument. Moreover, quantitative reverse transcription polymerase chain reaction showed that BmGC-I was strongly expressed in the fourth instar on day 2. Our results suggest that BmGC-I deficiency affects the pigment biosynthesis pathway, which supports the involvement of guanylyl cyclase in larval coloration.

Functional analysis of Bombyx Wnt1 during embryogenesis using the CRISPR/Cas9 system

Recently established, custom-designed nuclease technologies such as the clustered regularly interspaced short palindromic repeat (CRISPR)-associated system provide attractive genome editing tools. Targeted gene mutagenesis using the CRISPR/Cas9 system has been achieved in several orders of insects. However, outside of studies on Drosophila melanogaster and the lepidopteron model insect Bombyx mori, little success has been reported, which is largely due to a lack of effective genetic manipulation tools that can be used in other insect orders. To create a simple and effective method of gene knockout analysis, especially for dissecting gene functioning during insect embryogenesis, we performed a functional analysis of the Bombyx Wnt1 (BmWnt1) gene using Cas9/sgRNA-mediated gene mutagenesis. The Wnt1 gene is required for embryonic patterning in various organisms, and its crucial roles during embryogenesis have been demonstrated in several insect orders. Direct injection of Cas9 mRNA and BmWnt1-specific sgRNA into Bombyx embryos induced a typical Wnt-deficient phenotype: injected embryos could not hatch and exhibited severe defects in body segmentation and pigmentation in a dose-dependent manner. Quantitative real-time PCR (qRT-PCR) analysis revealed that Hox genes were down-regulated after BmWnt1 depletion. Furthermore, large deletion, up to 18Kb, ware generated. The current study demonstrates that using the CRISPR/Cas9 system is a promising approach to achieve targeted gene mutagenesis during insect embryogenesis.

Knockout silkworms reveal a dispensable role for juvenile hormones in holometabolous life cycle

Insect juvenile hormones (JHs) prevent precocious metamorphosis and allow larvae to undergo multiple rounds of status quo molts. However, the roles of JHs during the embryonic and very early larval stages have not been fully understood. We generated and characterized knockout silkworms (Bombyx mori) with null mutations in JH biosynthesis or JH receptor genes using genome-editing tools. We found that embryonic growth and morphogenesis are largely independent of JHs in Bombyx and that, even in the absence of JHs or JH signaling, pupal characters are not formed in first- or second-instar larvae, and precocious metamorphosis is induced after the second instar at the earliest. We also show by mosaic analysis that a pupal specifier gene broad, which is dramatically up-regulated in the late stage of the last larval instar, is essential for pupal commitment in the epidermis. Importantly, the mRNA expression level of broad, which is thought to be repressed by JHs, remained at very low basal levels during the early larval instars of JH-deficient or JH signaling-deficient knockouts. Therefore, our study suggests that the long-accepted paradigm that JHs maintain the juvenile status throughout larval life should be revised because the larval status can be maintained by a JH-independent mechanism in very early larval instars. We propose that the lack of competence for metamorphosis during the early larval stages may result from the absence of an unidentified broad-inducing factor, i.e., a competence factor.

Protruding structures on caterpillars are controlled by ectopic Wnt1 expression

Spine-like or protruding structures, which may be aposematic for predators, are often observed in multiple segments of lepidopteran larvae (caterpillars). For example, the larvae of the Chinese wheel butterfly, Byasa alcinous, display many protrusions on their backs as a warning that they are toxic. Although these protrusions are formed by an integument lined with single-layered epidermal cells, the molecular mechanisms underlying their formation have remained unclear. In this study, we focused on a spontaneous mutant of the silkworm, Bombyx mori, Knobbed, which shows similar protrusions to B. alcinous and demonstrates that Wnt1 plays a crucial role in the formation of protrusion structures. Using both transgene expression and RNAi-based knockdown approaches, we showed that Wnt1 designates the position where epidermal cells excessively proliferate, leading to the generation of knobbed structures. Furthermore, in the B. alcinous larvae, Wnt1 was also specifically expressed in association with the protrusions. Our results suggest that Wnt1 plays a role in the formation of protrusions on the larval body, and is conserved broadly among diverse species in Lepidoptera.

CYP18A1 regulates tissue-specific steroid hormone inactivation in Bombyx mori

Insect development and metamorphosis are regulated by two major hormones, juvenile hormone and ecdysteroids. Despite being the key regulator of insect developmental transitions, the metabolic pathway of the primary steroid hormone, 20-hydroxyecdysone (20E), especially its inactivation pathway, is still not completely elucidated. A cytochrome P450 enzyme, CYP18A1, has been shown to play key roles in insect steroid hormone inactivation through 26-hydroxylation. Here, we identified two CYP18 (BmCYP18A1 and BmCYP18B1) orthologs in the lepidopteran model insect, Bombyx mori. Interestingly, BmCYP18A1 gene is predominantly expressed in the middle silk gland (MSG) while BmCYP18B1 expresses ubiquitously in B. mori. BmCYP18A1 is induced by 20E in vitro, suggesting its role in 20E metabolism. Using the binary Gal4/UAS transgenic system, we ectopically overexpressed BmCYP18A1 in a MSG-specific manner with a Sericin1-Gal4 (Ser-Gal4) driver or in a ubiquitous manner with an Actin3-Gal4 (A3-Gal4) driver. Ectopic overexpression of BmCYP18A1 in MSG or in all tissues resulted in developmental arrestment of transgenic animals during the final instar larval stage. The 20E titers in the transgenic animals expressing BmCYP18A1 were lower compared to the levels in the control animals. Although the biological significance of MSG-specific expression of BmCYP18A1 is unclear, our results provide the first evidence that BmCYP18A1, which is conserved in most arthropods, is involved in a tissue-specific steroid hormone inactivation in B. mori.

Ancient homology underlies adaptive mimetic diversity across butterflies

Convergent evolution provides a rare, natural experiment with which to test the predictability of adaptation at the molecular level. Little is known about the molecular basis of convergence over macro-evolutionary timescales. Here we use a combination of positional cloning, population genomic resequencing, association mapping and developmental data to demonstrate that positionally orthologous nucleotide variants in the upstream region of the same gene, WntA, are responsible for parallel mimetic variation in two butterfly lineages that diverged >65 million years ago. Furthermore, characterization of spatial patterns of WntA expression during development suggests that alternative regulatory mechanisms underlie wing pattern variation in each system. Taken together, our results reveal a strikingly predictable molecular basis for phenotypic convergence over deep evolutionary time.

Little is known about the genetic basis of convergent evolution in deeply diverged species. Here, the authors show that variation in the WntA gene is associated with parallel wing pattern variation in two butterflies that diverged more than 65 million years ago.