SilkMeta: a comprehensive platform for sharing and exploiting pan-genomic and multi-omic silkworm data.

The silkworm Bombyx mori is a domesticated insect that serves as an animal model for research and agriculture. The silkworm super-pan-genome dataset, which we published last year, is a unique resource for the study of global genomic diversity and phenotype-genotype association. Here we present SilkMeta (http://silkmeta.org.cn), a comprehensive database covering the available silkworm pan-genome and multi-omics data. The database contains 1082 short-read genomes, 546 long-read assembled genomes, 1168 transcriptomes, 294 phenotype characterizations (phenome), tens of millions of variations (variome), 7253 long non-coding RNAs (lncRNAs), 18 717 full length transcripts and a set of population statistics. We have compiled publications on functional genomics research and genetic stock deciphering (mutant map). A range of bioinformatics tools is also provided for data visualization and retrieval. The large batch of omics data and tools were integrated in twelve functional modules that provide useful strategies and data for comparative and functional genomics research. The interactive bioinformatics platform SilkMeta will benefit not only the silkworm but also the insect biology communities.

Deciphering the Genetic Basis of Silkworm Cocoon Colors Provides New Insights into Biological Coloration and Phenotypic Diversification.

The genetic basis of phenotypic variation is a long-standing concern of evolutionary biology. Coloration has proven to be a visual, easily quantifiable, and highly tractable system for genetic analysis and is an ever-evolving focus of biological research. Compared with the homogenized brown-yellow cocoons of wild silkworms, the cocoons of domestic silkworms are spectacularly diverse in color, such as white, green, and yellow-red; this provides an outstanding model for exploring the phenotypic diversification and biological coloration. Herein, the molecular mechanism underlying silkworm green cocoon formation was investigated, which was not fully understood. We demonstrated that five of the seven members of a sugar transporter gene cluster were specifically duplicated in the Bombycidae and evolved new spatial expression patterns predominantly expressed in silk glands, accompanying complementary temporal expression; they synergistically facilitate the uptake of flavonoids, thus determining the green cocoon. Subsequently, polymorphic cocoon coloring landscape involving multiple loci and the evolution of cocoon color from wild to domestic silkworms were analyzed based on the pan-genome sequencing data. It was found that cocoon coloration involved epistatic interaction between loci; all the identified cocoon color-related loci existed in wild silkworms; the genetic segregation, recombination, and variation of these loci shaped the multicolored cocoons of domestic silkworms. This study revealed a new mechanism for flavonoids-based biological coloration that highlights the crucial role of gene duplication followed by functional diversification in acquiring new genetic functions; furthermore, the results in this work provide insight into phenotypic innovation during domestication.

Unusual tertiary pairs in eukaryotic tRNAAla.

tRNA molecules have well-defined sequence conservations that reflect the conserved tertiary pairs maintaining their architecture and functions during the translation processes. An analysis of aligned tRNA sequences present in the GtRNAdb database (the Lowe Laboratory, University of California, Santa Cruz) led to surprising conservations on some cytosolic tRNAs specific for alanine compared to other tRNA species, including tRNAs specific for glycine. First, besides the well-known G3oU70 base pair in the amino acid stem, there is the frequent occurrence of a second wobble pair at G30oU40, a pair generally observed as a Watson-Crick pair throughout phylogeny. Second, the tertiary pair R15/Y48 occurs as a purine-purine R15/A48 pair. Finally, the conserved T54/A58 pair maintaining the fold of the T-loop is observed as a purine-purine A54/A58 pair. The R15/A48 and A54/A58 pairs always occur together. The G30oU40 pair occurs alone or together with these other two pairs. The pairing variations are observed to a variable extent depending on phylogeny. Among eukaryotes, insects display all variations simultaneously, whereas mammals present either the G30oU40 pair or both R15/A48 and A54/A58. tRNAs with the anticodon 34A(I)GC36 are the most prone to display all those pair variations in mammals and insects. tRNAs with anticodon Y34GC36 have preferentially G30oU40 only. These unusual pairs are not observed in bacterial, nor archaeal, tRNAs, probably because of the avoidance of A34-containing anticodons in four-codon boxes. Among eukaryotes, these unusual pairing features were not observed in fungi and nematodes. These unusual structural features may affect, besides aminoacylation, transcription rates (e.g., 54/58) or ribosomal translocation (30/40).

Identification, expression, and artificial selection of silkworm epigenetic modification enzymes.

BACKGROUND: Understanding the genetic basis of phenotype variations during domestication and breeding is of great interest. Epigenetics and epigenetic modification enzymes (EMEs) may play a role in phenotypic variations; however, no comprehensive study has been performed to date. Domesticated silkworm (Bombyx mori) may be utilized as a model in determining how EMEs influence domestication traits. RESULTS: We identified 44 EMEs in the genome of silkworm (Bombyx mori) using homology searching. Phylogenetic analysis showed that genes in a subfamily among different animals were well clustered, and the expression pattern of EMEs is constant among Bombyx mori, Drosophila melanogaster, and Mus musculus. These are most highly expressed in brain, early embryo, and internal genitalia. By gene-related selective sweeping, we identified five BmEMEs under artificial selection during the domestication and breeding of silkworm. Among these selected genes, BmSuv4-20 and BmDNMT2 harbor selective mutations in their upstream regions that alter transcription factor-binding sites. Furthermore, these two genes are expressed higher in the testis and ovary of domesticated silkworm compared to wild silkworms, and correlations between their expression pattern and meiosis of the sperm and ova were observed. CONCLUSIONS: The domestication of silkworm has induced artificial selection on epigenetic modification markers that may have led to phenotypic changes during domestication. We present a novel perspective to understand the genetic basis underlying animal domestication and breeding.

Bmmp influences wing morphology by regulating anterior-posterior and proximal-distal axes development.

Insect wings are subject to strong selective pressure, resulting in the evolution of remarkably diverse wing morphologies that largely determine flight capacity. However, the genetic basis and regulatory mechanisms underlying wing size and shape development are not well understood. The silkworm Bombyx mori micropterous (mp) mutant exhibits shortened wing length and enlarged vein spacings, albeit without changes in total wing area. Thus, the mp mutant comprises a valuable genetic resource for studying wing development. In this study, we used molecular mapping to identify the gene responsible for the mp phenotype and designated it Bmmp. Phenotype-causing mutations were identified as indels and single nucleotide polymorphisms in noncoding regions. These mutations resulted in decreased Bmmp messenger RNA levels and changes in transcript isoform composition. Bmmp null mutants were generated by clustered regularly interspaced short palindromic repeats (CRISPR) / CRISPR-associated protein 9 and exhibited changed wing shape, similar to mp mutants, and significantly smaller total wing area. By examining the expression of genes critical to wing development in wildtype and Bmmp null mutants, we found that Bmmp exerts its function by coordinately modulating anterior-posterior and proximal-distal axes development. We also studied a Drosophila mp mutant and found that Bmmp is functionally conserved in Drosophila. The Drosophila mp mutant strain exhibits curly wings of reduced size and a complete loss of flight capacity. Our results increase our understanding of the mechanisms underpinning insect wing development and reveal potential targets for pest control.

High-resolution silkworm pan-genome provides genetic insights into artificial selection and ecological adaptation.

The silkworm Bombyx mori is 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.

Flight Muscle and Wing Mechanical Properties are Involved in Flightlessness of the Domestic Silkmoth, Bombyx mori.

Flight loss has occurred in many winged insect taxa. The flightless silkmoth Bombyx mori, is domesticated from the wild silkmoth, Bombyx mandarina, which can fly. In this paper, we studied morphological characteristics attributed to flightlessness in silkmoths. Three domestic flightless B. mori strains and one B. mandarina population were used to compare morphological components of the flight apparatus, including wing characteristics (shape, forewing area, loading, and stiffness), flight muscle (weight, ratio, and microscopic detail) and body mass. Compared with B. mandarina, B. mori strains have a larger body, greater wing loading, more flexible wings and a lower flight muscle ratio. The arrangement in microscopy of dorsal longitudinal flight muscles (DLFMs) of B. mori was irregular. Comparative analysis of the sexes suggests that degeneration of flight muscles and reduction of wing mechanical properties (stiffness) are associated with silkmoth flightlessness. The findings provide important clues for further research of the molecular mechanisms of B. mori flight loss.

Effects of P27/Bmdacapo, in the CIP/KIP family, on cell proliferation, growth and development in the silkworm (Bombyx mori).

We investigated changes in expression of the CIP/KIP family-related genes and the cycle-dependent factors Pcna, Cdk4 and Cdk2 during the growth and development of mice, Drosophila and silkworms. When the organism was in a period of rapid development, the related genes of the CIP/KIP family had low expression level and the cell cycle-dependent genes were highly expressed. In mammals, the CIP/KIP family includes three genes, p21, p27/Dacapo and p57. However, only one gene, P27/Dacapo, exists in the CIP/KIP family in silkworm and the orthologous gene in the silkworm is named Bmdacapo. Down-regulation of Bmdacapo in silkworm embryos caused overdevelopment of the embryos and indicated that Bmdacapo can inhibit silkworm growth and development. Up-regulation of Bmdacapo in silkworm cells inhibited cell proliferation, whereas down-regulation of Bmdacapo promoted cell proliferation. In order to explore the mechanism of Bmdacapo regulated silkworm development and cell proliferation, the effect of Bmdacapo on cell cycle changes was examined. The results demonstrate that Bmdacapo was able to induce G1/S phase arrest in the cell cycle. In silkworm cells, Bmdacapo inhibits the expression of Pcna, CDK4 and CDK2, which affects the cell cycle and ultimately inhibits cell proliferation. This regulatory mechanism is particularly different from mammals.

Structure analysis of the spinneret from Bombyx mori and its influence on silk qualities.

Silk is an excellent natural fiber, which has been widely used in versatile fields. Silk spinning is a complex process involving the larval spinneret. The spinneret is essential for silk spinning, but the sectional morphology of the spinneret that determines the silk monofilament, the muscular activities around the silk press as well as the relationships between the spinneret and the properties of the resulting silk remain poorly understood. We studied these factors by dissecting the spinneret and analyzing silk from different Bombyx mori strains. The sectional morphology of silk monofilament was found to be largely determined by the spinneret, especially by the silk press. Moreover, contractile activity of the muscles around the silk press is high, and the contraction frequency of the muscles was estimated to range from 11.42 to 50 HZ. A comparison of the fibroin filaments before they entered the common tube indicated that the spinneret determines both silk shape and silk size. This study provides insight into the silk spinning process, which may help develop bionic spinning in further studies and also provides a rationale to study the effect of the spinneret on silk fineness at the molecular level.