Silkworm Genomics and Its Applications: A Comprehensive Review

BackgroundIn contrast to wild species, which have typically evolved phenotypes over long periods of natural selection, domesticates rapidly gained human-preferred agronomic traits in a relatively short-time frame via artificial selection. Under domesticated conditions, many traits can be observed that cannot only be due to environmental alteration. In the case of silkworms, aside from genetic divergence, whether epigenetic divergence played a role in domestication is an unanswered question. The silkworm is still an enigma in that it has two DNA methyltransferases (DNMT1 and DNMT2) but their functionality is unknown. Even in particular the functionality of the widely distributed DNMT1 remains unknown in insects in general.ResultsBy embryonic RNA interference, we reveal that knockdown of silkworm Dnmt1 caused decreased hatchability, providing the first direct experimental evidence of functional significance of insect Dnmt1. In the light of this fact and those that DNA methylation is correlated with gene expression in silkworms and some agronomic traits in domesticated organisms are not stable, we comprehensively compare silk gland methylomes of 3 domesticated (Bombyx mori) and 4 wild (Bombyx mandarina) silkworms to identify differentially methylated genes between the two. We observed 2-fold more differentiated methylated cytosinces (mCs) in domesticated silkworms as compared to their wild counterparts, suggesting a trend of increasing DNA methylation during domestication. Further study of more domesticated and wild silkworms narrowed down the domesticates’ epimutations, and we were able to identify a number of differential genes. One such gene showing demethyaltion in domesticates correspondently displays lower gene expression, and more interestingly, has experienced selective sweep. A methylation-increased gene seems to result in higher expression in domesticates and the function of its Drosophila homolog was previously found to be essential for cell volume regulation, indicating a possible correlation with the enlargement of silk glands in domesticated silkworms.ConclusionsOur results imply epigenetic influences at work during domestication, which gives insight into long time historical controversies regarding acquired inheritance.

Effects of amino acid substitutions on the biological activity of anti-CD20 monoclonal antibody produced by transgenic silkworms (Bombyx mori)

The context sequence motif surrounding the ATG initiation codon influences mRNA translation efficiency and affects protein production; however, the optimal sequence differs among species. To determine the optimal sequence for production of recombinant proteins in a transgenic silkworm, we compared 14-nucleotide context motifs around the ATG (ATG-context) in 50 silkworm genes and found the following consensus: (A/T)AN(A/T)ATCAAAatgN. We were also able to define the least-common motif: CCN(C/G)CGN(C/T/G)(G/C/T)(T/G)atgC, which served as a negative control. To examine the regulatory role of these motifs in protein expression, we constructed reporter plasmids containing different ATG-context motifs together with either the luciferase gene or an enhanced green fluorescent protein (EGFP) gene. These constructs were then used for comparison of luciferase reporter activity and EGFP production in BmN4 cells in vitro as well as in transgenic silkworms in vivo. We detected 10-fold higher luciferase activity in BmN4 cells transfected with the consensus ATG-context motif construct, compared to the negative control plasmid. ELISA measurements of EGFP translation products with the corresponding constructs in BmN4 cells showed consistently similar results. Interestingly, the translation efficiency of the novel consensus ATG-context motif did not show the highest activity in the transgenic silkworms in vivo, except for the fat body. The highest efficiency in the middle and posterior silk glands was produced by the sericin 1 context. Our results show that the ATG-context motifs differ among silkworm tissues. This result is important for the further improvement of the transgenic silkworm system for the production of recombinant proteins.

Although there are some documented examples on population dynamics of transposable elements (TEs) in model organisms, the evolutionary dynamics of TEs in domesticated species has not been systematically investigated. The objective of this study is to understand population dynamics of TEs during silkworm domestication. In this work, using transposon-display we examined the polymorphism of seven TE families [they represent about 59% of silkworm (Bombyx mori) total TE content] in four domesticated silkworm populations and one wild silkworm population. Maximum likelihood (ML) was used to estimate selection pressure. Population differentiation and structure were performed by using AMOVA analysis and program DISTRUCT, respectively. The results of transposon-display showed that significant differentiation occurred between the domesticated silkworm and wild silkworm. These TEs have experienced expansions and fixation in the domesticated silkworm but not in wild silkworm. Furthermore, the ML results indicated that purifying selection of TEs in the domesticated silkworm were significantly weaker than that in the wild silkworm. Interestingly, an adaptation insertion induced by BmMITE-2 was found, and this insertion can reduce the polymorphism of the flanking regions of its neighboring COQ7 gene. Our results suggested that TEs expanded and were fixed in the domesticated silkworm might result from demographic effects and artificial selection during domestication. We concluded that the data presented in this study have general implication in animal and crop improvements as well as in domestication of new species.

Identification of stage- and tissue-specific cis-regulatory elements will enable more precise genomic editing. In previous studies of the silkworm Bombyx mori, we identified and characterized several tissue- and sex-specific cis-regulatory elements using transgenic technology, including a female- and fat body-specific promoter, vitellogenin, testis-specific promoters, Radial spoke head 1 (BmR1) and beta-tubulin 4 (Bmβ4). Here we report a cis-regulatory element specific for a somatic and germ cell-expressed promoter, nanos (Bmnos). We investigated activities of three truncated promoter sequences upstream of the transcriptional initiation site sequences of Bmnos in vitro (nos-0.6kb, nos-1kb and nos-2kb) and in vivo (nos-2kb). In BmN cultured cells, all three lengths drove expression of the gene encoding enhanced green fluorescence protein (EGFP), although nos-2kb had the highest fluorescence activity. In transgenic silkworms, nos-2kb drove EGFP expression at the early embryonic stage, and fluorescence was concentrated in the gonads at later embryonic stages. In addition, this cis-regulatory element was not sex differentiated. The fluorescence intensity gradually weakened following the larval developmental stage in the gonads and were broadly expressed in the whole body. The nos-2kb promoter drove the Cas9 system with efficiency comparable to that of the broad-spectrum strong IE1 promoter. These results indicate that Bmnos is an effective endogenous cis-regulatory element in the early embryo and in the gonad that can be used in applications involving the clustered, regularly interspaced, short palindromic repeats (CRISPR)/Cas9 system.

Adiponectin is an adipocyte hormone involved in glucose and lipid metabolism. The aim of this study was to develop a human adiponectin expression system in transgenic silkworm using a human adiponectin expression vector. The silk gland of the silkworm is a highly specialized organ that has the wonderful ability to synthesize and secrete silk protein. To express human adiponectin in the silk gland of transgenic silkworm, targeting vectors pB-A3-adiponectin-IRES-RFP and pB-Ser1-adiponectin-IRES-RFP were constructed and then introduced into the silkworm pupa. The transgenic silkworms were verified by PCR and then generated. The level of adiponectin in the transgenic silkworm was 6-10 ng/50 mg of freeze-dried powder, and western blotting using an antibody against human adiponectin demonstrated a specific band with a molecular weight of 30 kDa in the silkworm. These results showed that human adiponectin introduced into the silkworm genome was expressed successfully on a large-scale.

The Janus kinase/signal transducer and activator of transcription (JAK/STAT) signaling pathway plays a crucial role in the viral immune processes of organisms, with STAT being a key transcription factor downstream in this pathway. The STAT gene of Bombyx mori has two splicing forms, named BmSTAT-S and BmSTAT-L. This study compared the effects of the two splicing forms on Bombyx mori nucleopolyhedrovirus (BmNPV) infection through cell-level interference and further explored whether BmSTAT participates in the immune response to BmNPV infection via transgenic intervention at the individual level. Our research results indicated that BmNPV upregulates the expression of BmSTAT-S and BmSTAT-L in Bombyx mori BmE cells and larvae. Furthermore, BmE cells with interference of BmSTAT-S or BmSTAT-L displayed significantly higher expression levels of the viral gene GP41 and increased viral fluorescence compared to the control group after 48 h of infection with BmNPV. Then, we constructed transgenic silkworms with genetic interference, and the results showed that both the transgenic silkworms with systemic interference and midgut-specific interference of the two splice forms of BmSTAT exhibited significantly reduced survival rates and increased viral replication numbers after infection with BmNPV. The above results indicated that the BmSTAT gene is involved in the immune response of Bombyx mori to BmNPV and these findings lay the foundation for further research on the mechanism of JAK/STAT signaling pathway involvement in BmNPV infection.

The domesticated silkworm (Bombyx mori) was domesticated from wild silkworm (Bombyx mandarina) more than 5,000 years ago. During domestication, body color between B. mandarina and B. mori changed dramatically. However, the molecular mechanism of the silkworm body color transition is not known. In the present study, we examined within- and between-species nucleotide diversity for eight silkworm melanin synthesis pathway genes, which play a key role in cuticular pigmentation of insects. Our results showed that the genetic diversity of B. mori was significantly lower than that of B. mandarina and 40.7% of the genetic diversity of wild silkworm was lost in domesticated silkworm. We also examined whether position effect exists among melanin synthesis pathway genes in B. mandarina and B. mori. We found that the upstream genes have significantly lower levels of genetic diversity than the downstream genes, supporting a functional constraint hypothesis (FCH) of metabolic pathway, that is, upstream enzymes are under greater selective constraint than downstream enzymes because upstream enzymes participate in biosynthesis of a number of metabolites. We also investigated whether some of the melanin synthesis pathway genes experienced selection during domestication. Neutrality test, coalescent simulation, as well as network and phylogenetic analyses showed that tyrosine hydroxylase (TH) gene was a domestication locus. Sequence analysis further suggested that a putative expression enhancer (Abd-B-binding site) in the intron of TH gene might be disrupted during domestication. TH is the rate-limiting enzyme of melanin synthesis pathway in insects. Real-time polymerase chain reaction assay did show that the relative expression levels of TH gene in B. mori were significantly lower than that in B. mandarina at three different developmental stages, which is consistent with light body color of domesticated silkworm relative to wild silkworm. Therefore, we speculated that expression change of TH gene may contribute to the body color transition from B. mandarina to B. mori. Our results emphasize the exceptional role of gene expression regulation in morphological transition of domesticated animals.

Introduction.- The White spruce genome.- Targeted and whole genome re-sequencing for population and evolutionary genomic inferences in Norway spruce.- Transposable elements in Spruce.- An intact, but dormant LTR retrotransposon defines a moderately-sized family in white spruce (Picea glauca).- The DNA methylome of Norway spruce.- Epigenomics in Norway spruce.- Comparative Genomics of Spruce.- Comparative genomics with the comparisons of BACs between loblolly pine and white spruce.- Wood formation and genomic selection in white spruce, data imputation strategies.- Ecotypes in Norway spruce (Genomics of adaptation to drought across the Central European range of Norway Local adaptation in the interior spruce hybrid complex.- Genomics of secondary metabolism - Terpenoids in spruce (Evolutionary dynamics of TPS gene family in Chinese spruce).- Lignin biosynthesis in Norway spruce as an anti-fungal defense.- Future prospects.

The hycu-ep32 gene of Hyphantria cunea NPV can inhibit Bombyx mori nucleopolyhedrovirus (BmNPV) multiplication in co-infected cells, but it is not known whether the overexpression of the hycu-ep32 gene has an antiviral effect in the silkworm, Bombyx mori. Thus, we constructed four transgenic vectors, which were under the control of the 39 K promoter of BmNPV (39 KP), Bombyx mori A4 promoter (A4P), hr3 enhancer of BmNPV combined with 39 KP, and hr3 combined with A4P. Transgenic lines were created via embryo microinjection using practical diapause silkworm. qPCR revealed that the expression level of hycu-ep32 could be induced effectively after BmNPV infection in transgenic lines where hycu-ep32 was controlled by hr3 combined with 39 KP (i.e., HEKG). After oral inoculation of BmNPV with 3 × 105 occlusion bodies per third instar, the mortality with HEKG-B was approximately 30% lower compared with the non-transgenic line. The economic characteristics of the transgenic lines remained unchanged. These results suggest that overexpression of an exogenous antiviral gene controlled by an inducible promoter and enhancer is a feasible method for breeding silkworms with a high antiviral capacity.

Bombyx mori, the domesticated silkworm, is of great importance as a silk producer and as a powerful experimental model for the basic and applied research. Similar to other animals, abundant microorganisms live inside the silkworm gut; however, surprisingly, the microbiota of this model insect has not been well characterized to date. Here, we comprehensively characterized the gut microbiota of the domesticated silkworm and its wild relatives. Comparative analyses with the mulberry-feeding moths Acronicta major and Diaphania pyloalis revealed a highly diverse but distinctive silkworm gut microbiota despite thousands of years of domestication, and stage-specific signatures in both total (DNA-based) and active (RNA-based) bacterial populations, dominated by the phyla Proteobacteria, Firmicutes, Actinobacteria, and Bacteroidetes. Most fungal sequences were assigned to the phyla Ascomycota and Basidiomycota. Environmental factors, including diet and human manipulation (egg production), likely influence the silkworm gut composition. Despite a lack of spatial variation along the gut, microbial community shifts were apparent between early instars and late instars, in concert with host developmental changes. Our results demonstrate that the gut microbiota of silkworms assembles into increasingly identical community throughout development, which differs greatly from those of other mulberry-feeding lepidopterans from the same niche, highlighting host-specific effects on microbial associations and the potential roles these communities play in host biology.

Comparative population genomics: power and principles for the inference of functionality

Domestic species provides a powerful model for examining genetic mechanisms in the evolution of yield traits. The domestic silkworm (Bombyx mori) is an important livestock species in sericulture. While the mechanisms controlling cocoon yield are largely unknown. Here, using B. mori and its wild relative B. mandarina as intercross parents, 100 BC1 individuals were sequenced by restriction site-associated DNA sequencing (RAD-Seq). The linkage map contained 9,632 markers was constructed. We performed high-resolution quantitative trait locus (QTL) mapping for four cocoon yield traits. A total of 11 QTLs were identified, including one yield-enhancing QTL from wild silkworm. By integrating population genomics and transcriptomic analysis with QTLs, some favourable genes were revealed, including 14 domestication-related genes and 71 differentially expressed genes (DEGs) in the fifth-instar larval silk gland transcriptome between B. mori and B. mandarina. The relationships between the expression of two important candidate genes (KWMTBOMO04917 and KWMTBOMO12906) and cocoon yield were supported by quantitative real-time PCR (qPCR). Our results provide some new insights into the molecular mechanisms of complex yield traits in silkworm. The combined method might be an efficient approach for identifying putative causal genes in domestic livestock and wild relatives.

Population genetics is a part of evolutionary studies. Now with human genome sequences, population genomics emerged, starting from the analysis of multiple human mitochondrial DNA genome sequences. It was extended to nuclear DNA of human individuals, and genome-wide SNP data comparison is now flourishing, rapidly followed by comparisons of personal genomes. As the genome sequencing cost is becoming drastically reduced, population genomics will definitely expand to many other organisms. We discuss both methods and examples of population genomics in this chapter focusing human genomes. However, contents of this chapter can be applicable to population genomics of any diploid organisms.

Comparative genomic and phylogenetic analysis of a toxigenic clinical isolate of Corynebacterium diphtheriae strain B-D-16-78 from Malaysia