Fatty Acid Desaturase Bmdesat5, Suppressed in the Salivary Glands by Domestication, is Involved in Regulation of Food Intake in Silkworms.

Understanding the evolutionary genetics of food intake regulation in domesticated animals has relevance to evolutionary biology, animal improvement, and obesity treatment. Here, we observed that the fatty acid desaturase gene (Bmdesat5), which regulates food intake, is suppressed in domesticated silkworms, but expressed in the salivary glands of the wild silkworm Bombyx mandarina. The content of its catalytic product, cis-vaccenic acid, was related to the expression levels of Bmdesat5 in the salivary glands of domesticated and wild silkworm strains. These two strains also showed significant differences in food intake. Using orally administering cis-vaccenic acid and transgenic-mediated overexpression, we verified that cis-vaccenic acid functions as a satiation signal, regulating food intake and growth in silkworms. Selection analysis showed that Bmdesat5 experienced selection, especially in the potential promoter, 5'-untranslated, and intron regions. This study highlights the importance of the decrement of satiety in silkworm domestication and provides new insights into the potential involvement of salivary glands in the regulation of satiety in animals, by acting as a supplement to gut-brain nutrient signaling.

Recessive embryonic lethal mutations uncovered in heterozygous condition in silkworm semiconsomic strains.

In this study, we found two embryonic lethal mutations, t04 lethal (l-t04) and m04 lethal (l-m04), in semiconsomic strains T04 and M04, respectively. In these semiconsomic strains, the entire diploid genome, except for one chromosome 4 of the wild silkworm Bombyx mandarina, is substituted with chromosomes of the domesticated silkworm B. mori, and l-t04 and l-m04 mutations are located on B. mandarina-derived chromosome 4. To clarify the cause of the lethalities and the genes responsible for these mutations, positional cloning and CRISPR/Cas9 mediated knockout screening were performed. Finally, genetic complementation tests identified the mutations responsible for the l-t04 and l-m04 as the Bombyx homolog of imaginal discs arrested (Bmida) and TATA box binding protein-associated factor 5 (BmTaf5), respectively. Lethal stages of each knockout mutant indicated the importance of these genes in B. mori late embryogenesis. The lethal mutations responsible for l-t04 and l-m04 were not found in parental strains or wild B. mandarina collected from 39 distinct locations in Japan, indicating that both mutations were independently introduced during or after the development of the semiconsomic strains. We conclude that the recessive embryonic lethality in the T04 and M04 strains is due to deleterious mutations produced in B. mandarina-derived chromosome 4.

Allele-specific knockouts reveal a role for apontic-like in the evolutionary loss of larval melanin pigmentation in the domesticated silkworm, Bombyx mori.

The domesticated silkworm, Bombyx mori, and its wild progenitor, B. mandarina, are extensively studied as a model case of the evolutionary process of domestication. A conspicuous difference between these species is the dramatic reduction in melanin pigmentation in both larval and adult B. mori. Here we evaluate the efficiency of CRISPR/Cas9-targeted knockouts of pigment-related genes as a tool to understand their potential contributions to domestication-associated melanin pigmentation loss in B. mori. To demonstrate the efficacy of targeted knockouts in B. mandarina, we generated a homozygous CRISPR/Cas9-targeted knockout of yellow-y. In yellow-y knockout mutants, black body colour became lighter throughout the larval, pupal and adult stages, confirming a role for this gene in melanin pigment formation. Further, we performed allele-specific CRISPR/Cas9-targeted knockouts of the pigment-related transcription factor, apontic-like (apt-like) in B. mori × B. mandarina F1 hybrid individuals which exhibit B. mandarina-like larval pigmentation. Knockout of the B. mandarina allele of apt-like in F1 embryos results in white patches on the dorsal integument of larvae, whereas corresponding knockouts of the B. mori allele consistently exhibit normal F1 larval pigmentation. These results demonstrate a contribution of apt-like to the evolution of reduced melanin pigmentation in B. mori. Together, our results demonstrate the feasibility of CRISPR/Cas9-targeted knockouts as a tool for understanding the genetic basis of traits associated with B. mori domestication.

Phylogeographic Relationships among Bombyx mandarina (Lepidoptera: Bombycidae) Populations and Their Relationships to B. mori Inferred from Mitochondrial Genomes.

We report 37 mitochondrial genome (mitogenome) sequences of Bombyx mori strains (Lepidoptera: Bombycidae) and four of B. mandarina individuals, each preserved and collected, respectively, in South Korea. These mitogenome sequences combined with 45 public data showed a substantial genetic reduction in B. mori strains compared to the presumed ancestor B. mandarina, with the highest diversity detected in the Chinese origin B. mori. Chinese B. mandarina were divided into northern and southern groups, concordant to the Qinling-Huaihe line, and the northern group was placed as an immediate progenitor of monophyletic B. mori strains in phylogenetic analyses, as has previously been detected. However, one individual that was in close proximity to the south Qinling-Huaihe line was exceptional, belonging to the northern group. The enigmatic South Korean population of B. mandarina, which has often been regarded as a closer genetic group to Japan, was most similar to the northern Chinese group, evidencing substantial gene flow between the two regions. Although a substantial genetic divergence is present between B. mandarina in southern China and Japan, a highly supported sister relationship between the two regional populations may suggest the potential origin of Japanese B. mandarina from southern China instead of the Korean peninsula.

Bombyx mori nucleopolyhedrovirus ptp and egt genes are dispensable for triggering enhanced locomotory activity and climbing behavior in Bombyx mandarina larvae.

Baculoviruses are classic pathogens that alter host behavior to enhance their dispersal and transmission. While viral protein tyrosine phosphatase (ptp) has been considered as a critical factor for inducing enhanced locomotory activity, preceding investigations have reported that viral ecdysteroid UDP-glucosyltransferase (egt) contributes to triggering climbing behavior in some virus and host species. Here we found that both egt and ptp were dispensable for these abnormal behaviors in Bombyx mandarina larvae induced by Bombyx mori nucleopolyhedrovirus, thus implying that there is an unknown core mechanism of baculovirus-induced alteration of host behaviors.

Cytochrome P450 monooxygenase genes in the wild silkworm, Bombyx mandarina.

Wild (Bombyx mandarina) and domestic silkworms (B. mori) are good models for investigating insect domestication, as 5000 years of artificial breeding and selection have resulted in significant differences between B. mandarina and B. mori. In this study, we improved the genome assemblies to the chromosome level and updated the protein-coding gene annotations for B. mandarina. Based on this updated genome, we identified 68 cytochrome P450 genes in B. mandarina. The cytochrome P450 repository in B. mandarina is smaller than in B. mori. Certain currently unknown key genes, rather than gene number, are critical for insecticide resistance in B. mandarina, which shows greater resistance to insecticides than B. mori. Based on the physical maps of B. mandarina, we located 66 cytochrome P450s on 18 different chromosomes, and 27 of the cytochrome P450 genes were concentrated into seven clusters. KEGG enrichment analysis of the P450 genes revealed the involvement of cytochrome P450 genes in hormone biosynthesis. Analyses of the silk gland transcriptome identified candidate cytochrome P450 genes (CYP306A) involved in ecdysteroidogenesis and insecticide metabolism in B. mandarina.

[Comparison of chemical constituents of wild silkworm cocoon and domestic silkworm cocoon by UHPLC-MS technology].

Identifying and comparing the chemical constituents of wild silkworm cocoon and silkworm cocoon is of great significance for understanding the domestication of silkworm. In this study, we used high temperature and high pressure and methanol-water system to extract cocoon chemical constituents. We used UHPLC-MS to identify and compare cocoon chemical constituents of wild silkworm and domestic silkworm Dazao and Haoyue strains. The cocoon metabolic fingerprints of wild silkworm and domestic silkworm Dazao and Haoyue strains were obtained by using the UHPLC-MS in the positive ion mode and negative ion mode. By annotation, we found that cocoon chemical compounds with high abundances contained amino acids, flavonoids, alkaloids, terpenes, organic acids, and lignans. PLS-DA showed that the cocoon components were significantly different among the wild silkworm and two domestic silkworm strains Dazao and Haoyue. Proline, leucine/isoleucine and phenylalanine showed significantly higher abundances in the cocoon of domestic silkworm Dazao strain than in those of wild silkworm and domestic silkworm Haoyue strain. The flavonoid secondary metabolites are abundant in the Dazao cocoon, including quercetin, isoquercetin, quercetin 3-O-sophoroside, quercetin-3-O-α-L-rhamnoside, quercetin-3-O- rutinoside, and kaempferol. The other secondary metabolites, alkaloids, terpenes and lignans, showed higher abundances in the wild silkworm cocoon than in the domestic silkworm cocoon, including neurine, candicine, pilocarpidine, artemisiifolin, eupassopin, and eudesobovatol. By exposing cocoons to UV light and observing the green fluorescence of flavonoids, we found that Dazao cocoon had the most flavonoids, and Haoyue cocoon had least flavonoids and wild silkworm cocoon had mediate flavonoids. Alkaloids and organic acids are good anti-insect and antimicrobial agents, which have high abundance in the wild silkworm cocoon and could enhance the defense ability of wild silkworm cocoon. Flavonoids are abundant in the cocoon of domestic silkworm Dazao strain, which the main factors are leading to the yellow-green cocoon of Dazao.

Comparative mitochondrial genomes provide new insights into the true wild progenitor and origin of domestic silkworm Bombyx mori.

The domestication of domestic silkworm Bombyx mori, the only truly domesticated insect, is a distinctive event in agricultural history. The domestication and origin of domestic silkworm remains unclear, although it has connected with human for ~5500 years. In the present study, we would like to highlight our evidence from whole mitochondrial genome for the presence of two genetically distinctive subtypes in Chinese B. mandarina populations, corresponding to northern Chinese B. mandarina and southern Chinese B. mandarina, respectively. The mitochondrial genomes and mitochondrial phylogenetic tree provide a solid molecular evidence that the true wild ancestor of domestic silkworm is northern Chinese B. mandarina, rather than southern Chinese B. mandarina, thus implying that the early domestication event may have occurred in northern China. Our finding provides new insights into the origin and evolution of domestic silkworm.

Morphological and electrophysiological differences in tarsal chemosensilla between the wild silkmoth Bombyx mandarina and the domesticated species Bombyx mori.

Gustatory and olfactory senses of phytophagous insects play important roles in the recognition of host plants. In the domestic silkmoth Bombyx mori and its wild species Bombyx mandarina, the morphologies and responses of adult olfactory organs (antennae) have been intensely investigated. However, little is known about these features of adult gustatory organs and the influence of domestication on the gustatory sense. Here we revealed that both species have two types of sensilla (thick [T] and slim [S] types) on the fifth tarsomeres of the adult legs. In both species, females have 3.6-6.9 times more T-sensilla than males. Therefore, T-sensilla seem to play more important roles in females than in males. Moreover, gustatory cells of T-sensilla of B. mandarina females responded intensely to mulberry leaf extract in electrophysiological experiments, while T-sensilla of B. mori females (N4 strain) hardly responded to mulberry leaf extract. These results suggest that T-sensilla of B. mandarina females are involved in the recognition of oviposition sites. We also observed that, in three B. mori strains (N4, p50T, and Kinshu × Showa), the densities of sensilla on the fifth tarsomeres were much lower than in B. mandarina. These results indicate that domestication has influenced the tarsal gustatory system of B. mori.

Comparative Proteomic Analysis of Posterior Silk Glands of Wild and Domesticated Silkworms Reveals Functional Evolution during Domestication.

The wild silkworm Bombyx mandarina was domesticated to produce silk in China approximately 5000 years ago. Silk production is greatly improved in the domesticated silkworm B. mori, but the molecular basis of the functional evolution of silk gland remains elusive. We performed shotgun proteomics with label-free quantification analysis and identified 1012 and 822 proteins from the posterior silk glands (PSGs) of wild silkworms on the third and fifth days of the fifth instar, respectively, with 128 of these differentially expressed. Bioinformatics analysis revealed that, with the development of the PSG, the up-regulated proteins were mainly involved in the ribosome pathway, similar to what we previously reported for B. mori. Additionally, we screened 50 proteins with differential expression between wild and domesticated silkworms that might be involved in domestication at the two stages. Interestingly, the up-regulated proteins in domesticated compared to wild silkworms were enriched in the ribosome pathway, which is closely related to cell size and translation capacity. Together, these results suggest that functional evolution of the PSG during domestication was driven by reinforcing the advantageous pathways to increase the synthesis efficiency of silk proteins in each cell and thereby improve silk yield.

A novel splice variant of the decapentaplegic (dpp) gene in the wild silkworm, Bombyx mandarina.

Decapentaplegic (dpp) is a member of the transforming growth factor-ß superfamily. Although the dpp gene and related pathways are known to play important roles in insect development, few studies have examined its function in Bombyx mori and Bombyx mandarina. To date, there have been no previous reports on novel splice variants of dpp in silkworm. In the present study, we conducted RT-PCR to examine dpp expression in the mid-gut tissue of B. mandarina and discovered a novel dpp isoform. The isoform sequence was confirmed using sequencing analysis and found to have 333 bp deletion compared to full-length cDNA encoding dpp. The deleted sequence encodes a region of the latency associated peptide (LAP) region of transforming growth factor-ß (TGF-ß), which may affect the activity and specificity of TGF-ß. Using variant calling analyses, we detected 7 candidate single nucleotide variants (SNVs) for different alternative splicing in dpp. This is the first report of a novel splice variant of the dpp gene in B. mandarina and these results provide insight about the domestication process and distinct phenotypic traits of B. mori and B. mandarina.

IDENTIFICATION AND EXPRESSION ANALYSIS OF VITELLOGENIN RECEPTOR FROM THE WILD SILKWORM, Bombyx mandarina.

The vitellogenin receptor (VgR) plays a key role on embryonic development in oviparous animals. Here, we cloned a VgR gene, which was identified from the wild silkworm Bombyx mandarina (BmaVgR) using reverse transcriptase polymerase chain reaction (RT-PCR) and rapid amplification of cDNA ends (RACE). Sequence analysis revealed that BmaVgR is 5,861 bp long with an open reading frame encoded by 1,811 amino acid residues. The predicted amino acid sequence has 99.7 and 98.2% identity with the VgRs of Actias selene and Bombyx mori, respectively. The class B domain sequence of BmaVgR was cloned and expressed in Escherichia coli, and purified by a Ni-NTA column. Polyclonal antibodies were produced against the purified recombinant protein, and titer of the antibody was about 1:12,800 measured by enzyme-linked immunosorbent assay (ELISA). Western blot and RT-qPCR showed that BmaVgR was expressed in the ovary and fat body of female larvae and the ovary of moth, and the expression level was highest at the third day and then declined from third day to seventh in fat body of pupa. After knockdown of the BmaVgR gene through RNA interference (RNAi), other three BmaVgR-related genes (Vg, egg-specific protein, and low molecular weight lipoprotein LP gene) were all downregulated significantly.

Anatomical and functional analysis of domestication effects on the olfactory system of the silkmoth Bombyx mori.

The silkmoth Bombyx mori is the main producer of silk worldwide and has furthermore become a model organism in biological research, especially concerning chemical communication. However, the impact domestication might have had on the silkmoth's olfactory sense has not yet been investigated. Here, we show that the pheromone detection system in B. mori males when compared with their wild ancestors Bombyx mandarina seems to have been preserved, while the perception of environmental odorants in both sexes of domesticated silkmoths has been degraded. In females, this physiological impairment was mirrored by a clear reduction in olfactory sensillum numbers. Neurophysiological experiments with hybrids between wild and domesticated silkmoths suggest that the female W sex chromosome, so far known to have the sole function of determining femaleness, might be involved in the detection of environmental odorants. Moreover, the coding of odorants in the brain, which is usually similar among closely related moths, differs strikingly between B. mori and B. mandarina females. These results indicate that domestication has had a strong impact on odour detection and processing in the olfactory model species B. mori.

Mitochondrial genome nucleotide substitution pattern between domesticated silkmoth, Bombyx mori, and its wild ancestors, Chinese Bombyx mandarina and Japanese Bombyx mandarina.

Bombyx mori and Bombyx mandarina are morphologically and physiologically similar. In this study, we compared the nucleotide variations in the complete mitochondrial (mt) genomes between the domesticated silkmoth, B. mori, and its wild ancestors, Chinese B. mandarina (ChBm) and Japanese B. mandarina (JaBm). The sequence divergence and transition mutation ratio between B. mori and ChBm are significantly smaller than those observed between B. mori and JaBm. The preference of transition by DNA strands between B. mori and ChBm is consistent with that between B. mori and JaBm, however, the regional variation in nucleotide substitution rate shows a different feature. These results suggest that the ChBm mt genome is not undergoing the same evolutionary process as JaBm, providing evidence for selection on mtDNA. Moreover, investigation of the nucleotide sequence divergence in the A+T-rich region of Bombyx mt genomes also provides evidence for the assumption that the A+T-rich region might not be the fastest evolving region of the mtDNA of insects.

Mitochondrial genome nucleotide substitution pattern between domesticated silkmoth, Bombyx mori, and its wild ancestors, Chinese Bombyx mandarina and Japanese Bombyx mandarina

Bombyx mori and Bombyx mandarina are morphologically and physiologically similar. In this study, we compared the nucleotide variations in the complete mitochondrial (mt) genomes between the domesticated silkmoth, B. mori, and its wild ancestors, Chinese B. mandarina (ChBm) and Japanese B. mandarina (JaBm). The sequence divergence and transition mutation ratio between B. mori and ChBm are significantly smaller than those observed between B. mori and JaBm. The preference of transition by DNA strands between B. mori and ChBm is consistent with that between B. mori and JaBm, however, the regional variation in nucleotide substitution rate shows a different feature. These results suggest that the ChBm mt genome is not undergoing the same evolutionary process as JaBm, providing evidence for selection on mtDNA. Moreover, investigation of the nucleotide sequence divergence in the A+T-rich region of Bombyx mt genomes also provides evidence for the assumption that the A+T-rich region might not be the fastest evolving region of the mtDNA of insects.