SMYD2-Methylated PPARγ Facilitates Hypoxia-Induced Pulmonary Hypertension by Activating Mitophagy.

BACKGROUND: Hyperproliferation of pulmonary arterial smooth muscle cells (PASMCs) and consequent pulmonary vascular remodeling are the crucial pathological features of pulmonary hypertension (PH). Protein methylation has been shown to be critically involved in PASMC proliferation and PH, but the underlying mechanism remains largely unknown. METHODS: PH animal models were generated by treating mice/rats with chronic hypoxia for 4 weeks. SMYD2-vTg mice (vascular smooth muscle cell-specific suppressor of variegation, enhancer of zeste, trithorax and myeloid Nervy DEAF-1 (deformed epidural auto-regulatory factor-1) domain-containing protein 2 transgenic) or wild-type rats and mice treated with LLY-507 (3-cyano-5-{2-[4-[2-(3-methylindol-1-yl)ethyl]piperazin-1-yl]-phenyl}-N-[(3-pyrrolidin-1-yl)propyl]benzamide) were used to investigate the function of SMYD2 (suppressor of variegation, enhancer of zeste, trithorax and myeloid Nervy DEAF-1 domain-containing protein 2) on PH development in vivo. Primary cultured rat PASMCs with SMYD2 knockdown or overexpression were used to explore the effects of SMYD2 on proliferation and to decipher the underlying mechanism. RESULTS: We demonstrated that the expression of the lysine methyltransferase SMYD2 was upregulated in the smooth muscle cells of pulmonary arteries from patients with PH and hypoxia-exposed rats/mice and in the cytoplasm of hypoxia-induced rat PASMCs. More importantly, targeted inhibition of SMYD2 by LLY-507 significantly attenuated hypoxia-induced pulmonary vascular remodeling and PH development in both male and female rats in vivo and reduced rat PASMC hyperproliferation in vitro. In contrast, SMYD2-vTg mice exhibited more severe PH phenotypes and related pathological changes than nontransgenic mice after 4 weeks of chronic hypoxia treatment. Furthermore, SMYD2 overexpression promoted, while SMYD2 knockdown suppressed, the proliferation of rat PASMCs by affecting the cell cycle checkpoint between S and G2 phases. Mechanistically, we revealed that SMYD2 directly interacted with and monomethylated PPARγ (peroxisome proliferator-activated receptor gamma) to inhibit the nuclear translocation and transcriptional activity of PPARγ, which further promoted mitophagy to facilitate PASMC proliferation and PH development. Furthermore, rosiglitazone, a PPARγ agonist, largely abolished the detrimental effects of SMYD2 overexpression on PASMC proliferation and PH. CONCLUSIONS: Our results demonstrated that SMYD2 monomethylates nonhistone PPARγ and inhibits its nuclear translocation and activation to accelerate PASMC proliferation and PH by triggering mitophagy, indicating that targeting SMYD2 or activating PPARγ are potential strategies for the prevention of PH.

Comprehensive analysis of differentially expressed mRNAs, circRNAs, and miRNAs and their ceRNA network in the testis of cattle-yak, yak, and cattle.

Cattle-yak is a hybrid offspring resulting from the crossbreeding of yak and cattle, and it exhibits substantial heterosis in production performance. However, male sterility in cattle-yak remains a concern. Reports suggest that noncoding RNAs are involved in the regulation of spermatogenesis. Therefore, in this study, we comprehensively compared testicular transcription profiles among cattle, yak, and cattle-yak. Numerous differentially expressed genes (DEGs), differentially expressed circRNAs (DECs), and differentially expressed miRNAs (DEMs) were identified in the intersection of two comparison groups, namely cattle versus cattle-yak and yak versus cattle-yak, with the number of DEGs, DECs, and DEMs being 4968, 360, and 59, respectively. The DEGs in cattle-yaks, cattle, and yaks were mainly associated with spermatogenesis, male gamete generation, and sexual reproduction. Concurrently, GO and KEGG analyses indicated that DEC host genes and DEM source genes were involved in the regulation of spermatogenesis. The construction of a potential competing endogenous RNA network revealed that some differentially expressed noncoding RNAs may be involved in regulating the expression of genes related to testicular spermatogenesis, including miR-423-5p, miR-449b, miR-34b/c, and miR-15b, as well as previously unreported miR-6123 and miR-1306, along with various miRNA-circRNA interaction pairs. This study serves as a valuable reference for further investigations into the mechanisms underlying male sterility in cattle-yaks.

Integrative analysis of Iso-Seq and RNA-seq data reveals transcriptome complexity and differential isoform in skin tissues of different hair length Yak.

BACKGROUND: The hair follicle development process is regulated by sophisticated genes and signaling networks, and the hair grows from the hair follicle. The Tianzhu white yak population exhibits differences in hair length, especially on the forehead and shoulder region. However, the genetic mechanism is still unclear. Isoform sequencing (Iso-seq) technology with advantages in long reads sequencing. Hence, we combined the Iso-seq and RNA-seq methods to investigate the transcript complexity and difference between long-haired yak (LHY) and normal-haired yak (NHY). RESULTS: The hair length measurement result showed a significant difference between LHY and NHY on the forehead and the shoulder (P-value < 0.001). The skin samples from the forehead and the shoulder of LHY and NHY were pooled for isoform sequencing (Iso-seq). We obtained numerous long transcripts, including novel isoforms, long non-coding RNA, alternative splicing events, and alternative polyadenylation events. Combined with RNA-seq data, we performed differential isoforms (DEIs) analysis between LHY and NHY. We found that some hair follicle and skin development-related DEIs, like BMP4, KRT2, IGF2R, and COL1A2 in the forehead skin; BMP1, KRT1, FGF5, COL2A1, and IGFBP5 in the shoulder skin. Enrichment analysis revealed that DEIs in both two comparable groups significantly participated in skin and hair follicle development-related pathways, like ECM-receptor interaction, focal adhesion, and PI3K-Akt signaling pathways. The results indicated that the hair follicle development of Tianzhu white yak may influence the hair length difference. Besides, the protein-protein interaction (PPI) network of DEIs showed COL2A1 and COL3A1 exhibited a high degree of centrality, and these two genes were suggested as potential candidates for the hair length growth of Tianzhu white yak. CONCLUSIONS: The results provided a comprehensive analysis of the transcriptome complexity and identified differential transcripts that enhance our understanding of the molecular mechanisms underlying the variation in hair length growth in Tianzhu white yak.

Genetic diversity, phylogeography, and maternal origin of yak (Bos grunniens).

BACKGROUND: There is no consensus as to the origin of the domestic yak (Bos grunniens). Previous studies on yak mitochondria mainly focused on mitochondrial displacement loop (D-loop), a region with low phylogenetic resolution. Here, we analyzed the entire mitochondrial genomes of 509 yaks to obtain greater phylogenetic resolution and a comprehensive picture of geographical diversity. RESULTS: A total of 278 haplotypes were defined in 509 yaks from 21 yak breeds. Among them, 28 haplotypes were shared by different varieties, and 250 haplotypes were unique to specific varieties. The overall haplotype diversity and nucleotide diversity of yak were 0.979 ± 0.0039 and 0.00237 ± 0.00076, respectively. Phylogenetic tree and network analysis showed that yak had three highly differentiated genetic branches with high support rate. The differentiation time of clades I and II were about 0.4328 Ma, and the differentiation time of clades (I and II) and III were 0.5654 Ma. Yushu yak is shared by all haplogroups. Most (94.70%) of the genetic variation occurred within populations, and only 5.30% of the genetic variation occurred between populations. The classification showed that yaks and wild yaks were first clustered together, and yaks were clustered with American bison as a whole. Altitude had the highest impact on the distribution of yaks. CONCLUSIONS: Yaks have high genetic diversity and yak populations have experienced population expansion and lack obvious phylogeographic structure. During the glacial period, yaks had at least three or more glacial refugia.

A single-cell transcriptomic atlas characterizes cell types and their molecular features in yak ovarian cortex.

The ovary as one of the most dynamic organs produces steroids to orchestrate female secondary sexual characteristics, harbors ovarian reserve for oocytes, releases mature oocytes for fertilization, and maintains pregnancy. Yak (Bos grunniens) is the only bovid animal that can adapt to the harsh climatic conditions on the Qinghai-Tibetan Plateau (altitudes of over 3000 m above sea level). However, the cellular atlas is composed of oocytes and other somatic cells, and their individual molecular characteristics remain to be elucidated in the yak ovary. Here, single-cell RNA sequencing (scRNA-seq) was performed to delineate the molecular signature of various cell types in the yak ovarian cortex. A cellular atlas of yak ovarian cortex was constructed successfully on the basis of the differentially expressed genes (DEGs) from the distinct cell types and their functional enrichment analysis, comprising endothelial cells, nature kill cells, stromal cells, smooth muscle cells, oocytes, macrophages, epithelial cells, and granulosa cells. Meanwhile, the signature genes were determined based on their expression specificity in each cell type. A cell-to-cell communication network was built in light of the differentially overexpressed ligand and receptor genes from each cell type. Further, the oocytes were subdivided into four subtypes based on their individual DEGs and the functional enrichment of the DEGs. FST and TOP2A were identified as maker genes for oocytes by immunostaining in the yak ovarian cortex. The cellular atlas reveals the biological characteristics of the ovarian cortex at the cellular molecular level and provides insights into female reproductive biology via cellular communications in the yak.

Identification of circRNA-associated ceRNA networks in the longissimus dorsi of yak under different feeding systems.

BACKGROUND: Yaks (Bos grunniens), prized for their ability to thrive in high-altitude environments, are indispensable livestock in the plateau region. Modifying their feeding systems holds significant promise for improving their growth and meat quality. Tenderness, a key determinant of yak meat quality and consumer appeal, is demonstrably influenced by dietary regimen. Indoor feeding regimes have been shown to enhance tenderness by lowering shear stress and optimizing pH values. CircRNAs, well-known modulators of circulatory function, also play a crucial role in skeletal muscle development across various animal species. However, their functional significance in yak skeletal muscle remains largely unexplored. RESULTS: In this study, we identified a total of 5,534 circRNAs within the longissimus dorsi muscle, and we found 51 differentially expressed circRNAs (20 up-regulated and 31 down-regulated) between the two feeding groups. Constructing a comprehensive ceRNA network illuminated intricate regulatory mechanisms, with PGP and circRNA_0617 converging on bta-miR-2285q, mirrored by KLF15/circRNA_0345/bta-miR-20b and CTSF/circRNA_0348/bta-miR-146a. These findings shed light on the potential of circRNAs to influence yak muscle development and meat quality, offering valuable insights for future research. CONCLUSIONS: This investigation unraveled a complex interaction network between circRNAs、mRNAs and miRNAs in yak skeletal muscle. We further elucidated the target genes regulated by these target genes within the network, offering valuable insights into the potential regulatory mechanisms governing muscle development and meat quality-related traits in yaks.

Age-dependent changes in the expression and localization of LYZL4, LYZL6 and PCNA during testicular development in the Ashidan yak.

Lysozyme like 4 (LYZL4), lysozyme like 6 (LYZL6) and proliferating cell nuclear antigen (PCNA) are implicated in the regulation of testicular function, but there was no research reported available on the expression patterns of LYZL4, LYZL6 and PCNA genes at different developmental stages of yak testes. In this study, we used the qRT-PCR, western blotting and immunohistochemistry estimated the LYZL4, LYZL6 and PCNA gene expression and protein lo-calization at different developmental stages of yak testes. The qPCR results showed that the mRNA expression of LYZL4, LYZL6 and PCNA genes significantly increased with age in the testes of yaks. Western blot results showed that the protein abundance of LYZL4, LYZL6 and PCNA in yak testes was significantly higher after puberty than before puberty. Furthermore, the results of immunohistochemistry indicated that LYZL4, LYZL6 and PCNA may be involved in the regulation of spermatogonia proliferation and Leydig cell function in immature testis. In adult yak testes, LYZL4, LYZL6 and PCNA may involve in the development of round spermatids and primary spermatocytes during testicular development. Our results indicated that LYZL4, LYZL6 and PCNA may be involved in the development of Sertoli cells, Leydig cells and gonocytes in yak testes.

A comprehensive study on the longissius dorsi muscle of Ashdan yaks under different feeding regimes based on transcriptomic and metabolomic analyses.

Yak is an important dominant livestock species at high altitude, and the growth performance of yak has obvious differences under different feeding methods. This experiment was conducted to compare the effects of different feeding practices on growth performance and meat quality of yaks through combined transcriptomic and metabolomic analyses. In terms of yak growth performance, compared with traditional grazing, in-house feeding can significantly improve the average daily weight gain, carcass weight and net meat weight of yaks; in terms of yak meat quality, in-house feeding can effectively improve the quality of yak meat. A combined transcriptomic and metabolomic analysis revealed 31 co-enriched pathways, among which arginine metabolism, proline metabolism and glycerophospholipid metabolism may be involved in the development of the longissimus dorsi muscle of yak and the regulation of meat quality-related traits. The experimental results increased our understanding of yak meat quality and provided data materials for subsequent deep excavation of the mechanism of yak meat quality.

Prophylactic Intra-abdominal Drainage is Associated With Lower Postoperative Complications in Patients With Crohn's Disease: A Randomized Controlled Trial.

BACKGROUND: Prophylactic intraoperative drains have been shown not superior for patients underwent intestinal surgery. However, for patients with Crohn's disease (CD), this needs further exploration. METHODS: In this pilot study, CD patients were randomly assigned to drain (n = 50) and no-drain (n = 50) groups. The primary endpoint was the rate of postoperative prolonged ileus (PPOI). The secondary endpoints were postoperative abdominal ascites, postoperative systemic inflammatory response syndrome (SIRS) and C-reactive protein (CRP) levels. RESULTS: The incidences of PPOI and postoperative abdominal ascites were significantly lower in the drain group (12% vs 44%; 0% vs 24%, both P < .05). Postoperative SIRS incidence and CRP levels were significantly increased in the no-drain group [36% vs 10%; 54.9 vs 34.3 mg/L, both P < .05]. In multivariate analysis, prophylactic drainage was the independent protective factor for PPOI and postoperative LOS. CONCLUSIONS: Prophylactic drainage may be associated with improved clinical outcomes in CD patients.

[Overview of new approaches to ß-thalassemia treatment].

Hemoglobinopathies are one of the most common single-gene genetic disorders globally, with approximately 1% to 5% of the global population carrying the mutated gene for thalassemia. Thalassemia are classified into transfusion-dependent thalassemia and non-transfusion-dependent thalassemia based on the need for blood transfusion. Traditional treatment modalities include blood transfusion, splenectomy, hydroxyurea therapy, and iron chelation therapy, which are now widely used for clinical treatment and constitute the main methods recommended in the ß-thalassemia treatment guidelines. However, there are multiple barriers and limitations to the application of these approaches, and there is an urgent need to explore new therapeutic approaches. With the in-depth study of the pathophysiological process of ß-thalassemia, a deeper understanding of the pathogenesis of the disease has been gained. It has been demonstrated that the pathogenesis of thalassemia is closely related to ineffective erythropoiesis (IE), imbalance in the ratio of α/ß-globin protein chains and iron overload. New therapeutic approaches are emerging for different pathogenic mechanisms. Among them, new drugs for the treatment of IE mainly include activin receptor II trap ligands, Janus kinase 2 inhibitors, pyruvate kinase activators, and glycine transporter protein 1 inhibitors. Correcting the imbalance in the hemoglobin chain is mainly due to emerging technologies such as bone marrow transplantation and gene editing. Measures in reducing iron overload are associated with inhibiting the activity of transferrin and hepcidin. These new approaches provide new ideas and options for the treatment and management of ß-thalassemia.

Oxaliplatin related lncRNAs prognostic models predict the prognosis of patients given oxaliplatin-based chemotherapy.

BACKGROUND: Oxaliplatin-based chemotherapy is the first-line treatment for colorectal cancer (CRC). Long noncoding RNAs (lncRNAs) have been implicated in chemotherapy sensitivity. This study aimed to identify lncRNAs related to oxaliplatin sensitivity and predict the prognosis of CRC patients underwent oxaliplatin-based chemotherapy. METHODS: Data from the Genomics of Drug Sensitivity in Cancer (GDSC) was used to screen for lncRNAs related to oxaliplatin sensitivity. Four machine learning algorithms (LASSO, Decision tree, Random-forest, and support vector machine) were applied to identify the key lncRNAs. A predictive model for oxaliplatin sensitivity and a prognostic model based on key lncRNAs were established. The published datasets, and cell experiments were used to verify the predictive value. RESULTS: A total of 805 tumor cell lines from GDSC were divided into oxaliplatin sensitive (top 1/3) and resistant (bottom 1/3) groups based on their IC50 values, and 113 lncRNAs, which were differentially expressed between the two groups, were selected and incorporated into four machine learning algorithms, and seven key lncRNAs were identified. The predictive model exhibited good predictions for oxaliplatin sensitivity. The prognostic model exhibited high performance in patients with CRC who underwent oxaliplatin-based chemotherapies. Four lncRNAs, including C20orf197, UCA1, MIR17HG, and MIR22HG, displayed consistent responses to oxaliplatin treatment in the validation analysis. CONCLUSION: Certain lncRNAs were associated with oxaliplatin sensitivity and predicted the response to oxaliplatin treatment. The prognostic models established based on the key lncRNAs could predict the prognosis of patients given oxaliplatin-based chemotherapy.

Tumor microenvironment characterization in colorectal cancer to identify prognostic and immunotherapy genes signature.

BACKGROUND: The tumor microenvironment (TME) plays a crucial role in tumorigenesis, progression, and therapeutic response in many cancers. This study aimed to comprehensively investigate the role of TME in colorectal cancer (CRC) by generating a TMEscore based on gene expression. METHODS: The TME patterns of CRC datasets were investigated, and the TMEscores were calculated. An unsupervised clustering method was used to divide samples into clusters. The associations between TMEscores and clinical features, prognosis, immune score, gene mutations, and immune checkpoint inhibitors were analyzed. A TME signature was constructed using the TMEscore-related genes. The results were validated using external and clinical cohorts. RESULTS: The TME pattern landscape was for CRC was examined using 960 samples, and then the TMEscore pattern of CRC datasets was evaluated. Two TMEscore clusters were identified, and the high TMEscore cluster was associated with early-stage CRC and better prognosis in patients with CRC when compared with the low TMEscore clusters. The high TMEscore cluster indicated elevated tumor cell scores and tumor gene mutation burden, and decreased tumor purity, when compared with the low TMEscore cluster. Patients with high TMEscore were more likely to respond to immune checkpoint therapy than those with low TMEscore. A TME signature was constructed using the TMEscore-related genes superimposing the results of two machine learning methods (LASSO and XGBoost algorithms), and a TMEscore-related four-gene signature was established, which had a high predictive value for discriminating patients from different TMEscore clusters. The prognostic value of the TMEscore was validated in two independent cohorts, and the expression of TME signature genes was verified in four external cohorts and clinical samples. CONCLUSION: Our study provides a comprehensive description of TME characteristics in CRC and demonstrates that the TMEscore is a reliable prognostic biomarker and predictive indicator for patients with CRC undergoing immunotherapy.

Application of artificial intelligence to digital-rapid on-site cytopathology evaluation during endoscopic ultrasound-guided fine needle aspiration: A proof-of-concept study.

BACKGROUND: During endoscopic ultrasound-guided fine needle aspiration (EUS-FNA), cytopathology with rapid on-site evaluation (ROSE) can improve diagnostic yield and accuracy. However, ROSE is unavailable in most Asian and European institutions because of the shortage of cytopathologists. Therefore, developing computer-assisted diagnostic tools to replace manual ROSE is crucial. Herein, we reported the validation of an artificial intelligence (AI)-based model (ROSE-AI model) to substitute manual ROSE during EUS-FNA. METHODS: A total of 467 digitized images from Diff-Quik (D&F)-stained EUS-FNA slides were divided into training (3642 tiles from 367 images) and internal validation (916 tiles from 100 images) datasets. The ROSE-AI model was trained and validated using training and internal validation datasets, respectively. The specificity was emphasized while developing the model. Then, we evaluated the AI model on a 693-image external dataset. We assessed the performance of the AI model to detect cancer cells (CCs) regarding the accuracy, sensitivity, specificity, positive predictive value (PPV), and negative predictive value (NPV). RESULTS: The ROSE-AI model achieved an accuracy of 83.4% in the internal validation dataset and 88.7% in the external test dataset. The sensitivity and PPV were 79.1% and 71.7% in internal validation dataset and 78.0% and 60.7% in external test dataset, respectively. CONCLUSION: We provided a proof of concept that AI can be used to replace manual ROSE during EUS-FNA. The ROSE-AI model can address the shortage of cytopathologists and make ROSE available in more institutes.

Proanthocyanidins suppress NLRP3 inflammasome and M1 macrophage polarization to alleviate severe acute pancreatitis in mice.

The role of reactive oxygen species (ROS) is crucial for the pathogenesis of acute pancreatitis (AP). Proanthocyanidins (PAs) have been confirmed to exert antioxidant activity. Our study aimed to determine whether PAs alleviated SAP via reducing ROS, suppressing NLRP3 inflammasome, and inhibiting M1 macrophage polarization. Our study investigated the protective effects of PAs on pancreatic histopathological injury using SAP mice. The effects of PAs on macrophages were investigated in inflammatory RAW 264.7 cells or mouse bone marrow-derived macrophages (BMDMs) induced by lipopolysaccharide (LPS). Immunofluorescence staining and/or western blot assay were employed to evaluate NLRP3 inflammasome in macrophages and pancreatic tissue. Cell counting kit-8 (CCK-8) was used to access effects of PAs on cell viability and cytometry flow was used to determine the effects of the PAs on the ROS levels of the RAW 264.7 cells. Then, we evaluated M1 macrophage polarization using flow cytometry or real-time quantitative polymerase chain reaction (RT-qPCR). PAs administration alleviated pancreatic inflammation in SAP mice. The PAs depressed NLRP3 inflammasome and inhibited M1 macrophage polarization in pancreatic tissue. We also found that the PAs showed no cellular toxicity but decreased ROS levels in RAW 264.7 cells, downregulated the NLRP3 inflammasome in the macrophages, and inhibited cell M1 macrophage polarization. Our study indicates the anti-inflammatory properties of the PAs on SAP mice by decreasing ROS levels, suppressing NLRP3 inflammasome, and M1 macrophage polarization.

Identification and profiling of microRNAs during yak's testicular development.

BACKGROUND: Normal testicular development is highly crucial for male reproduction and is a precondition for spermatogenesis that is the production of spermatozoa in the testes. MiRNAs have been implicated in several testicular biological processes, including cell proliferation, spermatogenesis, hormone secretion, metabolism and reproductive regulation. In the present study, we used deep sequencing data to study the functions of miRNAs during testicular development and spermatogenesis, by analyzing the expression patterns of small RNAs in 6-, 18- and 30-month-old yak testis tissues. RESULTS: A total of 737 known and 359 novel miRNAs were obtained from 6-, 18- and 30-month-old yak testes. In all, we obtained 12, 142 and 139 differentially expressed (DE) miRNAs in 30- vs. 18-, 18- vs. 6-, and 30- vs. 6-month-old testes, respectively. Gene Ontology (GO) annotation and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis of all DE miRNA target genes revealed BMP2, TGFB2, GDF6, SMAD6, TGFBR2 and other target genes as participants in different biological processes, including TGF-ß, GnRH, Wnt, PI3K-Akt, MAPK signaling pathways and several other reproductive pathways. In addition, quantitative reverse transcriptase-polymerase chain reaction (qRT-PCR) was used to detect the expression of seven randomly selected miRNAs in 6-, 18- and 30-month-old testes, and the results were consistent with the sequencing data. CONCLUSIONS: The differential expression of miRNAs in yak testes at different development stages was characterized and investigated using deep sequencing technology. We believe that the results will contribute to further understanding the functions of miRNAs in regulating the development of yak testes and improving the reproductive performance of male yaks.

Genome-Wide Landscape of mRNAs, lncRNAs, and circRNAs during Testicular Development of Yak.

Testicular development is a tightly regulated process in mammals. Understanding the molecular mechanisms of yak testicular development will benefit the yak breeding industry. However, the roles of different RNAs, such as mRNA, lncRNA, and circRNA in the testicular development of yak, are still largely unclear. In this study, transcriptome analyses were performed on the expression profiles of mRNAs, lncRNAs, and circRNAs in testis tissues of Ashidan yak at different developmental stages, including 6-months-old (M6), 18-months-old (M18), and 30-months-old (M30). A total of 30, 23, and 277 common differentially expressed (DE) mRNAs, lncRNAs, and circRNAs were identified in M6, M18, and M30, respectively. Furthermore, functional enrichment analysis showed that the common DE mRNAs during the entire developmental process were mainly involved in gonadal mesoderm development, cell differentiation, and spermatogenesis processes. Additionally, co-expression network analysis identified the potential lncRNAs related to spermatogenesis, e.g., TCONS_00087394 and TCONS_00012202. Our study provides new information about changes in RNA expression during yak testicular development, which greatly improves our understanding of the molecular mechanisms regulating testicular development in yaks.

Single-Cell Transcriptomics Analysis Reveals a Cell Atlas and Cell Communication in Yak Ovary.

Yaks (Bos grunniens) are the only bovine species that adapt well to the harsh high-altitude environment in the Qinghai-Tibetan plateau. However, the reproductive adaptation to the climate of the high elevation remains to be elucidated. Cell composition and molecular characteristics are the foundation of normal ovary function which determines reproductive performance. So, delineating ovarian characteristics at a cellular molecular level is conducive to elucidating the mechanism underlying the reproductive adaption of yaks. Here, the single-cell RNA-sequencing (scRNA-seq) was employed to depict an atlas containing different cell types with specific molecular signatures in the yak ovary. The cell types were identified on the basis of their specifically expressed genes and biological functions. As a result, a cellular atlas of yak ovary was established successfully containing theca cells, stromal cells, endothelial cells, smooth muscle cells, natural killer cells, macrophages, and proliferating cells. A cell-to-cell communication network between the distinct cell types was constructed. The theca cells were clustered into five subtypes based on their biological functions. Further, CYP11A1 was confirmed as a marker gene for the theca cells by immunofluorescence staining. Our work reveals an ovarian atlas at the cellular molecular level and contributes to providing insights into reproductive adaption in yaks.

The Effect of the Feeding System on Fat Deposition in Yak Subcutaneous Fat.

Fat deposition is very important to the growth and reproduction of yaks. In this study, the effect of the feeding system on fat deposition in yaks was explored by transcriptomics and lipidomics. The thickness of the subcutaneous fat in yaks under stall (SF) and graze feeding (GF) was evaluated. The transcriptomes and lipidomes of the subcutaneous fat in yaks under different feeding systems were detected by RNA-sequencing (RNA-Seq) and non-targeted lipidomics based on ultrahigh-phase liquid chromatography tandem mass spectrometry (UHPLC-MS), respectively. The differences in lipid metabolism were explored, and the function of differentially expressed genes (DEGs) was evaluated by gene ontology (GO) and Kyoto encyclopedia of genes and genome (KEGG) analysis. Compared with GF yaks, SF yaks possessed stronger fat deposition capacity. The abundance of 12 triglycerides (TGs), 3 phosphatidylethanolamines (PEs), 3 diglycerides (DGs), 2 sphingomyelins (SMs) and 1 phosphatidylcholine (PC) in the subcutaneous fat of SF and GF yaks was significantly different. Under the mediation of the cGMP-PKG signaling pathway, the blood volume of SF and GF yaks may be different, which resulted in the different concentrations of precursors for fat deposition, including non-esterified fatty acid (NEFA), glucose (GLU), TG and cholesterol (CH). The metabolism of C16:0, C16:1, C17:0, C18:0, C18:1, C18:2 and C18:3 in yak subcutaneous fat was mainly realized under the regulation of the INSIG1, ACACA, FASN, ELOVL6 and SCD genes, and TG synthesis was regulated by the AGPAT2 and DGAT2 genes. This study will provide a theoretical basis for yak genetic breeding and healthy feeding.

Single-Cell RNA Sequencing Reveals Atlas of Yak Testis Cells.

Spermatogenesis is a complex process that involves proliferation and differentiation of diploid male germ cells into haploid flagellated sperm and requires intricate interactions between testicular somatic cells and germ cells. The cellular heterogeneity of this process presents a challenge in analyzing the different cell types at various developmental stages. Single-cell RNA sequencing (scRNA-seq) provides a useful tool for exploring cellular heterogeneity. In this study, we performed a comprehensive and unbiased single-cell transcriptomic study of spermatogenesis in sexually mature 4-year-old yak using 10× Genomics scRNA-seq. Our scRNA-seq analysis identified six somatic cell types and various germ cells, including spermatogonial stem cells, spermatogonia, early-spermatocytes, late-spermatocytes, and spermatids in yak testis. Pseudo-timing analysis showed that Leydig and myoid cells originated from common progenitor cells in yaks. Moreover, functional enrichment analysis demonstrated that the top expressed genes in yak testicular somatic cells were significantly enriched in the cAMP signaling pathway, PI3K-Akt signaling pathway, MAPK signaling pathway, and ECM receptor interactions. Throughout the spermatogenesis process, genes related to spermatogenesis, cell differentiation, DNA binding, and ATP binding were expressed. Using immunohistochemical techniques, we identified candidate marker genes for spermatogonial stem cells and Sertoli cells. Our research provides new insights into yak spermatogenesis and the development of various types of cells in the testis, and presents more reliable marker proteins for in vitro culture and identification of yak spermatogonial stem cells in the later stage.

Species diversity and related ecology of chiggers on small mammals in a unique geographical area of Yunnan Province, southwest China.

Chiggers are common ectoparasites and the exclusive vector of scrub typhus. Based on previous investigations from a unique geographical area in Yunnan Province of southwest China, the Three Parallel Rivers Area, we retrospectively studied the species diversity and related ecology of chiggers on rodents and other small mammals. A very high species diversity of 120 chigger species was identified. Five dominant chigger species accounted for 59.4% (5238/8965) of total chiggers, and among them Leptotrombidium scutellare is the second major vector of scrub typhus in China. Species diversity of the chigger community fluctuates greatly in different altitudinal and latitudinal gradients. There are significant differences in species composition, species diversity and dominant species of chiggers among hosts with apparent community heterogeneity. Based on the species abundance distribution, the expected total number of chigger species was estimated to be 170, 50 more than the number of actually collected species; this further indicates a very high chigger species diversity in this area. The bipartite ecological network analysis revealed the intricate relationships between chigger and host species-positive and negative correlations existed among some species of dominant and vector chiggers.