In the genomes of diploid organisms, runs of homozygosity (ROH), consecutive segments of homozygosity, are extended. ROH can be applied to evaluate the inbreeding situation of individuals without pedigree data and to detect selective signatures via ROH islands. We sequenced and analyzed data derived from the whole-genome sequencing of 97 horses, investigated the distribution of genome-wide ROH patterns, and calculated ROH-based inbreeding coefficients for 16 representative horse varieties from around the world. Our findings indicated that both ancient and recent inbreeding occurrences had varying degrees of impact on various horse breeds. However, recent inbreeding events were uncommon, particularly among indigenous horse breeds. Consequently, the ROH-based genomic inbreeding coefficient could aid in monitoring the level of inbreeding. Using the Thoroughbred population as a case study, we discovered 24 ROH islands containing 72 candidate genes associated with artificial selection traits. We found that the candidate genes in Thoroughbreds were involved in neurotransmission (CHRNA6, PRKN, and GRM1), muscle development (ADAMTS15 and QKI), positive regulation of heart rate and heart contraction (HEY2 and TRDN), regulation of insulin secretion (CACNA1S, KCNMB2, and KCNMB3), and spermatogenesis (JAM3, PACRG, and SPATA6L). Our findings provide insight into horse breed characteristics and future breeding strategies.

The regulatory role of non-CpG methylation in mammals has been important in whole-genome bisulfite sequencing. It has also been suggested that non-CpG methylation regulates gene expression to affect the development and health of mammals. However, the dynamic regulatory mechanisms of genome-wide, non-CpG methylation during testicular development still require intensive study. In this study, we analyzed the dataset from the whole-genome bisulfite sequencing (WGBS) and the RNA-seq of precocious porcine testicular tissues across two developmental stages (1 and 75 days old) in order to explore the regulatory roles of non-CpG methylation. Our results showed that genes regulated by non-CpG methylation affect the development of testes in multiple pathways. Furthermore, several hub genes that are regulated by non-CpG methylation during testicular development—such as VEGFA, PECAM1, and FZD7—were also identified. We also found that the relative expression of FZD7 was downregulated by the zebularine-induced demethylation of the first exon of FZD7. This regulatory relationship was consistent with the results of the WGBS and RNA-seq analysis. The immature porcine Sertoli cells were transfected with RNAi to mimic the expression patterns of FZD7 during testicular development. The results of the simulation test showed that cell proliferation was significantly impeded and that cell cycle arrest at the G2 phase was caused by the siRNA-induced FZD7 inhibition. We also found that the percentage of early apoptotic Sertoli cells was decreased by transfecting them with the RNAi for FZD7. This indicates that FZD7 is an important factor in linking the proliferation and apoptosis of Sertoli cells. We further demonstrated that Sertoli cells that were treated with the medium collected from apoptotic cells could stimulate proliferation. These findings will contribute to the exploration of the regulatory mechanisms of non-CpG methylation in testicular development and of the relationship between the proliferation and apoptosis of normal somatic cells.

RNA N⁶-methyladenosine (m⁶A) modification is one of the principal post-transcriptional modifications and plays a dynamic role in testicular development and spermatogenesis. However, the role of m⁶A in porcine testis is understudied. Here, we performed a comprehensive analysis of the m⁶A transcriptome-wide profile in Shaziling pig testes at birth, puberty, and maturity. We analyzed the total transcriptome m⁶A profile and found that the m⁶A patterns were highly distinct in terms of the modification of the transcriptomes during porcine testis development. We found that key m⁶A methylated genes (AURKC, OVOL, SOX8, ACVR2A, and SPATA46) were highly enriched during spermatogenesis and identified in spermatogenesis-related KEGG pathways, including Wnt, cAMP, mTOR, AMPK, PI3K-Akt, and spliceosome. Our findings indicated that m⁶A methylations are involved in the complex yet well-organized post-transcriptional regulation of porcine testicular development and spermatogenesis. We found that the m⁶A eraser ALKBH5 negatively regulated the proliferation of immature porcine Sertoli cells. Furthermore, we proposed a novel mechanism of m⁶A modification during testicular development: ALKBH5 regulated the RNA methylation level and gene expression of SOX9 mRNA. In addition to serving as a potential target for improving boar reproduction, our findings contributed to the further understanding of the regulation of m⁶A modifications in male reproduction.