Background: Microsporidian Nosema bombycis has received much attention because the pébrine disease of domesticated silkworms results in great economic losses in the silkworm industry. So far, no effective treatment could be found for pébrine. Compared to other known Nosema parasites, N. bombycis can unusually parasitize a broad range of hosts. To gain some insights into the underlying genetic mechanism of pathological ability and host range expansion in this parasite, a comparative genomic approach is conducted. The genome of two Nosema parasites, N. bombycis and N. antheraeae (an obligatory parasite to undomesticated silkworms Antheraea pernyi), were sequenced and compared with their distantly related species, N. ceranae (an obligatory parasite to honey bees). Results: Our comparative genomics analysis show that the N. bombycis genome has greatly expanded due to the following three molecular mechanisms: 1) the proliferation of host-derived transposable elements, 2) the acquisition of many horizontally transferred genes from bacteria, and 3) the production of abundnant gene duplications. To our knowledge, duplicated genes derived not only from small-scale events (e.g., tandem duplications) but also from large-scale events (e.g., segmental duplications) have never been seen so abundant in any reported microsporidia genomes. Our relative dating analysis further indicated that these duplication events have arisen recently over very short evolutionary time. Furthermore, several duplicated genes involving in the cytotoxic metabolic pathway were found to undergo positive selection, suggestive of the role of duplicated genes on the adaptive evolution of pathogenic ability. Conclusions: Genome expansion is rarely considered as the evolutionary outcome acting on those highly reduced and compact parasitic microsporidian genomes. This study, for the first time, demonstrates that the parasitic genomes can expand, instead of shrink, through several common molecular mechanisms such as gene duplication, horizontal gene transfer, and transposable element expansion. We also showed that the duplicated genes can serve as raw materials for evolutionary innovations possibly contributing to the increase of pathologenic ability. Based on our research, we propose that duplicated genes of N. bombycis should be treated as primary targets for treatment designs against pébrine.

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

Background: Carboxylesterase is a multifunctional superfamily and ubiquitous in all living organisms, including animals, plants, insects, and microbes. It plays important roles in xenobiotic detoxification, and pheromone degradation, neurogenesis and regulating development. Previous studies mainly used Dipteran Drosophila and mosquitoes as model organisms to investigate the roles of the insect COEs in insecticide resistance. However, genome-wide characterization of COEs in phytophagous insects and comparative analysis remain to be performed. Results: Based on the newly assembled genome sequence, 76 putative COEs were identified in Bombyx mori. Relative to other Dipteran and Hymenopteran insects, alpha-esterases were significantly expanded in the silkworm. Genomics analysis suggested that BmCOEs showed chromosome preferable distribution and 55% of which were tandem arranged. Sixty-one BmCOEs were transcribed based on cDNA/ESTs and microarray data. Generally, most of the COEs showed tissue specific expressions and expression level between male and female did not display obvious differences. Three main patterns could be classified, i.e. midgut-, head and integument-, and silk gland-specific expressions. Midgut is the first barrier of xenobiotics peroral toxicity, in which COEs may be involved in eliminating secondary metabolites of mulberry leaves and contaminants of insecticides in diet. For head and integument-class, most of the members were homologous to odorant-degrading enzyme (ODE) and antennal esterase. RT-PCR verified that the ODE-like esterases were also highly expressed in larvae antenna and maxilla, and thus they may play important roles in degradation of plant volatiles or other xenobiotics. Conclusion: B. mori has the largest number of insect COE genes characterized to date. Comparative genomic analysis suggested that the gene expansion mainly occurred in silkworm alpha-esterases. Expression evidence indicated that the expanded genes were specifically expressed in midgut, integument and head, implying that these genes may have important roles in detoxifying secondary metabolites of mulberry leaves, contaminants in diet, and odorants. Our results provide some new insights into functions and evolutionary characteristics of COEs in phytophagous insects.

Background: Glucosidation plays a major role in the inactivation and excretion of a great variety of both endogenous and exogenous compounds. A class of UDP-glycosyltransferases (UGTs) is involved in this process. Insect UGTs play important roles in several processes, including detoxication of substrates such as plant allelochemicals, cuticle formation, pigmentation, and olfaction. Identification and characterization of Bombyx mori UGT genes could provide valuable basic information for this important family and explain the detoxication mechanism and other processes in insects. Results: Taking advantage of the newly assembled genome sequence, we performed a genome-wide analysis of the candidate UGT family in the silkworm, B. mori. Based on UGT signature and their similarity to UGT homologs from other organisms, we identified 42 putative silkworm UGT genes. Most of them are clustered on the silkworm chromosomes, with two major clusters on chromosomes 7 and 28, respectively. The phylogenetic analysis of these identified 42 UGT protein sequences revealed five major groups. A comparison of the silkworm UGTs with homologs from other sequenced insect genomes indicated that some UGTs are silkworm-specific genes. The expression patterns of these candidate genes were investigated with known expressed sequence tags (ESTs), microarray data, and RT-PCR method. In total, 36 genes were expressed in tissues examined and showed different patterns of expression profile, indicating that these UGT genes might have different functions. Conclusion: B. mori possesses a largest insect UGT gene family characterized to date, including 42 genes. Phylogenetic analysis, genomic organization and expression profiles provide an overview for the silkworm UGTs and facilitate their functional studies in future.

Cuticular proteins (CPs) are crucial components of the insect cuticle. Although numerous genes encoding cuticular proteins have been identified in known insect genomes to date, their functions in maintaining insect body shape and adaptability remain largely unknown. In the current study, positional cloning led to the identification of a gene encoding an RR1-type cuticular protein, BmorCPR2, highly expressed in larval chitin-rich tissues and at the mulberry leaf-eating stages, which is responsible for the silkworm stony mutant. In the Dazao-stony strain, the BmorCPR2 allele is a deletion mutation with significantly lower expression, compared to the wild-type Dazao strain. Dysfunctional BmorCPR2 in the stony mutant lost chitin binding ability, leading to reduced chitin content in larval cuticle, limitation of cuticle extension, abatement of cuticle tensile properties, and aberrant ratio between internodes and intersegmental folds. These variations induce a significant decrease in cuticle capacity to hold the growing internal organs in the larval development process, resulting in whole-body stiffness, tightness, and hardness, bulging intersegmental folds, and serious defects in larval adaptability. To our knowledge, this is the first study to report the corresponding phenotype of stony in insects caused by mutation of RR1-type cuticular protein. Our findings collectively shed light on the specific role of cuticular proteins in maintaining normal larval body shape and will aid in the development of pest control strategies for the management of Lepidoptera.

Coloration is one of the most variable characters in animals and provides rich material for studying the developmental genetic basis of pigment patterns. In the silkworm, more than 100 gene mutation systems are related to aberrant color patterns. The melanism (mln) is a rare body color mutant that exhibits an easily distinguishable phenotype in both larval and adult silkworms. By positional cloning, we identified the candidate gene of the mln locus, Bm-iAANAT, whose homologous gene (Dat) converts dopamine into N-acetyldopamine, a precursor for N-acetyldopamine sclerotin in Drosophila. In the mln mutant, two types of abnormal Bm-iAANAT transcripts were identified, whose expression levels are markedly lower than the wild type (WT). Moreover, dopamine content was approximately twice as high in the sclerified tissues (head, thoracic legs, and anal plate) of the mutant as in WT, resulting in phenotypic differences between the two. Quantitative reverse transcription PCR analyses showed that other genes involved in the melanin metabolism pathway were regulated by the aberrant Bm-iAANAT activity in mln mutant in different ways and degrees. We therefore propose that greater accumulation of dopamine results from the functional deficiency of Bm-iAANAT in the mutant, causing a darker pattern in the sclerified regions than in the WT. In summary, our results indicate that Bm-iAANAT is responsible for the color pattern of the silkworm mutant, mln. To our knowledge, this is the first report showing a role for arylalkylamine-N-acetyltransferases in color pattern mutation in Lepidoptera.