Testis development and spermatogenesis are strictly regulated by numbers of genes and non-coding genes. However, Long non-coding RNAs (lncRNAs) as key regulators in multitudinous biological processes have not been systematically identified in bovine testes during sexual maturation. In this study, we comprehensively analyzed lncRNA and mRNA expression profiling of six bovine testes at 3d after birth and 13 months by RNA sequencing. 23,735 lncRNAs and 22,118 mRNAs were identified, in which, 540 lncRNAs (P-value < 0.05) and 3,525 mRNAs (P-adjust < 0.05) were significantly differentially expressed (DE) between two stages. Correspondingly, the results of RT-qPCR analysis showed well correlation with the transcriptome data. Moreover, GO and KEGG enrichment analyses showed that DE genes and target genes of DE lncRNAs were enriched in spermatogenesis. Furthermore, we constructed lncRNA-gene interaction networks, consequently, 15 DE lncRNAs and 12 cis-target genes were involved. The target genes (SPATA16, TCF21, ZPBP, PACRG, ATP8B3, COMP, ACE and OSBP2) were found that association with bovine sexual maturation. In addition, the expression of lncRNAs and cis-target genes were detected in bovine Leydig cells, Sertoli cell, and spermatogonia. Our study identified and analyzed lncRNAs and mRNAs in testis tissues, suggesting that lncRNAs may regulate testis development and spermatogenesis. Our findings provided new insights for further investigation of biological function in bovine lncRNA.

Circular RNA (circRNA) is a new class of non-coding RNA that has recently attracted researchers' interest. Studies have demonstrated that circRNA can function as microRNA sponges or competing endogenous RNAs. Although circRNA has been explored in some species and tissues, the genetic basis of testis development and spermatogenesis in cattle remains unknown. We performed ribo-depleted total RNA-Seq to detect circRNA expression profiles of neonatal (one week old) and adult (4 years old) Qinchuan cattle testes. We obtained 91 112 596 and 80 485 868 clean reads and detected 21 753 circRNAs. A total of 4248 circRNAs were significantly differentially expressed between neonatal and adult cattle testes. Among these circRNAs, 2225 were upregulated, and 2023 were downregulated in adult cattle testis. Genomic feature, length distribution and other characteristics of the circRNAs in cattle testis were studied. Moreover, Gene Ontology and KEGG pathway analyses were performed for source genes of circRNAs. These source genes were mainly involved in tight junction, adherens junction, TGFβ signalling pathway and reproduction, such as PIWIL1, DPY19L2, SLC26A8, IFT81, SMC1B, IQCG and TTLL5 . CircRNA expression patterns were validated by RT-qPCR. Our discoveries provide a solid foundation for the identification and characterization of key circRNAs involved in testis development or spermatogenesis.

Skeletal muscle is a complex heterogeneous tissue and characterizing its cellular heterogeneity and transcriptional and epigenetic signatures are important for understanding the details of its ontogeny. In our study, we applied scRNA‐seq and scATAC‐seq to investigate the cell types, molecular features, transcriptional and epigenetic regulation, and patterns of developing bovine skeletal muscle from gestational, lactational and adult stages. Detailed molecular analyses were used to dissect cellular heterogeneity, and we deduced the differentiation trajectory of myogenic cells and uncovered their dynamic gene expression profiles. SCENIC analysis was performed to demonstrate key regulons during cell fate decisions. We explored the future expression states of these heterogeneous cells by RNA velocity analysis and found extensive networks of intercellular communication using the toolkit CellChat. Moreover, the transcriptomic and chromatin accessibility modalities were confirmed to be highly concordant, and integrative analysis of chromatin accessibility and gene expression revealed key transcriptional regulators acting during myogenesis. In bovine skeletal muscle, by scRNA‐seq and scATAC‐seq analysis, different cell types such as adipocytes, endothelial cells, fibroblasts, lymphocytes, monocytes, pericyte cells and eight skeletal myogenic subpopulations were identified at the three developmental stages. The pseudotime trajectory exhibited a distinct sequential ordering for these myogenic subpopulations and eight distinct gene clusters were observed according to their expression pattern. Moreover, specifically expressed TFs (such as MSC, MYF5, MYOD1, FOXP3, ESRRA, BACH1, SIX2 and ATF4) associated with muscle development were predicted, and likely future transcriptional states of individual cells and the developmental dynamics of differentiation among neighbouring cells were predicted. CellChat analysis on the scRNA‐seq data set then classified many ligand–receptor pairs among these cell clusters, which were further categorized into significant signalling pathways, including BMP, IGF, WNT, MSTN, ANGPTL, TGFB, TNF, VEGF and FGF. Finally, scRNA‐seq and scATAC‐seq results were successfully integrated to reveal a series of specifically expressed TFs that are likely to be candidates for the promotion of cell fate transition during bovine skeletal muscle development. Overall, our results outline a single‐cell dynamic chromatin/transcriptional landscape for normal bovine skeletal muscle development; these provide an important resource for understanding the structure and function of mammalian skeletal muscle, which will promote research into its biology.

Spermatogenesis is a complex biological process regulated by well-coordinated gene regulation, including MicroRNAs (miRNAs). miRNAs are endogenous non-coding ribonucleic acids (ncRNAs) that mainly regulate the gene expression at post-transcriptional levels. Several studies have reported miRNAs expression in bull sperm and the process of spermatogenic arrest in cattle and yak. However, studies for the identification of differential miRNA expression and its mechanisms during the developmental stages of testis still remain uncertain. In the current study, we comprehensively analyzed the expression of miRNA in bovine testes at neonatal (3 days after birth, n = 3) and mature (13 months, n = 3) stages by RNA-seq. Moreover, the role of bta-miR-146b was also investigated in regulating the proliferation and apoptosis of bovine male germline stem cells (mGSCs) followed by a series of experiments. A total of 652 miRNAs (566 known and 86 novel miRNAs) were identified, whereas 223 miRNAs were differentially expressed between the two stages. Moreover, an elevated expression level of bta-miR-146b was found in bovine testis among nine tissues, and the functional studies indicated that the overexpression of bta-miR-146b inhibited the proliferation of bovine mGSCs and promoted apoptosis. Conversely, regulation of bta-miR-146b inhibitor promoted bovine mGSCs proliferation. This study provides a basis for understanding the regulation roles of miRNAs in bovine testis development and spermatogenesis.

In the biological process of testicular spermatogenesis, the expression and interaction of many genes are regulated by microRNAs (miRNAs). However, comparisons of miRNA expression between descended testes (DTs) and undescended testes (UDTs) are rarely done in horses. In this study, we selected two UDTs (CKY2b and GU4b) from Chakouyi (CKY) and Guanzhong (GU) horses and eight DTs (GU1–3, CKY1, CKY3, CKY2a, GU4a, and GU5). Three groups were compared to evaluate expression patterns of testicular miRNA in stallion testes. Group 1 compared normal CKY horses and GU horses (CKY1 and CKY3 vs. GU1–3). Group 2 (CKY2a and GU4a (DTs) vs. CKY2b and GU4b (UDTs)) and group 3 (GU1–3, CKY1, CKY3 (DTs) vs. CKY2b and GU4b (UDTs)) compared the expression levels in unilateral retained testes to normal testes. The results show that 42 miRNAs (7 upregulated and 35 downregulated) had significantly different expression levels in both comparisons. The expression levels of eca-miR-545, eca-miR-9084, eca-miR-449a, eca-miR-9024, eca-miR-9121, eca-miR-8908e, eca-miR-136, eca-miR-329b, eca-miR-370, and eca-miR-181b were further confirmed by quantitative real-time PCR assay. The target genes of differentially expressed miRNAs in three comparisons were predicted, and the functions were annotated. The putative target genes of the 42 co-differentially expressed miRNAs were annotated to 15 functional terms, including metal ion binding, GTPase activator activity, zinc ion binding, intracellular, cytoplasm, and cancer pathways, and osteoclast differentiation. Our data indicate that the differentially expressed miRNAs in undescended testis suggests a potential role in male fertility and a relationship with cryptorchidism in horses. The discovery of miRNAs in stallion testes might contribute to a new direction in the search for biomarkers of stallion fertility.

Donkeys (Equus asinus) are an important farm animal. After long-term natural and artificial selection, donkeys now exhibit a variety of body sizes and production performance values. In this study, six donkey breeds, representing different regions and phenotypes, were used for second-generation resequencing. The sequencing results revealed more than seven million single nucleotide variants (SNVs), with an average of more than four million SNVs per species. We combined two methods, Z-transformed heterozygosity (ZHp) and unbiased estimates of pairwise fixation index (di) values, to analyze the signatures of selection. We mapped 11 selected regions and identified genes associated with coat color, body size, motion capacity, and high-altitude adaptation. These candidate genes included staining (ASIP and KITLG), body type (ACSL4, BCOR, CDKL5, LCOR, NCAPG, and TBX3), exercise (GABPA), and adaptation to low-oxygen environments (GLDC and HBB). We also analyzed the SNVs of the breed-specific genes for their potential functions and found that there are three varieties in the conserved regions with breed-specific mutation sites. Our results provide data to support the establishment of the donkey SNV chip and reference information for the utilization of the genetic resources of Chinese domestic donkeys.