Comprehensive microarray-based analysis for stage-specific larval camouflage pattern-associated genes in the swallowtail butterfly, Papilio xuthus

Intrapuparial stage aging and PMI estimation based on the developmental transcriptomes of forensically important Aldrichina grahami (Diptera: Calliphoridae) gene expression

Background

The expression profiles of differentially expressed genes (DEGs) during pupal development have been demonstrated to be vital in age estimation of forensic entomological study. Here, using forensically important Aldrichina grahami (Diptera: Calliphoridae), we aimed to explore the potential of intrapuparial stage aging and postmortem interval (PMI) estimation based on characterization of successive developmental transcriptomes and gene expression patterns.

Methods

We collected A. grahami pupae at 11 successive intrapuparial stages at 20 °C and used the RNA-seq technique to build the transcriptome profiles of their intrapuparial stages. The DEGs were identified during the different intrapuparial stages using comparative transcriptome analysis. The selected marker DEGs were classified and clustered for intrapuparial stage aging and PMI estimation and then further verified for transcriptome data validation. Ultimately, we categorized the overall gene expression levels as the dependent variable and the age of intrapuparial A. grahami as the independent variable to conduct nonlinear regression analysis.

Results

We redefined the intrapuparial stages of A. grahami into five key successive substages (I, II, III, IV, and V), based on the overall gene expression patterns during pupal development. We screened 99 specific time-dependent expressed genes (stage-specific DEGs) to determine the different intrapuparial stages based on comparison of the gene expression levels during the 11 developmental intrapuparial stages of A. grahami. We observed that 55 DEGs showed persistent upregulation during the development of intrapuparial A. grahami. We then selected four DEGs (act79b, act88f, up and ninac) which presented consistent upregulation using RT-qPCR (quantitative real-time PCR) analysis, along with consideration of the maximum fold changes during the pupal development. We conducted nonlinear regression analysis to simulate the calculations of the relationships between the expression levels of the four selected DEGs and the developmental time of intrapuparial A. grahami and constructed fitting curves. The curves demonstrated that act79b and ninac showed continuous relatively increasing levels.

Conclusions

This study redefined the intrapuparial stages of A. grahami based on expression profiles of developmental transcriptomes for the first time. The stage-specific DEGs and those with consistent tendencies of expression were found to have potential in age estimation of intrapuparial A. grahami and could be supplementary to a more accurate prediction of PMI.

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.

Ovariole-specific Yellow-g and Yellow-g2 proteins are required for fecundity and egg chorion rigidity in the red flour beetle, Tribolium castaneum

Most insects reproduce by laying eggs that have an eggshell/chorion secreted by follicle cells, which serves as a protective barrier for developing embryos. Thus, eggshell formation is vital for reproduction. Insect yellow family genes encode for secreted extracellular proteins that perform different, context-dependent functions in different tissues at various stages of development involving, for example, cuticle/eggshell coloration and morphology, molting, courtship behavior and embryo hatching. In this study we investigated the function of two of this family's genes, yellow-g (TcY-g) and yellow-g2 (TcY-g2), on the formation and morphology of the eggshell of the red flour beetle, Tribolium castaneum. Real-time PCR analysis revealed that both TcY-g and TcY-g2 were specifically expressed in the ovarioles of adult females. Loss of function produced by injection of double-stranded RNA (dsRNA) for either TcY-g or TcY-g2 gene resulted in failure of oviposition. There was no effect on maternal survival. Ovaries dissected from those dsRNA-treated females exhibited ovarioles containing not only developing oocytes but also mature eggs in their egg chambers. However, the ovulated eggs were collapsed and ruptured, resulting in swollen lateral oviducts and calyxes. TEM analysis showed that lateral oviducts were filled with electron-dense material, presumably from some cellular content leakage out of the collapsed eggs. In addition, morphological abnormalities in lateral oviduct epithelial cells and the tubular muscle sheath were evident. These results support the hypothesis that both TcY-g and TcY-g2 proteins are required for maintaining the rigidity and integrity of the chorion, which is critical for resistance to mechanical stress and/or rehydration during ovulation and egg activation in the oviducts of T. castaneum. Because Yellow-g and Yellow-g2 are highly conserved among insect species, both genes are potential targets for development of gene-based insect pest population control methods.

Outer fold is sole effective tissue among three mantle folds with regard to oyster shell colour

Molluscs constitute the second largest phylum of animals in the world, and shell colour is one of their most important phenotypic characteristics. In this study, we found among three folds on the mantle edge of oyster, only the outer fold had the same colour as the shell. Transcriptome and mantle cutting experiment indicated that the outer fold may be mainly reflected in chitin framework formation and biomineralisation. There were obvious differences in SEM structure and protein composition between the black and white shell periostraca. The black shell periostraca had more proteins related to melanin biosynthesis and chitin binding. Additionally, we identified an uncharacterized protein gene (named as CgCBP) ultra-highly expressed only in the black outer fold and confirmed its function of chitin-binding and CaCO₃ precipitation promoting. RNAi also indicated that CgCBP knockdown could change the structure of shell periostracum and reduce shell pigmentation. All these results suggest that the mantle outer fold plays multiple key roles in shell periostraca bioprocessing, and shell periostracum structure affected by chitin-binding protein is functionally correlated with shell pigmentation. The investigation of oyster shell periostracum structure and shell colour will provide a better understanding in pigmentation during biological mineralisation in molluscs.

Transcriptome of the pygmy grasshopper Formosatettix qinlingensis (Orthoptera: Tetrigidae)

Formosatettix qinlingensis (Zheng, 1982) is a tiny grasshopper endemic to Qinling in China. For further study of its transcriptomic features, we obtained RNA-Seq data by Illumina HiSeq X Ten sequencing platform. Firstly, transcriptomic analysis showed that transcriptome read numbers of two female and one male samples were 25,043,314, 24,429,905, and 25,034,457, respectively. We assembled 65,977 unigenes, their average length was 1,072.09 bp, and the length of N50 was 2,031 bp. The average lengths of F. qinlingensis female and male unigenes were 911.30 bp, and 941.82 bp, and the N50 lengths were 1,745 bp and 1,735 bp, respectively. Eight databases were used to annotate the functions of unigenes, and 23,268 functional unigenes were obtained. Besides, we also studied the body color, immunity and insecticide resistance of F. qinlingensis. Thirty-nine pigment-related genes were annotated. Some immunity genes and signaling pathways were found, such as JAK-STAT and Toll-LIKE receptor signaling pathways. There are also some insecticide resistance genes and signal pathways, like nAChR, GST and DDT. Further, some of these genes were differentially expressed in female and male samples, including pigment, immunity and insecticide resistance. The transcriptomic study of F. qinlingensis will provide data reference for gene prediction and molecular expression study of other Tetrigidae species in the future. Differential genetic screening of males and females provides a basis for studying sex and immune balance in insects.

The Plastic Larval Body Color of the Pale Grass Blue Butterfly Zizeeria maha (Lepidoptera: Lycaenidae) in Response to the Host Plant Color: The Maternal Effect on Crypsis

Many lepidopteran larvae show body color polyphenism, and their colors may be cryptic on the host plant leaves. To elucidate the effect of the host plant color on the plastic larval body color, we focused on the lycaenid butterfly Zizeeria maha, which shows various larval body colors ranging from green to red, even within a sibling group. We showed that oviposition was normally performed on both green and red leaves, despite a green preference and the fact that the larvae grew equally by consuming either green or red leaves. The number of red larvae decreased from the second instar stage to the fourth instar stage, demonstrating a stage-dependent variation. When the larvae were fed either green or red leaves across multiple generations of the lineages, the red larvae were significantly more abundant in the red leaf lineage than in the green leaf lineage. Moreover, the red-fed siblings showed a significantly higher red larval frequency than the green-fed siblings in the red-leaf lineage but not in the green-leaf lineage. These results suggest that, in this butterfly species, the plastic larval body color for crypsis may be affected not only by the color of the leaves that the larvae consume (single-generation effect) but also by the color of the leaves that their mothers consume (maternal effect), in addition to a stage-dependent color variation.

Spatiotemporal Patterning Enabled by Gene Regulatory Networks

Spatiotemporal pattern formation plays a key role in various biological phenomena including embryogenesis and neural network formation. Though the reaction-diffusion systems enabling pattern formation have been studied phenomenologically, the biomolecular mechanisms behind these processes have not been modeled in detail. Here, we study the emergence of spatiotemporal patterns due to simple, synthetic and commonly observed two- and three-node gene regulatory network motifs coupled with their molecular diffusion in one- and two-dimensional space. We investigate the patterns formed due to the coupling of inherent multistable and oscillatory behavior of the toggle switch, toggle switch with double self-activation, toggle triad, and repressilator with the effect of spatial diffusion of these molecules. We probe multiple parameter regimes corresponding to different regions of stability (monostable, multistable, oscillatory) and assess the impact of varying diffusion coefficients. This analysis offers valuable insights into the design principles of pattern formation facilitated by these network motifs, and it suggests the mechanistic underpinnings of biological pattern formation.

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

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

Integrated Analysis of Transcriptome and Proteome to Reveal Pupal Color Switch in Papilio xuthus Butterflies

Pupal color polyphenism in Papilio butterflies, including green, intermediate, or brown, is an excellent study system for understanding phenotypic plasticity. Previous studies suggested that development of brown pupae may be controlled by a hormone called pupal-cuticle-melanizing-hormone (PCMH) which is synthesized and secreted from brain-suboesophageal ganglion and prothoracic ganglion complexes (Br-SG-TG1) during the pre-pupa stage. However, detailed molecular mechanisms of neuroendocrine regulation in pupal color development remain unknown. In this study, we integrated the expression profiles of transcriptome and proteome at pre-pupa stages [2 h after gut purge (T1) and 3 h after forming the garter around the body (T2)] and pigmentation stages [10 h after ecdysis (T3) and 24 h after ecdysis (T4)] to identify important genes and pathways underlying the development of green and brown pupa in the swallowtail butterfly Papilio xuthus. Combined comparisons of each developmental stage and each tissue under green and brown conditions, a total of 1042 differentially expressed genes (DEGs) and 430 different abundance proteins (DAPs) were identified. Weighted gene co-expression network analysis (WGCNA) and enrichment analysis indicate that these DEGs were mainly related to oxidation-reduction, structural constituent of cuticle, and pigment binding. Soft clustering by Mfuzz and enrichment analysis indicate that these DAPs are mainly involved in tyrosine metabolism, insect hormone biosynthesis, and melanogenesis. By homologous alignment, we further identified those genes encoding neuropeptides (51), GPCRs (116), G-proteins (8), cuticular proteins (226), chitinases (16), and chitin deacetylases (8) in the whole genome of P. xuthus and analyzed their expression profiles. Although we identified no gene satisfying with hypothesized expression profile of PCMH, we found some genes in the neuropeptide cascade showed differentially expressed under two pupal color conditions. We also found that Toll signaling pathway genes, juvenile hormone (JH) related genes, and multiple cuticular proteins play important roles in the formation of selective pupal colors during the prepupal-pupal transition. Our data also suggest that both green and brown pupa include complex pigment system that is regulated by genes involved in black, blue, and yellow pigments. Our results provide important insights into the evolution of pupal protective colors among swallowtail butterflies.

Yellow-b, -c, -d, and -h are required for normal body coloration of Henosepilachna vigintioctopunctata

The involvement of yellow genes y-b, y-c, y-e, and y-h in cuticle tanning has poorly been clarified. In the present paper, six putative yellow (y-y, y-b, y-c, y-e y-f, and y-h) genes were identified in Henosepilachna vigintioctopunctata. Hvy-b, Hvy-c, Hvy-e, and Hvy-h were abundantly transcribed at early larval and late pupal stages, especially in the epidermis. Accordingly, RNA interference (RNAi) experiments were performed by an injection of dsy-b, dsy-c, dsy-e, or dsy-h into the second instar larvae and 1-day-old pupae. The head capsule, scoli and strumae, and legs in the fourth-instar larvae became blacker; the blackish spots in the pupae were darkened and widened after RNAi of Hvy-b, compared with those of dsegfp-treated controls. Depletion of Hvy-b at the 1-day-old pupal stage expanded two pair of black markings on the sternum of the metathorax, and darkened the black patched on the sterna of the abdomen segments I-VI in the resultant adults. Depletion of Hvy-e caused darker pigmented adult body and elytral cuticles than those of dsegfp-introduced controls. However, there was no obvious difference in pigmentation of the black markings. Hvy-h-deficient larvae displayed dark yellow body color, whereas the body color of the dsegfp-injected control was pale yellow. There was no obvious difference in coloration of larval specific-black markings or pupal cuticle between dsHvy-h- and dsegfp-treated animals. Moreover, silence of Hvy-c at the second instar larval stage lightened black markings in the resulting larvae and pupae, but had no influence on pale yellow body color. Our results demonstrated their different roles of the four yellow genes during body pigmentation: HvY-b and HvY-c, respectively, inhibit and facilitate the coloration within dark markings, whereas HvY-e and HvY-h, respectively, repress the pigmentation in adult and larval body cuticles outside the black patches in H. vigintioctopunctata.

Sexually dimorphic dorsal coloration in a jumping spider: testing a potential case of sex-specific mimicry

To avoid predation, many animals mimic behaviours and/or coloration of dangerous prey. Here we examine potential sex-specific mimicry in the jumping spider Habronattus pyrrithrix. Previous work proposed that males' conspicuous dorsal coloration paired with characteristic leg-waving (i.e. false antennation) imperfectly mimics hymenopteran insects (e.g. wasps and bees), affording protection to males during mate-searching and courtship. By contrast, less active females are cryptic and display less leg-waving. Here we test the hypothesis that sexually dimorphic dorsal colour patterns in H. pyrrithrix are most effective when paired with sex-specific behaviours. We manipulated spider dorsal coloration with makeup to model the opposite sex and exposed them to a larger salticid predator (Phidippus californicus). We predicted that males painted like females should suffer higher predation rates than sham-control males. Likewise, females painted like males should suffer higher predation rates than sham-control females. Contrary to expectations, spiders with male-like coloration were attacked more than those with female-like coloration, regardless of their actual sex. Moreover, males were more likely to be captured, and were captured sooner, than females (regardless of colour pattern). With these unexpected negative results, we discuss alternative functional hypotheses for H. pyrrithrix colours, as well as the evolution of defensive coloration generally.

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.

Characterization and analysis of the transcriptome in Opisina arenosella from different developmental stages using single-molecule real-time transcript sequencing and RNA-seq

Opisina arenosella is one of the main pests harming coconut trees. To date, there have been few studies on the molecular genetics, biochemistry and physiology of O. arenosella at the transcriptional level, and there are no available reference genomes. Here, Illumina RNA sequencing combined with PacBio single-molecule real-time analysis was applied to study the transcriptome of this pest at different developmental stages, providing reference data for transcript expression analysis. Twelve samples of O. arenosella from different stages of development were sequenced using Illumina RNA sequencing, and the pooled RNA samples were sequenced with PacBio technology (Iso-Seq). A full-length transcriptome with 41,938 transcripts was captured, and the N50 and N90 lengths were 3543 bp and 1646 bp, respectively. A total of 36,925 transcripts were annotated in public databases, 6493 of which were long noncoding RNAs, while 2510 represented alternative splicing events. There were significant differences in the gene expression profiles at different developmental stages, with high levels of differential gene expression associated with growth, development, carbohydrate metabolism and immunity. This work provides resources and information for the study of the transcriptome and gene function of O. arenosella and provides a valuable foundation for understanding the changes in gene expression during development.

Yellow-g and Yellow-g2 proteins are required for egg desiccation resistance and temporal pigmentation in the Asian tiger mosquito, Aedes albopictus

Eggs from Aedes mosquitoes exhibit desiccation resistance that helps them to survive and spread as human disease vectors throughout the world. Previous studies have suggested that eggshell/chorion melanization and/or serosal cuticle formation are important for desiccation resistance. In this study, using dsRNAs for target genes, we analyzed the functional importance of two ovary-specific yellow genes, AalY-g and AalY-g2, in the resistance to egg desiccation of the Asian tiger mosquito, Aedes albopictus, a species in which neither the timing of the melanization nor temporal development of the serosal cuticle is correlated with desiccation resistance. Injections of dsAalY-g, dsAalY-g2 or dsAalY-g/g2 (co-injection) into adult females have no effect on their fecundity. However, initial melanization is delayed by 1-2 h with the eggshells eventually becoming black similar to that observed in eggs from dsEGFP-injected control females. In addition, the shape of the eggs from dsAalY-g, -g2 and -g/g2-treated females is abnormally crescent-shaped and the outermost exochorion is more fragile and partially peeled off. dsEGFP control eggs, like those from the wild-type strain, acquire resistance to desiccation between 18 and 24 h after oviposition (HAO). In contrast, ~80% of the 24 HAO dsAalY-g and dsAalY-g2 eggs collapse when they are transferred to a low humidity environment. In addition, there is no electron-dense outer endochorion evident in either dsAalY-g or dsAalY-g2 eggs. These results support the hypothesis that AalY-g and AalY-g2 regulate the timing of eggshell darkening and are required for integrity of the exochorion as well as for rigidity, normal morphology and formation of the outer endochorion, a structure that apparently is critical for desiccation resistance of the Ae. albopictus egg.

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

Excess melanin precursors rescue defective cuticular traits in stony mutant silkworms probably by upregulating four genes encoding RR1-type larval cuticular proteins

Melanin and cuticular proteins are vital cuticle components in insects. Cuticular defects caused by mutations in cuticular protein-encoding genes can obstruct melanin deposition. The effects of changes in melanin on the expression of cuticular protein-encoding genes, the cuticular and morphological traits, and the origins of these effects are unknown. We found that the cuticular physical characteristics and the expression patterns of larval cuticular protein-encoding genes markedly differed between the melanic and non-melanic integument regions. By using four p multiple-allele color pattern mutants with increasing degrees of melanism (+^p, p^M, p^S, and p^B), we found that the degree of melanism and the expression of four RR1-type larval cuticular protein-encoding genes (BmCPR2, BmLcp18, BmLcp22, and BmLcp30) were positively correlated. By modulating the content of melanin precursors and the expression of cuticular protein-encoding genes in cells in tissues and in vivo, we showed that this positive correlation was due to the induction of melanin precursors. More importantly, the melanism trait introduced into the BmCPR2 deletion strain Dazao-stony induced up-regulation of three other similar chitin-binding characteristic larval cuticular protein-encoding genes, thus rescuing the cuticular, morphological and adaptability defects of the Dazao-stony strain. This rescue ability increased with increasing melanism levels. This is the first study reporting the induction of cuticular protein-encoding genes by melanin and the biological importance of this induction in affecting the physiological characteristics of the cuticle.

Gene functions in adult cuticle pigmentation of the yellow mealworm, Tenebrio molitor

In many arthropod species including insects, the cuticle tanning pathway for both pigmentation and sclerotization begins with tyrosine and is responsible for production of both melanin- and quinoid-type pigments, some of which are major pigments for body coloration. In this study we identified and cloned cDNAs of the yellow mealworm, Tenebrio molitor, encoding seven key enzymes involved in this pathway including tyrosine hydroxylase (TmTH), DOPA decarboxylase (TmDDC), laccase 2 (TmLac2), Yellow-y (TmY-y), arylalkylamine N-acetyltransferase (TmAANAT1), aspartate 1-decarboxylase (TmADC) and N-β-alanyldopamine synthase (Tmebony). Expression profiles of these genes during development were analyzed by real-time PCR, revealing development-specific patterns of expression. Loss of function mediated by RNAi of either 1) TmTH or TmLac2, 2) TmDDC or TmY-y, and 3) TmAANAT1, TmADC or Tmebony resulted in pale/white, light yellow/brown and dark/black adult body coloration, respectively. In addition, there are three distinct layer/regional pigmentation differences in rigid types of adult cuticle, a brownish outer exocuticle (EX), a dark pigmented middle mesocuticle (ME) and a transparent inner endocuticle (EN). Decreases in pigmentation of the EX and/or ME layers were observed after RNAi of TmDDC or TmY-y. In TmADC- or Tmebony-deficient adults, a darker pigmented EX layer was observed. In TmAANAT1-deficient adults, trabeculae formed between the dorsal and ventral elytral cuticles as well as the transparent EN layer became highly pigmented. These results demonstrate that knocking down the level of gene expression of specific enzymes of this tyrosine metabolic pathway leads to abnormal pigmentation in individual layers and substructure of the rigid adult exoskeleton of T. molitor.

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.

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.

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.

De Novo Assembly and Characterization of the Xenocatantops brachycerus Transcriptome

Grasshoppers are common pests but also have high nutritional and commercial potential. Xenocatantops brachycerus Willemse (Orthoptera: Acrididae) is an economically important grasshopper species that is reared in China. Using the IlluminaHiSeqTM 4000 platform, three transcriptomes of the adult male, adult female, and nymph of X. brachycerus were sequenced. A total of 133,194,848 clean reads were obtained and de novo assembled into 43,187 unigenes with an average length of 964 bp (N50 of 1799 bp); of these, 24,717 (57.23%) unigenes matched known proteins. Based on these annotations, many putative transcripts related to X. brachycerus growth, development, environmental adaptability, and metabolism of nutritional components and bioactive components were identified. In addition, the expression profiles of all three transcriptome datasets were analyzed, and many differentially expressed genes were detected using RSEM and PossionDis. Unigenes. Unigenes with functions associated with growth and development exhibited higher transcript levels at the nymph stage, and unigenes associated with environmental adaptability showed increased transcription in the adults. These comprehensive X. brachycerus transcriptomic data will provide a useful molecular resource for gene prediction, molecular marker development, and studies on signaling pathways in this species and will serve as a reference for the efficient use of other grasshoppers.

A Tug-of-War Mechanism for Pattern Formation in a Genetic Network

Synthesizing spatial patterns with genetic networks is an ongoing challenge in synthetic biology. A successful demonstration of pattern formation would imply a better understanding of systems in the natural world and advance applications in synthetic biology. In developmental systems, transient patterning may suffice in order to imprint instructions for long-term development. In this paper we show that transient but persistent patterns can emerge from a realizable synthetic gene network based on a toggle switch. We show that a bistable system incorporating diffusible molecules can generate patterns that resemble Turing patterns but are distinctly different in the underlying mechanism: diffusion of mutually inhibiting molecules creates a prolonged "tug-of-war" between patches of cells at opposing bistable states. The patterns are transient but longer wavelength patterns persist for extended periods of time. Analysis of a representative small scale model implies the eigenvalues of the persistent modes are just above the threshold of stability. The results are verified through simulation of biologically relevant models.

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.

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.

spalt is functionally conserved in Locusta and Drosophila to promote wing growth

Locusta has strong fly wings to ensure its long distance migration, but the molecular mechanism that regulates the Locusta wing development is poorly understood. To address the developmental mechanism of the Locusta flying wing, we cloned the Dpp target gene spalt (sal) and analyzed its function in wing growth in the Locusta. The Locusta wing size is apparently reduced with vein defects when sal is interfered by injection of dsRNA, indicating that sal is required for locust wing growth and vein formation. This function is conserved during the Drosophila wing development. To better understand sal’s function in wing growth, we then used Drosophila wing disc as a model for further study. We found that sal promotes cell proliferation in the whole wing disc via positive regulation of a microRNA bantam. Our results firstly unravel sal’s function in the Locusta wing growth and confirm a highly conserved function of sal in Locusta and Drosophila.

Genetic Basis of Melanin Pigmentation in Butterfly Wings

Despite the variety, prominence, and adaptive significance of butterfly wing patterns, surprisingly little is known about the genetic basis of wing color diversity. Even though there is intense interest in wing pattern evolution and development, the technical challenge of genetically manipulating butterflies has slowed efforts to functionally characterize color pattern development genes. To identify candidate wing pigmentation genes, we used RNA sequencing to characterize transcription across multiple stages of butterfly wing development, and between different color pattern elements, in the painted lady butterfly Vanessa cardui This allowed us to pinpoint genes specifically associated with red and black pigment patterns. To test the functions of a subset of genes associated with presumptive melanin pigmentation, we used clustered regularly interspaced short palindromic repeats (CRISPR)/Cas9 genome editing in four different butterfly genera. pale, Ddc, and yellow knockouts displayed reduction of melanin pigmentation, consistent with previous findings in other insects. Interestingly, however, yellow-d, ebony, and black knockouts revealed that these genes have localized effects on tuning the color of red, brown, and ochre pattern elements. These results point to previously undescribed mechanisms for modulating the color of specific wing pattern elements in butterflies, and provide an expanded portrait of the insect melanin pathway.

Cloning, Expression, and Characterization of Prophenoloxidases from Asian Corn Borer, Ostrinia furnacalis (Gunée)

Insect phenoloxidase (PO) belongs to the type 3 copper protein family and possesses oxidoreductase activities. PO is typically synthesized as a zymogen called prophenoloxidase (PPO) and requires the proteolytic activation to function. We here cloned full-length cDNA for 3 previously unidentified PPOs, which we named OfPPO1a, OfPPO1b, and OfPPO3, from Asian corn borer, Ostrinia furnacalis (Gunée), in addition to the previously known OfPPO2. These conceptual PPOs and OfPPO2 all contain two common copper-binding regions, two potential proteolytic activation sites, a plausible thiol-ester site, and a conserved C-terminal region but lack a secretion signal peptide sequence at the N-terminus. O. furnacalis PPOs were highly similar to other insect PPOs (42% to 79% identity) and clustered well with other lepidopteran PPOs. RT-PCR assay showed the transcripts of the 4 OfPPOs were all detected at the highest level in hemocytes and at the increased amounts after exposure to infection by bacteria and fungi. Additionally, we established an Escherichia coli (E. coli) expression system to produce recombinant O. furnacalis PPO proteins for future use in investigating their functions. These insights could provide valuable information for better understanding the activation and functioning mechanisms of O. furnacalis PPOs.

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.

Characterization and analysis of a de novo transcriptome from the pygmy grasshopper Tetrix japonica

The pygmy grasshopper Tetrix japonica is a common insect distributed throughout the world, and it has the potential for use in studies of body colour polymorphism, genomics and the biology of Tetrigoidea (Insecta: Orthoptera). However, limited biological information is available for this insect. Here, we conducted a de novo transcriptome study of adult and larval T. japonica to provide a better understanding of its gene expression and develop genomic resources for future work. We sequenced and explored the characteristics of the de novo transcriptome of T. japonica using Illumina HiSeq 2000 platform. A total of 107 608 206 paired-end clean reads were assembled into 61 141 unigenes using the trinity software; the mean unigene size was 771 bp, and the N50 length was 1238 bp. A total of 29 225 unigenes were functionally annotated to the NCBI nonredundant protein sequences (Nr), NCBI nonredundant nucleotide sequences (Nt), a manually annotated and reviewed protein sequence database (Swiss-Prot), Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) databases. A large number of putative genes that are potentially involved in pigment pathways, juvenile hormone (JH) metabolism and signalling pathways were identified in the T. japonica transcriptome. Additionally, 165 769 and 156 796 putative single nucleotide polymorphisms occurred in the adult and larvae transcriptomes, respectively, and a total of 3162 simple sequence repeats were detected in this assembly. This comprehensive transcriptomic data for T. japonica will provide a usable resource for gene predictions, signalling pathway investigations and molecular marker development for this species and other pygmy grasshoppers.

Comparative analysis of the integument transcriptomes of the black dilute mutant and the wild-type silkworm Bombyx mori

The insect cuticle is a critical protective shell that is composed predominantly of chitin and various cuticular proteins and pigments. Indeed, insects often change their surface pigment patterns in response to selective pressures, such as threats from predators, sexual selection and environmental changes. However, the molecular mechanisms underlying the construction of the epidermis and its pigmentation patterns are not fully understood. Among Lepidoptera, the silkworm is a favorable model for color pattern research. The black dilute (bd) mutant of silkworm is the result of a spontaneous mutation; the larval body color is notably melanized. We performed integument transcriptome sequencing of the wild-type strain Dazao and the mutant strains +/bd and bd/bd. In these experiments, during an early stage of the fourth molt, a stage at which approximately 51% of genes were expressed genome wide (RPKM ≥1) in each strain. A total of 254 novel transcripts were characterized using Cuffcompare and BLAST analyses. Comparison of the transcriptome data revealed 28 differentially expressed genes (DEGs) that may contribute to bd larval melanism, including 15 cuticular protein genes that were remarkably highly expressed in the bd/bd mutant. We suggest that these significantly up-regulated cuticular proteins may promote melanism in silkworm larvae.

Comparative Analysis of Transcriptomes among Bombyx mori Strains and Sexes Reveals the Genes Regulating Melanic Morph and the Related Phenotypes

As a source of insect polymorphism, melanism plays an important role in ecological adaption and usually endows advantageous phenotypic-effects on insects. However, due to the mechanistic diversity, there are knowledge gaps in the molecular mechanisms underlying melanism and the related phenotypes. In silk moths, a recessive melanic mutant (sex-controlled melanism, sml) strain exhibits extended adult longevity. We took a transcriptome approach to perform a comparative analysis between this sml strain and a wild-type strain (Dazao). Our analysis resulted in the identification of 59 unique differentially expressed genes in the melanic mutant. Two key genes (laccase2 and yellow) involved in melanin formation were significantly up-regulated in melanic individuals. The laccase2 B-type isoform (BGIBMGA006746) was found to likely participate in the silkworm cuticular melanism process at late pupal stage. Moreover, we discovered 22 cuticular protein encoding genes with the possible function in melanin transport and/or maintenance. Based on our findings, we presume that the longer survival of the melanic sml male moths might be associated with the enhanced antioxidant defense systems and a reduction in the insulin/IGF-1 signaling pathway (IIS). These findings will facilitate the understanding of the molecular basis underlying melanism and the derived phenotypic-effects.

Morphological characterization and molecular mapping of an irradiation-induced Speckled mutant in the silkworm, Bombyx mori

Speckled (Spc), an X-ray-induced lethal mutant of Bombyx mori, exhibits a mosaic dark-brown-spotted larval epidermis in both sexes and egg-laying problems only in females. Here, we report the morphological characterization and molecular mapping of the Spc mutant. Morphological investigations revealed that the epidermal ultrastructure of the small, dark-brown spots was more dense than that of the white regions in both Spc/+ mutants and wild type, and that the lethality of the Spc/Spc mutants occurred during early embryogenesis. Furthermore, the ovarioles and ovipositor were disconnected in approximately 85.5% of Spc/+ females, a further 2.5% had a connection between the ovarioles and ovipositor that was too narrow to lay eggs. The remaining females showed a normal connection similar to that of the wild type. We successfully narrowed down the location of the Spc mutation to a region on chromosome 4 that was ∼1041 kb long. Gene-prediction analysis identified 25 candidate genes in this region. Chromosome structure analysis indicated that a ∼305 kb deletion was included in the mapping region. Temporal and spatial reverse transcription PCR (RT-PCR) analysis showed that several genes in the mapped region are associated with the Spc mutant. Although the genes responsible for the Spc mutation were not definitively identified, our results further the current understanding of the complex mechanism underlying the multiple morphological defects in Spc mutants.

Rapid identification of Helicoverpa armigera and Helicoverpa zea (Lepidoptera: Noctuidae) using ribosomal RNA internal transcribed spacer 1

Rapid identification of invasive species is crucial for deploying management strategies to prevent establishment. Recent Helicoverpa armigera (Hübner) invasions and subsequent establishment in South America has increased the risk of this species invading North America. Morphological similarities make differentiation of H. armigera from the native Helicoverpa zea (Boddie) difficult. Characteristics of adult male genitalia and nucleotide sequence differences in mitochondrial DNA are two of the currently available methods to differentiate these two species. However, current methods are likely too slow to be employed as rapid detection methods. In this study, conserved differences in the internal transcribed spacer 1 (ITS1) of the ribosomal RNA genes were used to develop species-specific oligonucleotide primers that amplified ITS1 fragments of 147 and 334 bp from H. armigera and H. zea, respectively. An amplicon (83 bp) from a conserved region of 18S ribosomal RNA subunit served as a positive control. Melting temperature differences in ITS1 amplicons yielded species-specific dissociation curves that could be used in high resolution melt analysis to differentiate the two Helicoverpa species. In addition, a rapid and inexpensive procedure for obtaining amplifiable genomic DNA from a small amount of tissue was identified. Under optimal conditions, the process was able to detect DNA from one H. armigera leg in a pool of 25 legs. The high resolution melt analysis combined with rapid DNA extraction could be used as an inexpensive method to genetically differentiate large numbers of H. armigera and H. zea using readily available reagents.

Aspartate Decarboxylase is Required for a Normal Pupa Pigmentation Pattern in the Silkworm, Bombyx mori

The pigmentation pattern of Lepidoptera varies greatly in different development stages. To date, the effects of key genes in the melanin metabolism pathway on larval and adult body color are distinct, yet the effects on pupal pigmentation remains unclear. In the silkworm, Bombyx mori, the black pupa (bp) mutant is only specifically melanized at the pupal stage. Using positional cloning, we found that a mutation in the Aspartate decarboxylase gene (BmADC) is causative in the bp mutant. In the bp mutant, a SINE-like transposon with a length of 493 bp was detected ~2.2 kb upstream of the transcriptional start site of BmADC. This insertion causes a sharp reduction in BmADC transcript levels in bp mutants, leading to deficiency of β-alanine and N-β-alanyl dopamine (NBAD), but accumulation of dopamine. Following injection of β-alanine into bp mutants, the color pattern was reverted that of the wild-type silkworms. Additionally, melanic pupae resulting from knock-down of BmADC in the wild-type strain were obtained. These findings show that BmADC plays a crucial role in melanin metabolism and in the pigmentation pattern of the silkworm pupal stage. Finally, this study contributes to a better understanding of pupa pigmentation patterns in Lepidoptera.

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.

Extensive transcriptional response associated with seasonal plasticity of butterfly wing patterns

In the eastern United States, the buckeye butterfly, Junonia coenia, shows seasonal wing colour plasticity where adults emerging in the spring are tan, while those emerging in the autumn are dark red. This variation can be artificially induced in laboratory colonies, thus making J. coenia a useful model system to examine the mechanistic basis of plasticity. To better understand the developmental basis of seasonal plasticity, we used RNA-seq to quantify transcription profiles associated with development of alternative seasonal wing morphs. Depending on the developmental stage, between 547 and 1420 transfrags were significantly differentially expressed between morphs. These extensive differences in gene expression stand in contrast to the much smaller numbers of differentially expressed transcripts identified in previous studies of genetic wing pattern variation in other species and suggest that environmentally induced phenotypic shifts arise from very broad systemic processes. Analyses of candidate endocrine and pigmentation transcripts revealed notable genes upregulated in the red morph, including several ecdysone-associated genes, and cinnabar, an ommochrome pigmentation gene implicated in colour pattern variation in other butterflies. We also found multiple melanin-related transcripts strongly upregulated in the red morph, including tan and yellow-family genes, leading us to speculate that dark red pigmentation in autumn J. coenia may involve nonommochrome pigments. While we identified several endocrine and pigmentation genes as obvious candidates for seasonal colour morph differentiation, we speculate that the majority of observed expression differences were due to thermal stress response. The buckeye transcriptome provides a basis for further developmental studies of phenotypic plasticity.

Molecular characterization of Rhodnius prolixus' embryonic cuticle

The embryonic cuticle (EC) of Rhodnius prolixus envelopes the entire body of the embryo during hatching and provides physical protection, allowing the embryo to pass through a narrow chorionic border. Most of the knowledge about the EC of insects is derived from studies on ultrastructure and secretion processes during embryonic development, and little is known about the molecular composition of this structure. We performed a comprehensive molecular characterization of the major components extracted from the EC of R. prolixus, and we discuss the role of the different molecules that were identified during the eclosion process. The results showed that, similar to the post-embryonic cuticles of insects, the EC of R. prolixus is primarily composed of carbohydrates (57%), lipids (19%), and proteins (8%). Considering only the carbohydrates, chitin is by far the major component (approximately 70%), and it is found primarily along the body of the EC. It is scarce or absent in its prolongations, which are composed of glycosaminoglycans. In addition to chitin, we also identified amino (15%), neutral (12%) and acidic (3%) carbohydrates in the EC of R. prolixus. In addition carbohydrates, we also identified neutral lipids (64.12%) and phospholipids (35.88%). Proteomic analysis detected 68 proteins (55 were identified and 13 are hypothetical proteins) using the sequences in the R. prolixus genome (http://www.vectorbase.org). Among these proteins, 8 out of 15 are associated with cuticle metabolism. These proteins are unequivocally cuticle proteins, and they have been described in other insects. Approximately 35% of the total proteins identified were classified as having a structural function. Chitin-binding protein, amino peptidase, amino acid oxidase, oxidoreductase, catalase and peroxidase are all proteins associated with cuticle metabolism. Proteins known to be cuticle constituents may be related to the function of the EC in assisting the insect during eclosion. To our knowledge, this is the first study to describe the global molecular composition of an EC in insects.

Transcriptome analysis of integument differentially expressed genes in the pigment mutant (quail) during molting of silkworm, Bombyx mori

In the silkworm Bombyx mori, pigment mutants with diverse body colors have been maintained throughout domestication for about 5000 years. The silkworm larval body color is formed through the mutual interaction of melanin, ommochromes, pteridines and uric acid. These pigments/compounds are synthesized by the cooperative action of various genes and enzymes. Previous reports showed that melanin, ommochrome and pteridine are increased in silkworm quail (q) mutants. To understand the pigment increase and alterations in pigment synthesis in q mutant, transcriptome profiles of the silkworm integument were investigated at 16 h after head capsule slippage in the fourth molt in q mutants and wild-type (Dazao). Compared to the wild-type, 1161 genes were differentially expressed in the q mutant. Of these modulated genes, 62.4% (725 genes) were upregulated and 37.6% (436 genes) were downregulated in the q mutant. The molecular function of differently expressed genes was analyzed by Blast2GO. The results showed that upregulated genes were mainly involved in protein binding, small molecule binding, transferase activity, nucleic acid binding, specific DNA-binding transcription factor activity and chromatin binding, while exclusively down-expressed genes functioned in oxidoreductase activity, cofactor binding, tetrapyrrole binding, peroxidase activity and pigment binding. We focused on genes related to melanin, pteridine and ommochrome biosynthesis; transport of uric acid; and juvenile hormone metabolism because of their importance in integument coloration during molting. This study identified differently expressed genes implicated in silkworm integument formation and pigmentation using silkworm q mutant. The results estimated the number and types of genes that drive new integument formation.