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

eLife assessment

This important study identifies the gene mamo as a new regulator of pigmentation in the silkworm Bombyx mori, a function that was previously unsuspected based on extensive work on Drosophila where the mamo gene is involved in gamete production. The evidence supporting the role of Bm-nano in pigmentation is convincing, including high-resolution linkage mapping of two mutant strains, expression profiling, and reproduction of the mutant phenotypes with state-of-the-art RNAi and CRISPR knock-out assays. The work will be of interest to evolutionary biologists and geneticists studying color patterns and evolution of gene networks.

https://doi.org/10.7554/eLife.90795.4.sa0

Significance of the findings:

Important: Findings that have theoretical or practical implications beyond a single subfield

Landmark
Fundamental
Important
Valuable
Useful

Strength of evidence:

Convincing: Appropriate and validated methodology in line with current state-of-the-art

Exceptional
Compelling
Convincing
Solid
Incomplete
Inadequate

During the peer-review process the editor and reviewers write an eLife Assessment that summarises the significance of the findings reported in the article (on a scale ranging from landmark to useful) and the strength of the evidence (on a scale ranging from exceptional to inadequate). Learn more about eLife Assessments

Abstract
Introduction
Results
Discussion
Materials and methods
Data availability
References
Article and author information
Metrics

Abstract

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.

Introduction

Insects often display stunning colors, and the patterns of these colors have been shown to be involved in behavior (Sword et al., 2000), immunity (Barnes and Siva-Jothy, 2000), thermoregulation (Hadley et al., 1992), and UV protection (Hu et al., 2013); in particular, visual antagonism of predators via aposematism (Stevens and Ruxton, 2012), mimicry (Dasmahapatra, 2012), cryptic color patterns (Gaitonde et al., 2018), or some combination of the above (Tullberg et al., 2005). In addition, the color patterns are divergent and convergent among populations (Wittkopp et al., 2003b). Due to these striking visual features and highly active adaptive evolutionary phenotypes, the genetic basis and evolutionary mechanism of color patterns have long been a topic of interest.

Insect coloration can be pigmentary, structural, or bioluminescent. Pigments are synthesized by insects themselves and form solid particles that are deposited within the cuticle of the body surface and the scales of the wings (Futahashi et al., 2010; Matsuoka and Monteiro, 2018). Interestingly, recent studies have shown that bile pigments and carotenoid pigments synthesized through biological synthesis are incorporated into the body fluids and fill in the wing membranes of two butterflies (Siproeta stelenes and Philaethria diatonica) via hemolymph circulation, providing color in the form of liquid pigments (Finet et al., 2023). These pigments form colors by selective absorption and/or scattering of light depending on their physical properties (Gürses et al., 2016). However, structural color refers to colors, such as metallic colors and iridescence, generated by optical interference and grating diffraction of the microstructure/nanostructure of the body surface or appendages (such as scales) (Lloyd and Nadeau, 2021; Seago et al., 2009). Pigment color and structural color are widely distributed in insects and can only be observed by the naked eye in illuminated environments. However, some insects, such as fireflies, exhibit colors (green to orange) in the dark due to bioluminescence (Oba et al., 2020). Bioluminescence occurs when luciferase catalyzes the oxidation of small molecules of luciferin (Syed and Anderson, 2021). In conclusion, the color patterns of insects have evolved to become highly sophisticated and are closely related to their living environment. For example, cryptic color can deceive animals via high similarity to the surrounding environment. However, the molecular mechanism by which insects form precise color patterns to match their living environment has not been determined.

Recent research has identified the metabolic pathways associated with related pigments, such as melanins, pterins, and ommochromes, in Lepidoptera (Tong et al., 2021). A deficiency of enzymes in the pigment metabolism pathway can lead to changes in color patterns (Dai et al., 2010; Wittkopp et al., 2002). In addition to pigment synthesis, the microstructure/nanostructure of the body surface and wing scales are important factors that influence body color patterns. The body surface and wing scales of Lepidoptera are composed mainly of cuticular proteins (CPs) (Liu et al., 2021; Yan et al., 2022). There are multiple genes encoding CPs in Lepidopteran genomes. For example, in Bombyx mori, more than 220 genes encode CPs (Futahashi et al., 2008a). However, the functions of CPs and the molecular mechanisms underlying their fine-scale localization in the cuticle are still unclear.

In addition, several pleiotropic factors, such as wnt1 (Yamaguchi et al., 2013), Apontic-like (Yoda et al., 2014), clawless (Jin et al., 2019), abdominal-A (Rogers et al., 2014), abdominal-B (Jeong et al., 2006), engrailed (Dufour et al., 2020), antennapedia (Prakash et al., 2022), optix (Reed et al., 2011), bric à brac (bab) (Kopp et al., 2000), and Distal-less (Dll) (Arnoult et al., 2013), play important roles in the regulation of color patterns.

The molecular mechanism by which these factors participate in the regulation of body color patterns needs further study. In addition, there are many undiscovered factors in the gene regulatory network involved in the regulation of insect color patterns. The identification of new color pattern regulatory genes and the study of their molecular mechanisms would be helpful for further understanding color patterns.

Silkworms (B. mori) have been completely domesticated for more than 5000 years and are famous for their silk fiber production (Tong et al., 2022). Because of this long history of domestication, artificial selection, and genetic research, its genetic background has been well studied, and approximately 200 strains of mutants with variable color patterns have been preserved (Cheng. et al., 2003), which provides a good resource for color pattern research. The black dilute (bd) mutant, which exhibits recessive Mendelian inheritance, has a dark gray larval body color, and the female is sterile. Black dilute fertile (bd^f) is an allele that leads to a lighter body color than bd and fertile females (Figure 1). The two mutations were mapped to a single locus at 22.9 cM of linkage group 9 (Kawaguchi et al., 2007). No pigmentation-related genes have been reported at this locus. Thus, this research may reveal a new color pattern-related gene, which stimulated our interest.

Figure 1

Download asset Open asset

Phenotypes of bd, *bd^f*, and wild-type Dazao larvae.

The epidermis of bd is dark gray, and the epidermis of *bd^f* is light gray at the fifth instar (Day 3) of silkworm larvae. The scale bar indicates 1 cm.

Results

Candidate gene of the bd allele

To identify the genomic region responsible for the bd alleles, positional cloning of the bd locus was performed. Due to the female infertility of the bd mutant and the fertility of females with the bd^f allele, we used bd^f and the wild-type Dazao strain as parents for mapping analysis.

The 1162 back-crossed filial 1st (BC1M) generation individuals from bd^f and Dazao were subjected to fine mapping with molecular markers (Figure 2). A genomic region of approximately 390 kb was responsible for the bd phenotype (Supplementary file 1–Table S1). According to the SilkDB database (Duan et al., 2010), this region included five predicted genes (BGIBMGA012516, BGIBMGA012517, BGIBMGA012518, BGIBMGA012519, and BGIBMGA014089). In addition, we analyzed the predictive genes for this genomic region from the GenBank (Sayers et al., 2022) and SilkBase (Kawamoto et al., 2022) databases. The number of predicted genes varied among the different databases. We performed sequence alignment analysis of the predicted genes in the three databases to determine their correspondence. Real-time quantitative polymerase chain reaction (qPCR) was subsequently performed, which revealed that BGIBMGA012517 and BGIBMGA012518 were significantly downregulated on Day 3 of the fifth instar of the larvae in the bd phenotype individuals, while there was no difference in the expression levels of the other genes (Figure 2—figure supplement 1). These two genes were predicted to be associated with a single locus (LOC101738295) in GenBank. To determine the gene structures of BGIBMGA012517 and BGIBMGA012518, we used forward primers for the BGIBMGA012517 gene and reverse primers for the BGIBMGA012518 gene to amplify cDNA from the wild-type Dazao strain. By gene cloning, the predicted genes BGIBMGA012517 and BGIBMGA012518 were proven to be one such gene. For convenience, we temporarily called this the 12517/8 gene.

Figure 2 with 5 supplements see all

Download asset Open asset

Positional cloning of the bd locus.

(A) We used 532 BC1 individuals to map the bd locus between PCR markers c and i. The numbers above the DNA markers indicate recombination events. (B) A total of 1162 BC1 individuals were used to narrow the bd locus to an approximately 400 kb genomic region. (C) Partial enlarged view of the region responsible for bd. This region contains four predicted genes, *BGIBMGA012517*, *BGIBMGA012518*, *BGIBMGA012519*, and *BGIBMGA014089*. (D) Analysis of nucleotide differences in the region responsible for bd. The green block indicates the deletion of the genome in bd mutants. The black vertical lines indicate the single-nucleotide polymorphisms (SNPs) and indels of *bd^f* mutants.

The 12517/8 gene produces two transcripts; the open reading frame (ORF) of the long transcript is 2397 bp, and the ORF of the short transcript is 1824 bp, which is the same 5′-terminus, in the wild-type Dazao strain (Figure 2—figure supplement 2). The 12517/8 gene showed significantly lower expression in the bd^f mutant (Figure 2—figure supplement 3), and comparative genomic analysis revealed multiple variations in the region near this gene between bd^f and Dazao (Supplementary file 2–Table S2). In addition, 12517/8 was completely silenced due to the deletion of DNA fragments (approximately 168 kb) from the first upstream intron and an insertion of 3560 bp in the bd mutant (Figure 2—figure supplement 4).

To predict the function of the 12517/8 gene, we performed a BLAST search using its full-length sequence and found a transcription factor, the maternal gene required for meiosis (mamo) in Drosophila melanogaster, that had high sequence similarity to that of the 12517/8 gene (Figure 2—figure supplement 5). Therefore, we named the 12517/18 gene Bm-mamo; the long transcript was designated Bm-mamo-L, and the short transcript was designated Bm-mamo-S.

The mamo gene belongs to the Broad-complex, Tramtrack and bric à brac/poxvirus zinc finger protein (BTB-ZF) family. In the BTB-ZF family, the zinc finger serves as a recognition motif for DNA-specific sequences, and the BTB domain promotes oligomerization and the recruitment of other regulatory factors (Bardwell and Treisman, 1994). Most of these factors are transcriptional repressors, such as nuclear receptor corepressor (N-CoR) and silencing mediator for retinoid and thyroid hormone receptor (SMRT) (Huynh and Bardwell, 1998), but some are activators, such as p300 (Staller et al., 2001). Therefore, these features commonly serve as regulators of gene expression. mamo is enriched in embryonic primordial germ cells in D. melanogaster. Individuals deficient in mamo are able to undergo oogenesis but fail to execute meiosis properly, leading to female infertility in D. melanogaster (Mukai et al., 2007). Bm-mamo was identified as an important candidate gene for further analysis.

Expression pattern analysis of Bm-mamo

To analyze the expression profiles of Bm-mamo, we performed quantitative PCR. The expression levels of the Bm-mamo gene in the Dazao strain were investigated throughout the body at different developmental stages, from the embryonic stage to the adult stage. The gene was highly expressed in the molting stage of caterpillars, and its expression was upregulated in the later pupal and adult stages (Figure 3A). This finding suggested that the Bm-mamo gene responds to ecdysone and participates in the processes of molting and metamorphosis in silkworms. According to the investigation of tissue-specific Dazao strain expression in 5th-instar 3rd-day larvae, the midgut, head, and epidermis exhibited high expression levels; the trachea, nerves, silk glands, testis, ovary, muscle, wing disc, and Malpighian tubules exhibited moderate expression levels; and the blood and fat bodies exhibited low expression levels (Figure 3B). These findings suggested that Bm-mamo is involved in the regulation of the development of multiple silkworm tissues. Due to the melanism of the epidermis of the bd mutant and the high expression level of the Bm-mamo gene in the epidermis, we measured the expression level of this gene in the epidermis of the 4th to 5th instars of the Dazao strain. In the epidermis, the Bm-mamo gene was upregulated during the molting period, and the highest expression was observed at the beginning of molting (Figure 3C).

Figure 3

Download asset Open asset

Spatiotemporal expression of *Bm-mamo*.

(A) Temporal expression of *Bm-mamo*. In the molting and adult stages, this gene was significantly upregulated. M: molting stage, W: wandering stage, P: pupal stage, P1: Day 1 of the pupal stage, A: adult stage, 4th3d indicates the 4th instar on Day 3. 1st to 5th denote the first instar of larvae to fifth instar of larvae, respectively. (B) Tissue-specific expression of genes in 4th-instar molting larvae. *Bm-mamo* had relatively high expression levels in the midgut, head, and epidermis. AMSG: anterior division of silk gland and middle silk gland, PSG: posterior silk gland. (C) Detailed analysis of *Bm-mamo* at the 4th larval stage in the epidermis of the Dazao strain. *Bm-mamo* expression is upregulated during the molting stage. HCS indicates the head capsule stage. The ‘h’ indicates the hour, 90 h: at 90 h of the 4th instar. The mean ± SD, n=3.

Functional analyses of Bm-mamo

To study the function of the Bm-mamo gene, we carried out an RNA interference (RNAi) experiment. Short interfering RNA (siRNA) was injected into the hemolymph of silkworms, and an electroporation experiment (Ando and Fujiwara, 2013) was immediately conducted. We found significant melanin pigmentation in the epidermis of the newly molted 5th-instar larvae. These findings indicate that Bm-mamo deficiency can cause melanin pigmentation (Figure 4). The melanistic phenotype of the RNAi individuals was similar to that of the bd mutants.

Figure 4

Download asset Open asset

RNA interference (RNAi) of *Bm-mamo*.

(A) Short interfering RNA (siRNA) was introduced by microinjection followed by electroporation. ‘+’ and ‘−’ indicate the positive and negative poles of the electrical current, respectively. Scale bars, 1 cm. (B) Partial magnification of the siRNA experimental group and negative control group. Scale bars, 0.2 cm. (C, D) Relative expression levels of *Bm-mamo* in the negative control and RNAi groups were determined by quantitative polymerase chain reaction (qPCR) analysis. The mean ± SD. Numbers of samples are shown in the upper right in each graph. **p < 0.01, paired Student’s t test (ns, not significant). (E) Statistical analysis of the efficiency of the RNAi.

In addition, gene knockout was performed. Ribonucleoproteins (RNPs) generated from guide RNA (gRNA) and recombinant Cas9 protein were injected into 450 silkworm embryos. In the G0 generation, individuals with a mosaic melanization phenotype were found. These melanistic individuals were raised to moths and subsequently crossed. A homozygous line with a gene knockout was obtained through generations of screening. The gene-edited individuals had a significantly melanistic body color, and the female moths were sterile (Figure 5).

Figure 5

Download asset Open asset

*Bm-mamo* knockout.

(A) Larval phenotype of G3 fourth-instar larvae of *Bm-mamo*-targeted Dazao. Scale bars, 1 cm. (B) Partial magnification of *Bm-mamo* knockout and control individuals. Scale bars, 0.2 cm. (C) After 48 hr of egg laying, the pigmentation of the eggs indicates that those produced by homozygous knockout females cannot undergo normal pigmentation and development. (D) Genomic structure of *Bm-mamo*. The open reading frame (blue) and untranslated region (black) are shown. The gRNA 1 and gRNA 2 sequences are shown. (E) Sequences of the *Bm-mamo* knockout individuals. Lines 1 and 2 indicate deletions of 15 and 71 bp, respectively.

These results indicated that the Bm-mamo gene negatively regulates melanin pigmentation in caterpillars and participates in the reproductive regulation of female moths.

Downstream target gene analysis

Bm-mamo belongs to the zinc finger protein family, which specifically recognizes downstream DNA sequences according to their zinc fingers. We identified homologous genes of mamo in multiple species and conducted a phylogenetic analysis (Figure 6). Because the sequences of zinc finger motifs are key for downstream target gene recognition, we compared the zinc finger motifs of these orthologous proteins. Sequence alignment revealed that the amino acid residues of the zinc finger motif of the mamo-S protein were highly conserved among 57 species. This may also be due to incomplete prediction of alternative splicing in GenBank, resulting in fewer homologous proteins in mamo-L than in mamo-S (Figure 6—figure supplement 1). As the mamo protein contains a tandem Cys₂His₂ zinc finger (C2H2-ZF) motif, it can directly bind to DNA sequences. Previous research has suggested that the ZF-DNA-binding interface can be understood as a ‘canonical binding model’, in which each finger contacts DNA in an antiparallel manner. The binding sequence of the C2H2-ZF motif is determined by the amino acid residue sequence of its α-helical component. The first amino acid residue in the α-helical region of the C2H2-ZF domain is at position 1, and positions −1, 2, 3, and 6 are key amino acids involved in the recognition and binding of DNA. The residues at positions −1, 3, and 6 specifically interact with base 3, base 2, and base 1 of the DNA sense sequence, respectively, while the residue at position 2 interacts with the complementary DNA strand (Pabo et al., 2001; Wolfe et al., 2000). To analyze the downstream target genes of mamo, we first predicted the DNA-binding motifs of these genes using online software (http://zf.princeton.edu) based on the canonical binding model (Persikov and Singh, 2014). In addition, the DNA-binding sequence of mamo (TGCGT) in Drosophila was confirmed by electrophoretic mobility shift assay (EMSA) (Hira et al., 2013), which has a consensus sequence with the predicted binding site sequence of Bm-mamo-S (GTGCGTGGC), and the predicted sequence was longer. This indicates that the predicted results for the DNA-binding site have good reliability. Furthermore, the predicted DNA-binding sites of Bm-mamo-L and Bm-mamo-S were highly consistent with those of mamo orthologs in different species (Figure 6—figure supplement 2). This finding suggested that the protein may regulate similar target genes between species.

Figure 6 with 2 supplements see all

Download asset Open asset

Phylogenetic analysis.

The phylogenetic tree was subsequently constructed using the neighbor-joining method with the MEGA7 program (Pearson model). The confidence levels for various phylogenetic lineages were estimated by bootstrap analysis (2000 replicates). The magenta color indicates the homologous Bm-mamo protein, which has three zinc finger motifs. The blue color indicates the homologous Bm-mamo protein, which has five zinc finger motifs.

Figure 6—source data 1 The amino acid sequences of the mamo orthologs among multiple species.: https://cdn.elifesciences.org/articles/90795/elife-90795-fig6-data1-v1.xlsx
Download elife-90795-fig6-data1-v1.xlsx

C2H2-ZF transcription factors function by recognizing and binding to cis-regulatory sequences in the genome, which harbor cis-regulatory elements (CREs) (Wittkopp and Kalay, 2012). CREs are broadly classified as promoters, enhancers, silencers, or insulators (Preissl et al., 2023). CREs are often near their target genes, such as enhancers, which are typically located upstream (5′) or downstream (3′) of the gene they regulate or in introns, but approximately 12% of CREs are located far from their target gene (Kvon et al., 2014). Therefore, we first investigated the 2 kb upstream and downstream regions of the predicted genes in silkworms.

The predicted position weight matrices (PWMs) of the recognized sequences of Bm-mamo protein and the Find Individual Motif Occurrences (FIMO) software of MEME were used to perform silkworm whole-genome scanning for possible downstream target genes. The 2 kb upstream and downstream regions of the 14,623 predicted genes in silkworms were investigated. A total of 10,622 genes contained recognition sites within 2 kb of the upstream/downstream region of the Bm-mamo protein in the silkworm genome (Figure 7, Supplementary file 3–Table S3, and Supplementary file 4–Table S4).

Figure 7

Download asset Open asset

The DNA-binding site sequence logo of the Bm-mamo protein and the analysis of downstream target genes.

(A) The DNA-binding site sequence logo of the Bm-mamo-S protein. (B) The DNA-binding site sequence logo of the Bm-mamo-L protein. (C) Potential downstream target genes of the Bm-mamo protein in the silkworm genome were identified via the MEME FIMO program.

Moreover, we compared the transcriptome data of integument tissue between homozygotes and heterozygotes of the bd mutant at the 4th instar/beginning molting stage (Wu et al., 2016). In the integument tissue, 10,072 genes (~69% of the total predicted genes of silkworm) were expressed in heterozygotes, and 9853 genes (~67% of the total predicted genes) were expressed in homozygotes of the bd mutant. In addition, there were 191 genes whose expression significantly differed between homozygotes (bd/bd) and heterozygotes (+/bd) according to comparative transcriptome analysis (Supplementary file 5–Table S5; Wu et al., 2016). Protein functional annotation was performed, and 19 CP genes were significantly differentially expressed between heterozygotes and homozygotes of bd. In addition, the orthologs of these CPs were analyzed in Danaus plexippus, Papilio xuthu, and D. melanogaster (Supplementary file 6–Table S6). Furthermore, we identified 53 enzyme-encoding genes, 17 antimicrobial peptide genes, 6 transporter genes, 5 transcription factor genes, 5 cytochrome genes, and others. Among the differentially expressed genes (DEGs), CPs were significantly enriched, and previous studies have shown that CPs can affect pigmentation (Xiong et al., 2017). Therefore, we first investigated the expression of the CP genes. Among them, 18 CP genes had Bm-mamo binding sites in the upstream and downstream 2 kb genomic regions. In addition, we investigated the expression levels of the 18 CP genes in the integument from the 4th instar (Day 1) to the beginning of the 5th instar in the Bm-mamo knockout lines (Figure 8). Compared with those in heterozygous (mamo⁻/+) Bm-mamo knockout individuals, CP genes were significantly upregulated at one or several time points in homozygous (mamo⁻/mamo⁻) individuals. Interestingly, the expression of the CP gene BmorCPH24 was significantly upregulated at the feeding stage in Bm-mamo knockout homozygotes. Previous studies have shown that BmorCPH24 deficiency can lead to a marked decrease in pigmentation in silkworm larvae (Xiong et al., 2017). Therefore, the expression of some CP genes may be necessary for determining the color patterns of caterpillars.

Figure 8

Download asset Open asset

The expression levels of 18 cuticular protein genes were measured via quantitative polymerase chain reaction (qPCR).

The red line indicates the homozygous *Bm-mamo* knockout individuals, and the blue line indicates the heterozygous *Bm-mamo* knockout individuals. The mean ± SD, n=3. *p < 0.05 ( Student’s t test). 4th1d indicates 4-instar Day 1; 0 h indicates the beginning of molting, h indicates hours, and 5th0h indicates the beginning of the 5th instar.

In addition, the synthesis of pigment substances is an important driver of color patterns (Wittkopp et al., 2003a). Eight key genes (TH, DDC, aaNAT, ebony, black, tan, yellow, and laccase2) involved in melanin synthesis (Matsuoka and Monteiro, 2018) were investigated in the heterozygous and homozygous Bm-mamo gene knockout lines. Among them, yellow, tan, and DDC were significantly upregulated during the molting period in the homozygous Bm-mamo knockout individuals (Figure 9). The upstream and/or downstream 2 kb genomic regions of yellow, tan, and DDC contain the binding site of the Bm-mamo protein. In addition, the expression of yellow, DDC, and tan can promote the generation of melanin (Gibert, 2020).

Figure 9 with 3 supplements see all

Download asset Open asset

Melanin metabolism pathway and quantitative polymerase chain reaction (qPCR) of related genes.

(A) The melanin metabolism pathway. Blue indicates amino acids and catecholamines, red indicates the names of genes, and black indicates the names of enzymes. (B) Relative expression levels of eight genes in the heterozygous *Bm-mamo* knockout group (blue) and homozygous *Bm-mamo* knockout group (red) determined via qPCR analysis. The mean ± SD, n=3. *p < 0.05 (Student’s t test). 4th3d indicates the 4th instar on Day 3, M indicates molting, and h indicates hours.

To explore the interaction between the Bm-mamo protein and its binding sequence, EMSA was conducted. The binding site of Bm-mamo-S (CTGCGTGGT) was located approximately 70 bp upstream of the transcription initiation site of the Bm-yellow gene. The EMSA results showed that the Bm-mamo-S protein expressed in prokaryotes can bind to the CTGCGTGGT sequence in vitro (Figure 9—figure supplement 1). This finding suggested that the Bm-mamo-S protein can bind to the upstream region of the Bm-yellow gene and regulate its transcription. Therefore, the Bm-mamo gene may control the color pattern of caterpillars by regulating key melanin synthesis genes and 18 CP genes.

Discussion

Insects have evolved many important phenotypes during the process of adapting to the environment. Among these traits, the color pattern is one of the most interesting. The biochemical metabolic pathways of pigments, such as melanin, ommochromes, and pteridines, have been identified in insects (Tong et al., 2021). However, the regulation of pigment metabolism-related genes and the processes involved in the transport and deposition of pigment substances are unclear. In this study, we discovered that the Bm-mamo gene negatively regulates melanin pigmentation in caterpillars. When this gene is deficient, the body color of silkworms exhibits substantial melanism, and changes in the expression of some melanin synthesis genes and some CP genes are also significant.

Structural genes involved in melanin synthesis

The deficiencies of some genes in the melanin synthesis pathway, such as TH, yellow, ebony, tan, and aaNAT, can lead to variations in color patterns (Matsuoka and Monteiro, 2018). However, they often lead to localized pigmentation changes in later-instar caterpillars. For example, in yellow-mutant silkworms, the eye spot, lunar spot, star spot, spiracle plate, head, and some sclerotized areas appear reddish brown, and the other cuticle is consistent with that in the wild type in later instars of larvae (third to fifth instars) (Futahashi et al., 2008b). This situation is highly similar to that of the tan mutant in silkworms (Noh et al., 2016). Why do the yellow and tan phenotypes appear only in a limited cuticle region in later instars of silkworm larvae? One possible reason is that the expression of yellow and tan is limited to a certain region by transcription regulators. Alternatively, other factors, such as CPs in the cuticle, may limit pigmentation by interacting with pigment substances in the cuticle.

On the one hand, we investigated the expression levels of yellow and tan in the pigmented region (lunar spot) and nonpigmented region of the epidermis during the 4th molting of the wild-type Dazao strain. The expression level of yellow was significantly upregulated in the lunar spot (the epidermis on the dorsum of the fifth body segment) compared with the nonpigmented region (the epidermis on the dorsum of the sixth body segment) at 6 and 12 hr after molting. Moreover, tan was significantly upregulated 18 hr after molting in the lunar spot region (Figure 9—figure supplement 2). It is suggested that the upregulated expression of pigment synthesis genes at key time points may be important for pigmentation. However, yellow and tan were still moderately expressed in the nonpigmented epidermis, although they did not cause significant melanin pigmentation. This finding indicates that pigment synthesis alone cannot determine the predominant color pattern of the cuticle in caterpillars.

CPs participate in pigmentation

On the other hand, synthesized pigment substances need to be transported from epidermal cells and embedded into the cuticle to allow pigmentation. Therefore, the correct cuticle structure and location of CPs in the cuticle may be important factors affecting pigmentation. Previous studies have shown that a lack of BmorCPH24, which encodes an important component of the endocuticle, can lead to dramatic changes in body shape and a significant reduction in the pigmentation of caterpillars (Xiong et al., 2017). We crossed Bo (BmorCPH24 null mutation) and bd to obtain F₁ (Bo/+^Bo, bd/+) and then self-crossed F₁ to observe the phenotype of F₂. The area of lunar spots and star spots decreased, and light-colored stripes appeared on the body segments; however, the other areas still exhibited significant melanin pigmentation in double-mutation (Bo, bd) individuals (Figure 10). However, in previous studies, the introduction of Bo into L (the ectopic expression of wnt1 results in lunar stripes generated on each body segment) (Yamaguchi et al., 2013) and U (the overexpression of SoxD results in excessive melanin pigmentation of the epidermis) (Wang et al., 2022) strains by genetic crosses markedly reduced the pigmentation of L and U (Xiong et al., 2017). Interestingly, there was a more significant decrease in pigmentation in the double mutants (Bo, L) and (Bo, U) than in (Bo, bd). These findings suggested that Bm-mamo has a greater ability than does wnt1 and SoxD to regulate pigmentation. On the one hand, mamo may be a stronger regulator of the melanin metabolic pathway, and on the other hand, mamo may regulate other CP genes to reduce the impact of BmorCPH24 deficiency.

Figure 10

Download asset Open asset

Phenotypes of offspring from hybridization between Bo and bd.

(A) The phenotypes of 5th instar and Day 3. the. The bars indicate 1 cm. (B) Magnification of body segments. The bars indicate 0.5 cm. The area where the pigment is significantly reduced is indicated by the arrow. The triangle shape indicates that the number of lunar spots is reduced.

How do CPs affect pigmentation? One study showed that some CPs can form ‘pore canals’ to transport macromolecules (Kramer et al., 2001). In addition, some CPs can be crosslinked with catecholamines, which are synthesized via the melanin metabolism pathway (Noh et al., 2015). Because there are no live cells in the cuticle, melanin precursor substances may be transported by the pore canals formed by some CPs and fixed to specific positions through cross-linking with CPs. The CP TcCPR4 is needed for the formation of pore canals in the cuticle of Tribolium castaneum (Noh et al., 2014). In contrast, the vertical pore canal is lacking in the less pigmented cuticles of T. castaneum (Mun et al., 2015). This finding suggested that the pore canals constructed by TcCPR4 may transport pigments and contribute to cuticle pigmentation in T. castaneum. Moreover, the melanin metabolites N-acetyldopamine and N-β-alanyldopamine can target and sclerotize the cuticle by cross-linking with specific CPs (Noh et al., 2016). This finding suggested that pigments interact with specific CPs, thereby affecting pigmentation, hardening properties, and the structure of the cuticle. Interestingly, a study showed that in addition to absorbing specific wavelengths, pigments can affect cuticle polymerization, density, and the refractive index, which in turn affects the reflected wavelengths that produce structural color in butterfly wing scales (Prakash et al., 2023). This implies that the interaction between pigments and CPs can be very subtle, resulting in the formation of unique nanostructures, such as those on wing scales, that produce brilliant structural colors.

Consequently, to maintain the accuracy of the color pattern, the localization of CPs in the cuticle may be very important. In a previous study employing microarray analysis, different CPs were found in differently colored areas of the epidermis in P. xuthus larvae (Futahashi et al., 2012). We investigated whether the CP genes were highly expressed in the black region of P. xuthus caterpillars. Thirteen orthologous genes were found in silkworms (Supplementary file 7–Table S7). Among them, 11 genes (BmorCPR67, BmorCPR71, BmorCPR76, BmorCPR79, BmorCPR99, BmorCPR107, BmorCPT4, BmorCPH5, BmorCPFL4, BmorCPG27, and BmorCPG4) were significantly upregulated in homozygous (mamo⁻/mamo⁻) knockout individuals at some time points from the 3rd day of the 4th instar to the 5th instar; BmorCPG4 was also among the 18 previously detected CP genes, and two genes (BmorCPG38 and BmorCPR129) were not differentially expressed between homozygous (mamo⁻/mamo⁻) and heterozygous (mamo⁻/+) individuals (Figure 11).

Figure 11 with 1 supplement see all

Download asset Open asset

The expression levels of 12 cuticular protein-encoding genes were detected.

The expression levels of the orthologous genes of the cuticular protein genes highly expressed in the black marking of *Papilio xuthus* caterpillars in the silkworm caterpillars were detected. Blue indicates heterozygous *Bm-mamo* gene knockout individuals, and red indicates homozygous *Bm-mamo* gene knockout individuals. The mean ± SD, n=3. *p < 0.05 (Student’s t test). 4th3d indicates the 4th instar on Day 3, M indicates molting, and h indicates hours.

The expression of the 28 CP genes mentioned above was significantly upregulated in homozygous (mamo⁻/mamo⁻) gene knockout individuals at several stages. These CPs may be involved in the transport or cross-linking of melanin in the cuticle. However, there were no differences in the expression levels of these genes during some periods compared with those in the control group, and the expression of some genes was significantly downregulated at some time points in the melanic individuals (Figure 11—figure supplement 1). This finding suggested that the regulation of CP genes is complex and may involve other transcription factors and feedback effects. CPs are essential components of the insect cuticle and are involved in cuticular microstructure construction (Noh et al., 2017), body shape development (Tajiri et al., 2017), wing morphogenesis (Zhao et al., 2019), and pigmentation (Xiong et al., 2017). CP genes usually account for more than 1% of the total genes in an insect genome and can be categorized into several families, including CPR, CPG, CPH, CPAP1, CPAP3, CPT, CPF, and CPFL (Willis, 2010). The CPR family is the largest group of CPs and contains a chitin-binding domain called the Rebers and Riddiford motif (R&R) (Rebers and Riddiford, 1988). The R&R consensus sequences can be divided into three subfamilies (RR-1, RR-2, and RR-3) (Karouzou et al., 2007). Among the 28 CPs, 11 RR-1 genes, 6 RR-2 genes, 4 hypothetical cuticular protein (CPH) genes, 3 glycine-rich cuticular protein (CPG) genes, 3 cuticular protein Tweedle motif (CPT) genes, and 1 CPFL (like the CPFs in a conserved C-terminal region) gene were identified. The RR-1 consensus among species is usually more variable than the RR-2 consensus, which suggests that RR-1 may have a species-specific function. RR-2 often clustered into several branches, which may be due to gene duplication events in co-orthologous groups and may result in conserved functions between species (Chen and Hou, 2021). The classification of CPH was based on the lack of known motifs. In the epidermis of Lepidoptera, CPH genes often exhibit high expression levels. For example, BmorCPH24 has the highest expression level in the silkworm larval epidermis (Fu et al., 2023). The CPG protein is rich in glycine. The CPH and CPG genes are less commonly found in insects outside the order Lepidoptera (Yang et al., 2017). This finding suggested that these genes may provide species-specific functions for Lepidoptera. CPT contains a Tweedle motif, and the TweedleD1 mutation has a dramatic effect on body shape in D. melanogaster (Guan et al., 2006). The CPFL family members are relatively conserved among species and may be involved in the synthesis of larval cuticles (Togawa et al., 2007). CPT and CPFL may have relatively conserved functions among insects. The CP genes are a group of rapidly evolving genes, and their copy numbers may undergo significant changes in different species. In addition, RNAi experiments on 135 CP genes in the brown planthopper (Nilaparvata lugens) revealed that deficiency of 32 CP genes leads to significant defects in phenotypes, such as lethality and developmental retardation. These findings suggested that the 32 CP genes are indispensable and that other CP genes may have redundant and complementary functions (Pan et al., 2018). Previous studies revealed that the construction of the larval cuticle of silkworms requires the precise expression of more than two hundred CP genes (Yan et al., 2022). The production, interaction, and deposition of CPs and pigments are complex and precise processes, and our research showed that Bm-mamo plays an important regulatory role in this process in silkworm caterpillars. To further understand the role of CPs, future work should aim to identify the functions of important CP genes and the deposition mechanism in the cuticle.

Maturation of pigment granules

In addition, among the 191 DEGs found in the comparative transcriptome data, we also discovered several interesting genes. For example, BGIBMGA013242, which encodes a major facilitator superfamily protein (MFS) named BmMFS and is responsible for the cheek and tail spot (cts) mutant, was significantly upregulated in the bd mutant. A deficiency of this gene (BmMFS) results in chocolate-colored head and anal plates on the silkworm caterpillar (Ito et al., 2012). In the ommochrome metabolic pathway, Bm-re, which encodes an MFS protein, may function in the transportation of several amino acids, such as cysteine or methionine, into pigment granules (Osanai-Futahashi et al., 2012; Luo et al., 2021). Therefore, the encoded product of BmMFS may participate in pigmentation by promoting the maturation of pigment granules. Moreover, BGIBMGA013576 and BGIBMGA013656, which encode MFS domain-containing proteins belonging to the solute carrier family 22 (SLC22) and solute carrier family 2 (SLC2) families, respectively, were significantly upregulated in bd. These two genes may participate in pigmentation in a similar manner to that of the BmMFS gene. In addition, the red Malpighian tubule (red) gene encodes a LysM domain-containing protein. Red deficiency results in a significant decrease in orange pterin in the wing of Colias butterflies. Research suggests that the product of red may interact with V-ATPase to modulate vacuolar pH across a variety of endosomal organelles, thereby affecting the maturation of pigment granules in cells (Hanly et al., 2023). These findings indicate that the maturation of pigment granules plays an important role in the coloring of Lepidoptera and that Bm-mamo may be involved in regulating the maturation of pigment granules in the epidermal cells of silkworms.

Bm-mamo may affect the synthesis of melanin in epidermal cells by regulating yellow, DDC, and tan, regulating the maturation of melanin granules in epidermal cells through BmMFS, and affecting the deposition of melanin granules in the cuticle by regulating CP genes, thereby comprehensively regulating the color pattern of caterpillars.

The upstream sequence of Bm-mamo

Moreover, in D. melanogaster, mamo is needed for functional gamete production. In the silkworm, Bm-mamo has a conserved function in female reproduction. The DNA-binding motifs of the mamo proteins in silkworms and fruit flies are highly conserved, suggesting that they may regulate the same downstream target genes. However, mamo has developed new functions in color pattern regulation in silkworm caterpillars. We found that in D. melanogaster, the mamo gene is expressed mainly during the adult stage; it is not expressed during the 1st instar or 2nd instar larval stage and has very low expression at the 3rd instar (Supplementary file 8–Table S8). In addition, several binding sites of mamo were found near the TSS of yellow in D. melanogaster (Figure 9—figure supplement 3). The yellow gene is a key melanin metabolic gene with upstream and intron sequences that have been identified as multiple CREs, and it has been considered a research model for CREs (Kalay et al., 2019; Xin et al., 2020). This difference may be due to a change in the expression pattern of mamo, which led to the development of a new function for this gene in regulating coloration in silkworm caterpillars.

Changes in gene expression patterns are generally believed to be the result of the evolution of CREs. Because CREs play important roles in the spatiotemporal expression pattern regulation of genes, many CREs are found in noncoding regions, which are relatively prone to sequence variation compared with the sequences of coding genes (Wittkopp and Kalay, 2012). However, the molecular mechanism underlying the sequence changes in CREs is unclear.

Transcription factors (TFs), because they recognize relatively short sequences, generally between 4 and 20 bp in length, can have many binding sites in the genome (Jolma and Taipale, 2011). Therefore, one member of the TF family has the potential to regulate many genes in one genome. Single regulation of specific TFs can lead to patterned changes in the expression of multiple downstream genes, enabling organisms to adapt to the environment by altering a given type of trait. For example, the marine form of the three-spined stickleback (Gasterosteus aculeatus) has thick armor, whereas the lake population (which was recently derived from the marine form) does not. Research has shown that pelvic loss in different natural populations of three-spined stickleback fishes occurs via regulatory mutations resulting from the deletion of a tissue-specific enhancer (Pel) of the pituitary homeobox transcription factor 1 (Pitx1) gene. The researchers genotyped 13 pelvic-reduced populations of three-spined sticklebacks from disparate geographic locations. Nine of the 13 pelvic-reduced stickleback populations had sequence deletions of varying lengths, all of which were located at the Pel enhancer (Chan et al., 2010).

We investigated nucleotide diversity in 51 wild silkworms and 171 domesticated silkworms (Tong et al., 2022; Lu et al., 2024). The nucleotide diversity of introns and upstream sequences of Bm-mamo in wild silkworms was significantly greater than that in domestic silkworms. In addition, the approximately 1 kb genomic region upstream had high fixation indices (F_ST) (Figure 12). This indicates a significant degree of differentiation between the wild strains in this genomic region. Multiple sequence alignment analysis of this region was performed between 12 wild silkworms and 12 domestic silkworms (Figure 12—figure supplement 1 and Supplementary file 9–Table S9). The results showed that the wild silkworm group had a greater degree of nucleotide variation, while the sequences in domestic silkworms were highly conserved. The domestic silkworm group had two different forms: one had a long interspersed nuclear element (LINE) inserted at approximately 4.5 kb, and the other did not have this transposon sequence. It is suggested that this genomic region is prone to variation in wild silkworm and is relatively fixed in domestic silkworms. The larvae of domesticated strains often have a lighter or even pale body color, but the wild-type strains have a darker color. The Bm-mamo gene may be involved in the domestication of silkworms.

Figure 12 with 1 supplement see all

Download asset Open asset

Investigation of nucleotide diversity and fixation index in the *Bm-mamo* genomic region using 51 wild and 171 domestic silkworm strains.

(A) Nucleotide diversity analysis. The greater the value of pi (π) is, the greater the diversity of nucleotide sequences within the population. The dashed line represents the average nucleotide diversity level of the silkworm. The X-axis represents the position of chromosome 9. (B) Fixation index analysis. The greater the value of *F_ST* is, the greater the degree of differentiation of alleles between populations. The dashed line represents the top 5% fixation level of silkworms. The X-axis represents the position of chromosome 9.

TFs function by precisely regulating the temporal and spatial expression of target genes. Therefore, the expression of TFs requires strict regulation. There are long intergenic regions upstream of many important TFs, dozens of kilobase pairs (kb) to hundreds of kb, which may contain many CREs for better control of their expression pattern. It is often believed that changes in CREs are caused by random mutations. However, recent studies have shown that there is a mutation bias in the genome; compared with that in the intergenic region, the mutation frequency is reduced by half inside gene bodies and by two-thirds in essential genes. In addition, they compared the mutation rates of genes with different functions. The mutation rate in the coding region of essential genes (such as translation) is the lowest, and the mutation rates in the coding region of specialized functional genes (such as environmental response) are the highest. These patterns are mainly affected by the features of the epigenome (Monroe et al., 2022). Due to the plasticity of epigenomic features, mutation bias associated with epigenomes may even lead to environmental influences on mutations (Belfield et al., 2021). This finding suggested that different functional regions in the genome are subject to distinct controls and that some sequences can undergo procedural changes under environmental changes. Therefore, some variations in the CREs of TFs in response to environmental changes may be mutation bias.

Materials and methods

Silkworm strains

Request a detailed protocol

The bd and bd^f mutant strains and the wild-type Dazao and N4 strains were obtained from the bank of genetic resources of Southwest University. Silkworms were reared on mulberry leaves or artificial diets at 25°C and 73% relative humidity in the dark for 12 hr and light for 12 hr.

Positional cloning of bd

Request a detailed protocol

For mapping of the bd locus, F1 heterozygous individuals were obtained from a cross between a bd^f strain and a Dazao strain. Then, an F1 female was crossed with a bd^f male (BC1F), and an F1 male was backcrossed with a bd^f female (BC1M). A total of 1162 BC1M individuals were used for recombination analysis. Genomic DNA was extracted from the parents (Dazao, bd^f, and F1) and each BC1 individual using the phenol chloroform extraction method. Available DNA molecular markers were identified through polymorphism screening and linkage analysis. The primers used for mapping are listed in Supplementary file 10–Table S10.

Phylogenetic analysis

Request a detailed protocol

To determine whether Bm-mamo orthologs were present in other species, the BlastX program of the National Center for Biotechnology Information (NCBI) (http://www.ncbi.nlm.nih.gov/BLAST/) was used. The Bm-mamo-L and Bm-mamo-S sequences were subjected to BLAST searches against the nonredundant protein sequence (nr) database. Sequences with a maximum score and an E-value ≤10⁻⁴ were downloaded. The sequences of multiple species were subjected to multiple sequence alignment of the predicted amino acid sequences by MUSCLE (Edgar, 2004). A phylogenetic tree was subsequently constructed using the neighbor-joining method with the MEGA7 program (Pearson model). The confidence levels for various phylogenetic lineages were estimated by bootstrap analysis (2000 replicates).

Quantitative PCR

Request a detailed protocol

Total RNA was isolated from the whole body and integument of the silkworms using TRIzol reagent (Invitrogen, California, USA), purified by phenol chloroform extraction and then reverse transcribed with a PrimeScript RT Reagent Kit (TAKARA, Dalian, China) according to the manufacturer’s protocol.

qPCR was performed using a CFX96 Real-Time PCR Detection System (Bio-Rad, Hercules, CA) with a qPCR system and an iTaq Universal SYBR Green Supermix System (Bio-Rad). The cycling parameters were as follows: 95°C for 3 min, followed by 40 cycles of 95°C for 10 s and annealing for 30 s. The primers used for the target genes are listed in Supplementary file 10–Table S10. The expression levels of the genes in each sample were determined with biological replicates, and each sample was analyzed in triplicate. The gene expression levels were normalized against the expression levels of the ribosomal protein L3 (RpL3). The relative expression levels were analyzed using the classical R = 2^−ΔΔCt method.

siRNA for gene knockdown

Request a detailed protocol

siRNAs for Bm-mamo were designed with the siDirect program (http://sidirect2.rnai.jp). The target siRNAs and negative controls were synthesized by Tsingke Biotechnology Company Limited. siRNA (5 μl, 1 μl/μg) was injected from the abdominal spiracle into the hemolymph at the fourth-instar (Day 3) larval stage. Immediately after injection, phosphate-buffered saline (pH 7.3) droplets were placed nearby, and a 20-V pulse for 1 s and pause for 1 s were applied three times. The phenotype was observed for fifth-instar larvae. The left and right epidermis were separately dissected from the injected larvae, after which RNA was extracted. Then, cDNA was synthesized, and the expression level of the gene was detected via qPCR.

Single-guide RNA synthesis and RNP complex assembly

Request a detailed protocol

CRISPRdirect (http://crispr.dbcls.jp/doc/) online software was used to screen appropriate single-guide RNA (sgRNA) target sequences. The gRNAs were synthetized by the Beijing Genomics Institute. The sgRNA templates were transcribed using T7 polymerase with RiboMAX Large-Scale RNA Production Systems (Promega, Beijing, China) according to the manufacturer’s instructions. The RNA transcripts were purified using 3 M sodium acetate (pH 5.2) and anhydrous ethanol (1:30) precipitation, washed with 75% ethanol, and eluted in nuclease-free water. All injection mixtures contained 300 ng/µl Cas9 nuclease (Invitrogen, California, USA) and 300 ng/µl purified sgRNA. Before injection, mixtures of Cas9 nuclease and gRNA were incubated for 15 min at 37°C to reconstitute active RNPs (Zou et al., 2022).

Microinjection of embryos

Request a detailed protocol

For embryo microinjection, microcapillary needles were manufactured using a PC-10 model micropipette puller (Narishige, Tokyo, Japan). Microinjection was performed using an Eppendorf TransferMan NK 2 and a FemtoJet 4i system (Eppendorf, Hamburg, Germany). The eggs used for microinjection were generated by mating female and male wild-type Dazao moths. Within 4 hr of culture, the eggs were allowed to adhere to a clean glass slide. CRISPR/Cas9-messenger RNP mixtures with volumes of approximately 1 nl were injected into the middle of the eggs, and the wound was sealed with glue. All the injected embryos were allowed to develop in a climate chamber at 25°C and 80% humidity (Long et al., 2021).

Comparative genomics

Request a detailed protocol

The reference genome of Dazao was downloaded from Silkbase (https://silkbase.ab.a.u-tokyo.ac.jp/cgi-bin/index.cgi). The bd^f genome was obtained from the silkworm pangenome project (Tong et al., 2022).

The short reads of the bd^f strains were mapped to the silkworm reference genome by BWA55 v0.7.17 mem with default parameters. The SAMtools v1.11(https://github.com/samtools/samtools; Danecek et al., 2024) and Picard v2.23.5 (https://broadinstitute.github.io/picard/) programs were used to filter the unmapped and duplicated reads. A GVCF file of the samples was obtained using GATK57 v4.1.8.1 HaplotypeCaller with the parameter -ERC=GVCF. The VCF files of insertions/deletions (indels) and single-nucleotide polymorphisms were used for further analysis via eGPS software.

Downstream target gene screening

Request a detailed protocol

Online software (http://zf.princeton.edu) was used to predict the DNA-binding site for Cys2His2 zinc finger proteins. The confident ZF domains with scores higher than 17.7 were chosen. RF regression on the B1H model was used to predict the DNA-binding sites. Then, the sequence logo and PWMs of the DNA-binding sites were obtained. The sequences 2 kb upstream and downstream of the predicted genes were extracted by a Perl script. The FIMO package of the MEME suite was used to search for binding sites in the silkworm genome.

Analysis of EMSA

Request a detailed protocol

Primers with binding sites and flanking sequences were designed according to the FIMO results. A biotin label was added to one end of the upstream primer (probe), and the downstream primer was used as a reference. Primers with the same sequence as the labeled probes were used as competitive probes. The EMSA was conducted with an EMSA reagent kit (Beyotime, Shanghai, China) according to the manufacturer’s protocol.

Data availability

Sequencing data have been deposited in GenBank under accession numbers PP426589 and PP426590. In additional, all data generated or analyzed during this study are included in the manuscript and supporting files; source data files have been provided for Figures 2, 6 , 9, and 11.

The following data sets were generated

1. Wu S
(2024) NCBI GenBank
ID PP426589. Bombyx mori strain Dazao mamo-L transcript 2 (mamo) mRNA, complete cds.

https://www.ncbi.nlm.nih.gov/nuccore/PP426589.1/
1. Wu S
(2024) NCBI GenBank
ID PP426590. Bombyx mori strain Dazao mamo-S transcript 1 (mamo) mRNA, complete cds.

https://www.ncbi.nlm.nih.gov/nuccore/PP426590.1/

References

1. Ando T
2. Fujiwara H
(2013) Electroporation-mediated somatic transgenesis for rapid functional analysis in insects
Development 140:454–458.

https://doi.org/10.1242/dev.085241
- PubMed
- Google Scholar
1. Arnoult L
2. Su KFY
3. Manoel D
4. Minervino C
5. Magriña J
6. Gompel N
7. Prud’homme B
(2013) Emergence and diversification of fly pigmentation through evolution of a gene regulatory module
Science 339:1423–1426.

https://doi.org/10.1126/science.1233749
- PubMed
- Google Scholar
1. Bardwell VJ
2. Treisman R
(1994) The POZ domain: a conserved protein-protein interaction motif
Genes & Development 8:1664–1677.

https://doi.org/10.1101/gad.8.14.1664
- PubMed
- Google Scholar
1. Barnes AI
2. Siva-Jothy MT
(2000) Density-dependent prophylaxis in the mealworm beetle Tenebrio molitor L. (Coleoptera: Tenebrionidae): cuticular melanization is an indicator of investment in immunity
Proceedings. Biological Sciences 267:177–182.

https://doi.org/10.1098/rspb.2000.0984
- PubMed
- Google Scholar
1. Belfield EJ
2. Brown C
3. Ding ZJ
4. Chapman L
5. Luo M
6. Hinde E
7. van Es SW
8. Johnson S
9. Ning Y
10. Zheng SJ
11. Mithani A
12. Harberd NP
(2021) Thermal stress accelerates Arabidopsis thaliana mutation rate
Genome Research 31:40–50.

https://doi.org/10.1101/gr.259853.119
- PubMed
- Google Scholar
1. Chan YF
2. Marks ME
3. Jones FC
4. Villarreal G
5. Shapiro MD
6. Brady SD
7. Southwick AM
8. Absher DM
9. Grimwood J
10. Schmutz J
11. Myers RM
12. Petrov D
13. Jónsson B
14. Schluter D
15. Bell MA
16. Kingsley DM
(2010) Adaptive evolution of pelvic reduction in sticklebacks by recurrent deletion of a Pitx1 enhancer
Science 327:302–305.

https://doi.org/10.1126/science.1182213
- PubMed
- Google Scholar
1. Chen EH
2. Hou QL
(2021) Identification and expression analysis of cuticular protein genes in the diamondback moth, Plutella xylostella (Lepidoptera: Plutellidae)
Pesticide Biochemistry and Physiology 178:104943.

https://doi.org/10.1016/j.pestbp.2021.104943
- PubMed
- Google Scholar
(2003)
Studies on the mutant systems of the Bombyx mori gene bank

Scientia Agricultura Sinica 36:968–975.
- Google Scholar
1. Dai F
2. Qiao L
3. Tong X
4. Cao C
5. Chen P
6. Chen J
7. Lu C
8. Xiang Z
(2010) Mutations of an arylalkylamine-N-acetyltransferase, Bm-iAANAT, are responsible for silkworm melanism mutant
The Journal of Biological Chemistry 285:19553–19560.

https://doi.org/10.1074/jbc.M109.096743
- PubMed
- Google Scholar
Software
1. Danecek P
2. Bonfield JK
3. Liddle J
4. Marshall J
5. Ohan V
6. Pollard MO
7. Whitwham A
8. Keane T
9. McCarthy SA
10. Davies RM
11. Li H
(2024) Samtools, version 4db9a36
GitHub.

https://github.com/samtools/samtools
1. Dasmahapatra KK
(2012) Butterfly genome reveals promiscuous exchange of mimicry adaptations among species
Nature 487:94–98.

https://doi.org/10.1038/nature11041
- Google Scholar
1. Duan J
2. Li R
3. Cheng D
4. Fan W
5. Zha X
6. Cheng T
7. Wu Y
8. Wang J
9. Mita K
10. Xiang Z
11. Xia Q
(2010) SilkDB v2.0: a platform for silkworm (Bombyx mori ) genome biology
Nucleic Acids Research 38:D453–D456.

https://doi.org/10.1093/nar/gkp801
- PubMed
- Google Scholar
(2020) Temporal flexibility of gene regulatory network underlies a novel wing pattern in flies
PNAS 117:11589–11596.

https://doi.org/10.1073/pnas.2002092117
- Google Scholar
1. Edgar RC
(2004) Muscle: Multiple sequence alignment with high accuracy and high throughput
Nucleic Acids Research 32:1792–1797.

https://doi.org/10.1093/nar/gkh340
- Google Scholar
1. Finet C
2. Ruan Q
3. Bei YY
4. You En Chan J
5. Saranathan V
6. Yang JKW
7. Monteiro A
(2023) Multi-scale dissection of wing transparency in the clearwing butterfly Phanus vitreus
Journal of the Royal Society, Interface 20:20230135.

https://doi.org/10.1098/rsif.2023.0135
- PubMed
- Google Scholar
1. Fu X
2. Chen M
3. Xia R
4. Li X
5. Li Q
6. Li Y
7. Cao H
8. Liu Y
(2023) Genome-wide identification and transcriptome-based expression profile of cuticular protein genes in Antheraea pernyi
International Journal of Molecular Sciences 24:6991.

https://doi.org/10.3390/ijms24086991
- PubMed
- Google Scholar
1. Futahashi R
2. Okamoto S
3. Kawasaki H
4. Zhong Y-S
5. Iwanaga M
6. Mita K
7. Fujiwara H
(2008a) Genome-wide identification of cuticular protein genes in the silkworm, Bombyx mori
Insect Biochemistry and Molecular Biology 38:1138–1146.

https://doi.org/10.1016/j.ibmb.2008.05.007
- Google Scholar
1. Futahashi R
2. Sato J
3. Meng Y
4. Okamoto S
5. Daimon T
6. Yamamoto K
7. Suetsugu Y
8. Narukawa J
9. Takahashi H
10. Banno Y
11. Katsuma S
12. Shimada T
13. Mita K
14. Fujiwara H
(2008b) yellow and ebony are the responsible genes for the larval color mutants of the silkworm Bombyx mori
Genetics 180:1995–2005.

https://doi.org/10.1534/genetics.108.096388
- PubMed
- Google Scholar
(2010) Caterpillar color patterns are determined by a two-phase melanin gene prepatterning process: new evidence from tan and laccase2
Evolution & Development 12:157–167.

https://doi.org/10.1111/j.1525-142X.2010.00401.x
- PubMed
- Google Scholar
(2012) Comprehensive microarray-based analysis for stage-specific larval camouflage pattern-associated genes in the swallowtail butterfly, Papilio xuthus
BMC Biology 10:46.

https://doi.org/10.1186/1741-7007-10-46
- PubMed
- Google Scholar
(2018) Evolution of ontogenic change in color defenses of swallowtail butterflies
Ecology and Evolution 8:9751–9763.

https://doi.org/10.1002/ece3.4426
- PubMed
- Google Scholar
1. Gibert JM
(2020) Phenotypic plasticity in insects
Biologie Aujourd’hui 214:33–44.

https://doi.org/10.1051/jbio/2020005
- PubMed
- Google Scholar
(2006) Mutation of TweedleD, a member of an unconventional cuticle protein family, alters body shape in Drosophila
PNAS 103:16794–16799.

https://doi.org/10.1073/pnas.0607616103
- PubMed
- Google Scholar
Book
(2016) Dyes and Pigments
Springer International Publishing.

https://doi.org/10.1007/978-3-319-33892-7
- Google Scholar
(1992) Coloration and its thermal consequences in the New Zealand tiger beetle Neocicindela perhispida
Journal of Thermal Biology 17:55–61.

https://doi.org/10.1016/0306-4565(92)90020-G
- Google Scholar
1. Hanly JJ
2. Francescutti CM
3. Loh LS
4. Corning O
5. Long DJ
6. Nakatani MA
7. Porter AH
8. Martin A
(2023) Genetics of yellow-orange color variation in a pair of sympatric sulphur butterflies
Cell Reports 42:112820.

https://doi.org/10.1016/j.celrep.2023.112820
- PubMed
- Google Scholar
1. Hira S
2. Okamoto T
3. Fujiwara M
4. Kita H
5. Kobayashi S
6. Mukai M
(2013) Binding of Drosophila maternal Mamo protein to chromatin and specific DNA sequences
Biochemical and Biophysical Research Communications 438:156–160.

https://doi.org/10.1016/j.bbrc.2013.07.045
- Google Scholar
1. Hu YG
2. Shen YH
3. Zhang Z
4. Shi GQ
(2013) Melanin and urate act to prevent ultraviolet damage in the integument of the silkworm, Bombyx mori
Archives of Insect Biochemistry and Physiology 83:41–55.

https://doi.org/10.1002/arch.21096
- PubMed
- Google Scholar
1. Huynh KD
2. Bardwell VJ
(1998) The BCL-6 POZ domain and other POZ domains interact with the co-repressors N-CoR and SMRT
Oncogene 17:2473–2484.

https://doi.org/10.1038/sj.onc.1202197
- PubMed
- Google Scholar
1. Ito K
2. Kidokoro K
3. Katsuma S
4. Shimada T
5. Yamamoto K
6. Mita K
7. Kadono-Okuda K
(2012) Positional cloning of a gene responsible for the cts mutation of the silkworm, Bombyx mori
Genome 55:493–504.

https://doi.org/10.1139/g2012-033
- PubMed
- Google Scholar
(2006) Regulation of body pigmentation by the Abdominal-B Hox protein and its gain and loss in Drosophila evolution
Cell 125:1387–1399.

https://doi.org/10.1016/j.cell.2006.04.043
- PubMed
- Google Scholar
1. Jin H
2. Seki T
3. Yamaguchi J
4. Fujiwara H
(2019) Prepatterning of Papilio xuthus caterpillar camouflage is controlled by three homeobox genes: clawless, abdominal-A, and Abdominal-B
Science Advances 05:eaav7569.

https://doi.org/10.1126/sciadv.aav7569
- Google Scholar
1. Jolma A
2. Taipale J
(2011) Methods for analysis of transcription factor DNA-binding specificity in vitro
Sub-Cellular Biochemistry 52:155–173.

https://doi.org/10.1007/978-90-481-9069-0_7
- PubMed
- Google Scholar
1. Kalay G
2. Lachowiec J
3. Rosas U
4. Dome MR
5. Wittkopp P
(2019) Redundant and cryptic enhancer activities of the Drosophila yellow gene
Genetics 212:343–360.

https://doi.org/10.1534/genetics.119.301985
- Google Scholar
(2007) Drosophila cuticular proteins with the R&R Consensus: Annotation and classification with a new tool for discriminating RR-1 and RR-2 sequences
Insect Biochemistry and Molecular Biology 37:754–760.

https://doi.org/10.1016/j.ibmb.2007.03.007
- Google Scholar
1. Kawaguchi Y
2. Kusakabe T
3. Lee JM
4. Koga K
(2007) Manifestation of the Dilute Black (bd) Mutation and Constitution of the bd Locus in Bombyx mori
Journal of the Faculty of Agriculture, Kyushu University 52:355–359.

https://doi.org/10.5109/9343
- Google Scholar
(2022) SilkBase: an integrated transcriptomic and genomic database for Bombyx mori and related species
Database 2022:baac040.

https://doi.org/10.1093/database/baac040
- PubMed
- Google Scholar
1. Kopp A
2. Duncan I
3. Godt D
4. Carroll SB
(2000) Genetic control and evolution of sexually dimorphic characters in Drosophila
Nature 408:553–559.

https://doi.org/10.1038/35046017
- PubMed
- Google Scholar
1. Kramer KJ
2. Kanost MR
3. Hopkins TL
4. Jiang H
5. Zhu YC
6. Xu R
7. Kerwin JL
8. Turecek F
(2001) Oxidative conjugation of catechols with proteins in insect skeletal systems
Tetrahedron 57:385–392.

https://doi.org/10.1016/S0040-4020(00)00949-2
- Google Scholar
1. Kvon EZ
2. Kazmar T
3. Stampfel G
4. Yáñez-Cuna JO
5. Pagani M
6. Schernhuber K
7. Dickson BJ
8. Stark A
(2014) Genome-scale functional characterization of Drosophila developmental enhancers in vivo
Nature 512:91–95.

https://doi.org/10.1038/nature13395
- PubMed
- Google Scholar
1. Liu J
2. Chen Z
3. Xiao Y
4. Asano T
5. Li S
6. Peng L
7. Chen E
8. Zhang J
9. Li W
10. Zhang Y
11. Tong X
12. Kadono-Okuda K
13. Zhao P
14. He N
15. Arunkumar KP
16. Gopinathan KP
17. Xia Q
18. Willis JH
19. Goldsmith MR
20. Mita K
(2021) Lepidopteran wing scales contain abundant cross-linked film-forming histidine-rich cuticular proteins
Communications Biology 4:491.

https://doi.org/10.1038/s42003-021-01996-4
- PubMed
- Google Scholar
1. Lloyd VJ
2. Nadeau NJ
(2021) The evolution of structural colour in butterflies
Current Opinion in Genetics & Development 69:28–34.

https://doi.org/10.1016/j.gde.2021.01.004
- PubMed
- Google Scholar
1. Long D
2. Cheng X
3. Hao Z
4. Sun J
5. Umuhoza D
6. Liu Y
7. Chen L
8. Xiang Z
9. Yang F
10. Zhao A
(2021) Genetic hybridization of highly active exogenous functional proteins into silk-based materials using “light-clothing” strategy
Matter 4:2039–2058.

https://doi.org/10.1016/j.matt.2021.03.020
- Google Scholar
1. Lu K
2. Pan Y
3. Shen J
4. Yang L
5. Zhan C
6. Liang S
7. Tai S
8. Wan L
9. Li T
10. Cheng T
11. Ma B
12. Pan G
13. He N
14. Lu C
15. Westhof E
16. Xiang Z
17. Han M-J
18. Tong X
19. Dai F
(2024) SilkMeta: a comprehensive platform for sharing and exploiting pan-genomic and multi-omic silkworm data
Nucleic Acids Research 52:D1024–D1032.

https://doi.org/10.1093/nar/gkad956
- PubMed
- Google Scholar
1. Luo J-W
2. An E-X
3. Lu Y-R
4. Yang L
5. Gai T-T
6. He S-Z
7. Wu S-Y
8. Hu H
9. Li C-L
10. Lu C
11. Tong X-L
12. Dai F-Y
(2021) Molecular basis of the silkworm mutant re^l causing red egg color and embryonic death
Insect Science 28:1290–1299.

https://doi.org/10.1111/1744-7917.12871
- PubMed
- Google Scholar
1. Matsuoka Y
2. Monteiro A
(2018) Melanin pathway genes regulate color and morphology of butterfly wing scales
Cell Reports 24:56–65.

https://doi.org/10.1016/j.celrep.2018.05.092
- PubMed
- Google Scholar
1. Monroe JG
2. Srikant T
3. Carbonell-Bejerano P
4. Becker C
5. Lensink M
6. Exposito-Alonso M
7. Klein M
8. Hildebrandt J
9. Neumann M
10. Kliebenstein D
11. Weng M-L
12. Imbert E
13. Ågren J
14. Rutter MT
15. Fenster CB
16. Weigel D
(2022) Mutation bias reflects natural selection in Arabidopsis thaliana
Nature 602:101–105.

https://doi.org/10.1038/s41586-021-04269-6
- Google Scholar
1. Mukai M
2. Hayashi Y
3. Kitadate Y
4. Shigenobu S
5. Arita K
6. Kobayashi S
(2007) MAMO, a maternal BTB/POZ-Zn-finger protein enriched in germline progenitors is required for the production of functional eggs in Drosophila
Mechanisms of Development 124:570–583.

https://doi.org/10.1016/j.mod.2007.05.001
- Google Scholar
1. Mun S
2. Noh MY
3. Dittmer NT
4. Muthukrishnan S
5. Kramer KJ
6. Kanost MR
7. Arakane Y
(2015) Cuticular protein with a low complexity sequence becomes cross-linked during insect cuticle sclerotization and is required for the adult molt
Scientific Reports 5:10484.

https://doi.org/10.1038/srep10484
- PubMed
- Google Scholar
(2014) Two major cuticular proteins are required for assembly of horizontal laminae and vertical pore canals in rigid cuticle of Tribolium castaneum
Insect Biochemistry and Molecular Biology 53:22–29.

https://doi.org/10.1016/j.ibmb.2014.07.005
- Google Scholar
(2015) Tribolium castaneum RR-1 cuticular protein TcCPR4 is required for formation of pore canals in rigid cuticle
PLOS Genetics 11:e1004963.

https://doi.org/10.1371/journal.pgen.1004963
- PubMed
- Google Scholar
(2016) Cuticle formation and pigmentation in beetles
Current Opinion in Insect Science 17:1–9.

https://doi.org/10.1016/j.cois.2016.05.004
- PubMed
- Google Scholar
(2017) Development and ultrastructure of the rigid dorsal and flexible ventral cuticles of the elytron of the red flour beetle, Tribolium castaneum
Insect Biochemistry and Molecular Biology 91:21–33.

https://doi.org/10.1016/j.ibmb.2017.11.003
- PubMed
- Google Scholar
1. Oba Y
2. Konishi K
3. Yano D
4. Shibata H
5. Kato D
6. Shirai T
(2020) Resurrecting the ancient glow of the fireflies
Science Advances 6:eabc5705.

https://doi.org/10.1126/sciadv.abc5705
- PubMed
- Google Scholar
1. Osanai-Futahashi M
2. Tatematsu K-I
3. Yamamoto K
4. Narukawa J
5. Uchino K
6. Kayukawa T
7. Shinoda T
8. Banno Y
9. Tamura T
10. Sezutsu H
(2012) Identification of the Bombyx red egg gene reveals involvement of a novel transporter family gene in late steps of the insect ommochrome biosynthesis pathway
The Journal of Biological Chemistry 287:17706–17714.

https://doi.org/10.1074/jbc.M111.321331
- PubMed
- Google Scholar
(2001) Design and selection of novel Cys2His2 zinc finger proteins
Annual Review of Biochemistry 70:313–340.

https://doi.org/10.1146/annurev.biochem.70.1.313
- PubMed
- Google Scholar
1. Pan PL
2. Ye YX
3. Lou YH
4. Lu JB
5. Cheng C
6. Shen Y
7. Moussian B
8. Zhang CX
(2018) A comprehensive omics analysis and functional survey of cuticular proteins in the brown planthopper
PNAS 115:5175–5180.

https://doi.org/10.1073/pnas.1716951115
- PubMed
- Google Scholar
1. Persikov AV
2. Singh M
(2014) De novo prediction of DNA-binding specificities for Cys2His2 zinc finger proteins
Nucleic Acids Research 42:97–108.

https://doi.org/10.1093/nar/gkt890
- PubMed
- Google Scholar
(2022) Antennapedia and optix regulate metallic silver wing scale development and cell shape in Bicyclus anynana butterflies
Cell Reports 40:111052.

https://doi.org/10.1016/j.celrep.2022.111052
- PubMed
- Google Scholar
Preprint
1. Prakash A
2. Balakrishnan D
3. Daurer BJ
4. Finet C
5. Lim YC
(2023) Nanoscale Cuticle Density Variations Correlate with Pigmentation and Color in Butterfly Wing Scales
arXiv.

https://arxiv.org/abs/2305.16628
- Google Scholar
(2023) Characterizing cis-regulatory elements using single-cell epigenomics
Nature Reviews. Genetics 24:21–43.

https://doi.org/10.1038/s41576-022-00509-1
- PubMed
- Google Scholar
1. Rebers JE
2. Riddiford LM
(1988) Structure and expression of a Manduca sexta larval cuticle gene homologous to Drosophila cuticle genes
Journal of Molecular Biology 203:411–423.

https://doi.org/10.1016/0022-2836(88)90009-5
- PubMed
- Google Scholar
1. Reed RD
2. Papa R
3. Martin A
4. Hines HM
5. Counterman BA
6. Pardo-Diaz C
7. Jiggins CD
8. Chamberlain NL
9. Kronforst MR
10. Chen R
11. Halder G
12. Nijhout HF
13. McMillan WO
(2011) optix drives the repeated convergent evolution of butterfly wing pattern mimicry
Science 333:1137–1141.

https://doi.org/10.1126/science.1208227
- PubMed
- Google Scholar
1. Rogers WA
2. Grover S
3. Stringer SJ
4. Parks J
5. Rebeiz M
6. Williams TM
(2014) A survey of the trans-regulatory landscape for Drosophila melanogaster abdominal pigmentation
Developmental Biology 385:417–432.

https://doi.org/10.1016/j.ydbio.2013.11.013
- PubMed
- Google Scholar
(2022) GenBank
Nucleic Acids Research 50:D161–D164.

https://doi.org/10.1093/nar/gkab1135
- PubMed
- Google Scholar
(2009) Gold bugs and beyond: a review of iridescence and structural colour mechanisms in beetles (Coleoptera)
Journal of the Royal Society, Interface 6 Suppl 2:S165–S184.

https://doi.org/10.1098/rsif.2008.0354.focus
- PubMed
- Google Scholar
1. Staller P
2. Peukert K
3. Kiermaier A
4. Seoane J
5. Lukas J
6. Karsunky H
7. Möröy T
8. Bartek J
9. Massagué J
10. Hänel F
11. Eilers M
(2001) Repression of p15INK4b expression by Myc through association with Miz-1
Nature Cell Biology 3:392–399.

https://doi.org/10.1038/35070076
- PubMed
- Google Scholar
1. Stevens M
2. Ruxton GD
(2012) Linking the evolution and form of warning coloration in nature
Proceedings. Biological Sciences 279:417–426.

https://doi.org/10.1098/rspb.2011.1932
- PubMed
- Google Scholar
(2000) Density-dependent aposematism in the desert locust
Proceedings. Biological Sciences 267:63–68.

https://doi.org/10.1098/rspb.2000.0967
- PubMed
- Google Scholar
1. Syed AJ
2. Anderson JC
(2021) Applications of bioluminescence in biotechnology and beyond
Chemical Society Reviews 50:5668–5705.

https://doi.org/10.1039/d0cs01492c
- PubMed
- Google Scholar
1. Tajiri R
2. Ogawa N
3. Fujiwara H
4. Kojima T
(2017) Mechanical control of whole body shape by a single cuticular protein obstructor-E in Drosophila melanogaster
PLOS Genetics 13:e1006548.

https://doi.org/10.1371/journal.pgen.1006548
- PubMed
- Google Scholar
(2007) CPF and CPFL, two related gene families encoding cuticular proteins of Anopheles gambiae and other insects
Insect Biochemistry and Molecular Biology 37:675–688.

https://doi.org/10.1016/j.ibmb.2007.03.011
- Google Scholar
1. Tong XL
2. Qiao L
3. Luo JW
4. Ding X
5. Wu SY
(2021) The evolution and genetics of lepidopteran egg and caterpillar coloration
Current Opinion in Genetics & Development 69:140–146.

https://doi.org/10.1016/j.gde.2021.04.006
- Google Scholar
1. Tong X
2. Han MJ
3. Lu K
4. Tai S
5. Liang S
6. Liu Y
7. Hu H
8. Shen J
9. Long A
10. Zhan C
11. Ding X
12. Liu S
13. Gao Q
14. Zhang B
15. Zhou L
16. Tan D
17. Yuan Y
18. Guo N
19. Li YH
20. Wu Z
21. Liu L
22. Li C
23. Lu Y
24. Gai T
25. Zhang Y
26. Yang R
27. Qian H
28. Liu Y
29. Luo J
30. Zheng L
31. Lou J
32. Peng Y
33. Zuo W
34. Song J
35. He S
36. Wu S
37. Zou Y
38. Zhou L
39. Cheng L
40. Tang Y
41. Cheng G
42. Yuan L
43. He W
44. Xu J
45. Fu T
46. Xiao Y
47. Lei T
48. Xu A
49. Yin Y
50. Wang J
51. Monteiro A
52. Westhof E
53. Lu C
54. Tian Z
55. Wang W
56. Xiang Z
57. Dai F
(2022) High-resolution silkworm pan-genome provides genetic insights into artificial selection and ecological adaptation
Nature Communications 13:33366.

https://doi.org/10.1038/s41467-022-33366-x
- Google Scholar
(2005) Aposematism and crypsis combined as a result of distance dependence: functional versatility of the colour pattern in the swallowtail butterfly larva
Proceedings of the Royal Society B 272:1315–1321.

https://doi.org/10.1098/rspb.2005.3079
- Google Scholar
1. Wang N
2. Zhang Y
3. Li W
4. Peng Z
5. Pan H
6. Li S
7. Cheng T
8. Liu C
(2022) Abnormal overexpression of SoxD enhances melanin synthesis in the Ursa mutant of Bombyx mori
Insect Biochemistry and Molecular Biology 149:103832.

https://doi.org/10.1016/j.ibmb.2022.103832
- PubMed
- Google Scholar
1. Willis JH
(2010) Structural cuticular proteins from arthropods: Annotation, nomenclature, and sequence characteristics in the genomics era
Insect Biochemistry and Molecular Biology 40:189–204.

https://doi.org/10.1016/j.ibmb.2010.02.001
- Google Scholar
(2002) Reciprocal functions of the Drosophila yellow and ebony proteins in the development and evolution of pigment patterns
Development 129:1849–1858.

https://doi.org/10.1242/dev.129.8.1849
- PubMed
- Google Scholar
(2003a) Evolution in black and white: genetic control of pigment patterns in Drosophila
Trends in Genetics 19:495–504.

https://doi.org/10.1016/S0168-9525(03)00194-X
- PubMed
- Google Scholar
(2003b) Drosophila pigmentation evolution: divergent genotypes underlying convergent phenotypes
PNAS 100:1808–1813.

https://doi.org/10.1073/pnas.0336368100
- PubMed
- Google Scholar
1. Wittkopp PJ
2. Kalay G
(2012) Cis-regulatory elements: molecular mechanisms and evolutionary processes underlying divergence
Nature Reviews Genetics 13:59–69.

https://doi.org/10.1038/nrg3095
- Google Scholar
(2000) DNA recognition by Cys2His2 zinc finger proteins
Annual Review of Biophysics and Biomolecular Structure 29:183–212.

https://doi.org/10.1146/annurev.biophys.29.1.183
- PubMed
- Google Scholar
1. Wu S
2. Tong X
3. Peng C
4. Xiong G
5. Lu K
6. Hu H
7. Tan D
8. Li C
9. Han M
10. Lu C
11. Dai F
(2016) Comparative analysis of the integument transcriptomes of the black dilute mutant and the wild-type silkworm Bombyx mori
Scientific Reports 6:26114.

https://doi.org/10.1038/srep26114
- PubMed
- Google Scholar
1. Xin Y
2. Le Poul Y
3. Ling L
4. Museridze M
5. Mühling B
6. Jaenichen R
7. Osipova E
8. Gompel N
(2020) Enhancer evolutionary co-option through shared chromatin accessibility input
PNAS 117:20636–20644.

https://doi.org/10.1073/pnas.2004003117
- PubMed
- Google Scholar
1. Xiong G
2. Tong X
3. Gai T
4. Li C
5. Qiao L
6. Monteiro A
7. Hu H
8. Han M
9. Ding X
10. Wu S
11. Xiang Z
12. Lu C
13. Dai F
(2017) Body shape and coloration of silkworm larvae are influenced by a novel cuticular protein
Genetics 207:1053–1066.

https://doi.org/10.1534/genetics.117.300300
- PubMed
- Google Scholar
1. Yamaguchi J
2. Banno Y
3. Mita K
4. Yamamoto K
5. Ando T
6. Fujiwara H
(2013) Periodic Wnt1 expression in response to ecdysteroid generates twin-spot markings on caterpillars
Nature Communications 4:1857.

https://doi.org/10.1038/ncomms2778
- PubMed
- Google Scholar
1. Yan Z
2. Tong X
3. Xiong G
4. Yang W
5. Lu K
6. Yuan Y
7. Han M
8. Hu H
9. Wei W
10. Dai F
(2022) A blueprint of microstructures and stage-specific transcriptome dynamics of cuticle formation in Bombyx mori
International Journal of Molecular Sciences 23:5155.

https://doi.org/10.3390/ijms23095155
- PubMed
- Google Scholar
1. Yang C-H
2. Yang P-C
3. Zhang S-F
4. Shi Z-Y
5. Kang L
6. Zhang A-B
(2017) Identification, expression pattern, and feature analysis of cuticular protein genes in the pine moth Dendrolimus punctatus (Lepidoptera: Lasiocampidae)
Insect Biochemistry and Molecular Biology 83:94–106.

https://doi.org/10.1016/j.ibmb.2017.03.003
- Google Scholar
1. Yoda S
2. Yamaguchi J
3. Mita K
4. Yamamoto K
5. Banno Y
6. Ando T
7. Daimon T
8. Fujiwara H
(2014) The transcription factor Apontic-like controls diverse colouration pattern in caterpillars
Nature Communications 5:4936.

https://doi.org/10.1038/ncomms5936
- PubMed
- Google Scholar
1. Zhao X
2. Gou X
3. Liu W
4. Ma E
5. Moussian B
6. Li S
7. Zhu K
8. Zhang J
(2019) The wing-specific cuticular protein LmACP7 is essential for normal wing morphogenesis in the migratory locust
Insect Biochemistry and Molecular Biology 112:103206.

https://doi.org/10.1016/j.ibmb.2019.103206
- PubMed
- Google Scholar
1. Zou Y-L
2. Ding X
3. Zhang L
4. Xu L-F
5. Liang S-B
6. Hu H
7. Dai F-Y
8. Tong X-L
(2022) Bmmp influences wing morphology by regulating anterior-posterior and proximal-distal axes development
Insect Science 29:1569–1582.

https://doi.org/10.1111/1744-7917.12998
- PubMed
- Google Scholar

Article and author information

Author details

Songyuan Wu

State Key Laboratory of Resource Insects, Key Laboratory of Sericultural Biology and Genetic Breeding, Ministry of Agriculture and Rural Affairs, College of Sericulture, Textile and Biomass Sciences, Southwest University, Chongqing, China

Contribution
Data curation, Funding acquisition, Validation, Investigation, Visualization, Methodology, Writing - original draft, Writing - review and editing

Competing interests
No competing interests declared

"This ORCID iD identifies the author of this article:" 0000-0002-5732-8539
Xiaoling Tong

State Key Laboratory of Resource Insects, Key Laboratory of Sericultural Biology and Genetic Breeding, Ministry of Agriculture and Rural Affairs, College of Sericulture, Textile and Biomass Sciences, Southwest University, Chongqing, China

Contribution
Supervision, Writing - review and editing, Discussion

Competing interests
No competing interests declared

"This ORCID iD identifies the author of this article:" 0000-0002-2649-899X
Chenxing Peng

State Key Laboratory of Resource Insects, Key Laboratory of Sericultural Biology and Genetic Breeding, Ministry of Agriculture and Rural Affairs, College of Sericulture, Textile and Biomass Sciences, Southwest University, Chongqing, China

Contribution
Investigation

Competing interests
No competing interests declared
Jiangwen Luo

State Key Laboratory of Resource Insects, Key Laboratory of Sericultural Biology and Genetic Breeding, Ministry of Agriculture and Rural Affairs, College of Sericulture, Textile and Biomass Sciences, Southwest University, Chongqing, China

Contribution
Investigation

Competing interests
No competing interests declared

"This ORCID iD identifies the author of this article:" 0009-0004-2084-0725
Chenghao Zhang

State Key Laboratory of Resource Insects, Key Laboratory of Sericultural Biology and Genetic Breeding, Ministry of Agriculture and Rural Affairs, College of Sericulture, Textile and Biomass Sciences, Southwest University, Chongqing, China

Contribution
Investigation

Competing interests
No competing interests declared
Kunpeng Lu

State Key Laboratory of Resource Insects, Key Laboratory of Sericultural Biology and Genetic Breeding, Ministry of Agriculture and Rural Affairs, College of Sericulture, Textile and Biomass Sciences, Southwest University, Chongqing, China

Contribution
Investigation

Competing interests
No competing interests declared
Chunlin Li

State Key Laboratory of Resource Insects, Key Laboratory of Sericultural Biology and Genetic Breeding, Ministry of Agriculture and Rural Affairs, College of Sericulture, Textile and Biomass Sciences, Southwest University, Chongqing, China

Contribution
Writing - review and editing, Discussion

Competing interests
No competing interests declared
Xin Ding

State Key Laboratory of Resource Insects, Key Laboratory of Sericultural Biology and Genetic Breeding, Ministry of Agriculture and Rural Affairs, College of Sericulture, Textile and Biomass Sciences, Southwest University, Chongqing, China

Contribution
Investigation

Competing interests
No competing interests declared
Xiaohui Duan

State Key Laboratory of Resource Insects, Key Laboratory of Sericultural Biology and Genetic Breeding, Ministry of Agriculture and Rural Affairs, College of Sericulture, Textile and Biomass Sciences, Southwest University, Chongqing, China

Contribution
Investigation

Competing interests
No competing interests declared
Yaru Lu

State Key Laboratory of Resource Insects, Key Laboratory of Sericultural Biology and Genetic Breeding, Ministry of Agriculture and Rural Affairs, College of Sericulture, Textile and Biomass Sciences, Southwest University, Chongqing, China

Contribution
Investigation

Competing interests
No competing interests declared
Hai Hu

State Key Laboratory of Resource Insects, Key Laboratory of Sericultural Biology and Genetic Breeding, Ministry of Agriculture and Rural Affairs, College of Sericulture, Textile and Biomass Sciences, Southwest University, Chongqing, China

Contribution
Investigation

Competing interests
No competing interests declared
Duan Tan

State Key Laboratory of Resource Insects, Key Laboratory of Sericultural Biology and Genetic Breeding, Ministry of Agriculture and Rural Affairs, College of Sericulture, Textile and Biomass Sciences, Southwest University, Chongqing, China

Contribution
Investigation

Competing interests
No competing interests declared
Fangyin Dai

State Key Laboratory of Resource Insects, Key Laboratory of Sericultural Biology and Genetic Breeding, Ministry of Agriculture and Rural Affairs, College of Sericulture, Textile and Biomass Sciences, Southwest University, Chongqing, China

Contribution
Resources, Supervision, Funding acquisition, Writing - review and editing

For correspondence
fydai@swu.edu.cn

Competing interests
No competing interests declared

"This ORCID iD identifies the author of this article:" 0000-0002-0215-2177

Funding

National Natural Science Foundation of China (32002230)

Songyuan Wu

National Natural Science Foundation of China (32330102)

Fangyin Dai

National Natural Science Foundation of China (U20A2058)

Fangyin Dai

Fundamental Research Funds for the Central Universities (SWU120024)

Songyuan Wu

National Key Research and Development Program of China (2023YFD1600901)

Fangyin Dai

National Key Research and Development Program of China (2023YFF1103801)

Fangyin Dai

Natural Science Foundation of Chongqing (cstc2021jcyj-cxtt0005)

Fangyin Dai

High-level Talents Program of Southwest University (SWURC2021001)

Fangyin Dai

The funders had no role in study design, data collection, and interpretation, or the decision to submit the work for publication.

Acknowledgements

This work was supported by the National Natural Science Foundation of China (Nos. 32002230, 32330102, and U20A2058), the Fundamental Research Funds for the Central Universities in China (No. SWU120024), National Key Research and Development Program (Project Nos. 2023YFD1600901 and 2023YFF1103801), the Natural Science Foundation of Chongqing, China (No. cstc2021jcyj-cxtt0005), and the High-level Talents Program of Southwest University (No. SWURC2021001).

Version history

Preprint posted: April 7, 2023
Sent for peer review: July 13, 2023
Reviewed Preprint version 1: September 1, 2023
Reviewed Preprint version 2: January 25, 2024
Reviewed Preprint version 3: March 19, 2024
Version of Record published: April 8, 2024

Cite all versions

You can cite all versions using the DOI https://doi.org/10.7554/eLife.90795. This DOI represents all versions, and will always resolve to the latest one.

Copyright

This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.