A comparative analysis of differentially expressed mRNAs, miRNAs and circRNAs provides insights into the key genes involved in the high-altitude adaptation of yaks

Pangenome and multi-tissue gene atlas provide new insights into the domestication and highland adaptation of yaks

BACKGROUND

The genetic diversity of yak, a key domestic animal on the Qinghai-Tibetan Plateau (QTP), is a vital resource for domestication and breeding efforts. This study presents the first yak pangenome obtained through the de novo assembly of 16 yak genomes.

RESULTS

We discovered 290 Mb of nonreference sequences and 504 new genes. Our pangenome-wide presence and absence variation (PAV) analysis revealed 5,120 PAV-related genes, highlighting a wide range of variety-specific genes and genes with varying frequencies across yak populations. Principal component analysis (PCA) based on binary gene PAV data classified yaks into three new groups: wild, domestic, and Jinchuan. Moreover, we proposed a 'two-haplotype genomic hybridization model' for understanding the hybridization patterns among breeds by integrating gene frequency, heterozygosity, and gene PAV data. A gene PAV-GWAS identified a novel gene (BosGru3G009179) that may be associated with the multirib trait in Jinchuan yaks. Furthermore, an integrated transcriptome and pangenome analysis highlighted the significant differences in the expression of core genes and the mutational burden of differentially expressed genes between yaks from high and low altitudes. Transcriptome analysis across multiple species revealed that yaks have the most unique differentially expressed mRNAs and lncRNAs (between high- and low-altitude regions), especially in the heart and lungs, when comparing high- and low-altitude adaptations.

CONCLUSIONS

The yak pangenome offers a comprehensive resource and new insights for functional genomic studies, supporting future biological research and breeding strategies.

Effect of Different Feeds on the Fungi Microbiome of Suffolk Crossed with Tibetan Sheep

The gut microbiome plays an important role in the metabolism, nutrient absorption and immunocompetency of animals. The dynamics of the microbiota can be influenced by modulatory factors that involve nutrition, environment, health, diseases, etc. Few reports have been documented regarding the effects of different feeds on the fungi microbiome of Suffolk crossed with Tibetan sheep. A total of 30 Suffolk crossed with Tibetan sheep (ST sheep) were selected for the study and randomly divided into five equal groups (n = 6): AZ, BZ, CZ, DZ and EZ. Group AZ was fed with alfalfa and oat grass, whereas group BZ was fed with mixture of concentrated feed, alfalfa and oat grass. Groups CZ, DZ and EZ were fed with concentrated feed #1, #2 and #3, respectively. All experimental animals were fed twice a day for four months, and rectum samples were collected for microbiota analysis. Results revealed that 2,781,461 raw reads and 2,333,239 clean reads were achieved in the ST sheep. When compared with the sheep of groups AZ and BZ (164), the shared amplicon sequence variants (ASVs) between AZ and CZ (109), AZ (113) and DZ (118) as well as AZ along with EZ were fewer. Conspicuous different phyla (8) and genera (56) were examined and compared with free-range sheep in AZ. Genera including Xeromyces, Kazachstania, Cordyceps, Rhodotorula, Pichia, Spor, etc. were found higher in animals in the CZ, DZ and EZ groups. The results of this study provide new insights regarding the effects of different feeds on the fungi microbiome of sheep farmed on the plateau. We concluded that the differences in feed in Suffolk crossed with Tibetan sheep altered their gut microbiota.

Study of Transcriptomic Analysis of Yak (Bos grunniens) and Cattle (Bos taurus) Pulmonary Artery Smooth Muscle Cells under Oxygen Concentration Gradients and Differences in Their Lung Histology and Expression of Pyruvate Dehydrogenase Kinase 1-Related Factors

The aim of this study was to investigate the molecular mechanisms by which hypoxia affects the biological behavior of yak PASMCs, the changes in the histological structure of yak and cattle lungs, and the relationships and regulatory roles that exist regarding the differences in the distribution and expression of PDK1 and its hypoxia-associated factors screened for their role in the adaptation of yak lungs to the plateau hypoxic environment. The results showed that, at the level of transcriptome sequencing, the molecular regulatory mechanisms of the HIF-1 signaling pathway, glucose metabolism pathway, and related factors (HK2/PGK1/ENO1/ENO3/ALDOC/ALDOA) may be closely related to the adaptation of yaks to the hypoxic environment of the plateau; at the tissue level, the presence of filled alveoli and semi-filled alveoli, thicker alveolar septa and basement membranes, a large number of erythrocytes, capillary distribution, and collagen fibers accounted for all levels of fine bronchioles in the lungs of yaks as compared to cattle. A higher percentage of goblet cells was found in the fine bronchioles of yaks, and PDK1, HIF-1α, and VEGF were predominantly distributed and expressed in the monolayers of ciliated columnar epithelium in the branches of the terminal fine bronchioles of yak and cattle lungs, with a small amount of it distributed in the alveolar septa; at the molecular level, the differences in PDK1 mRNA relative expression in the lungs of adult yaks and cattle were not significant (p > 0.05), the differences in HIF-1α and VEGF mRNA relative expression were significant (p < 0.05), and the expression of PDK1 and HIF-1α proteins in adult yaks was stronger than that in adult cattle. PDK1 and HIF-1α proteins were more strongly expressed in adult yaks than in adult cattle, and the difference was highly significant (p < 0.01); the relative expression of VEGF proteins was not significantly different between adult yaks and cattle (p > 0.05). The possible regulatory relationship between the above results and the adaptation of yak lungs to the plateau hypoxic environment paves the way for the regulatory mechanisms of PDK1, HIF-1α, and VEGF, and provides basic information for studying the mechanism of hypoxic adaptation of yaks in the plateau. At the same time, it provides a reference for human hypoxia adaptation and a target for the prevention and treatment of plateau diseases in humans and plateau animals.

Gut flora alterations among aquatic firefly Aquatica leii inhabiting various dissolved oxygen in fresh water

Knowledge about the impact of different dissolved oxygen (DO) on the composition and function of gut bacteria of aquatic insects is largely unknown. Herein, we constructed freshwater environments with different DOs (hypoxia: 2.50 ± 0.50, normoxia: 7.00 ± 0.50, and hyperoxia: 13.00 ± 0.50 mg/L) where aquatic firefly Aquatica leii larvae lived for three months. Their gut flora was analyzed using the combination of 16S rRNA amplicon sequencing and metagenomics. The results showed no difference in alpha diversity of the gut flora between A. leii inhabiting various DOs. However, the relative abundance of several bacterial lineages presented significant changes, such as Pseudomonas. In addition, bacterial genes with an altered relative abundance in response to various DOs were primarily related to metabolism. The alteration of these functions correlated with the DO change. This is the first to uncover structure of gut flora under various DOs in aquatic insect larvae.

Regulatory effects of circular RNA on hypoxia adaptation in chicken embryos

The Tibetan chicken, a native breed of the Tibetan plateau, is adapted to the high-altitude and hypoxic environment of the plateau. As endogenous molecules, circular RNAs (circRNAs) have been shown to play an important role in the adaptation to hypoxic environments and regulation of angiogenesis. In this study, highland Tibetan and lowland Chahua chicken eggs were incubated in a hypoxic environment and the chorionic allantoic membrane was collected for Ribo-Zero RNA sequencing. A total of 1,414 circRNAs, mostly derived from exons, were identified. Of these, 93 differentially expressed circRNAs were detected between Tibetan and Chahua chickens. Combined with the differentially expressed miRNAs and mRNAs identified in our previous study, we identified four circRNAs (circBRD1, circPRDM2, circPTPRS, and circDENND4C). These circRNAs may act as competing endogenous RNA (ceRNA) to upregulate APOA1 expression by absorbing novel_miR_589, thereby regulating angiogenesis and affecting hypoxia adaptation in chicken embryos. The regulatory circRNAs/novel_miR_589/APOA1 axis provides valuable evidence for a better understanding of the specific functions and molecular mechanisms of circRNAs in plateau hypoxia adaptation in Tibetan chickens.

Combined transcriptome and proteome analysis of yak PASMCs under hypoxic and normoxic conditions

Background

Yaks are animals that have lived in plateau environments for generations. Yaks can adapt to the hypoxic plateau environment and also pass this adaptability on to the next generation. The lungs are the most important respiratory organs for mammals to adapt to their environment. Pulmonary artery smooth muscle cells play an important role in vascular remodeling under hypoxia, but the genetic mechanism underpinning the yak's ability to adapt to challenging plateau conditions is still unknown.

Methods

A tandem mass tag (TMT) proteomics study together with an RNA-seq transcriptome analysis were carried out on pulmonary artery smooth muscle cells (PASMCs) that had been grown for 72 hours in both normoxic (20% O2) and hypoxic (1% O2) environments. RNA and TP (total protein) were collected from the hypoxic and normoxic groups for RNA-seq transcriptome sequencing and TMT marker protein quantification, and RT-qPCR validation was performed.

Results

A total of 17,711 genes and 6,859 proteins were identified. Further, 5,969 differentially expressed genes (DEGs) and 531 differentially expressed proteins (DEPs) were identified in the comparison group, including 2,924 and 186 upregulated genes and proteins and 3,045 and 345 down-regulated genes and proteins, respectively. The transcriptomic and proteomic analyses revealed that 109 DEGs and DEPs were highly positively correlated, with 77 genes showing the same expression trend. Nine overlapping genes were identified in the HIF-1 signaling pathway, glycolysis / gluconeogenesis, central carbon metabolism in cancer, PPAR signaling pathway, AMPK signaling pathway, and cholesterol metabolism (PGAM1, PGK1, TPI1, HMOX1, IGF1R, OLR1, SCD, FABP4 and LDLR), suggesting that these differentially expressed genes and protein functional classifications are related to the hypoxia-adaptive pathways. Overall, our study offers abundant data for further analysis of the molecular mechanisms in yak PASMCs and their adaptability to different oxygen concentrations.

Circular RNA Profiling Identifies Novel circPPARA that Promotes Intramuscular Fat Deposition in Pigs

Intramuscular fat (IMF) content plays an important role in pork quality. Circular RNAs (circRNAs) implicate various biological processes; however, the regulatory mechanisms and functions of circRNAs in porcine IMF remains elusive. Hence, the study assessed the circRNA expression profiling in the longissimus dorsi muscle of pigs with high (H) and low (L) IMF content to unravel their regulatory functions in improving meat quality. The RNA sequencing analysis identified 29,732 circRNAs from six sampled pigs, most of which were exon-derived. In the muscle, 336 were differentially expressed (DE) between the H and L IMF groups; 196 circRNAs were upregulated, and 140 were downregulated. Subsequent qRT-PCR validation of 10 DE circRNAs revealed expression patterns consistent with the RNA-seq data. Gene ontology and KEGG enrichment analysis revealed that most significantly enriched DE circRNAs' host genes were linked to lipid metabolism and adipogenesis processes. The circRNA-miRNA regulatory network analysis found several circRNAs targeting miRNAs associated with adipogenesis. Finally, a novel circRNA, circPPARA, was identified with the expression positively correlated with the IMF content. Detailed analysis revealed that circPPARA was formed via head-to-tail splicing and was more stable than the linear PPARA, predominantly located in the cytoplasm. Functional studies using overexpression and siRNA constructs demonstrated that circPPARA promotes differentiation and hinders the proliferation of porcine intramuscular preadipocytes. Moreover, the dual-luciferase assay revealed that circPPARA adsorbed miR-429 and miR-200b, thereby promoting intramuscular adipogenesis in pigs. Our results identified a candidate circRNA, circPPARA, that affects porcine IMF content. The study provides knowledge of the regulatory functions of circRNAs in intramuscular adipogenesis and abundant resource for future research on circRNAs in pigs.

Integrative Analysis of the lncRNA-Associated ceRNA Regulatory Network Response to Hypoxia in Alveolar Type II Epithelial Cells of Tibetan Pigs

The function of alveolar type II epithelial (ATII) cells is severely hampered by oxygen deficiency, and understanding the regulatory mechanisms controlling responses to hypoxia may assist in relieving injury induced by hypoxia. In this study, we cultured ATII cells from Tibetan pigs and Landrace pigs under hypoxic and normoxic environments to screen for differentially expressed (DE) lncRNAs, DEmiRNAs, and construct their associated ceRNA regulatory networks in response to hypoxia. Enrichment analysis revealed that target genes of DElncRNAs of Tibetan pigs and Landrace pig between the normoxic (TN, LN) and hypoxic (TL, LL) groups significantly enriched in the proteoglycans in cancer, renal cell carcinoma, and erbB signaling pathways, while the target genes of DEmiRNAs were significantly enriched in the axon guidance, focal adhesion, and mitogen-activated protein kinase (MAPK) signaling pathways. Hypoxia induction was shown to potentially promote apoptosis by activating the focal adhesion/PI3K-Akt/glycolysis pathway. The ssc-miR-20b/MSTRG.57127.1/ssc-miR-7-5p axis potentially played a vital role in alleviating hypoxic injury by regulating ATII cell autophagy under normoxic and hypoxic conditions. MSTRG.14861.4-miR-11971-z-CCDC12, the most affected axis, regulated numerous RNAs and may thus regulate ATII cell growth in Tibetan pigs under hypoxic conditions. The ACTA1/ssc-miR-30c-3p/MSTRG.23871.1 axis is key for limiting ATII cell injury and improving dysfunction and fibrosis mediated by oxidative stress in Landrace pigs. Our findings provide a deeper understanding of the lncRNA/miRNA/mRNA regulatory mechanisms of Tibetan pigs under hypoxic conditions.