Exogenous substances regulate silkworm fat body protein synthesis through MAPK and PI3K/Akt signaling pathways

Gene expression differentials driven by mass rearing and artificial selection in black soldier fly colonies

The micro-evolutionary forces that shape genetic diversity during domestication have been assessed in many plant and animal systems. However, the impact of these processes on gene expression, and consequent functional adaptation to artificial environments, remains under-investigated. In this study, whole-transcriptome dynamics associated with the early stages of domestication of the black soldier fly (BSF), Hermetia illucens, were assessed. Differential gene expression (DGE) was evaluated in relation to (i) generational time within the cultured environment (F2 vs. F3), and (ii) two selection strategies [no artificial selective pressure (NS); and selection for greater larval mass (SEL)]. RNA-seq was conducted on 5th instar BSF larvae (n = 36), representing equal proportions of the NS (F2 = 9; F3 = 9) and SEL (F2 = 9; F3 = 9) groups. A multidimensional scaling plot revealed greater gene expression variability within the NS and F2 subgroups, while the SEL group clustered separately with lower levels of variation. Comparisons between generations revealed 898 differentially expressed genes (DEGs; FDR-corrected p < 0.05), while between selection strategies, 213 DEGs were observed (FDR-corrected p < 0.05). Enrichment analyses revealed that metabolic, developmental, and defence response processes were over-expressed in the comparison between F2 and F3 larvae, while metabolic processes were the main differentiating factor between NS and SEL lines. This illustrates the functional adaptations that occur in BSF colonies across generations due to mass rearing; as well as highlighting genic dynamics associated with artificial selection for production traits that might inform future selective breeding strategies.

Evaluation of tolerance to λ-cyhalothrin and response of detoxification enzymes in silkworms reared on artificial diet

A representative silkworm rearing mode of Ⅰ-Ⅲ instars reared on artificial diet and Ⅳ-Ⅴ instars reared on fresh mulberry leaves has been recognized in some sericultural areas of China. Under this rearing mode, silkworms are prone to be poisoned by pesticide residues on mulberry leaves at the Ⅳ and Ⅴ instar stages. As one of the most widely applied insecticides, λ-cyhalothrin was used to study the insecticide tolerance of silkworm reared on artificial diet (referred as the AD group). Our results showed that the newly ecdysized Ⅳ instar larvae in the AD group were less tolerant to λ-cyhalothrin compared to the mulberry leaves reared group (referred as the ML group). After continuous exposure to trace λ-cyhalothrin, the weight gain and the survival rate of silkworms were significantly lower in the AD group than those in the ML group, even though compensatory growth was observed in the control of the AD group. Histopathology and ultrastructure of fat body showed that λ-cyhalothrin induced more severe cell injuries in the AD group, such as shrunken nucleus, dilatation of endoplasmic reticulum, and mitochondrial swelling. The transcription levels of detoxification related genes (CYP4M5, CYP6AB4, CarE2, CarE5, GSTe1 and GSTe3) and the enzyme activities of P450s, CarEs and GSTs were inducible by trace λ-cyhalothrin in a time-specific manner, and the data showed that the response of P450 enzyme activity was retarded in the AD group, indicating a potential reason for a higher sensitivity to λ-cyhalothrin. Our results provided a new clue for the study of the relationship between feed nutrition and detoxification ability, and also provided an important reference for the development of modern silkworm rearing mode.

Transcriptome Analysis Reveals the Gene Expression Changes in the Silkworm (Bombyx mori) in Response to Hydrogen Sulfide Exposure

Simple Summary

The fat body is one of the most important tissues in the body of insects due to its number of functions. Nowadays the new physiological function of H₂S has gained attention as a novel signaling molecule. H₂S performs crucial regulatory functions involving growth, the cardiovascular system, oxidative stress, and inflammation in many organisms. In this study, RNA-seq technology was used to investigate the fat body of the silkworm at the transcriptional level after H₂S exposure during the 5th larvae stage. A total of 1200 (DEGs) was identified after 7.5 µM H₂S treatment, of which 977 DEGs were up-regulated and 223 DEGs were down-regulated. DEGs were mainly involved in the transport pathway, cellular community, carbohydrate metabolism, and immune-associated signal transduction. Present research provides new insights on the gene expression changes in the fat body of silkworms after H₂S exposure.

Abstract

Hydrogen sulfide (H₂S) has been recognized for its beneficial influence on physiological alterations. The development (body weight) and economic characteristics (cocoon weight, cocoon shell ratio, and cocoon shell weight) of silkworms were increased after continuous 7.5 µM H₂S treatment. In the present study, gene expression changes in the fat body of silkworms at the 5th instar larvae in response to the H₂S were investigated through comparative transcriptome analysis. Moreover, the expression pattern of significant differentially expressed genes (DEGs) at the 5th instar larvae was confirmed by quantitative real-time PCR (qRT-PCR) after H₂S exposure. A total of 1200 (DEGs) was identified, of which 977 DEGs were up-regulated and 223 DEGs were down-regulated. Most of the DEGs were involved in the transport pathway, cellular community, carbohydrate metabolism, and immune-associated signal transduction. The up regulated genes under H₂S exposure were involved in endocytosis, glycolysis/gluconeogenesis, the citrate cycle (TCA cycle), and the synthesis of fibroin, while genes related to inflammation were down-regulated, indicating that H₂S could promote energy metabolism, the transport pathway, silk synthesis, and inhibit inflammation in the silkworm. In addition, the expression levels of these genes were increased or decreased in a time-dependent manner during the 5th instar larvae. These results provided insight into the effects of H₂S on silkworms at the transcriptional level and a substantial foundation for understanding H₂S function.

Zearalenone promotes apoptosis of mouse Leydig cells by targeting phosphatase and tensin homolog and thus inhibiting the PI3K/AKT signal pathway

Zearalenone (ZEA) is a mycotoxin with estrogenic activity whose main effect is to impair the reproductive systems of animals. It leads to reproductive disorders in livestock and thus causes serious losses to agriculture and animal husbandry. This study aims to examine whether ZEA induces toxicity in Leydig cells through the PI3K/AKT signaling pathway and also to investigate the role played by the upstream phosphatase and tensin homolog (PTEN) gene. An adenovirus vector model was constructed to interfere with the PTEN gene to investigate whether ZEA promotes the apoptosis of TM3 cells through the PI3K/AKT pathway. Apoptosis was detected cytometrically and the protein expression levels of PTEN, AKT, p-AKT, Bax, and Bcl-2 were evaluated via western blot analysis. The results show that ZEA induces apoptosis of TM3 cells. PTEN expression is significantly increased (P < 0.01), Bax expression is increased (P < 0.05), AKT and p-AKT expression of anti-apoptotic protein is significantly decreased (P < 0.01), and Bcl-2 protein expression is decreased (P < 0.05) in the ZEA group compared with the control group. In the shRNA+ZEA group, the expression levels of PTEN and Bax proteins are significantly decreased (P < 0.01), AKT protein is significantly increased (P < 0.01), and p-AKT protein is increased (P < 0.05) compared with the ZEA group. This study thus demonstrates that ZEA promotes apoptosis of TM3 cells by targeting PTEN and thus inhibiting the PI3K/AKT signal pathway.

Organophosphorus Compounds and MAPK Signaling Pathways

The molecular signaling pathways that lead to cell survival/death after exposure to organophosphate compounds (OPCs) are not yet fully understood. Mitogen-activated protein kinases (MAPKs) including the extracellular signal-regulated protein kinase (ERK), the c-Jun NH2-terminal kinase (JNK), and the p38-MAPK play the leading roles in the transmission of extracellular signals into the cell nucleus, leading to cell differentiation, cell growth, and apoptosis. Moreover, exposure to OPCs induces ERK, JNK, and p38-MAPK activation, which leads to oxidative stress and apoptosis in various tissues. However, the activation of MAPK signaling pathways may differ depending on the type of OPCs and the type of cell exposed. Finally, different cell responses can be induced by different types of MAPK signaling pathways after exposure to OPCs.

Identification of the nucleotide exchange factor BmGrpE and its role under high-temperature stress in silkworm, Bombyx mori

The high-temperature stress gene GrpE plays an important role in coping with high-temperature stress. The mutation of key sites of this gene can improve the high-temperature resistance of organisms. In the present study, using complementary DNAs from the silkworm fat body as the template, the open reading frame sequence of the GrpE gene (BmGrpE) was amplified and was found to be 644 bp in length and encode a protein with a predicted molecular weight of 24.1 kDa. The presence of a binding site for the heat shock transcription factor (Hsf1) at -1440 bp upstream of its coding region indicates that BmGrpE may respond to high-temperature stress. BmGrpE was constitutively expressed throughout developmental stages, with the highest level observed in the 5th instar larvae stage. Moreover, in 5th instar larvae (the 3th day), BmGrpE was expressed in all tissues examined, with the highest levels in the fat body, silk gland, and midgut. Interestingly, under high-temperature stress, TiO₂ nanoparticle treatment increased the messenger RNA levels of BmGrpE in the fat body and silk gland. After treatment with dsRNA of BmGrpE, the cell viability of BmN cells was significantly decreased under 34°C and H₂ O₂ stress (p < .05). Mutation of BmGrpE (H163L) enhanced the resistance of BmN cells under high-temperature stress. These results provide new clues for the study of molecular mechanisms of insect resistance to high temperatures.

Identification of genetic effects and potential causal polymorphisms of CPM gene impacting milk fatty acid traits in Chinese Holstein

Our previous GWAS revealed 83 significant SNPs and 20 promising candidate genes associated with milk fatty acid traits in dairy cattle. Out of them, the carboxypeptidase M (CPM) gene contains a genome-wide significant SNP, Hapmap49848-BTA-106779, which is strongly associated with myristic acid (C14:0; P = 0.0064). Herein, we aimed to confirm the genetic effects of CPM on milk fatty acids in Chinese Holstein. Seven SNPs were detected by re-sequencing the sequences of entire exons and 3000 bp of up-/downstream flanking regions of the CPM gene, of which three were in 5' flanking region, one in the 3' UTR and three were in the 3' flanking region. Using the Haploview 4.1, we estimated the LD among the identified SNPs and found two haplotype blocks. With the animal model, we performed the SNP- and haplotype-based association analyses, and observed that these SNPs and haplotype blocks mainly had strong genetic associations with medium-chain saturated fatty acids (caproic acid, C6:0; caprylic acid, C8:0; capric acid, C10:0; and lauric acid, C12:0) (P < 0.0001-0.0257). In addition, using the Genomatix software, we predicted that three SNPs in the 5' flanking region of CPM (g.45079507A>G, g.45080228C>A and g.45080335C>G) changed the transcription factor binding sites for PREF (progesterone receptor biding site), ZBRK1 (transcription factor with eight central zinc fingers and an N-terminal KRAB domain), SOX9 (sex-determining region Y-box 9, dimeric binding sites), SOX6 (sex-determining region Y-box 6) and FOXP1-ES (alternative splicing variant of FOXP1, activated in ESCs). Further, the dual-luciferase reporter assay showed these three SNPs altered the transcriptional activity of CPM gene (P ≤ 0.0006). In summary, using the post-GWAS strategy, we first confirmed the significant genetic effects of CPM with milk fatty acids in dairy cattle, and identified three potential causal mutations.

Transcriptome analysis reveals gene expression changes of the fat body of silkworm (Bombyx mori L.) in response to selenium treatment

A comparative transcriptome analysis was conducted to investigate the gene expression changes in the fat body of silkworm after treatment with different concentrations (50 μM and 200 μM) of selenium (Se). 912 differential expression genes (DEGs) (371 up-regulated and 541 down-regulated) and 1420 DEGs (1078 up-regulated and 342 down-regulated) were identified in silkworm fat body treated with 50 μM and 200 μM of Se, respectively. In case of 50 μM group, DEGs were mainly enriched in the peroxisome pathway and fatty acid metabolism pathway, and later were associated with antioxidant defense and nutrition regulation. After 200 μM Se-treatment, DEGs were mainly located in the glycerolipid metabolism and arachidonic acid metabolism pathways, which further encoded detoxification related genes. Furthermore, 32 candidate DEGs from these pathways had been selected to confirm the RNA-seq data. Among these DEGs, 14 genes were up-regulated in the 50 μM Se-treated group (only three genes in the 200 μM Se-treated group) which were involved in lipid metabolism and antioxidant defense, and 13 up-regulated genes (only two genes were up-regulated in the 50 μM Se-treated group) were involved in detoxification of the 200 μM Se-treated group. These changes showed that lower concentration of Se could regulate the nutrition and promote antioxidation pathways; whereas, high levels of Se promoted the detoxification of silkworm. These findings can be helpful to understand the possible mechanisms of Se action and detoxification in silkworm and other insects.

Melatonin suppresses milk fat synthesis by inhibiting the mTOR signaling pathway via the MT1 receptor in bovine mammary epithelial cells

Milk fat content is an important criterion for assessing milk quality and is one of the main target traits of dairy cattle breeding. Recent studies have shown the importance of melatonin in regulating lipid metabolism, but the potential effects of melatonin on milk fat synthesis in bovine mammary epithelial cells (BMECs) remain unclear. Here, we showed that melatonin supplementation at 10 μmol/L significantly downregulated the mRNA expression of lipid metabolism-related genes and resulted in lower lipid droplet formation and triglyceride accumulation. Moreover, melatonin significantly upregulated melatonin receptor subtype melatonin receptor 1a (MT1) gene expression, and the negative effects of melatonin on milk fat synthesis were reversed by treatment with the nonselective MT1/melatonin receptor subtype melatonin receptor 1b (MT2) antagonist. However, a selective MT2 antagonist did not modify the negative effects of melatonin on milk fat synthesis. In addition, KEGG analysis revealed that melatonin inhibition of milk fat synthesis may occur via the mTOR signaling pathway. Further analysis revealed that melatonin significantly suppressed the activation of the mTOR pathway by restricting the phosphorylation of mTOR, 4E-BP1, and p70S6K, and the inhibition of melatonin on milk fat synthesis was reversed by mTOR activator MHY1485 in BMECs. Furthermore, in vivo experiments in Holstein dairy cows showed that exogenous melatonin significantly decreased milk fat concentration. Our data from in vitro and in vivo studies revealed that melatonin suppresses milk fat synthesis by inhibiting the mTOR signaling pathway via the MT1 receptor in BMECs. These findings lay a foundation to identify a new potential means for melatonin to modulate the fat content of raw milk in Holstein dairy cows.

Transcriptomic analysis of boron hyperaccumulation mechanisms in Puccinellia distans

Puccinellia distans, common alkali grass, is found throughout the world and can survive in soils with boron concentrations that are lethal for other plant species. Indeed, P. distans accumulates very high levels of this element. Despite these interesting features, very little research has been performed to elucidate the boron tolerance mechanism in this species. In this study, P. distans samples were treated for three weeks with normal (0.5 mg L^-1) and elevated (500 mg L^-1) boron levels in hydroponic solution. Expressed sequence tags (ESTs) derived from shoot tissue were analyzed by RNA sequencing to identify genes up and down-regulated under boron stress. In this way, 3312 differentially expressed transcripts were detected, 67.7% of which were up-regulated and 32.3% of which were down-regulated in boron-treated plants. To partially confirm the RNA sequencing results, 32 randomly selected transcripts were analyzed for their expression levels in boron-treated plants. The results agreed with the expected direction of change (up or down-regulation). A total of 1652 transcripts had homologs in A. thaliana and/or O. sativa and mapped to 1107 different proteins. Functional annotation of these proteins indicated that the boron tolerance and hyperaccumulation mechanisms of P. distans involve many transcriptomic changes including: alterations in the malate pathway, changes in cell wall components that may allow sequestration of excess boron without toxic effects, and increased expression of at least one putative boron transporter and two putative aquaporins. Elucidation of the boron accumulation mechanism is important in developing approaches for bioremediation of boron contaminated soils.

Cloning and Functional Analysis of CncC and Keap1 Genes in Silkworm

CncC/keap1-ARE is an important signaling pathway for detoxification and antioxidation in Diptera and Coleoptera insects. However, such a signaling pathway has not been studied in Bombyx mori. In this study, BmCncC and Bmkeap1 genes were cloned, their amino acid sequences were analyzed, and each functional domain was mapped. Through phylogenetic analysis and sequence comparison among multiple species, we found that the Neh1 motif of CncC was highly conserved and the DLG motif was replaced by the DMG motif in Neh2. Conformational analysis showed that Neh1 of BmCncC forms a hairpin structure to bind DNA. The DGR region of Bmkeap1 contained abundant β sheets, which was involved in the recognition of Neh2. The transcription and expression analyses showed that both BmCncC and Bmkeap1 were highly expressed in the first instar larvae, and these two genes were expressed at a high level in the reproductive gland, fat body, and head. The transcriptional and expression levels of Akt and BmCncC in the fat body were significantly upregulated, and the expression of Bmkeap1 was downregulated after the phoxim treatment in silkworm. The transcriptional levels of CncC-regulated detoxification enzymes GST, cyp4M5, cyp6AE2, and cyp9G3 were increased by 4.026-, 5.246-, 3.821-, and 9.787-fold, respectively, while the activities of GST and CYP450 were increased by 1.521- and 1.231-fold, respectively, after phoxim treatment. These results indicated that the BmCncC/Bmkeap1 signaling pathway was activated by phoxim, leading to the expression of downstream detoxifying enzymes and detoxification of phoxim in silkworm.