Integrated Analysis of Transcriptomic, miRNA and Proteomic Changes of a Novel Hybrid Yellow Catfish Uncovers Key Roles for miRNAs in Heterosis

Synergy of gut microbiota and host genome in driving heterosis expression of chickens

Heterosis has been widely utilized in agricultural production. Despite over a century of extensive research, the underlying mechanisms of heterosis remain elusive. Most hypotheses and research have focused on the genetic basis of heterosis. However, the potential role of gut microbiota in heterosis has been largely ignored. Here, we carefully design a crossbreeding experiment with two distinct broiler breeds and conduct 16S rRNA amplicon and transcriptome sequencing to investigate the synergistic role of gut microbiota and host genes in driving heterosis. We find that the breast muscle weight of the hybrids exhibits a high heterosis, 6.28% higher than mid-parent value. A notable difference is observed in the composition and potential function of cecal microbiota between hybrids and their parents. Over 90% of the differentially colonized microbiota and differentially expressed genes exhibit nonadditive patterns. Integrative analyses uncover associations between nonadditive genes and nonadditive microbiota, including a connection between the expression of cellular signaling pathway and metabolism-related genes and the abundance of Odoribacter, Oscillibacter, and Alistipes in hybrids. Moreover, higher abundances of these microbiota are related to better meat yield. In summary, these findings highlight the importance of gut microbiota in heterosis, serving as crucial factors that modulate heterosis expression in chickens.

Multiple trait comparison and global intestine transcriptional provide new insights into bases of heterosis in hybrid tilapia (Oreochromis niloticus × Oreochromis aureus)

Heterosis has been utilized in aquaculture for many years, yet its molecular basis remains elusive. Therefore, a comprehensive analysis of heterosis was conducted by comparing growth, digestion and biochemistry indices, as well as the intestinal gene expression profiles of Nile tilapia, blue tilapia and their hybrids. The results revealed that hybrid tilapia demonstrated an enhanced growth traits and elevated digestive enzyme activity compared to Nile and blue tilapia. Additionally, the hybrid tilapia displayed superior antioxidants and non-specific immune levels, with increased levels of catalase (CAT), alkaline phosphatase (AKP), acid phosphatase (ACP), glutathione (GSH), superoxide dismutase (SOD), total antioxidant capacity (TAOC), lysozyme, and immunoglobulin M (IgM) relative to Nile and blue tilapia. Moreover, 3392, 2470 and 1261 differentially expressed genes (DEGs) were identified in the intestinal tissues when comparing Nile tilapia to blue tilapia, hybrid tilapia to blue tilapia, and hybrid tilapia to Nile tilapia. Upon classifying the differentially expressed genes (DEGs), non-additively expressed DEGs accounted for 68.1 % of the total DEGs, with dominant and over-dominant expressed DEGs comprising 63.7 % and 4.4 % in the intestines, respectively. These non-additively expressed DEGs were primarily associated with metabolic, digestive, growth, and developmental pathways. This enrichment enhances our comprehension of the molecular underpinnings of growth heterosis in aquatic species.

Integrated miRNA and mRNA Sequencing Reveals the Sterility Mechanism in Hybrid Yellow Catfish Resulting from Pelteobagrus fulvidraco (♀) × Pelteobagrus vachelli (♂)

The hybrid yellow catfish exhibits advantages over pure yellow catfish in terms of fast growth, fast development, a high feeding rate, and strong immunity; additionally, it is almost sterile, thus ensuring the conservation of the genetic stock of fish populations. To investigate the sterility mechanism in hybrid yellow catfish (P. fulvidraco (♀) × P. vachelli (♂)), the mRNA and miRNA of the gonads of P. fulvidraco, P. vachelli, and a hybrid yellow catfish were analyzed to characterize the differentially expressed genes; this was carried out to help establish gene expression datasets to assist in the further determination of the mechanisms of genetic sterility in hybrid yellow catfish. In total, 1709 DEGs were identified between the hybrid and two pure yellow catfishes. A KEGG pathway analysis indicated that several genes related to reproductive functions were upregulated, including those involved in the cell cycle, progesterone-mediated oocyte maturation, and oocyte meiosis, and genes associated with ECM-receptor interaction were downregulated. The spermatogenesis-related GO genes CFAP70, RSPH6A, and TSGA10 were identified as being downregulated DEGs in the hybrid yellow catfish. Sixty-three DEmiRNAs were identified between the hybrid and the two pure yellow catfish species. The upregulated DEmiRNAs ipu-miR-194a and ipu-miR-499 were found to target the spermatogenesis-related genes CFAP70 and RSPH6A, respectively, playing a negative regulatory role, which may underscore the miRNA-mRNA regulatory mechanism of sterility in hybrid yellow catfish. The differential expression of ipu-miR-196d, ipu-miR-125b, and ipu-miR-150 and their target genes spidr, cep85, and kcnn4, implicated in reproductive processes, was verified via qRT-PCR, consistent with the transcriptome sequencing expression trends. This study provides deep insights into the mechanism of hybrid sterility in vertebrate groups, thereby contributing to achieving a better understanding and management of fish conservation related to hybrid sterility.

Strand-Specific RNA Sequencing Reveals Gene Expression Patterns in F1 Chick Breast Muscle and Liver after Hatching

Heterosis refers to the phenomenon where hybrids exhibit superior performance compared to the parental phenotypes and has been widely utilized in crossbreeding programs for animals and crops, yet the molecular mechanisms underlying this phenomenon remain enigmatic. A better understanding of the gene expression patterns in post-hatch chickens is very important for exploring the genetic basis underlying economically important traits in the crossbreeding of chickens. In this study, breast muscle and liver tissues (n = 36) from full-sib F1 birds and their parental pure lines were selected to identify gene expression patterns and differentially expressed genes (DEGs) at 28 days of age by strand-specific RNA sequencing (ssRNA-seq). This study indicates that additivity is the predominant gene expression pattern in the F1 chicken post-hatch breast muscle (80.6% genes with additivity) and liver (94.2% genes with additivity). In breast muscle, Gene Ontology (GO) enrichment analysis revealed that a total of 11 biological process (BP) terms closely associated with growth and development were annotated in the identified DEG sets and non-additive gene sets, including STAT5A, TGFB2, FGF1, IGF2, DMA, FGF16, FGF12, STAC3, GSK3A, and GRB2. Kyoto Encyclopedia of Genes and Genomes (KEGG) annotation presented that a total of six growth- and development-related pathways were identified, involving key genes such as SLC27A4, GLUL, TGFB2, COX17, and GSK3A, including the PPAR signaling pathway, TGF-beta signaling pathway, and mTOR signaling pathway. Our results may provide a theoretical basis for crossbreeding in domestic animals.

Non-Additive and Asymmetric Allelic Expression of p38 mapk in Hybrid Tilapia (Oreochromis niloticus ♀ × O. aureus ♂)

Hybridization is a widely used breeding technique in fish species that enhances desirable traits in cultured species through heterosis. However, the mechanism by which hybrids alter gene expression to form heterosis remains unclear. In this study, a group of hybrid tilapia was used to elucidate heterosis through interspecies crossing. Specifically, p38 was analyzed to describe the regulation of gene expression variation in hybrid tilapia. Transcripts from the Nile tilapia allele were found to be significantly higher than those from the blue tilapia allele in hybrid individuals, indicating that the expression of p38 was dominated by Nile tilapia sub-genomic alleles. The study also found a compensatory interaction of cis- and trans-acting elements of the Nile tilapia and blue tilapia sub-genomes, inducing a non-additive expression of p38 in hybrids. Eight specific SNPs were identified in the p38 promoter regions of Nile tilapia and blue tilapia, and were found to be promoter differences leading to differences in gene expression efficiencies between parental alleles using a dual-luciferase reporter system. This study provides insights into the non-additive expression patterns of key functional genes in fish hybrids related to growth and immunity response.

Integrated Network Analysis of microRNAs, mRNAs, and Proteins Reveals the Regulatory Interaction between hsa-mir-200b and CFL2 Associated with Advanced Stage and Poor Prognosis in Patients with Intestinal Gastric Cancer

Intestinal gastric cancer (IGC) carcinogenesis results from a complex interplay between environmental and molecular factors, ultimately contributing to disease development. We used integrative bioinformatic analysis to investigate IGC high-throughput molecular data to uncover interactions among differentially expressed genes, microRNAs, and proteins and their roles in IGC. An integrated network was generated based on experimentally validated microRNA-gene/protein interaction data, with three regulatory circuits involved in a complex network contributing to IGC progression. Key regulators were determined, including 23 microRNA and 15 gene/protein hubs. The regulatory circuit networks were associated with hallmarks of cancer, e.g., cell death, apoptosis and the cell cycle, the immune response, and epithelial-to-mesenchymal transition, indicating that different mechanisms of gene regulation impact similar biological functions. Altered expression of hubs was related to the clinicopathological characteristics of IGC patients and showed good performance in discriminating tumors from adjacent nontumor tissues and in relation to T stage and overall survival (OS). Interestingly, expression of upregulated hub hsa-mir-200b and its downregulated target hub gene/protein CFL2 were related not only to pathological T staging and OS but also to changes during IGC carcinogenesis. Our study suggests that regulation of CFL2 by hsa-miR-200b is a dynamic process during tumor progression and that this control plays essential roles in IGC development. Overall, the results indicate that this regulatory interaction is an important component in IGC pathogenesis. Also, we identified a novel molecular interplay between microRNAs, proteins, and genes associated with IGC in a complex biological network and the hubs closely related to IGC carcinogenesis as potential biomarkers.

The mRNA-lncRNA landscape of multiple tissues uncovers key regulators and molecular pathways that underlie heterosis for feed intake and efficiency in laying chickens

BACKGROUND

Heterosis is routinely exploited to improve animal performance. However, heterosis and its underlying molecular mechanism for feed intake and efficiency have been rarely explored in chickens. Feed efficiency continues to be an important breeding goal trait since feed accounts for 60 to 70% of the total production costs in poultry. Here, we profiled the mRNA-lncRNA landscape of 96 samples of the hypothalamus, liver and duodenum mucosa from White Leghorn (WL), Beijing-You chicken (YY), and their reciprocal crosses (WY and YW) to elucidate the regulatory mechanisms of heterosis.

RESULTS

We observed negative heterosis for both feed intake and residual feed intake (RFI) in YW during the laying period from 43 to 46 weeks of age. Analysis of the global expression pattern showed that non-additivity was a major component of the inheritance of gene expression in the three tissues for YW but not for WY. The YW-specific non-additively expressed genes (YWG) and lncRNA (YWL) dominated the total number of non-additively expressed genes and lncRNA in the hypothalamus and duodenum mucosa. Enrichment analysis of YWG showed that mitochondria components and oxidation phosphorylation (OXPHOS) pathways were shared among the three tissues. The OXPHOS pathway was enriched by target genes for YWL with non-additive inheritance of expression in the liver and duodenum mucosa. Weighted gene co-expression network analysis revealed divergent co-expression modules associated with feed intake and RFI in the three tissues from WL, YW, and YY. Among the negatively related modules, the OXPHOS pathway was enriched by hub genes in the three tissues, which supports the critical role of oxidative phosphorylation. Furthermore, protein quantification of ATP5I was highly consistent with ATP5I expression in the liver, which suggests that, in crossbred YW, non-additive gene expression is down-regulated and decreases ATP production through oxidative phosphorylation, resulting in negative heterosis for feed intake and efficiency.

CONCLUSIONS

Our results demonstrate that non-additively expressed genes and lncRNA involved in oxidative phosphorylation in the hypothalamus, liver, and duodenum mucosa are key regulators of the negative heterosis for feed intake and RFI in layer chickens. These findings should facilitate the rational choice of suitable parents for producing crossbred chickens.

Multi-omics analysis revealed the brain dysfunction induced by energy metabolism in Pelteobagrus vachelli under hypoxia stress

Hypoxia in water environment has become increasingly frequent and serious due to global warming and environmental pollution. Revealing the molecular mechanism of fish hypoxia adaptation will help to develop markers of environmental pollution caused by hypoxia. Here, we used a multi-omics method to identify the hypoxia-associated mRNA, miRNA, protein, and metabolite involved in various biological processes in Pelteobagrus vachelli brain. The results showed that hypoxia stress caused brain dysfunction by inhibiting energy metabolism. Specifically, the biological processes involved in energy synthesis and energy consumption are inhibited in P. vachelli brain under hypoxia, such as oxidative phosphorylation, carbohydrate metabolism and protein metabolism. Brain dysfunction is mainly manifested as blood-brain barrier injury accompanied by neurodegenerative diseases and autoimmune diseases. In addition, compared with previous studies, we found that P. vachelli has tissue specificity in response to hypoxia stress and the muscle suffers more damage than the brain. This is the first report to the integrated analysis of the transcriptome, miRNAome, proteome, and metabolome in fish brain. Our findings could provide insights into the molecular mechanisms of hypoxia, and the approach could also be applied to other fish species. DATA AVAILABILITY: The raw data of transcriptome has been uploaded to NCBI database (ID: SUB7714154 and SUB7765255). The raw data of proteome has been uploaded to ProteomeXchange database (PXD020425). The raw data of metabolome has been uploaded to Metabolight (ID: MTBLS1888).

Transcriptome Analysis of Reciprocal Hybrids Between Crassostrea gigas and C. angulata Reveals the Potential Mechanisms Underlying Thermo-Resistant Heterosis

Heterosis, also known as hybrid vigor, is widely used in aquaculture, but the molecular causes for this phenomenon remain obscure. Here, we conducted a transcriptome analysis to unveil the gene expression patterns and molecular bases underlying thermo-resistant heterosis in Crassostrea gigas ♀ × Crassostrea angulata ♂ (GA) and C. angulata ♀ × C. gigas ♂ (AG). About 505 million clean reads were obtained, and 38,210 genes were identified, of which 3779 genes were differentially expressed between the reciprocal hybrids and purebreds. The global gene expression levels were toward the C. gigas genome in the reciprocal hybrids. In GA and AG, 95.69% and 92.00% of the differentially expressed genes (DEGs) exhibited a non-additive expression pattern, respectively. We observed all gene expression modes, including additive, partial dominance, high and low dominance, and under- and over-dominance. Of these, 77.52% and 50.00% of the DEGs exhibited under- or over-dominance in GA and AG, respectively. The over-dominance DEGs common to reciprocal hybrids were significantly enriched in protein folding, protein refolding, and intrinsic apoptotic signaling pathway, while the under-dominance DEGs were significantly enriched in cell cycle. As possible candidate genes for thermo-resistant heterosis, GRP78, major egg antigen, BAG, Hsp70, and Hsp27 were over-dominantly expressed, while MCM6 and ANAPC4 were under-dominantly expressed. This study extends our understanding of the thermo-resistant heterosis in oysters.

Ubiquitylome study reveals the regulatory effect of α-lipoic acid on ubiquitination of key proteins in tryptophan metabolism pathway of pig liver

The decline in antioxidant defenses make it easily for human and animals to suffer from liver damage and diseases induced by oxidative stress, causing enormous losses to human health and livestock production. As one of the canonical protein post-translational modifications (PTMs), ubiquitination is widely involved in cell proliferation, apoptosis and damage/repair response, and is proven to be involved in the ability of mammals to resist oxidative stress. To explore whether α-lipoic acid (LA), a safe and efficient antioxidant, plays a role in regulating liver antioxidant status by PTMs, proteins in livers of pigs fed with LA were analyzed at the level of proteome and ubiquitylome. Based on proteome-wide enrichment of ubiquitination, a total of 7274 proteins were identified and 5326 were quantified, we also identified 1564 ubiquitination sites in 580 ubiquitinated proteins, among which there were 136 differentially ubiquitinated sites in 103 differentially ubiquitinated proteins upon LA. Further bioinformatics analysis showed that these differential proteins were mainly enriched in tryptophan metabolic pathway, and accompanied by significantly improvement of liver antioxidant capacity. We revealed the regulatory effect of LA on ubiquitination of kynurenine 3-monooxygenase (KMO) and other key proteins in tryptophan metabolism pathway of pig liver for the first time.

Integrated Analysis of Transcriptome and Metabolome Reveals Distinct Responses of Pelteobagrus fulvidraco against Aeromonas veronii Infection at Invaded and Recovering Stage

Yellow catfish (Pelteobagrus fulvidraco) is an important aquaculture fish susceptible to Aeromonas veronii infection, which causes acute death resulting in huge economic losses. Understanding the molecular processes of host immune defense is indispensable to disease control. Here, we conducted the integrated and comparative analyses of the transcriptome and metabolome of yellow catfish in response to A. veronii infection at the invaded stage and recovering stage. The crosstalk between A. veronii-induced genes and metabolites uncovered the key biomarkers and pathways that strongest contribute to different response strategies used by yellow catfish at corresponding defense stages. We found that at the A. veronii invading stage, the immune defense was strengthened by synthesizing lipids with energy consumption to repair the skin defense line and accumulate lipid droplets promoting intracellular defense line; triggering an inflammatory response by elevating cytokine IL-6, IL-10 and IL-1β following PAMP-elicited mitochondrial signaling, which was enhanced by ROS produced by impaired mitochondria; and activating apoptosis by up-regulating caspase 3, 7 and 8 and Prostaglandin F1α, meanwhile down-regulating FoxO3 and BCL6. Apoptosis was further potentiated via oxidative stress caused by mitochondrial dysfunction and exceeding inflammatory response. Additionally, cell cycle arrest was observed. At the fish recovering stage, survival strategies including sugar catabolism with D-mannose decreasing; energy generation through the TCA cycle and Oxidative phosphorylation pathways; antioxidant protection by enhancing Glutathione (oxidized), Anserine, and α-ketoglutarate; cell proliferation by inducing Cyclin G2 and CDKN1B; and autophagy initiated by FoxO3, ATG8 and ATP6V1A were highlighted. This study provides a comprehensive picture of yellow catfish coping with A. veronii infection, which adds new insights for deciphering molecular mechanisms underlying fish immunity and developing stage-specific disease control techniques in aquaculture.

Research Progress of the Gut Microbiome in Hybrid Fish

Fish, including hybrid species, are essential components of aquaculture, and the gut microbiome plays a vital role in fish growth, behavior, digestion, and immune health. The gut microbiome can be affected by various internal and/or external factors, such as host development, diet, and environment. We reviewed the effects of diet and dietary supplements on intestinal microorganisms in hybrid fish and the difference in the gut microbiome between the hybrid and their hybrids that originate. Then, we summarized the role of the gut microbiome in the speciation and ecological invasion of hybrid fish. Finally, we discussed possible future studies on the gut microbiome in hybrid fish, including the potential interaction with environmental microbiomes, the effects of the gut microbiome on population expansion, and fish conservation and management.

Integration of mRNA and miRNA Profiling Reveals Heterosis in Oreochromis niloticus × O. aureus Hybrid Tilapia

Heterosis is a widespread biological phenomenon in fishes, in which hybrids have superior traits to parents. However, the underlying molecular basis for heterosis remains uncertain. Heterosis in growth and survival rates is apparent in hybrid tilapia (Oreochromis niloticus ♀ × O. aureus ♂). Comparisons of growth and hematological biochemical characteristics and mRNA and miRNA transcriptional analyses were performed in hybrid and parents tilapia stocks to investigate the underlying molecular basis for heterosis. Growth characteristics and hematological glucose and cholesterol parameters were significantly improved in hybrids. Of 3097 differentially expressed genes (DEGs) and 120 differentially expressed miRNAs (DEMs) identified among three stocks (O. niloticus, O. aureus, and hybrids), 1598 DEGs and 62 DEMs were non-additively expressed in hybrids. Both expression level dominance and overdominance patterns occurred for DEGs and DEMs, indicating that dominance and overdominance models are widespread in the transcriptional and post-transcriptional regulation of genes involved in growth, metabolism, immunity, and antioxidant capacity in hybrid tilapia. Moreover, potential negative regulation networks between DEMs and predicted target DEGs revealed that most DEGs from miRNA-mRNA pairs are up-regulated. Dominance and overdominance models in levels of transcriptome and miRNAome facilitate the integration of advantageous parental alleles into hybrids, contributing to heterosis of growth and improved survival. The present study provides new insights into molecular heterosis in hybrid tilapia, advancing our understanding of the complex mechanisms involved in this phenomenon in aquatic animals.

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.

Integrated application of multi-omics strategies provides insights into the environmental hypoxia response in Pelteobagrus vachelli muscle

Increasing pressures on aquatic ecosystems because of pollutants, nutrient enrichment, and global warming have severely depleted oxygen concentrations. This sudden and significant lack of oxygen has resulted in persistent increases in fish mortality rates. Revealing the molecular mechanism of fish hypoxia adaptation will help researchers to find markers for hypoxia induced by environmental stress. Here, we used a multiomics approach to identify several hypoxia-associated miRNAs, mRNAs, proteins, and metabolites involved in diverse biological pathways in the muscles of Pelteobagrus vachelli. Our findings revealed significant hypoxia-associated changes in muscles over 4 h of hypoxia exposure and discrete tissue-specific patterns. We have previously reported that P. vachelli livers exhibit increased anaerobic glycolysis, heme synthesis, erythropoiesis, and inhibit apoptosis when exposed to hypoxia for 4 h. However, the opposite was observed in muscles. According to our comprehensive analysis, fishes show an acute response to hypoxia, including activation of catabolic pathways to generate more energy, reduction of biosynthesis to decrease energy consumption, and shifting from aerobic to anaerobic metabolic contributions. Also, we found that hypoxia induced muscle dysfunction by impairing mitochondrial function, activating inflammasomes, and apoptosis. The hypoxia-induced mitochondrial dysfunction enhanced oxidative stress, apoptosis, and further triggered interleukin-1β production via inflammasome activation. In turn, interleukin-1β further impaired mitochondrial function or apoptosis by suppressing downstream mitochondrial biosynthesis–related proteins, thus resulting in a vicious cycle of inflammasome activation and mitochondrial dysfunction. Our findings contribute meaningful insights into the molecular mechanisms of hypoxia, and the methods and study design can be utilized across different fish species.

•

First multiomics analysis of mRNA, miRNA, protein, and metabolite in fishes.

•

Liver and muscle were tissue-specific induced by hypoxia.

•

About 70 genes and 16 miRNAs related to hypoxia adaptation were detected.

•

Hypoxia affects muscle function by mediating energy metabolism via HIF pathway.

Proteomic Analysis of Human Serum for Patients at Different Pathological Stages of Hepatic Fibrosis

Background

Hepatic fibrosis is a severe liver disease that has threatened human health for a long time. In order to undergo timely and adequate therapy, it is important for patients to obtain an accurate diagnosis of fibrosis. Laboratory inspection methods have been efficient in distinguishing between advanced hepatic fibrosis stages (F3, F4), but the identification of early stages of fibrosis has not been achieved. The development of proteomics may provide us with a new direction to identify the stages of fibrosis.

Methods

We established serum proteomic maps for patients with hepatic fibrosis at different stages and identified differential expression of proteins between fibrosis stages through ultra-high-performance liquid chromatography tandem mass spectrometry proteomic analysis.

Results

From the proteomic profiles of the serum of patients with different stages of liver fibrosis, a total of 1,338 proteins were identified. Among three early fibrosis stages (control, F1, and F2), 55 differential proteins were identified, but no proteins simultaneously exhibited differential expression between control, F1, and F2. Differential proteins were detected in the comparison between different fibrosis stages. Significant differences were found between advanced fibrosis stages (F2-vs.-F3 and F3-vs.-F4) through a series of statistical analysis, including hierarchical clustering, Gene Ontology (GO) functional annotation, Kyoto Encyclopedia of Genes and Genomes pathway, and protein-protein interaction network analysis. The differential proteins identified by GO annotation were associated with biological processes (mainly platelet degranulation and cell adhesion), molecular functions, and cellular components.

Conclusions

All potential biomarkers identified between the stages of fibrosis could be key points in determining the fibrosis staging. The differences between early stages may provide a useful reference in addressing the challenge of early fibrosis staging.

Novel Corneal Protein Biomarker Candidates Reveal Iron Metabolic Disturbance in High Myopia Eyes

Myopia is a major public health concern with increasing global prevalence and is the leading cause of vision loss and complications. The potential role of the cornea, a substantial component of refractive power and the protective fortress of the eye, has been underestimated in the development of myopia. Our study acquired corneal stroma tissues from myopic patients undergoing femtosecond laser-assisted small incision lenticule extraction (SMILE) surgery and investigated the differential expression of circulating proteins between subjects with low and high myopia by means of high-throughput proteomic approaches—the quantitative tandem mass tag (TMT) labeling method and parallel reaction monitoring (PRM) validation. Across all corneal stroma tissue samples, a total of 2,455 proteins were identified qualitatively and quantitatively, 103 of which were differentially expressed between those with low and high myopia. The differentially abundant proteins (DAPs) between the groups of stroma samples mostly demonstrated catalytic activity and molecular function regulator and transporter activity and participated in metabolic processes, biological regulation, response to stimulus, and so forth. Pathway enrichment showed that mineral absorption, ferroptosis, and HIF-1 signaling pathways were activated in the human myopic cornea. Furthermore, TMT analysis and PRM validation revealed that the expression of ferritin light chain (FTL, P02792) and ferritin heavy chain (FTH1, P02794) was negatively associated with myopia development, while the expression of serotransferrin (TF, P02787) was positively related to myopia status. Overall, our results indicated that subjects with low and high myopia could have different proteomic profiles or signatures in the cornea. These findings revealed disturbances in iron metabolism and corneal oxidative stress in the more myopic eyes. Iron metabolic proteins could serve as an essential modulator in the pathogenesis of myopia.

Integrated Analysis of Liver Transcriptome, miRNA, and Proteome of Chinese Indigenous Breed Ningxiang Pig in Three Developmental Stages Uncovers Significant miRNA-mRNA-Protein Networks in Lipid Metabolism

Liver is an important metabolic organ of mammals. During each transitional period of life, liver metabolism is programmed by a complex molecular regulatory system for multiple physiological functions, many pathways of which are regulated by hormones and cytokines, nuclear receptors, and transcription factors. To gain a comprehensive and unbiased molecular understanding of liver growth and development in Ningxiang pigs, we analyzed the mRNA, microRNA (miRNA), and proteomes of the livers of Ningxiang pigs during lactation, nursery, and fattening periods. A total of 22,411 genes (19,653 known mRNAs and 2758 novel mRNAs), 1122 miRNAs (384 known miRNAs and 738 novel miRNAs), and 1123 unique proteins with medium and high abundance were identified by high-throughput sequencing and mass spectrometry. We show that the differences in transcriptional, post-transcriptional, or protein levels were readily identified by comparing different time periods, providing evidence that functional changes that may occur during liver development are widespread. In addition, we found many overlapping differentially expressed genes (DEGs)/differentially expressed miRNAs (DEMs)/differentially expressed proteins (DEPs) related to glycolipid metabolism in any group comparison. These overlapping DEGs/DEMs/DGPs may play an important role in functional transformation during liver development. Short Time-series Expression Miner (STEM) analysis revealed multiple expression patterns of mRNA, miRNA, and protein in the liver. Furthermore, several diverse key Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways, including immune defense, glycolipid metabolism, protein transport and uptake, and cell proliferation and development, were identified by combined analysis of DEGs and DGPs. A number of predicted miRNA–mRNA–protein pairs were found and validated by qRT-PCR and parallel reaction monitoring (PRM) assays. The results provide new and important information about the genetic breeding of Ningxiang pigs, which represents a foundation for further understanding the molecular regulatory mechanisms of dynamic development of liver tissue, functional transformation, and lipid metabolism.

Transcriptome analysis reveals the molecular mechanisms of heterosis on thermal resistance in hybrid abalone

BACKGROUND

Heterosis has been exploited for decades in different animals and crops due to it resulting in dramatic increases in yield and adaptability. Hybridization is a classical breeding method that can effectively improve the genetic characteristics of organisms through heterosis. Abalone has become an increasingly economically important aquaculture resource with high commercial value. However, due to changing climate, abalone is now facing serious threats of high temperature in summer. Interspecific hybrid abalone (Haliotis gigantea ♀ × H. discus hannai ♂, SD) has been cultured at large scale in southern China and has been shown high survival rates under heat stress in summer. Therefore, SD has become a good model material for heterosis research, but the molecular basis of heterosis remains elusive.

RESULTS

Heterosis in thermal tolerance of SD was verified through Arrhenius break temperatures (ABT) of cardiac performance in this study. Then RNA-Sequencing was conducted to obtain gene expression patterns and alternative splicing events at control temperature (20 °C) and heat stress temperature (30 °C). A total of 356 (317 genes), 476 (435genes), and 876 (726 genes) significantly diverged alternative splicing events were identified in H. discus hannai (DD), H. gigantea (SS), and SD in response to heat stress, respectively. In the heat stress groups, 93.37% (20,512 of 21,969) of the expressed genes showed non-additive expression patterns, and over-dominance expression patterns of genes account for the highest proportion (40.15%). KEGG pathway enrichment analysis showed that the overlapping genes among common DEGs and NAGs were significantly enriched in protein processing in the endoplasmic reticulum, mitophagy, and NF-κB signaling pathway. In addition, we found that among these overlap genes, 39 genes had undergone alternative splicing events in SD. These pathways and genes may play an important role in the thermal resistance of hybrid abalone.

CONCLUSION

More alternative splicing events and non-additive expressed genes were detected in hybrid under heat stress and this may contribute to its thermal heterosis. These results might provide clues as to how hybrid abalone has a better physiological regulation ability than its parents under heat stress, to increase our understanding of heterosis in abalone.

De Novo Transcriptome Assembly, Functional Annotation and SSR Marker Discovery of Qinling Takin (Budorcas taxicolor bedfordi)

The takin (Budorcas taxicolor) is an endemic ruminant species belonging to the bovine family. The International Union for Conservation of Nature (IUCN) has listed it as an endangered and vulnerable species. However, little is known about its molecular characterization since it lacks a reference genome. This study used RNA sequencing followed by de novo assembly, annotation and simple sequence repeats (SSRs) prediction to assess the transcriptome of Qinling takin (Budorcas taxicolor bedfordi) muscles. In total, 21,648 unigenes with an N50 and mean length of 1388 bp and 817 bp, respectively, were successfully detected and annotated against the public databases (NR, GO, KEGG, and EggNOG). Furthermore, 6222 SSRs were identified using the MIcroSAtellite (MISA) identification tool software. Taken together, these findings will provide valuable information for genetic, genomic, and evolutionary studies on takin.

Multiomics Analysis Reveals New Insights into the Apple Fruit Quality Decline under High Nitrogen Conditions

Excessive application of nitrogen (N) fertilizer is common in Chinese apple production. High N reduced the contents of soluble sugar and total flavonoids by 16.05 and 19.01%, respectively, resulting in poor fruit quality. Moreover, high N increased the total N and decreased the total C and C/N ratio of apple fruits. On the basis of the transcriptomic, proteomic, and metabolomic analyses, the global network was revealed. High N inhibited the accumulation of carbohydrates (sucrose, glucose, and trehalose) and flavonoids (rhamnetin-3-O-rutinoside, rutin, and trihydroxyisoflavone-7-O-galactoside) in fruits, and more C skeletons were used to synthesize amino acids and their derivatives (especially low C/N ratio, e.g., arginine) to be transferred to N metabolism. This study revealed new insights into the decline in soluble sugar and flavonoids caused by high N, and hub genes (MD07G1172700, MD05G1222800, MD16G1227200, MD01G1174400, and MD02G1207200) and hub proteins (PFK, gapN, and HK) were obtained.

Proteomic profiling analysis of postmenopausal osteoporosis and osteopenia identifies potential proteins associated with low bone mineral density

Postmenopausal osteoporosis (PMOP) is a major global public health concern and older women are more susceptible to experiencing fragility fractures. Our study investigated the associations between circulating proteins with bone mineral density (BMD) in postmenopausal women with or without low BMD (osteoporosis and osteopenia) using a tandem mass tag (TMT) labeling proteomic experiment and parallel reaction monitoring testing. Across all plasma samples, we quantitatively measured 1,092 proteins, and the OP and normal control (NC) samples were differentiated by principal component analysis and a partial least squares-discrimination analysis model based on the protein profiling data. The differentially abundant proteins between the low BMD and NC samples mostly exhibited binding, molecular function regulator, transporter and molecular transducer activity, and were involved in metabolic and cellular processes, stimulus response, biological regulation, immune system processes and so forth. TMT analysis and RRM validation indicated that the expression of protein Lysozyme C (P61626) was negatively related to BMD, while the expression of proteins Glucosidase (A0A024R592) and Protein disulfideisomerase A5 (Q14554) was positively related to BMD values. Collectively, our results suggest that postmenopausal women with low BMD have a different proteomic profile or signature. Protein alterations may play an important role in the pathogenesis of PMOP, and they may act as novel biomarkers and targets of therapeutic agents for this disease.