Transcriptomics and proteomics-based analysis of heterosis on main economic traits of silkworm, Bombyx mori

Comprehensive review on glucose 6 phosphate dehydrogenase: A critical immunometabolic and redox switch in insects

Mounting an active immune response is energy intensive and demands the reallocation of nutrients to maintain the body's resistance and tolerance against infections. Central to this metabolic adaptation is Glucose-6-phosphate dehydrogenase (G6PDH), a housekeeping enzyme involve in pentose phosphate pathway (PPP). PPP play an essential role in generating ribose, which is critical for nicotinamide adenine dinucleotide phosphate (NADPH). It is vital for physiological and cellular processes such as generating nucleotides, fatty acids and reducing oxidative stress. The G6PDH is extremely conserved enzyme across species in PP shunt. The deficiency of enzymes leads to serious consequences on organism, particularly on adaptation and development. Acute deficiency can lead to impaired cell development, halted embryonic growth, reduce sensitivity to insulin, hypertension and increase inflammation. Historically, research focusing on G6PDH and PPP have primarily targeted diseases on mammalian. However, our review has investigated the unique functions of the G6PDH enzyme in insects and greatly improved mechanistic understanding of its operations. This review explore how G6PDH in insects plays a crucial role in managing the redox balance and immune related metabolism. This study aims to investigate the enzyme's role in different metabolic adaptations.

Schistosoma haematobium and Schistosoma bovis first generation hybrids undergo gene expressions changes consistent with species compatibility and heterosis

When two species hybridize, the two parental genomes are brought together and some alleles might interact for the first time. To date, the extent of the transcriptomic changes in first hybrid generations, along with their functional outcome constitute an important knowledge gap, especially in parasite species. Here we explored the molecular and functional outcomes of hybridization in first-generation hybrids between the blood fluke parasites Schistosoma haematobium and S. bovis. Through a transcriptomic approach, we measured gene expression in both parental species and hybrids. We described and quantified expression profiles encountered in hybrids along with the main biological processes impacted. Up to 7,100 genes fell into a particular hybrid expression profile (intermediate between the parental expression levels, over-expressed, under-expressed, or expressed like one of the parental lines). Most of these genes were different depending on the direction of the parental cross (S. bovis mother and S. haematobium father or the reverse) and depending on the sex. For a given sex and cross direction, the vast majority of genes were hence unassigned to a hybrid expression profile: either they were differentially expressed genes but not typical of any hybrid expression profiles or they were not differentially expressed neither between hybrids and parental lines nor between parental lines. The most prevalent profile of gene expression in hybrids was the intermediate one (24% of investigated genes). These results suggest that transcriptomic compatibility between S. haematobium and S. bovis remains quite high. We also found support for an over-dominance model (over- and under-expressed genes in hybrids compared to parental lines) potentially associated with heterosis. In females in particular, processes such as reproductive processes, metabolism and cell interactions as well as signaling pathways were indeed affected. Our study hence provides new insight on the biology of Schistosoma hybrids with evidences supporting compatibility and heterosis.

DIA-Based Proteomic Analysis Reveals MYOZ2 as a Key Protein Affecting Muscle Growth and Development in Hybrid Sheep

Hybridization of livestock can be used to improve varieties, and different hybrid combinations produce unique breeding effects. In this study, male Southdown and Suffolk sheep were selected to hybridize with female Hu sheep to explore the effects of male parentage on muscle growth and the development of offspring. Using data-independent acquisition technology, we identified 119, 187, and 26 differentially abundant proteins (DAPs) between Hu × Hu (HH) versus Southdown × Hu (NH), HH versus Suffolk × Hu (SH), and NH versus SH crosses. Two DAPs, MYOZ2 and MYOM3, were common to the three hybrid groups and were mainly enriched in muscle growth and development-related pathways. At the myoblast proliferation stage, MYOZ2 expression decreased cell viability and inhibited proliferation. At the myoblast differentiation stage, MYOZ2 expression promoted myoblast fusion and enhanced the level of cell fusion. These findings provide new insights into the key proteins and metabolic pathways involved in the effect of male parentage on muscle growth and the development of hybrid offspring in sheep.

Transcriptome analysis of ovarian tissues highlights genes controlling energy homeostasis and oxidative stress as potential drivers of heterosis for egg number and clutch size in crossbred laying hens

Heterosis is the major benefit of crossbreeding and has been exploited in laying hens breeding for a long time. This genetic phenomenon has been linked to various modes of nonadditive gene action. However, the molecular mechanism of heterosis for egg production in laying hens has not been fully elucidated. To fill this research gap, we sequenced mRNAs and lncRNAs of the ovary stroma containing prehierarchical follicles in White Leghorn, Rhode Island Red chickens as well as their reciprocal crossbreds that demonstrated heterosis for egg number and clutch size. We further delineated the modes of mRNAs and lncRNAs expression to identify their potential functions in the observed heterosis. Results showed that dominance was the principal mode of nonadditive expression exhibited by mRNAs and lncRNAs in the prehierarchical follicles of crossbred hens. Specifically, low-parent dominance was the main mode of mRNA expression, while high-parent dominance was the predominant mode of lncRNA expression. Important pathways enriched by genes that showed higher expression in crossbreds compared to either one or both parental lines were cell adhesion molecules, tyrosine and purine metabolism. In contrast, ECM-receptor interaction, focal adhesion, PPAR signaling, and ferroptosis were enriched in genes with lower expression in the crossbred. Protein network interaction identified nonadditively expressed genes including apolipoprotein B (APOB), transferrin, acyl-CoA synthetase medium-chain family member (APOBEC) 3, APOBEC1 complementation factor, and cathepsin S as hub genes. Among these potential hub genes, APOB was the only gene with underdominance expression common to the 2 reciprocal crossbred lines, and has been linked to oxidative stress. LncRNAs with nonadditive expression in the crossbred hens targeted natriuretic peptide receptor 1, epidermal differentiation protein beta, spermatogenesis-associated gene 22, sperm-associated antigen 16, melanocortin 2 receptor, dolichol kinase, glycine amiinotransferase, and prolactin releasing hormone receptor. In conclusion, genes with nonadditive expression in the crossbred may play crucial roles in follicle growth and atresia by improving follicle competence and increasing oxidative stress, respectively. These 2 phenomena could underpin heterosis for egg production in crossbred laying hens.

G6PDH as a key immunometabolic and redox trigger in arthropods

The enzyme glucose-6-phosphate dehydrogenase (G6PDH) plays crucial roles in glucose homeostasis and the pentose phosphate pathway (PPP), being also involved in redox metabolism. The PPP is an important metabolic pathway that produces ribose and nicotinamide adenine dinucleotide phosphate (NADPH), which are essential for several physiologic and biochemical processes, such as the synthesis of fatty acids and nucleic acids. As a rate-limiting step in PPP, G6PDH is a highly conserved enzyme and its deficiency can lead to severe consequences for the organism, in particular for cell growth. Insufficient G6PDH activity can lead to cell growth arrest, impaired embryonic development, as well as a reduction in insulin sensitivity, inflammation, diabetes, and hypertension. While research on G6PDH and PPP has historically focused on mammalian models, particularly human disorders, recent studies have shed light on the regulation of this enzyme in arthropods, where new functions were discovered. This review will discuss the role of arthropod G6PDH in regulating redox homeostasis and immunometabolism and explore potential avenues for further research on this enzyme in various metabolic adaptations.

Transcriptomic analysis of the fat body of resistant and susceptible silkworm strains, Bombyx mori (Lepidoptera), after oral treatment with fenpropathrin

The widespread use of pyrethroid pesticides has brought serious economic losses in sericulture, but there is still no viable solution. The key to solving the problem is to improve silkworm resistance to pesticides, which depends on understanding the resistance mechanism of silkworms to pesticides. This study aimed to use transcriptomes to understand the underlying mechanism of silkworm resistance to fenpropathrin, which will provide a theoretical molecular reference for breeding pesticide-resistant silkworm varieties. In this study, the fat bodies of two strains with differential resistance after 12 h of fenpropathrin feeding were analyzed using RNA-Seq. After feeding fenpropathrin, 760 differentially expressed genes (DEGs) were obtained in the p50(r) strain and 671 DEGs in the 8y strain. The DEGs involved in resistance to fenpropathrin were further identified by comparing the two strains, including 207 upregulated DEGs in p50(r) and 175 downregulated DEGs in 8y. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis showed that these fenpropathrin-related DEGs are mainly enriched in the metabolism and transporter pathways. Moreover, 28 DEGs involved in the metabolic pathway and 18 in the transporter pathway were identified. Furthermore, organic cation transporter protein 6 (BmOCT6), a transporter pathway member, was crucial in enhancing the tolerance of BmN cells to fenpropathrin. Finally, the knockdown of the expression of the homologs of BmOCT6 in Glyphodes pyloalis (G. pyloalis) significantly decreased the resistant level of larvae to fenpropathrin. The findings showed that the metabolism and transporter pathways are associated with resistance to fenpropathrin in silkworm, and OCT6 is an effective and potential target not only for silkworm breeding but also for pest biocontrol.

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.

Multi-tissue metabolomic profiling reveals potential mechanisms of cocoon yield in silkworms (Bombyx mori) fed formula feed versus mulberry leaves

Use of formula feed (FF) for silkworms for all instars, has promoted transformation and progress in traditional sericulture. However, the cocoon yield of FF silkworms has failed to reach that of silkworms fed mulberry leaves (ML). The biological mechanisms underlying this phenomenon have not been well described. This study aimed to identify metabolic mechanisms and potential biomarkers relating to the poor cocoon yield of FF silkworms. In this study, silkworms received treatments of either ML (ML group) or FF (FF group) for all instars. At the 3rd day of the 5th instar, the midgut (MG), hemolymph (HL) and posterior silk gland (PSG) were collected for the metabolome profiles detection. The remaining silkworms were fed ML or FF until cocooning for investigation. The whole cocoon yield (WCY) was significantly higher in the FF group than the ML group (p < 0.05), whereas the cocoon shell weight (CSW) and cocoon shell rate (CSR) were significantly lower in the FF group (p < 0.05). A total of 845, 867 and 831 metabolites were qualified and quantified in the MG, HL and PSG of the FF silkworms, respectively. Correspondingly, 789, 833 and 730 metabolites were quantified in above three tissues of the ML group. Further, 230, 249 and 304 significantly different metabolites (SDMs) were identified in the MG, HL and PSG between the FF and ML group, respectively. Eleven metabolic pathways enriched by the SDMs were mutual among the three tissues. Among them, cysteine and methionine metabolism, arginine biosynthesis, and arginine and proline metabolism were the top three pathways with the highest impact value in the PSG. Six biomarkers were obtained through biomarker analysis and Pearson correlation calculation. Among them, homocitrulline, glycitein, valyl-threonine, propyl gallate and 3-amino-2,3-dihydrobenzoic acid were positively correlated with WCY, but negatively correlated with CSW and CSR (p < 0.05). An opposite correlation pattern was observed between 3-dimethylallyl-4-hydroxyphenylpyruvate and the three cocoon performance traits. Overall, three key metabolic pathways and six biomarkers associated with cocoon yield were interpreted, and should provide directions for formula feed optimization in factory-raised silkworms.

Transcriptome Profiling Revealed Basis for Growth Heterosis in Hybrid Tilapia (Oreochromis niloticus ♀ × O. aureus ♂)

Hybrid tilapia were produced from hybridization of Nile tilapia (Oreochromis niloticus) and blue tilapia (O. aureus). Comparative transcriptome analysis was carried out on the liver of hybrid tilapia and their parents by RNA sequencing. A total of 2319 differentially expressed genes (DEGs) were identified. Trend co-expression analysis showed that non-additive gene expression accounted for 67.1% of all DEGs. Gene Ontology and KEGG enrichment analyses classified the respective DEGs. Gene functional enrichment analysis indicated that most up-regulated genes, such as FASN, ACSL1, ACSL3, ACSL6, ACACA, ELOVL6, G6PD, ENO1, GATM, and ME3, were involved in metabolism, including fatty acid biosynthesis, unsaturated fatty acid biosynthesis, glycolysis, pentose phosphate pathway, amino acid metabolism, pyruvate metabolism, and the tricarboxylic acid cycle. The expression levels of a gene related to ribosomal biosynthesis in eukaryotes, GSH-Px, and those associated with heat shock proteins (HSPs), such as HSPA5 and HSP70, were significantly down-regulated compared with the parent tilapia lineages. The results revealed that the metabolic pathway in hybrid tilapia was up-regulated, with significantly improved fatty acid metabolism and carbon metabolism, whereas ribosome biosynthesis in eukaryotes and basal defense response were significantly down-regulated. These findings provide new insights into our understanding of growth heterosis in hybrid tilapia.

Time-Course Transcriptome Analysis Reveals Global Gene Expression Profiling and Dynamic Developmental Signatures across Complete Life Cycle of Bombyx mori

Background: The silkworm (Bombyx mori) is an important lepidopteran model insect worldwide which undergoes a complete metamorphosis developmental process. Although genome sequencing has been long performed, no transcriptome data covering the complete life cycle are available. Methods: Herein, a total of 10 samples were collected consecutively at four different developmental stages, including eggs of 24 h after oviposition (Ed) and eggs of 24 h after artificial egg-hatching (E); larvae from fist to fifth instar (L1–L5); early and late pupa (P4 and P8); and adult moth (M), were subjected to Illumina RNA-Seq and time-course analysis. Results: The summations of the gene expression of the silkworm ten developmental stages show: at Ed stage, eggs develop towards diapause status, the total gene expression level is relatively low; at E stage, after artificial egg-hatching, the expression level improves rapidly; during larval stages from L1–L5, the expression level rises gradually and reaches a peak at L5 stage; during pupae and moth stages, the total gene expression decline stage by stage. The results revealed a dynamical gene expression profile exhibiting significant differential expressions throughout the silkworm life cycle. Moreover, stage-specific key genes were identified at different developmental stages, suggesting their functions mainly characterized in maintaining insect development and immunity homeostasis or driving metamorphosis. GO annotation and KEGG enrichment analysis further revealed the most significantly enriched and fundamentally biological processes during silkworm growth. Conclusion: Collectively, our omics data depicted the first comprehensive landscape of dynamic transcriptome throughout complete developmental processes of B. mori. Our findings also provide valuable references and novel insights into understanding the molecular developmental remodeling events for other Lepidoptera species.