Identification of Genes that Control Silk Yield by RNA Sequencing Analysis of Silkworm (Bombyx mori) Strains of Variable Silk Yield

Comparative transcriptome and proteome reveal synergistic functions of differentially expressed genes and proteins implicated in an over-dominant silkworm heterosis of increased silk yield

We previously observed an over-dominant silkworm heterosis of increased yield in a cross of Bombyx mori nuclear polyhydrosis virus-resistant strain NB with a susceptible strain 306. In the present study, we found that heterosis also exists in crosses of NB with other susceptible strains, indicating it is a more general phenomenon. We performed comparative transcriptome and proteome and identified 1624 differentially expressed genes (DEGs) and 298 differentially expressed proteins (DEPs) in silk glands between parents and F1 hybrids, of which 24 DEGs/DEPs showed consistent expression at mRNA and protein levels revealed by Venn joint analysis. Their expressions are completely non-additive, mainly transgressive and under low-parent, suggesting recombination of parental genomes may be the major genetic mechanism for the heterosis. GO and KEGG analyses revealed that they may function in generally similar but distinctive aspects of metabolisms and processes with signal transduction and translation being most affected. Notably, they may not only up-regulate biosynthesis and transport of silk proteins but also down-regulate other unrelated processes, synergistically and globally remodelling the silk gland to increase yield and cause the heterosis. Our findings contribute insights into the understanding of silkworm heterosis and silk gland development and provide targets for transgenic manipulation to further increase the silk yield.

Bombyx mori Vps13d is a key gene affecting silk yield

Bombyx mori is an important economic insect, its economic value mainly reflected in the silk yield. The major functional genes affecting the silk yield of B. mori have not been determined yet. Bombyx mori vacuolar protein sorting-associated protein 13d (BmVps13d) has been identified, but its function is not reported. In this study, BmVps13d protein shared 30.84% and 34.35% identity with that of in Drosophila melanogaster and Homo. sapiens, respectively. The expressions of BmVps13d were significantly higher in the midgut and silk gland of JS (high silk yield) than in that of L10 (low silk yield). An insertion of 9 bp nucleotides and two deficiencies of adenine ribonucleotides in the putative promoter region of BmVps13d gene in L10 resulted in the decline of promoter activity was confirmed using dual luciferase assay. Finally, the functions of BmVps13d in B. mori were studied using the CRISPR/Cas9 system, and the mutation of BmVps13d resulted in a 24.7% decline in weight of larvae, as well as a 27.1% (female) decline and a 11.8% (male) decline in the silk yield. This study provides a foundation for studying the molecular mechanism of silk yield and breeding the silkworm with high silk yield.

Perfluorooctanoic acid (PFOA) exposure induces splenic atrophy via overactivation of macrophages in male mice

Perfluorooctanoic acid (PFOA), a synthetic and widely used chemical, has aroused wide public concern due to its persistence, bioaccumulation, and potential toxicity. To investigate splenic atrophy induced by PFOA, male mice were exposed to 0, 0.4, 2, or 10 mg/kg/d PFOA for 28 d. Results demonstrated that spleen weight and relative spleen weight (RSW) decreased in the 2 and 10 mg/kg/d PFOA exposure groups. Iron levels in the spleen and serum were also reduced in all PFOA exposure groups. Weighted gene co-expression network analysis (WGCNA) of 7 043 genes highlighted enrichment in cell cycle, autoimmunity, and anemia in the spleen. In addition, changes in the levels of hemoglobin, platelets, bilirubin, and heme oxygenase-1 were consistent with anemia. The ratio of total macrophages to M1 macrophages in the spleen, phagocytic ability of macrophages, and levels of cytokines such as TNF-α, IL-1β, and IL-6 all increased, thus suggesting the occurrence of autoimmune disorder. Therefore, we concluded that overactivation of macrophages may be an important reason for splenic atrophy induced by PFOA exposure.

Gene mapping reveals the association between tyrosine protein kinase Abl1 and the silk yield of Bombyx mori

The Z chromosome of the silkworm contains a major gene that influences silk yield. This major locus on chromosome Z accounts for 35.10% of the phenotypic variance. The location and identification of the gene have been a focus of silkworm genetics research. Unfortunately, identification of this gene has been difficult. We used extreme phenotype subpopulations and selected from a backcross population, BC₁ M, which was obtained using the high-yield strain 872B and the low-yield strain IS-Dazao as parents, for mapping the gene on the chromosome Z. The candidate region was narrowed down to 134 kb at the tip of the chromosome. BmAbl1 in this region correlated with silk gland development by spatiotemporal expression analysis. This gene was differentially expressed in the posterior silk glands of the high- and low-yield strains. In BmAbl1, an insertion-deletion (indel) within the 10th exonic region and an SNP within the 6th intronic region were detected and shown to be associated with cocoon shell weight in 84 Bombyx mori strains with different yields. Nucleotide diversity analysis of BmAbl1 and its 50 kb flanking regions indicated that BmAbl1 has experienced strong artificial selection during silkworm domestication. This study is the first to identify the genes controlling silk yield in the major QTL of the Z chromosome using forward genetics.

Comparative transcriptome analysis of the rice leaf folder (Cnaphalocrocis medinalis) to heat acclimation

Background

The rice leaf folder Cnaphalocrocis medinalis Güenée is a serious insect pest of rice in Asia. This pest occurs in summer, and it is sensitive to high temperature. However, the larvae exhibit heat acclimation/adaptation. To understand the underlying mechanisms, we established a heat-acclimated strain via multigenerational selection at 39 °C. After heat shock at 41 °C for 1 h, the transcriptomes of the heat-acclimated (S-39) and unacclimated (S-27) larvae were sequenced, using the unacclimated larvae without exposure to 41 °C as the control.

Results

Five generations of selection at 39 °C led larvae to acclimate to this heat stress. Exposure to 41 °C induced 1160 differentially expressed genes (DEGs) between the heat-acclimated and unacclimated larvae. Both the heat-acclimated and unacclimated larvae responded to heat stress via upregulating genes related to sensory organ development and structural constituent of eye lens, whereas the unacclimated larvae also upregulated genes related to structural constituent of cuticle. Compared to unacclimated larvae, heat-acclimated larvae downregulated oxidoreductase activity-related genes when encountering heat shock. Both the acclimated and unacclimated larvae adjusted the longevity regulating, protein processing in endoplasmic reticulum, antigen processing and presentation, MAPK and estrogen signaling pathway to responsed to heat stress. Additionally, the unacclimated larvae also adjusted the spliceosome pathway, whereas the heat-acclimated larvae adjusted the biosynthesis of unsaturated fatty acids pathway when encountering heat stress. Although the heat-acclimated and unacclimated larvae upregulated expression of heat shock protein genes under heat stress including HSP70, HSP27 and CRYAB, their biosynthesis, metabolism and detoxification-related genes expressed differentially.

Conclusions

The rice leaf folder larvae could acclimate to a high temperature via multigenerational heat selection. The heat-acclimated larvae induced more DEGs to response to heat shock than the unacclimated larvae. The changes in transcript level of genes were related to heat acclimation of larvae, especially these genes in sensory organ development, structural constituent of eye lens, and oxidoreductase activity. The DEGs between heat-acclimated and unacclimated larvae after heat shock were enriched in the biosynthesis and metabolism pathways. These results are helpful to understand the molecular mechanism underlying heat acclimation of insects.

Effects of Certain cis-Regulatory Elements on Stage-Specific vitellogenin Expression in the Bombyx mori (Lepidoptera: Bombycidae)

Bombyx mori vitellogenin (BmVg) is highly upregulated during pupation, and the 20-hydroxyecdysone and amino acids may regulate stage-specific BmVg expression. However, previous studies showed that other factors may also affect stage-specific BmVg expression. Here, we characterized effective BmVg transcription factors by identifying the corresponding cis-regulatory elements (CREs). We prepared transgenic B. mori, in which DsRed was driven by various lengths of BmVg promoter. qRT-PCR analysis showed that DsRed expression driven by a 1.0-kb BmVg promoter (VgP1.0K) was consistent with endogenous BmVg. VgP1.0K specificity was closer to the endogenous BmVg promoter than that of VgP0.8K. These results suggest that CREs affecting stage-specific BmVg expression were localized to the 1.0-kb BmVg promoter. We investigated the effects of certain CREs that could influence the stage specificity of BmVg promoter on BmVg expression in transgenic B. mori. The relative DsRed expression was significantly reduced in transgenic female B. mori and the peak in DsRed expression was delayed after E-box CRE mutation. These results demonstrate that the E-box element enhanced BmVg expression and also affected stage-specific BmVg expression. Moreover, the relative DsRed expression was significantly increased in transgenic female of B. mori after 3×BD CRE mutation in BmVg promoter. However, the stage specificity of the mutated promoter was consistent with that of the endogenous BmVg promoter. The 3×BD element downregulated BmVg but had no effect on stage-specific BmVg expression. The present study promoted the process of elucidating the regulatory network for stage-specific BmVg expression and furnished a theoretical basis for the application of BmVg promoter.

A Comparison of Co-expression Networks in Silk Gland Reveals the Causes of Silk Yield Increase During Silkworm Domestication

Long-term domestication and selective breeding have increased the silk yield of the domestic silkworm (Bombyx mori) by several times the amount of the silk yield of its wild ancestor (Bombyx mandarina). However, little is known about the molecular mechanisms behind the increase in silk yield during domestication. Based on dynamic patterns of functional divergence in the silk gland between domestic and wild silkworms, we found that at early and intermediate stages of silk gland development, the up-regulated genes of the domestic silkworm were mainly involved in DNA integration, nucleic acid binding, and transporter activity, which are related to the division and growth of cells. This has led to the posterior silk gland (PSG) of the domestic silkworm having significantly more cells (“factories” of fibroin protein synthesis) than that of the wild silkworm. At the late stage of silk gland development, the up-regulated genes in the domestic silkworm was enriched in protein processing and ribosome pathways, suggesting protein synthesis efficiency is greatly improved during silkworm domestication. While there was an increase in fibroin protein synthesis, the production of sericin protein was simultaneously reduced in the silk gland of the domestic silkworm. This reflects that domestic and wild silkworms have been under different selection pressures. Importantly, we found that the network co-expressed with the silk-coding genes of the domestic silkworm was larger than that of the wild silkworm. Furthermore, many more genes co-expressed with silk-coding genes in the domestic silkworm were subjected to artificial selection than those in the wild silkworm. Our results revealed that the increase of silk yield during silkworm domestication is involved in improvement of a biological system which includes not only expansion of “factories” (cells of PSG) of protein synthesis, but also a high expression of silk-coding genes and silk production-related genes such as biological energy, transport, and ribosome pathway genes.