Numerous genes involved in male reproduction regulate testis development and spermatogenesis. In this study, the testis tissue transcriptome was used to identify candidate genes and key pathways associated with fecundity in sheep. Histological analysis of testis tissue using hematoxylin–eosin (HE) routine staining was performed for two sheep breeds. Overall, 466 differentially expressed genes (DEGs) were identified between Hu sheep (HS) and Tibetan sheep (TS) through RNA sequencing technology (RNA-Seq), including 226 upregulated and 240 downregulated genes. Functional analysis showed that several terms and pathways, such as “protein digestion and absorption”, “cAMP signaling pathway”, “focal adhesion”, and “p53 signaling pathway” were closely related to testis development and spermatogenesis. Several genes (including COL1A1, COL1A2, COL3A1, SOX9, BCL2, HDC, and GGT5) were significantly enriched in these terms and pathways and might affect the reproduction of sheep by regulating the migration of spermatogenic cells, apoptosis of spermatogenic cells, and secretion of sterol hormones via testicular interstitial cells. Our results provide a theoretical basis for better understanding the molecular mechanisms of reproduction in sheep.

Alternative splicing (AS) allows the generation of multiple transcript variants from a single gene and affects biological processes by generating protein diversity in organisms. In total, 41,642 AS events corresponding to 9,924 genes were identified, and SE is the most abundant alternatively spliced type. The analysis of functional categories demonstrates that alternatively spliced differentially expressed genes (DEGs) were enriched in the MAPK signaling pathway and hypoxia-inducible factor 1 (HIF-1) signaling pathway. Proteoglycans in cancer between the normoxic (21% O₂, TN and LN) and hypoxic (2% O₂, TL and LL) groups, such as SLC2A1, HK1, HK2, ENO3, and PFKFB3, have the potential to rapidly proliferate alveolar type II epithelial (ATII) cells by increasing the intracellular levels of glucose and quickly divert to anabolic pathways by glycolysis intermediates under hypoxia. ACADL, EHHADH, and CPT1A undergo one or two AS types with different frequencies in ATII cells between TN and TL groups (excluding alternatively spliced DEGs shared between normoxic and hypoxic groups), and a constant supply of lipids might be obtained either from the circulation or de novo synthesis for better growth of ATII cells under hypoxia condition. MCM7 and MCM3 undergo different AS types between LN and LL groups (excluding alternatively spliced DEGs shared between normoxic and hypoxic groups), which may bind to the amino-terminal PER-SIM-ARNT domain and the carboxyl terminus of HIF-1α to maintain their stability. Overall, AS and expression levels of candidate mRNAs between Tibetan pigs and Landrace pigs revealed by RNA-seq suggest their potential involvement in the ATII cells grown under hypoxia conditions.

The Tibetan pig is an endemic economic animal in the plateau region of China, and has a unique adaptation mechanism to the plateau hypoxic environment. Research into microRNAs (miRNAs) involved in the mechanism underlying hypoxia adaptation of Tibetan pig is very limited. Therefore, we isolated alveolar type II epithelial (ATII) cells from the lungs of the Tibetan pig, cultured them in normoxia/hypoxia (21% O₂; 2% O₂) for 48 h, and performed high-throughput sequencing analysis. We identified a hypoxic stress-related ssc-miR-141 and predicted its target genes. The target genes of ssc-miR-141 were mainly enriched in mitogen-activated protein kinase (MAPK), autophagy-animal, and Ras signaling pathways. Further, we confirmed that PDCD4 may serve as the target gene of ssc-miR-141. Real-time quantitative polymerase chain reaction (RT-qPCR) analysis was performed to confirm the expression levels of ssc-miR-141 and PDCD4, and a dual-luciferase gene reporter system was used to verify the targeted linkage of ssc-miR-141 to PDCD4. The results showed that the expression level of ssc-miR-141 in the hypoxia group was higher than that in the normoxia group, while the expression level of PDCD4 tended to show the opposite trend and significantly decreased under hypoxia. These findings suggest that ssc-miR-141 is associated with hypoxia adaptation and provide a new insight into the role of miRNAs from ATII cells of Tibetan pig in hypoxia adaptation.