Sequence-specific retrotransposition of 28S rDNA-specific LINE R2Ol in human cells.

R2 is a long interspersed element (LINE) found in a specific sequence of the 28S rDNA among a wide variety of animals. Recently, we observed that R2Ol isolated from medaka fish, Oryzias latipes, retrotransposes sequence specifically into the target sequence of zebrafish. Because the 28S target and flanking regions are widely conserved among vertebrates, we examined whether R2Ol can also integrate in a sequence-specific manner in human cells. Using adenovirus-mediated expression of R2Ol constructs, we confirmed an accurate insertion of R2Ol into the 28S target of human 293T cells. However, the R2Ol mutant devoid of endonuclease (EN) activity showed no retrotransposition ability, suggesting that the sequence-specific integration of R2Ol into 28S rDNA occurs via the cleavage activity of EN. By introducing both R2Ol helper virus and donor plasmid in human cells, we succeeded in retrotransposing an exogenous EGFP gene into the 28S target site by the trans-complementation system, which enabled simplification of specific gene knock-in in a time-efficient manner. We believe that R2Ol may provide an alternative targeted gene knock-in method for practical applications such as gene therapy in future.

Mechanical Control of Whole Body Shape by a Single Cuticular Protein Obstructor-E in Drosophila melanogaster.

Body shapes are much more variable than body plans. One way to alter body shapes independently of body plans would be to mechanically deform bodies. To what extent body shapes are regulated physically, or molecules involved in physical control of morphogenesis, remain elusive. During fly metamorphosis, the cuticle (exoskeleton) covering the larval body contracts longitudinally and expands laterally to become the ellipsoidal pupal case (puparium). Here we show that Drosophila melanogaster Obstructor-E (Obst-E) is a protein constituent of the larval cuticle that confers the oriented contractility/expandability. In the absence of obst-E function, the larval cuticle fails to undergo metamorphic shape change and finally becomes a twiggy puparium. We present results indicating that Obst-E regulates the arrangement of chitin, a long-chain polysaccharide and a central component of the insect cuticle, and directs the formation of supracellular ridges on the larval cuticle. We further show that Obst-E is locally required for the oriented shape change of the cuticle during metamorphosis, which is associated with changes in the morphology of those ridges. Thus, Obst-E dramatically affects the body shape in a direct, physical manner by controlling the mechanical property of the exoskeleton.

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 (pS ), 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.

The pivotal role of aristaless in development and evolution of diverse antennal morphologies in moths and butterflies.

BACKGROUND: Antennae are multi-segmented appendages and main odor-sensing organs in insects. In Lepidoptera (moths and butterflies), antennal morphologies have diversified according to their ecological requirements. While diurnal butterflies have simple, rod-shaped antennae, nocturnal moths have antennae with protrusions or lateral branches on each antennal segment for high-sensitive pheromone detection. A previous study on the Bombyx mori (silk moth) antenna, forming two lateral branches per segment, during metamorphosis has revealed the dramatic change in expression of antennal patterning genes to segmentally reiterated, branch-associated pattern and abundant proliferation of cells contributing almost all the dorsal half of the lateral branch. Thus, localized cell proliferation possibly controlled by the branch-associated expression of antennal patterning genes is implicated in lateral branch formation. Yet, actual gene function in lateral branch formation in Bombyx mori and evolutionary mechanism of various antennal morphologies in Lepidoptera remain elusive. RESULTS: We investigated the function of several genes and signaling specifically in lateral branch formation in Bombyx mori by the electroporation-mediated incorporation of siRNAs or morpholino oligomers. Knock down of aristaless, a homeobox gene expressed specifically in the region of abundant cell proliferation within each antennal segment, during metamorphosis resulted in missing or substantial shortening of lateral branches, indicating its importance for lateral branch formation. aristaless expression during metamorphosis was lost by knock down of Distal-less and WNT signaling but derepressed by knock down of Notch signaling, suggesting the strict determination of the aristaless expression domain within each antennal segment by the combinatorial action of them. In addition, analyses of pupal aristaless expression in antennae with various morphologies of several lepidopteran species revealed that the aristaless expression pattern has a striking correlation with antennal shapes, whereas the segmentally reiterated expression pattern was observed irrespective of antennal morphologies. CONCLUSIONS: Our results presented here indicate the significance of aristaless function in lateral branch formation in B. mori and imply that the diversification in the aristaless expression pattern within each antennal segment during metamorphosis is one of the significant determinants of antennal morphologies. According to these findings, we propose a mechanism underlying development and evolution of lepidopteran antennae with various morphologies.

Electroporation-mediated somatic transgenesis for rapid functional analysis in insects.

Transgenesis is a powerful technique for determining gene function; however, it is time-consuming. It is virtually impossible to carry out in non-model insects in which egg manipulation and screening are difficult. We have established a rapid genetic functional analysis system for non-model insects using a low-cost electroporator (costing under US$200) designed for somatic transformation with the piggyBac transposon. Using this system, we successfully generated somatic transgenic cell clones in various target tissues (e.g. olfactory neurons, wing epidermis, larval epidermis, muscle, fat body and trachea) of the silkworm Bombyx mori during development. We also induced stable and transient RNA interference (RNAi) using short hairpin RNA (shRNA)-mediating DNA vectors and direct transfer of small interfering RNAs (siRNAs), respectively. We found that these electroporation-mediated approaches could also be applied to the swallowtail butterfly Papilio xuthus and the red flour beetle Tribolium castaneum. Thus, this method could be a powerful genetic tool for elucidating various developmental phenomena in non-model insects.

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

BACKGROUND: Body coloration is an ecologically important trait that is often involved in prey-predator interactions through mimicry and crypsis. Although this subject has attracted the interest of biologists and the general public, our scientific knowledge on the subject remains fragmentary. In the caterpillar of the swallowtail butterfly Papilio xuthus, spectacular changes in the color pattern are observed; the insect mimics bird droppings (mimetic pattern) as a young larva, and switches to a green camouflage coloration (cryptic pattern) in the final instar. Despite the wide variety and significance of larval color patterns, few studies have been conducted at a molecular level compared with the number of studies on adult butterfly wing patterns. RESULTS: To obtain a catalog of genes involved in larval mimetic and cryptic pattern formation, we constructed expressed sequence tag (EST) libraries of larval epidermis for P. xuthus, and P. polytes that contained 20,736 and 5,376 clones, respectively, representing one of the largest collections available in butterflies. A comparison with silkworm epidermal EST information revealed the high expression of putative blue and yellow pigment-binding proteins in Papilio species. We also designed a microarray from the EST dataset information, analyzed more than five stages each for six markings, and confirmed spatial expression patterns by whole-mount in situ hybridization. Hence, we succeeded in elucidating many novel marking-specific genes for mimetic and cryptic pattern formation, including pigment-binding protein genes, the melanin-associated gene yellow-h3, the ecdysteroid synthesis enzyme gene 3-dehydroecdysone 3b-reductase, and Papilio-specific genes. We also found many cuticular protein genes with marking specificity that may be associated with the unique surface nanostructure of the markings. Furthermore, we identified two transcription factors, spalt and ecdysteroid signal-related E75, as genes expressed in larval eyespot markings. This finding suggests that E75 is a strong candidate mediator of the hormone-dependent coordination of larval pattern formation. CONCLUSIONS: This study is one of the most comprehensive molecular analyses of complicated morphological features, and it will serve as a new resource for studying insect mimetic and cryptic pattern formation in general. The wide variety of marking-associated genes (both regulatory and structural genes) identified by our screening indicates that a similar strategy will be effective for understanding other complex traits.

Functional unit of supergene in female-limited Batesian mimicry of Papilio polytes.

Supergenes are sets of genes and genetic elements that are inherited like a single gene and control complex adaptive traits, but their functional roles and units are poorly understood. In Papilio polytes, female-limited Batesian mimicry is thought to be regulated by a ∼130 kb inversion region (highly diversified region: HDR) containing 3 genes, UXT, U3X, and doublesex (dsx) which switches non-mimetic and mimetic types. To determine the functional unit, we here performed electroporation-mediated RNAi analyses (and further Crispr/Cas9 for UXT) of genes within and flanking the HDR in pupal hindwings. We first clarified that non-mimetic dsx-h had a function to form the non-mimetic trait in female and only dsx-H isoform 3 had an important function in the formation of mimetic traits. Next, we found that UXT was involved in making mimetic-type pale-yellow spots and adjacent gene sir2 in making red spots in hindwings, both of which refine more elaborate mimicry. Furthermore, downstream gene networks of dsx, U3X, and UXT screened by RNA sequencing showed that U3X upregulated dsx-H expression and repressed UXT expression. These findings demonstrate that a set of multiple genes, not only inside but also flanking HDR, can function as supergene members, which extends the definition of supergene unit than we considered before. Also, our results indicate that dsx functions as the switching gene and some other genes such as UXT and sir2 within the supergene unit work as the modifier gene.

Dramatic changes in patterning gene expression during metamorphosis are associated with the formation of a feather-like antenna by the silk moth, Bombyx mori.

Many moths use sex pheromones to find their mates in the dark. Their antennae are well developed with lateral branches to receive the pheromone efficiently. However, how these structures have evolved remains elusive, because the mechanism of development of these antennae has not been studied at a molecular level. To elucidate the developmental mechanism of this type of antenna, we observed morphogenesis, cell proliferation, cell death and antennal patterning gene expression in the branched antenna of the silk moth, Bombyx mori. Region-specific cell proliferation and almost ubiquitous apoptosis occur during early pupal stages and appear to shape the lateral branch cooperatively. Antennal patterning genes are expressed in a pattern largely conserved among insects with branchless antennae until the late 5th larval instar but most of them change their expression dramatically to a pattern prefiguring the lateral branch during metamorphosis. These findings imply that although antennal primordium is patterned by conserved mechanisms before metamorphosis, most of the antennal patterning genes are reused to form the lateral branch during metamorphosis. We propose that the acquisition of a new regulatory circuit of antennal patterning genes may have been an important event during evolution of the sensory antenna with lateral branches in the Lepidoptera.

Coevolution of telomeric repeats and telomeric repeat-specific non-LTR retrotransposons in insects.

In the telomeres of the silkworm Bombyx mori, telomeric repeat-specific non-long terminal repeat (LTR) retrotransposon SARTBm1 is accumulated in the TTAGG telomeric repeats. Here, we identify novel telomeric repeat-specific non-LTR retrotransposons, SARTTc family, from the red flour beetle Tribolium castaneum in the unconventional TCAGG telomeric repeats. To compare the sequence specificity of SARTBm1 and SARTTc1, we developed a comparable ex vivo retrotransposition assay. Both SARTBm1 and SARTTc1 preferred the telomeric sequence of their hosts, suggesting that the target specificity of these retrotransposons coevolved with their host's telomeric repeats. Swapping experiment indicated that the endonuclease domain is involved in recognizing the target sequence. Moreover, SARTBm1 proteins could retrotranspose 3'untranslated region (UTR) sequence of SARTTc1 as well as their own 3'UTR, whereas SARTTc1 proteins could only retrotranspose their own 3'UTRs. These results provide insights to the mechanism and divergence of sequence specificity and 3'UTR recognition in non-LTR retrotransposons.

Genomic architecture and functional unit of mimicry supergene in female limited Batesian mimic Papilio butterflies.

It has long been suggested that dimorphic female-limited Batesian mimicry of two closely related Papilio butterflies, Papilio memnon and Papilio polytes, is controlled by supergenes. Whole-genome sequencing, genome-wide association studies and functional analyses have recently identified mimicry supergenes, including the doublesex (dsx) gene. Although supergenes of both the species are composed of highly divergent regions between mimetic and non-mimetic alleles and are located at the same chromosomal locus, they show critical differences in genomic architecture, particularly with or without an inversion: P. polytes has an inversion, but P. memnon does not. This review introduces and compares the detailed genomic structure of mimicry supergenes in two Papilio species, including gene composition, repetitive sequence composition, breakpoint/boundary site structure, chromosomal inversion and linkage disequilibrium. Expression patterns and functional analyses of the respective genes within or flanking the supergene suggest that dsx and other genes are involved in mimetic traits. In addition, structural comparison of the corresponding region for the mimicry supergene among further Papilio species suggests three scenarios for the evolution of the mimicry supergene between the two Papilio species. The structural features revealed in the Papilio mimicry supergene provide insight into the formation, maintenance and evolution of supergenes. This article is part of the theme issue 'Genomic architecture of supergenes: causes and evolutionary consequences'.

doublesex Controls Both Hindwing and Abdominal Mimicry Traits in the Female-Limited Batesian Mimicry of Papilio memnon.

Papilio butterflies are known to possess female-limited Batesian mimicry polymorphisms. In Papilio memnon, females have mimetic and non-mimetic forms, whereas males are monomorphic and non-mimetic. Mimetic females are characterized by color patterns and tails in the hindwing and yellow abdomens. Recently, an analysis of whole-genome sequences has shown that an approximately 160 kb region of chromosome 25 is responsible for mimicry and has high diversity between mimetic (A) and non-mimetic (a) alleles (highly diversified region: HDR). The HDR includes three genes, UXT, doublesex (dsx), and Nach-like, but the functions of these genes are unknown. Here, we investigated the function of dsx, a gene involved in sexual differentiation, which is expected to be functionally important for hindwing and abdominal mimetic traits in P. memnon. Expression analysis by reverse transcription quantitative PCR (RT-qPCR) and RNA sequencing showed that mimetic dsx (dsx-A) was highly expressed in the hindwings in the early pupal stage. In the abdomen, both dsx-A and dsx-a were highly expressed during the early pupal stage. When dsx was knocked down using small interfering RNAs (siRNAs) designed in the common region of dsx-A and dsx-a, a male-like pattern appeared on the hindwings of mimetic and non-mimetic females. Similarly, when dsx was knocked down in the abdomen, the yellow scales characteristic of mimetic females changed to black. Furthermore, when dsx-a was specifically knocked down, the color pattern of the hindwings changed, as in the case of dsx knockdown in non-mimetic females but not mimetic females. These results suggest that dsx-a is involved in color pattern formation on the hindwings of non-mimetic females, whereas dsx-A is involved in hindwing and abdominal mimetic traits. dsx was involved in abdominal and hindwing mimetic traits, but dsx expression patterns in the hindwing and abdomen were different, suggesting that different regulatory mechanisms may exist. Our study is the first to show that the same gene (dsx) regulates both the hindwing and abdominal mimetic traits. This is the first functional analysis of abdominal mimicry in butterflies.

A corset function of exoskeletal ECM promotes body elongation in Drosophila.

Body elongation is a general feature of development. Postembryonically, the body needs to be framed and protected by extracellular materials, such as the skeleton, the skin and the shell, which have greater strength than cells. Thus, body elongation after embryogenesis must be reconciled with those rigid extracellular materials. Here we show that the exoskeleton (cuticle) coating the Drosophila larval body has a mechanical property to expand less efficiently along the body circumference than along the anteroposterior axis. This "corset" property of the cuticle directs a change in body shape during body growth from a relatively round shape to an elongated one. Furthermore, the corset property depends on the functions of Cuticular protein 11 A and Tubby, protein components of a sub-surface layer of the larval cuticle. Thus, constructing a stretchable cuticle and supplying it with components that confer circumferential stiffness is the fly's strategy for executing postembryonic body elongation.

Genetic switch in UV response of mimicry-related pale-yellow colors in Batesian mimic butterfly, Papilio polytes.

In a Batesian mimic butterfly Papilio polytes, mimetic females resemble an unpalatable model, Pachliopta aristolochiae, but exhibit a different color pattern from nonmimetic females and males. In particular, the pale-yellow region on hind wings, which correspondingly sends important putative signals for mimicry and mate preference, is different in shape and chemical features between nonmimetic and mimetic morphs. Recently, we found that mimetic-type doublesex [dsx (H)] causes mimetic traits; however, the control of dimorphic pale-yellow colors remains unclear. Here, we revealed that dsx (H) switched the pale-yellow colors from UV-excited fluorescent type (nonmimetic) to UV-reflecting type (mimetic), by repressing the papiliochrome II synthesis genes and nanostructural changes in wing scales. Photoreceptor reactivities showed that some birds and butterflies could effectively recognize mimetic and nonmimetic pale-yellow colors, suggesting that a genetic switch in the UV response of pale-yellow colors may play essential roles in establishing the dimorphic female-limited Batesian mimicry.

Species-specific coordinated gene expression and trans-regulation of larval color pattern in three swallowtail butterflies.

The diversity of butterfly larval color pattern has been attracted to people since Darwin's time; however, its molecular mechanisms still remain largely unknown. Larval body markings often differ completely between closely related species under natural selection. The final instar larvae of the swallowtail butterflies Papilio xuthus and Papilio polytes show a green camouflage pattern, whereas those of Papilio machaon show a warning color pattern, although P. xuthus and P. machaon are closely related species. To identify the genes that contribute to species divergence, we compared the expression pattern of eight pigment-associated genes between three Papilio species. The spatial expression pattern of melanin-related genes coincided with the species-specific cuticular markings. We newly found that the combination of bilin-binding protein and yellow-related gene (YRG) correlated perfectly with larval blue, yellow, and green coloration. To distinguish whether the interspecific differences in pigment-associated genes are caused by cis-regulatory changes or distribution differences in trans-regulatory proteins, we compared species-specific mRNA expression in an F1 hybrid specimen. Px-YRG and Pp-YRG showed a similar expression pattern, suggesting that the change in expression of YRG is caused mainly by changes in the distribution of trans-regulatory proteins. Our findings shed light on the gene regulatory networks for butterfly larval color pattern.

Caterpillar color patterns are determined by a two-phase melanin gene prepatterning process: new evidence from tan and laccase2.

The larval color patterns in Lepidoptera exhibit splendid diversity, and identifying the genes responsible for pigment distribution is essential to understanding color-pattern evolution. The swallowtail butterfly, Papilio xuthus, is a good candidate for analyzing marking-associated genes because its body markings change dramatically at the final molt. Moreover, the silkworm Bombyx mori is most suitable for identification of lab-generated color mutants because genome information and many color mutants are available. Here, we analyzed the expression pattern of 10 melanin-related genes in P. xuthus, and analyzed whether these genes were responsible for Bombyx larval color mutants. We found that seven genes correlated strongly with the stage-specific larval cuticular markings of P. xuthus, suggesting that, compared with Drosophila, more genes showed marking specificity in lepidopteran larvae. We newly found that the expression of both tan and laccase2 is strongly correlated with the larval black markings in both P. xuthus and B. mori. The results of F2 linkage analysis and mutant analysis strongly suggest that tan is the responsible gene for Bombyx larval color mutant rouge, and that tan is important in emphasizing black markings of lepidopteran larvae. Detailed comparison of temporal and spatial expression patterns showed that larval cuticular markings were regulated at two different phases. Marking-specific expression of oxidizing enzymes preceded the marking-specific expression of melanin synthesis enzymes at mRNA level, which is the reverse of the melanin synthesis step.

Essential factors involved in the precise targeting and insertion of telomere-specific non-LTR retrotransposon, SART1Bm.

Telomere length maintenance is essential for most eukaryotes to ensure genome stability and integrity. A non-long terminal repeat (LTR) retrotransposon, SART1Bm, targets telomeric repeats (TTAGG)n of the silkworm Bombyx mori and is presumably involved in telomere length maintenance. However, how many telomeric repeats are required for its retrotransposition and how reverse transcription is initiated at the target site are not well understood. Here, using an ex vivo and trans-in vivo recombinant baculovirus retrotransposition system, we demonstrated that SART1Bm requires at least three (TTAGG) telomeric repeats and a longer poly(A) tail for its accurate retrotransposition. We found that SART1Bm retrotransposed only in the third (TTAGG) tract of three repeats and that the A residue of the (TTAGG) unit was essential for its retrotransposition. Interestingly, SART1Bm also retrotransposed into telomeric repeats of other species, such as human (TTAGGG)n repeats, albeit with low retrotransposition efficiency. We further showed that the reverse transcription of SART1Bm occurred inaccurately at the internal site of the 3' untranslated region (UTR) when using a short poly(A) tail but at the accurate site when using a longer poly(A) tail. These findings promote our understanding of the general mechanisms of site-specific retrotransposition and aid the development of a site-specific gene knock-in tool.

Batesian mimicry has evolved with deleterious effects of the pleiotropic gene doublesex.

Dimorphic female-limited Batesian mimicry in the swallowtail butterfly Papilio polytes is regulated by the supergene locus H, harbouring the mimetic (H) and non-mimetic (h) doublesex (dsx) gene. In the present study, we demonstrated that dsx-H negatively affects the number of eggs laid, hatching rate, larval survival rate, and adult lifespan. When crossed with hh males, the number of eggs laid of mimetic females (genotype HH) was lower than that of non-mimetic females (hh). Moreover, hh and Hh females laid fewer eggs when crossed with HH males. The hatching and larval survival rates were lower when both female and male parents harboured dsx-H. The adult lifespan of HH females was shorter than that of hh females, while it was similar in males regardless of the genotype. These findings suggest the presence of a cost-benefit balance of Batesian mimicry, which is evolved to avoid predation but is accompanied by physiological deficits, in this species.

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.

Positional cloning of a Bombyx wingless locus flugellos (fl) reveals a crucial role for fringe that is specific for wing morphogenesis.

Mutations at the flügellos (fl) locus in Bombyx mori produce wingless pupae and moths because of the repressed response of wing discs to ecdysteroid. Four recessive fl alleles occurred spontaneously and were mapped at 13.0 of the silkworm genetic linkage group 10. By positional cloning, we confirmed that the gene responsible for fl is fringe (fng) encoding Fng glycosyltransferase, which is involved in regulating the Notch signaling pathway. In four different fl alleles, we detected a large deletion of the fng gene in fl(k) and nonsense mutations in fl, fl(o), and fl(n). In the wild-type (WT) silkworm, fng is expressed actively in the wing discs, brain, and reproductive organs from the fourth to final instars but barely in the other tissues tested. In situ hybridization showed that fng mRNA is expressed in the dorsal layer of the WT wing discs. The wingless (wg) mRNA, a downstream marker of Fng-mediated Notch signaling, is localized at the dorsoventral boundary in the WT wing discs but repressed markedly in the fl wing discs. Although null mutants of Drosophila fng result in postembryonic lethality, loss of fng function in Bombyx affects only wing morphogenesis, suggesting different essential roles for fng in tissue differentiation among insects.

yellow and ebony are the responsible genes for the larval color mutants of the silkworm Bombyx mori.

Many larval color mutants have been obtained in the silkworm Bombyx mori. Mapping of melanin-synthesis genes on the Bombyx linkage map revealed that yellow and ebony genes were located near the chocolate (ch) and sooty (so) loci, respectively. In the ch mutants, body color of neonate larvae and the body markings of elder instar larvae are reddish brown instead of normal black. Mutations at the so locus produce smoky larvae and black pupae. F(2) linkage analyses showed that sequence polymorphisms of yellow and ebony genes perfectly cosegregated with the ch and so mutant phenotypes, respectively. Both yellow and ebony were expressed in the epidermis during the molting period when cuticular pigmentation occurred. The spatial expression pattern of yellow transcripts coincided with the larval black markings. In the ch mutants, nonsense mutations of the yellow gene were detected, whereas large deletions of the ebony ORF were detected in the so mutants. These results indicate that yellow and ebony are the responsible genes for the ch and so loci, respectively. Our findings suggest that Yellow promotes melanization, whereas Ebony inhibits melanization in Lepidoptera and that melanin-synthesis enzymes play a critical role in the lepidopteran larval color pattern.