Characterization of N6-methyladenosine in cattle-yak testis tissue

Introduction

The yak (Bos grunniens) is a highly recognizable symbol of the Qinghai-Tibet Plateau (QTP) and nearby high-altitude regions where it is endemic (1, 2), primarily residing at altitudes of 2,500–6,000 m (3). Yaks are well-suited to surviving in the harsh climate of the QTP, and can readily tolerate environmental stressors including hypoxia, cold, and local diseases (3), serving as the primary form of livestock for humans living in this area (4). Highlanders use yaks both as pack animals and as a source of meat, milk, fuel, and fur (5). They are also central components of efforts to conserve agrobiodiversity, promote socio-economic development, and maintain rangeland ecosystems in these plateau regions (6). Relative to other cattle species, however, yaks exhibit lower rates of growth, reduced production performance (7), lower levels of fertility, and late sexual maturity (8). In most cases, yaks present with short-lived, relatively inconspicuous estrus behaviors beginning at 3–4 years of age, producing just once every 2 years or twice every 3 years (9).

Interspecific hybridization can play a central role in the adaptation and evolution of particular species in nature, with heterosis often resulting from the hybridization of genetically distinct animals (10). In an effort to improve yak production performance, researchers and breeders have sought to combine the excellent productivity of other cattle species with the QTP-adapted traits of yaks, such as cold tolerance, through interspecific hybridization (11). The resultant cattle-yaks, which are a hybrid of yaks and conventional cattle (Bos taurus), possess the high productivity of cattle and the adaptability of yaks to high-altitude conditions (12). Cattle-yaks exhibit excellent heterosis with respect to growth speed, disease resistance, drought tolerance, meat quality, and adaptability to the plateau environment (13, 14). Meat from these cattle-yak hybrids, which are larger than yaks, contains higher levels of protein and lower fat content as compared to yak meat, thus fulfilling a need for the production of healthier high-quality foods fit for public consumption (15). However, this heterosis cannot be effectively leveraged owing to the fact that F1 male cattle-yaks exhibit an inability to produce sperm, resulting in sterility (16). Despite the normal development of external cattle-yak reproductive organs, spermatogenesis in these hybrids is blocked after the primary spermatocyte stage, with only a relatively small number of autosomes from a few spermatocytes being evident in the meiotic synaptonemal complex (SC) in these animals (17). Prior research efforts have leveraged multifunctional strategies including transcriptomic, proteomic, genetic, physiological, and endocrinological approaches in an effort to clarify the mechanistic basis for male hybrid sterility (12). However, the mechanisms that ultimately drive male cattle-yak infertility have yet to be firmly established.

Spermatogenesis is a complex process wherein diploid spermatogonial stem cells (SSCs) undergo differentiation to produce haploid spermatozoa through tightly regulated mitotic, meiotic, and spermatogenic processes that are controlled at the transcriptional, post-transcriptional, and translational levels (18). Epigenetic regulation can also shape the development of the male reproductive system, with both developmental and environmental factors influencing the establishment of epigenetic marks that, in turn, govern both the early stages of embryonic development and gametogenic processes (19). Several studies have shown that impaired epigenetic dysregulation can disrupt human spermatogenesis, contributing to male infertility and associated spermatogenesis disorders (20). Over 100 different chemical RNA modifications have been detected to date, including the N¹-methyladenosine (m^lA), m⁶A, and 5-methylcytosine (m⁵C) modifications (21). In eukaryotic cells, m⁶A is the most common form of RNA modification (22), controlling all stages of the RNA metabolism process including the folding, maturation, stabilization, and translation of modified mRNAs (23). Notably, m⁶A RNA modification has been shown to be crucial for male germ line development, particularly in the context of mammalian spermatogenesis (24, 25).

Analyses of testes tissue samples from a range of species have highlighted relationships between the process of spermatogenesis and RNA m⁶A modification (26–28). Male mice in which the ALKBH5 gene has been knocked out exhibit increased m⁶A levels and impaired fertility as compared to wild-type (WT) littermates, suggesting a role for this m⁶A demethylase in this physiological setting (29, 30). Notably, METTL3-deficient murine embryonic stem cells do not undergo normal differentiation, and the mutation of METTL3 is associated with embryonic lethality in mice (29). Two FTO mutations have been linked to reductions in semen quality, and the dysfunction of the FTO protein has been tied to decreased male fertility (31). Single-cell sequencing data from human testis samples have revealed the expression of RNA m⁶A regulatory genes in almost all testis cell types, including both somatic and spermatogenic cells (32). Consistently, knocking down FTO, ALKBH5, METTL3, METTL14, YTHDF2, or YTHDC2 impairs normal gametogenesis and fertility (33). These prior results thus strongly suggest that altered RNA m⁶A modification may play a causative role in the regulatory processes underlying the molecular pathogenesis of male infertility.

The present study was developed with the goal of studying the mechanisms underlying m⁶A modification and associated regulatory processes in samples of testis tissue from sterile male cattle-yaks. To that end, a MeRIP-seq approach was used to establish whole transcriptomic m⁶A profiles for samples of testis tissue from normal sexually matured sterile male cattle-yaks. Using this approach, differentially methylated peaks were identified by comparing 5-years-old cattle-yaks (T group) with pre-sexually mature 18-months-old yaks (Y group) and post-sexually mature 5-years-old cattle-yaks (M group) that have been published previously (34), thus offering insight into the regulatory importance of m⁶A methylation in the context of cattle-yak sterility.

Materials and methods

Ethics statement

All animal-related procedures were consistent with guidelines established by the China Council on Animal Care and the Ministry of Agriculture of the People's Republic of China. The Animal Care and Use Committee of the Lanzhou Institute of Husbandry and Pharmaceutical Sciences Chinese Academy of Agricultural Sciences approved all yak handling procedures for this study (Permit No: SYXK-2014-0002).

Tissue sample isolation

Samples of testicular tissue were harvested following castration from three normal 5-years-old cattle-yaks in Xiahe County, Gannan Tibetan Autonomous Prefecture (N34°51 ', E102°26 '). Prior to tissue sample collection, iodophor was used to disinfect the samples. After sample isolation, a surgical needle was used to suture the wound site, with penicillin/streptomycin then being administered to protect against infection. The white testicular membrane was removed, and tissues were rinsed using 1 × PBS. Samples were then minced into 5 cm³ segments and transferred to a cryotube, after which they were snap-frozen using liquid nitrogen. A subsample of these testicular tissue isolates was also immobilized in Bouin's Fluid (SolarBio, Beijing, China). Samples were transferred to Lanzhou Institute of Husbandry and Pharmaceutical Sciences, Chinese Academy of Agricultural Sciences for subsequent use immediately after collection.

Hematoxylin and eosin staining

After fixation, testicular tissue samples were dehydrated using 75% ethanol, paraffin-embedded, and cut into 6 mm sections that were then stained with an improved H&E staining kit (SolarBio, Beijing, China) based on provided direction. Sections were then sealed using neutral gum and imaged with a Pannoramic 250 digital section scanner (Drnjier, Jinan, China).

RNA isolation and CDNA preparation

TRIzol (Invitrogen, CA, USA) was used to extract RNA from tissue samples, after which the purity (OD260/280 ratio) and concentration of RNA in these samples was measured with a NanoDrop 2,000 instrument (ThermoFisher Scientific, MA, United States). RNA at a concentration of 500–5,000 ng/mL with an OD260/280 ratio of 1.9–2.1 was selected and diluted to 500 ng/mL. Then, cDNA was prepared using a Transcriptor First Strand cDNA Synthesis Kit (Takara Bio Inc., Dalian, China), after which it was stored at−80°C for subsequent use.

qPCR

Analyses of cattle-yak testis m⁶A status were performed using appropriate primers (Supplementary Table S1), a LightCycler^® 96 Instrument (Roche, Beijing, China), and reaction conditions (Supplementary material) reported previously (34). Analyzed RNA methylation-associated genes detected via qPCR included METTL3, METTL14, WTAP, FTO, ALKBH5, YTHDF1/2/3, YTHDC1/2, RBM15, VIRMA, and ZC3H13. Relative gene expression levels were assessed via the 2⁻ΔΔCT method (35), with GAPDH serving as a normalization control. Analyses were repeated in triplicate, and differences in gene expression were compared via analyses of variance (ANOVAs).

m⁶A content analyses

An EpiQuik RNA Methylation Quantification Kit (Epigentek, P-9005, NY, United States) was used to detect mRNA m⁶A levels in cattle-yak testis samples based on provided directions.

MeRIP-seq and RNA-seq analyses

Sequencing libraries (Supplementary material) were prepared with the Illumina TrueSeq Stranded mRNA platform as in prior reports (34). The Illumina HiSeq X10 System at OE Biotech Co., Ltd. (Shanghai, China) was used to perform paired-end sequencing. All data were submitted to the Gene Expression Omnibus (GEO) database (Accession number: GSE205649). RNA-seq data were subjected to quality control and statistical analyses (Supplementary material) as in prior reports (34).

Statistical analysis

SPSS v 21.0 was used for all statistical analyses (36). Data were compared via one-way ANOVAs, with P < 0.05 or P < 0.01 as the threshold of significance (37).

Results

Histopathological analyses of cattle-yak testicular tissue samples

Initially, histopathological examination of cattle-yak testicular tissue samples was performed under a light microscope (10x, Figure 1A). Visually, the seminiferous tubules were loosely arranged with air bubbles between them. When examined under further magnification (40×, Figures 1B,C) germ cells in these testis tissues were found to primarily be arranged in a spermatogonia monolayer attached to the seminiferous tubule basement membrane with only limited numbers of visible spermatocytes.

FIGURE 1

Figure 1. Cattle-yak testicular tissue sections. (A) 10× magnification; (B,C) 40× magnification.

qPCR and global m⁶A quantification analyses

Next, overall methylation levels were compared between testicular tissue samples collected from 5-years-old cattle-yaks, 18-months-old yaks (pre-sexually mature yaks), and 5-years-old yaks (post-sexually mature yaks). While the levels of methylation in 5-year-old yaks were higher than in 5-years-old cattle-yaks, these differences were not significant. However, methylation levels were significantly higher in testicular samples from both 5-years-old yaks and cattle-yaks in comparison to 18-months-old yaks (Figure 2A). The mRNA level expression of several methylation-related genes was also assessed in these samples (Figure 2B), revealing that the majority of these genes were downregulated in cattle-yak testicular samples relative to those in samples from both 5-years-old and 18-months-old yaks. Specifically, testicular METTL14, YTHDF1, YTHDF2, RBM15, ZC3H13, and VIRMA expression was significantly reduced in cattle-yaks relative to both analyzed yak samples, whereas testicular WTAP, FTO, ALKBH5, YTHDC1, and YTHDC2 mRNA levels were significantly lower in cattle-yaks relative to 5-years-old yaks, although they did not differ significantly relative to 18-months-old yaks. Moreover, METTL3 and YTHDF3 expression was significantly increased in cattle-yaks relative to 18-months-old yaks, whereas they did not differ significantly from levels in 5-years-old yaks. The observed alterations in methylation-associated enzyme expression may underlie profound differences in the m⁶A methylation status of cattle-yak testicular tissues.

FIGURE 2

Figure 2. Analyses of m⁶A levels in yak and cattle-yak testicular tissues. (A) Overall levels of testis tissue methylation. Y: 18-months-old yak, T: cattle-yak, M: 5-years-old yak, **P < 0.01, *P < 0.05. (B) Methylation-related enzyme expression in the indicated yak groups. Different lowercase letters correspond to significant differences among groups (P < 0.05), whereas identical letters indicate a lack of any significant differences (P > 0.05).

Sequencing data quality control and reference genome comparisons

Next, the MeRIP-seq analysis of 6 samples was conducted (Supplementary Table S2), yielding 43.31 Gb of clean data. Each sample yielded 6.31–8.15 Gb of data, with Q30 base distributions of 94.83%-96.30%, with an average GC content of 56.70%. These reads were successfully aligned to a reference genome (LU_Bosgru_v3.0: http://ftp://ftp.ensembl.org/pub/release-99/fasta/bos_grunniens/dna/Bos_grunniens.LU_Bosgru_v3.0.dna_sm.toplevel.Fa.gz, with alignment rates of 76.01%-89.07% (Supplementary Table S3).

Methylation peak detection and annotation

Next, the whole transcriptome m⁶A profile for cattle-yak testicular tissues were obtained by high-throughput sequencing. A total of 16,186 peaks (Supplementary Table S4) were detected, and the widths of most peaks was distributed between 1 bp and 1,000 bp (n = 8,059), for details of the numbers of methylation peaks in each width range (see Table 1). For further details regarding methylation peak width distributions, see Supplementary Figure S1A. As the predicted methylation sites varied among samples, they were next separated into five reliability-based subcategories: Non-reliable, Low, Moderate, High, and Very High (Supplementary Figure S1B). CHIPseeker-based annotation results revealed 1–18 m⁶A peaks per gene, with 54% of genes exhibiting only a single peak (Supplementary Table S5), while the SYNE2 gene encoded on chromosome 11 exhibited the highest number of peaks (18 m⁶A peaks).

TABLE 1

Table 1. The numbers of methylation peaks in each width range.

Most peaks were concentrated in the exonic regions of associated genes, with progressively fewer peaks in the 3'-UTR and 5'-UTR regions (Figure 3A). To demonstrate representative patterns of m⁶A methylation, two model genes were selected for depiction. For the VGLL3 mRNA, peaks were present in the CDS and 3'-UTR regions, while in the DONSON mRNA they were located in the 3'-UTR, CDS, and 5'-UTR regions (Figure 3B). As m⁶A modification is generally associated with the 5-RRACH-3 sequence (R = A or G; H = A, C or U) (38), the enrichment of m⁶A peaks in the common RRACH sequence was assessed, revealing comparable enrichment in these cattle-yak samples to previously reported data (Figure 3C). These results thus offer further credibility to these m⁶A peak data, supporting the regulatory role of a ubiquitous mechanism governing patterns of RNA methylation.

FIGURE 3

Figure 3. Cattle-yak methylation peak analyses. (A) mRNA methylation peak distributions. (B) IGV plots demonstrating peaks in the VGLL3 and DONSON genes. (C) The motifs most commonly subjected to m⁶A modification in cattle-yak testicular samples.

Detection and annotation of differentially methylated peaks

Differentially methylated peaks (DMPs) were next analyzed with the MeTDiff software (screening criteria: diff. p ≤ 0.05; diff.fc ≥ 1.5), with 5,948 DMPs being identified for the T vs. Y comparison, of which 3,953 and 2,843 were upregulated and downregulated, respectively, in cattle-yak samples (Figure 4A). Similarly, 5,946 DMPs were identified for the T vs. M comparison, of which 3,103 and 2,843 were respectively upregulated and downregulated in cattle-yak samples (Figure 4B). DMP width distributions are presented in Supplementary Figure S1C. Distributions of detected DMPs on gene functional elements were further annotated, revealing the majority of these peaks to be concentrated in exonic regions, followed by the 3'-UTR and 5'-UTR regions (Figure 4C).

FIGURE 4

Figure 4. Identification and annotation of differentially methylated peaks. (A,B) Volcano plots demonstrating differentially methylated peak distributions. (C) Pie charts representing the annotation of DMPs associated with different gene regions.

Enrichment analyses of differentially methylated peak-associated genes

To gain potential insight into the effects of m⁶A methylation on gene expression and function in the testes of sterile male cattle-yaks, GO, and KEGG functional enrichment analyses of genes bearing DMPs were next conducted. GO analyses of the DMPs identified for the T vs. Y comparison revealed these genes to be enriched for molecular functions including cytoskeletal protein binding, tubulin binding, and ubiquitin-like protein ligase binding, cellular components including the microtubule organizing center, cytoskeleton, and Golgi apparatus, and biological processes including the cell cycle, cytoskeleton organization, and cell differentiation (Figure 5A). KEGG analyses for the T vs. Y comparison further indicated that DMP-associated genes were significantly enriched in the regulation of actin cytoskeleton, Notch signaling, cell cycle, and TGF-β signaling pathways (Figure 5B). GO analyses for the T vs. M group comparison similarly revealed DMPs to be enriched for molecular function terms including microtubule binding, kinesin binding, and protein serine/threonine kinase activity, cellular components including microtubule-organizing centers, actin filaments, and the ubiquitin ligase complex, and biological processes including meiotic spindle organization, ubiquitin-dependent protein catabolic process, and endocytosis (Figure 5C). KEGG analyses for the T vs. M comparison additionally indicated that these DMPs were associated with the homologous recombination, apoptosis, steroid hormone biosynthesis, and Fanconi anemia pathways (Figure 5B).

FIGURE 5

Figure 5. Enrichment analyses. (A) GO analyses for the T vs. Y group. (B) KEGG analyses for the T vs. Y group. (C) GO analyses for the T vs. M group. (D) KEGG analyses for the T vs. M group.

RNA-seq-based differential gene enrichment analyses

Next, RNA-seq analyses were used to compare changes in mRNA expression in cattle-yak testis tissues (Supplementary Table S6). Relative to testis samples from 18-months-old yaks, cattle-yak samples exhibited 5,895 and 6,577 genes that were respectively up- and downregulated (P < 0.05, log2FC > 0.58). Relative to 5-years-old yak tissue samples, cattle-yak samples exhibited 5,385 and 6,146 genes that were up- and down-regulated, respectively (P < 0.05, log2FC > 0.58) (Figure 6C). GO analyses for differentially expressed genes (DEGs) identified for the T vs. Y group comparison revealed these genes to be enriched for molecular functions including actin filament binding, protein serine/threonine kinase activity, and ATP binding, cell components including the cytoskeleton, cytoplasm, and centrosome, and biological processes including spermatogenesis, actin filament organization, and intracellular signal transduction (Figure 6A). KEGG analyses for the T vs. Y comparison revealed these DEGs to be primarily enriched in the cellular senescence, regulation of actin cytoskeleton, and MAPK signaling pathways (Figure 6D). GO analyses for the T vs. M comparison revealed associated DEGs to be enriched for molecular function terms including ATP binding, actin filament binding, and protein serine/threonine kinase activity, cellular component terms including the cytoplasm, centrosome, and microtubules, and biological process terms including spermatid development, protein ubiquitination, and DNA replication-dependent nucleosome assembly (Figure 6B). Similarly, KEGG analyses for the T vs. M comparison revealed associated DEGs to primarily be enriched in the apoptosis, regulation of actin cytoskeleton, and Fanconi anemia pathways (Figure 6E).

FIGURE 6

Figure 6. Functional enrichment analyses for DEGs. (A) GO analyses for the T vs. Y group. (B) KEGG analyses for the T vs. Y group. (C) DEG statistics. (D) GO analyses for the T vs. M group. (E) KEGG analyses for the T vs. M group.

Combined analysis of genes differentially expressed and differentially methylated in yak and cattle-yak testis samples

To gain additional insight into the link between m⁶A methylation and mRNA expression, we further analyzed genes that were both differentially methylated (Log₂ FC > 1, P < 0.05) and differentially expressed (log₂ FC > 1, P < 0.05). In total, 1,554 genes for the T vs. Y comparison were significantly differentially methylated and expressed, of which 812 exhibited increased m⁶A peaks and increased mRNA expression, 371 exhibited increased m⁶A peaks and decreased mRNA expression, 113 exhibited reduced m⁶A peaks and increased mRNA expression, and 258 exhibited decreased m⁶A peaks and decreased mRNA expression (Supplementary Table S7). Moreover, 1,364 genes for the T vs. M comparison were differentially expressed and differentially methylated, of which 541 exhibited increased m⁶A peaks and increased mRNA expression, 301 exhibited increased m⁶A peaks and decreased mRNA expression, 226 exhibited decreased m⁶A peaks and increased mRNA expression, and 296 exhibited decreased m⁶A peaks and decreased mRNA expression (Supplementary Table S8). For further details regarding the link between m⁶A methylation status and mRNA expression levels (see Figure 7A).

FIGURE 7

Figure 7. Combined analyses of differentially expressed and differentially methylated genes. (A) A combined four-quadrant scatterplot highlighting genes exhibiting significantly altered methylation and expression patterns when comparing the indicated samples, with gray dots corresponding to genes that did not exhibit any significant differences. (B) qPCR-based validation of DEGs identified via RNA-seq analyses. **P < 0.01, *P<0.05.

Lastly, qPCR analyses were performed to confirm the expression of differentially methylated genes in testicular tissue samples from yaks and cattle-yaks. The detected expression patterns were consistent with RNA-seq results (Figure 7B), thus validating these transcriptomic analyses.

Discussion

Distant hybridization is a term used to refer to hybridization between relatively distantly related species or genera (39), ultimately promoting the exchange of genes between different biotypes and thereby contributing to altered genotypic and phenotypic characteristics in the resultant offspring. Genotypic changes associated with distant hybridization can arise at the DNA and chromosomal levels, while at the phenotypic level heterosis can arise as a result of the integration of the relative advantages of the parents of a given offspring (40). Distant hybridization can support the development of new species (41–43), but male sterility is common among these hybrid offspring as a consequence of reproductive isolation (44, 45).

Research to date on hybridization-related male sterility has focused on crosses between chickens and quails, silver-black foxes and blue foxes, horses and donkeys, Muscovy ducks and domestic ducks, and yaks and cattle. Wodsedalek et al. (46) posited that male mule infertility following the hybridization of donkeys and horses may be attributable to differences in the numbers of chromosomes in these parental species, contributing to meiotic block and impaired spermatogenesis. Sun et al. (47) and Feng et al. (48), in contrast, found that hybrid offspring of quails and chickens exhibited the same number of chromosomes as both parental species (2n = 78), but that the chromosome morphology in these hybrid offspring differed significantly from that of either parent. Imbalances in gonadal enzyme ratios and abnormal hormone offspring in the hybrid offspring of Muscovy ducks and domestic ducks reportedly contribute to abnormalities in testicular anatomy and germ cell division in these animals, with an incomplete blood-testis barrier also potentially contributing to male sterility in this case (49). The prophase arrest of germ cell meiosis in the hybrid offspring of silver-black foxes and blue foxes, together with lower testosterone concentrations and increased prolactin and LH concentrations, can adversely impact spermatogenesis (50).

In prior studies, cattle-yak seminiferous tubules were reported to only harbor Sertoli cells and a limited number of spermatogonia, whereas they were devoid of spermatocytes or other identifiable germ cells. Moreover, these hybrid testis samples exhibited a spermatogenic cell monolayer that was loosely attached to the basement membrane, with increased spermatogonia apoptosis and the disruption of meiotic progression at the mid-pachytene (51). Consistently, H&E staining of cattle-yak testis samples in the present study revealed that only Sertoli cells and spermatogonia were present within the seminiferous tubules, confirming the accuracy of these sampling strategies. The anterior pituitary of cattle-yaks was largely devoid of basophilic cells, with severe segmentation of the nucleus of follicle-stimulating cells and limited numbers of secretory granules, thereby contributing to decreased FSH secretion, limiting seminiferous tubule development (52). Shah et al. (6) employed a STA-PUT approach to isolate spermatocytes and spermatogonia from cattle-yak testis, revealing the diameters of these cells to be significantly reduced relative to those of corresponding cells from yaks or cattle. Cai et al. (53) conducted an RNA-seq analysis assessing genes differentially expressed between yak and cattle-yak testicular tissues, leading to the identification of a link between NLRP14 and STRA8 upregulation and undifferentiated spermatogonia accumulation and apoptosis in cattle-yak testes. In contrast, the downregulation of the SPP1, SPIN2B, and PIWIL1 genes was linked to impaired spermatogonia genomic integrity and cell cycle progression, with several other meiosis-assocaited genes also exhibiting some level of downregulation in this context. While many comprehensive analyses of these cattle-yak hybrids have thus been conducted to date, the specific mechanisms underlying male sterility in these animals has yet to be established.

Further qPCR-based quantitative analyses of methylation-associated enzymes in these animals revealed the majority to be expressed at lower levels in cattle-yak testes relative to testes from both 18-month-old and 5-year-old yaks. For example, METTL14, YTHDF1, and RBM15 expression levels were significantly reduced in cattle-yak testicular tissues relative to tissues from 18-month-old and 5-year-old yaks, while YTHDC2, FTO, and ALKBH5 levels were significantly reduced in cattle-yak samples relative to those from 5-year-old yaks. ALKBH5-mediated m⁶A modification can influence the stability and splicing of mRNAs with long 3'-UTRs in spermatocytes and round spermatids (54). As such, the m⁶A modification of ALKBH5 is critical to the meiotic and haploid stages of the spermatogenic process (55). In mice, the knockout of ALKBH5 can lead to the impairment of spermatogenesis and associated male sterility (56). YTHDC2-KO mice, in contrast, exhibit germ cells that fail to mature beyond the zygotic stage (57). The dual knockout of both METTL3 and METTL14 can suppress the translation of critical spermatogenesis-associated transcripts linked to m⁶A modification, contributing to further anaphase spermatogenesis abnormalities (18). Mutations in FTO are positively correlated with reductions in semen quality and may be linked to decreased male fertility (31). As such, enzymes associated with the process of m⁶A methylation are likely to be critical to the maintenance of normal male reproductive function, with decreases in methylation-related enzyme expression thus representing one potential cause of male cattle-yak sterility.

Methylation peaks in cattle-yak testis tissues were primarily concentrated in exonic regions, with lower levels in the 3'-UTR and 5'-UTR regions. Partial exon methylation peaks were largely concentrated at sites proximal to stop codons, in line with prior data from human and murine studies (58–60). Ke et al. (61) reported that the majority of m⁶A-modified residues are located near the final exon in a given transcript, providing a mechanism for 3'-UTR regulation. As such, overall m⁶A site distributions seem to be similar across different mammalian species. In prior studies, m⁶A peaks were shown to be enriched in regions harboring conserved RRACH motifs (58, 59), and the same was true for many of the m⁶A-associated sequences in the present study, in line with data from other studies of yeast, amphibians, plants, and mammals (26). Overall, these results provide support for the conservation of m⁶A methylation in mammals.

Functional enrichment analyses were conducted in an effort to explore the biological roles of DEGs and DMPs identified when comparing testis samples from cattle-yaks to those from 18-months-old or 5-year-old yaks. GO analyses for the T vs. Y group comparison revealed genes that were differentially expressed and differentially methylated to be associated with spermatogenesis-related GO terms including cytoskeleton and actin binding, and with the actin cytoskeleton regulation and MAPK signaling KEGG pathways. Moreover, GO analyses for the T vs. M group comparison revealed differentially expressed and differentially methylated genes to be associated with GO terms including protein ubiquitination, ubiquitin ligase complexes, ubiquitin-dependent protein catabolism, and endocytosis, and with the apoptosis and Fanconi anemia KEGG pathways. Fanconi anemia (FA) is a rare autosomal recessive genetic disorder in which patients experience progressive bone marrow failure associated with a range of congenital abnormalities, reduced fertility, and a higher risk of developing leukemia and head and neck squamous cell carcinoma (62, 63). Roughly 50% of women affected by FA are infertile, while fertility is rare among males diagnosed with FA (64, 65). In male FA patients, reduced fertility-related clinical symptoms include major reductions in spermatozoa with concomitant spermatozoa abnormalities, with Sertoli cell syndrome and non-obstructive azoospermia often being diagnosed in these individuals (66).

While methylation-associated enzymes are known to play a key role in the spermatogenic process and cattle-yak m⁶A transcriptional mapping has been conducted, how altered methylation levels ultimately contribute to cattle-yak male sterility has yet to be established. Here, analyses of cattle-yak testicular methylation profiles were performed at the whole-tissue level. However, the process of spermatogenesis is highly complex and relies upon synergistic interactions among multiple cell types. As such, further work will be necessary to establish the cell type-specific effects of RNA methylation in the context of male cattle-yak infertility. Future advances in single-cell sequencing technologies have the potential to offer new insight into the methylation of individual cells within cattle-yak testes, providing a foundation for future functional research exploring the link between m⁶A methylation and male cattle-yak sterility.

Conclusion

In summary, analyses m⁶A methylation levels in cattle-yak testis tissues and associated gene enrichment analyses were performed in this study. The expression of related methylases in cattle-yak samples was significantly reduced relative to that in yaks both before and after sexual maturity. Enrichment analyses of differentially expressed genes and genes associated with differentially methylated peaks indicated that the differentially expressed genes identified when comparing cattle-yaks and pre-sexually mature yaks were mainly associated with spermatogenesis, including cytoskeleton and actin binding, and actin cytoskeleton regulation. In contrast, differentially expressed genes identified when comparing cattle-yaks and sexually mature yaks were primarily associated with protein ubiquitination, ubiquitin ligase complexes, ubiquitin-dependent protein catabolism, and endocytosis, as well as with the apoptosis and Fanconi anemia KEGG pathways.

Data availability statement

The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found in the article/Supplementary material.

Ethics statement

All animal-related procedures were consistent with guidelines established by the China Council on Animal Care and the Ministry of Agriculture of the People's Republic of China. The Animal Care and Use Committee of the Lanzhou Institute of Husbandry and Pharmaceutical Sciences Chinese Academy of Agricultural Sciences approved all yak handling procedures for this study (Permit No: SYXK-2014-0002). Written informed consent was obtained from the owners for the participation of their animals in this study.

Author contributions

Conceptualization: XG and JP. Data curation and writing—original draft: XW. Investigation: SG, MC, and YK. Software: YL and CL. Writing—Review and editing: LX, PB, and PY. Funding acquisition: XG. All authors contributed to the interpretation of the results, writing of the article, read, and agreed to the published version of the manuscript. All authors contributed to the article and approved the submitted version.

Funding

This work was supported by the China Agriculture Research System of MOF and MARA (CARS-37) and the Innovation Project of Chinese Academy of Agricultural Sciences (25-LZIHPS-01).

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher's note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Supplementary material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fvets.2022.971515/full#supplementary-material

References

1. Wang K, Hu Q, Ma H, Wang L, Yang Y, Luo W, et al. Genome-wide variation within and between wild and domestic yak. Mol Ecol Resour. (2014) 14:794–801. doi: 10.1111/1755-0998.12226

PubMed Abstract | CrossRef Full Text | Google Scholar

2. Zhao J, Xu Z, Chen A, You X, Zhao Y, He W, et al. Identification of meat from yak and cattle using snp markers with integrated allele-specific polymerase chain reaction-capillary electrophoresis method. Meat Sci. (2019) 148:120–6. doi: 10.1016/j.meatsci.2018.08.019

PubMed Abstract | CrossRef Full Text | Google Scholar

3. Wang H, Chai Z, Hu D, Ji Q, Xin J, Zhang C, et al. A global analysis of Cnvs in diverse yak populations using whole-genome resequencing. BMC Genom. (2019) 20:61. doi: 10.1186/s12864-019-5451-5

PubMed Abstract | CrossRef Full Text | Google Scholar

4. Das P, Krishnan G, Doley J, Biswas T, Paul V, Chakravarty P, et al. Identification and expression profiling of Msy genes of yak for bull fertility. J Genet. (2019) 98:1–10. doi: 10.1007/s12041-019-1091-4

PubMed Abstract | CrossRef Full Text | Google Scholar

5. Das P, Krishnan G, Doley J, Bhattacharya D, Deb S, Chakravarty P, et al. Establishing gene Amelogenin as sex-specific marker in yak by genomic approach. J Genet. (2019) 98:7. doi: 10.1007/s12041-019-1061-x

PubMed Abstract | CrossRef Full Text | Google Scholar

6. Shah M, Xu C, Wu S, Zhao W, Luo H, Yi C, et al. Isolation and characterization of spermatogenic cells from cattle, yak and cattleyak. Anim Reprod Sci. (2018) 193:182–90. doi: 10.1016/j.anireprosci.2018.04.067

PubMed Abstract | CrossRef Full Text | Google Scholar

7. Ding LM, Wang YP, Brosh A, Chen JQ, Gibb MJ, Shang ZH, et al. Seasonal heat production and energy balance of grazing yaks on the Qinghai-Tibetan Plateau. Anim Feed Sci Tech. (2014) 198:83–93. doi: 10.1016/j.anifeedsci.2014.09.022

CrossRef Full Text | Google Scholar

8. Lan D, Xiong X, Wei Y, Xu T, Zhong J, Zhi X, et al. Rna-seq analysis of yak ovary: improving yak gene structure information and mining reproduction-related genes. Sci China Life Sci. (2014) 57:925–35. doi: 10.1007/s11427-014-4678-2

PubMed Abstract | CrossRef Full Text | Google Scholar

9. Liu P, Dong Q, Liu S, Degen A, Zhang J, Qiu Q, et al. Postpartum oestrous cycling resumption of yak cows following different calf weaning strategies under range conditions. Anim Sci J. (2018) 89:1492–503. doi: 10.1111/asj.13097

PubMed Abstract | CrossRef Full Text | Google Scholar

10. Li Y, Wang G, Xu S, Zhang X, Yang Q. The expression of histone methyltransferases and distribution of selected histone methylations in testes of yak and cattle-yak hybrid. Theriogenology. (2020) 144:164–73. doi: 10.1016/j.theriogenology.2020.01.001

PubMed Abstract | CrossRef Full Text | Google Scholar

11. Jiang Y, Dai P, Dai Q, Ma J, Wang Z, Hu R, et al. Effects of the higher concentrate ratio on the production performance, ruminal fermentation, and morphological structure in male cattle-yaks. Vet Med Sci. (2022) 8:771–80. doi: 10.1002/vms3.678

PubMed Abstract | CrossRef Full Text | Google Scholar

12. Niayale R, Cui Y, Adzitey F. Male hybrid sterility in the cattle-yak and other bovines: a review. BioL Reprod. (2021) 104:495–507. doi: 10.1093/biolre/ioaa207

PubMed Abstract | CrossRef Full Text | Google Scholar

13. Robert N, Yan C, Si-Jiu Y, Bo L, He H, Pengfei Z, et al. Expression of Rad51 and the histo-morphological evaluation of testis of the sterile male cattle-yak. Theriogenology. (2021) 172:239–54. doi: 10.1016/j.theriogenology.2021.06.018

PubMed Abstract | CrossRef Full Text | Google Scholar

14. Zhao S, Chen T, Luo X, Chen S, Wang J, Lai S, et al. Identification of novel Lncrna and differentially expressed genes (Degs) of testicular tissues among cattle, yak, and cattle-yak associated with male infertility. Animals (Basel). (2021) 11:2420. doi: 10.3390/ani11082420

PubMed Abstract | CrossRef Full Text | Google Scholar

15. Huang C, Ge F, Ma X, Dai R, Dingkao R, Zhaxi Z, et al. Longissimus dorsicomprehensive analysis of Mrna, Lncrna, Circrna, and Mirna expression profiles and their Cerna networks in the muscle of cattle-yak and yak. Front Genet. (2021) 12:772557. doi: 10.3389/fgene.2021.772557

PubMed Abstract | CrossRef Full Text | Google Scholar

16. Xu C, Shah M, Mipam T, Wu S, Yi C, Luo H, et al. Bovid micrornas involved in the process of spermatogonia differentiation into spermatocytes. Int J Biol Sci. (2020) 16:239–50. doi: 10.7150/ijbs.38232

PubMed Abstract | CrossRef Full Text | Google Scholar

17. Wang S, Pan Z, Zhang Q, Xie Z, Liu H, Li Q. Differential Mrna expression and promoter methylation status of Sycp3 gene in testes of yaks and cattle-yaks. Reprod Domest Anim. (2012) 47:455–62. doi: 10.1111/j.1439-0531.2011.01902.x

PubMed Abstract | CrossRef Full Text | Google Scholar

18. Lin Z, Hsu P, Xing X, Fang J, Lu Z, Zou Q, et al. Mettl3-/Mettl14-Mediated Mrna N-Methyladenosine modulates murine spermatogenesis. Cell Res. (2017) 27:1216–30. doi: 10.1038/cr.2017.117

PubMed Abstract | CrossRef Full Text | Google Scholar

19. McGowan M, Holland M, Boe-Hansen G. Review: ontology and endocrinology of the reproductive system of bulls from fetus to maturity. Animal. (2018) 12:s19–26. doi: 10.1017/S1751731118000460

PubMed Abstract | CrossRef Full Text | Google Scholar

20. Cai Z, Niu Y, Li H. Rna N6-methyladenosine modification, spermatogenesis, and human male infertility. Mol Hum Reprod. (2021) 27:gaab020. doi: 10.1093/molehr/gaab020

PubMed Abstract | CrossRef Full Text | Google Scholar

21. Liu Z, Li L, Sun H, Liu S. Link between M6a Modification and Cancers. Front Bioeng Biotech. (2018) 6:89. doi: 10.3389/fbioe.2018.00089

PubMed Abstract | CrossRef Full Text | Google Scholar

22. Sun T, Wu R, Ming L. The role of M6a Rna methylation in cancer. Biomed Pharmacother. (2019) 112:108613. doi: 10.1016/j.biopha.2019.108613

PubMed Abstract | CrossRef Full Text | Google Scholar

23. Yu F, Wei J, Cui X, Yu C, Ni W, Bungert J, et al. Post-translational modification of Rna M6a demethylase Alkbh5 regulates Ros-induced DNA damage response. Nucleic Acids Res. (2021) 49:5779–97. doi: 10.1093/nar/gkab415

PubMed Abstract | CrossRef Full Text | Google Scholar

24. Lin Z, Tong M. Ma Mrna modification regulates mammalian spermatogenesis. BBA-Gene Regul Mech. (2019) 1862:403–11. doi: 10.1016/j.bbagrm.2018.10.016

PubMed Abstract | CrossRef Full Text | Google Scholar

25. Xu K, Yang Y, Feng G, Sun B, Chen J, Li Y, et al. Mettl3-mediated ma regulates spermatogonial differentiation and meiosis initiation. Cell Res. (2017) 27:1100–14. doi: 10.1038/cr.2017.100

PubMed Abstract | CrossRef Full Text | Google Scholar

26. Sai L, Li Y, Zhang Y, Zhang J, Qu B, Guo Q, et al. Distinct Ma Methylome profiles in Poly(a) Rna from Xenopus Laevis Testis and that treated with atrazine. Chemosphere. (2020) 245:125631. doi: 10.1016/j.chemosphere.2019.125631

PubMed Abstract | CrossRef Full Text | Google Scholar

27. Xia H, Zhong C, Wu X, Chen J, Tao B, Xia X, et al. Mettl3 mutation disrupts gamete maturation and reduces fertility in zebrafish. Genetics. (2018) 208:729–43. doi: 10.1534/genetics.117.300574

PubMed Abstract | CrossRef Full Text | Google Scholar

28. Zhao T, Wang J, Shen L, Long C, Liu B, Wei Y, et al. Increased M6a Rna modification is related to the inhibition of the nrf2-mediated antioxidant response in Di-(2-Ethylhexyl) phthalate-induced prepubertal testicular injury. Environ Pollut. (2020) 259:113911. doi: 10.1016/j.envpol.2020.113911

PubMed Abstract | CrossRef Full Text | Google Scholar

29. Wojtas M, Pandey R, Mendel M, Homolka D, Sachidanandam R, Pillai R. Regulation of Ma transcripts by the 3' → 5' Rna helicase Ythdc2 is essential for a successful meiotic program in the mammalian germline. Mol Cell. (2017) 68:374–87.e12. doi: 10.1016/j.molcel.2017.09.021

PubMed Abstract | CrossRef Full Text | Google Scholar

30. Zheng G, Dahl J, Niu Y, Fu Y, Klungland A, Yang Y, et al. Sprouts of Rna epigenetics: the discovery of mammalian Rna demethylases. RNA Biol. (2013) 10:915–8. doi: 10.4161/rna.24711

PubMed Abstract | CrossRef Full Text | Google Scholar

31. Landfors M, Nakken S, Fusser M, Dahl J, Klungland A, Fedorcsak P. Sequencing of Fto and Alkbh5 in men undergoing infertility work-up identifies an infertility-associated variant and two missense mutations. Fertil Steril. (2016) 105:1170–9.e5. doi: 10.1016/j.fertnstert.2016.01.002

PubMed Abstract | CrossRef Full Text | Google Scholar

32. Wang M, Liu X, Chang G, Chen Y, An G, Yan L, et al. Single-cell Rna sequencing analysis reveals sequential cell fate transition during human spermatogenesis. Cell Stem Cell. (2018) 23:599–614.e4. doi: 10.1016/j.stem.2018.08.007

PubMed Abstract | CrossRef Full Text | Google Scholar

33. Wang L, Wu Z, Zou C, Liang S, Zou Y, Liu Y, et al. N6sex-dependent Rna editing and -adenosine rna methylation profiling in the gonads of a fish, the olive flounder. Front Cell Dev Biol. (2020) 8:751. doi: 10.3389/fcell.2020.00751

PubMed Abstract | CrossRef Full Text | Google Scholar

34. Wang X, Pei J, Guo S, Cao M, Bao P, Xiong L, et al. Characterization of N-Methyladenosine in domesticated yak testes before and after sexual maturity. Front Cell Dev Biol. (2021) 9:755670. doi: 10.3389/fcell.2021.755670

PubMed Abstract | CrossRef Full Text | Google Scholar

35. Jozefczuk J, Adjaye J. Quantitative real-time Pcr-based analysis of gene expression. Method Enzymol. (2011) 500:99–109. doi: 10.1016/B978-0-12-385118-5.00006-2

PubMed Abstract | CrossRef Full Text | Google Scholar

36. Qiu H, Jin G, Jin R, Zhao W. Analysis of Variance of Repeated Data Measured by Water Maze with Spss. Zhong xi yi jie he xue bao. (2007) 5:101–5. doi: 10.3736/jcim20070121

PubMed Abstract | CrossRef Full Text | Google Scholar

37. Lan X, Cretney E, Kropp J, Khateeb K, Berg M, Peñagaricano F, et al. Maternal diet during pregnancy induces gene expression and DNA methylation changes in fetal tissues in sheep. Front Genet. (2013) 4:49. doi: 10.3389/fgene.2013.00049

PubMed Abstract | CrossRef Full Text | Google Scholar

38. Yi D, Wang Q, Zhao Y, Song Y, You H, Wang J, et al. Nalteration of -Methyladenosine Mrna Methylation in a rat model of cerebral ischemia-reperfusion injury. Front Neurosci-Switz. (2021) 15:605654. doi: 10.3389/fnins.2021.605654

PubMed Abstract | CrossRef Full Text | Google Scholar

39. Zhang ZH, Chen J, Li L, Tao M, Zhang C, Tan QB, et al. Research advances in animal distant hybridization. Sci China Life Sci. (2014) 44:161–74.

PubMed Abstract | Google Scholar

40. Liu SJ. Distant hybridization leads to different ploidy fishes. Sci China Life Sci. (2010) 40:104–14.

PubMed Abstract | Google Scholar

41. Bierne N, Welch J, Loire E, Bonhomme F, David P. The coupling hypothesis: why genome scans may fail to map local adaptation genes. Mol Ecol. (2011) 20:2044–72. doi: 10.1111/j.1365-294X.2011.05080.x

PubMed Abstract | CrossRef Full Text | Google Scholar

42. Cheng J, Czypionka T, Nolte A. The genomics of incompatibility factors and sex determination in hybridizing species of cottus (Pisces). Heredity. (2013) 111:520–9. doi: 10.1038/hdy.2013.76

PubMed Abstract | CrossRef Full Text | Google Scholar

43. Mallet J. Hybrid speciation. Nature. (2007) 446:279–83. doi: 10.1038/nature05706

PubMed Abstract | CrossRef Full Text | Google Scholar

44. Schilthuizen M, Giesbers M, Beukeboom L. Haldane's rule in the 21st Century. Heredity. (2011) 107:95–102. doi: 10.1038/hdy.2010.170

PubMed Abstract | CrossRef Full Text | Google Scholar

45. Wünsch L, Pfennig K. Failed Sperm development as a reproductive isolating barrier between species. Evol Dev. (2013) 15:458–65. doi: 10.1111/ede.12054

PubMed Abstract | CrossRef Full Text | Google Scholar

46. Wodsedalek J E. Cause of the sterility in the mule. Biol Bull-US. (1966) 30:1–57. doi: 10.2307/1536434

CrossRef Full Text

47. Sun YQ, Zhao ZS, Meng QM, Li DQ. Karyotype analysis of chicken-quail hybrids. J Anhui Agri Sci. (2007) 14:4152–3.

48. Feng XL, Liu XX, Yao SX, Yao XJ, Zhao ZS. Analysis of possible mechanisms of early embryonic death in chicken-quail hybrids. China Poultry. (2009) 31:36–9.

Google Scholar

49. Rigdon R. Gonads in hermaphroditic ducks: a pathologic study. Am J Vet Res. (1967) 28:1125–34.

PubMed Abstract | Google Scholar

50. Yang T, Yang Y, Peng Y, Cong B, Diao Y, Bao K, et al. Comparative studies on testicular and epididymal morphology, and serum hormone concentrations in foxes and the hybrids during the breeding season. Anim Reprod Sci. (2016) 168:66–72. doi: 10.1016/j.anireprosci.2016.02.027

PubMed Abstract | CrossRef Full Text | Google Scholar

51. Yan P, Xiang L, Guo X, Bao P, Jin S, Wu X. The low expression of Dmrt7 is associated with spermatogenic arrest in cattle-yak. Mol Biol Rep. (2014) 41:7255–63. doi: 10.1007/s11033-014-3611-x

PubMed Abstract | CrossRef Full Text | Google Scholar

52. Min XY, Xiong XR, Xiong Y, Zhang H, He ZY, Li J. Research progress on sterility mechanisms of male cattle-yak. Chinese J Anim Sci. (2021) 57:46–52.

53. Cai X, Yu S, Mipam T, Yang F, Zhao W, Liu W, et al. Comparative analysis of testis transcriptomes associated with male infertility in Cattleyak. Theriogenol. (2017) 88:28–42. doi: 10.1016/j.theriogenology.2016.09.047

PubMed Abstract | CrossRef Full Text | Google Scholar

54. Tang C, Klukovich R, Peng H, Wang Z, Yu T, Zhang Y, et al. Alkbh5-dependent M6a demethylation controls splicing and stability of Long 3'-Utr Mrnas in male germ cells. Proc Natl Acad Sci USA. (2018) 115:E325–33. doi: 10.1073/pnas.1717794115

PubMed Abstract | CrossRef Full Text | Google Scholar

55. Chen H, Zhang J, Yan Y, Zhu C, Wang L, Fu S, et al. N6-methyladenosine Rna demethylase Alkbh5 is testis-specifically downregulated in hybrid male sterile Dzo and is a target gene of Bta-Mir-200a. Theriogenology. (2022) 187:51–7. doi: 10.1016/j.theriogenology.2022.04.022

PubMed Abstract | CrossRef Full Text | Google Scholar

56. Zheng G, Dahl J, Niu Y, Fedorcsak P, Huang C, Li C, et al. Alkbh5 Is a mammalian Rna demethylase that impacts Rna metabolism and mouse fertility. Mol Cell. (2013) 49:18–29. doi: 10.1016/j.molcel.2012.10.015

PubMed Abstract | CrossRef Full Text | Google Scholar

57. Hsu P, Zhu Y, Ma H, Guo Y, Shi X, Liu Y, et al. Ythdc2 is an N-Methyladenosine binding protein that regulates mammalian spermatogenesis. Cell Res. (2017) 27:1115–27. doi: 10.1038/cr.2017.99

PubMed Abstract | CrossRef Full Text | Google Scholar

58. Dominissini D, Moshitch-Moshkovitz S, Schwartz S, Salmon-Divon M, Ungar L, Osenberg S, et al. Topology of the human and mouse M6a Rna Methylomes revealed by M6a-Seq. Nature. (2012) 485:201–6. doi: 10.1038/nature11112

PubMed Abstract | CrossRef Full Text | Google Scholar

59. Meyer K, Saletore Y, Zumbo P, Elemento O, Mason C, Jaffrey S. Comprehensive analysis of Mrna methylation reveals Enrichment in 3' Utrs and near Stop Codons. Cell. (2012) 149:1635–46. doi: 10.1016/j.cell.2012.05.003

PubMed Abstract | CrossRef Full Text | Google Scholar

60. Schwartz S, Mumbach M, Jovanovic M, Wang T, Maciag K, Bushkin G, et al. Perturbation of M6a writers reveals two distinct classes of Mrna methylation at internal and 5' sites. Cell Rep. (2014) 8:284–96. doi: 10.1016/j.celrep.2014.05.048

PubMed Abstract | CrossRef Full Text | Google Scholar

61. Ke S, Alemu E, Mertens C, Gantman E, Fak J, Mele A, et al. A majority of M6a residues are in the last exons, allowing the potential for 3' Utr regulation. Gene Dev. (2015) 29:2037–53. doi: 10.1101/gad.269415.115

PubMed Abstract | CrossRef Full Text | Google Scholar

62. Ceccaldi R, Sarangi P, D'Andrea A. The fanconi anaemia pathway: new players and new functions. Nat Rev Mol Cell Bio. (2016) 17:337–49. doi: 10.1038/nrm.2016.48

PubMed Abstract | CrossRef Full Text | Google Scholar

63. Tsui V, Crismani W. The fanconi anemia pathway and fertility. Trends Genet. (2019) 35:199–214. doi: 10.1016/j.tig.2018.12.007

PubMed Abstract | CrossRef Full Text | Google Scholar

64. Alter B, Frissora C, Halpérin D, Freedman M, Chitkara U, Alvarez E, et al. Fanconi's anaemia and pregnancy. Brit J Haematol. (1991) 77:410–8. doi: 10.1111/j.1365-2141.1991.tb08593.x

PubMed Abstract | CrossRef Full Text | Google Scholar

65. Auerbach A. Fanconi anemia and its diagnosis. Mutat Res. (2009) 668:4–10. doi: 10.1016/j.mrfmmm.2009.01.013

PubMed Abstract | CrossRef Full Text | Google Scholar

66. Krausz C, Riera-Escamilla A, Chianese C, Moreno-Mendoza D, Ars E, Rajmil O, et al. From exome analysis in idiopathic azoospermia to the identification of a high-risk subgroup for occult fanconi anemia. Genet Med. (2019) 21:189–94. doi: 10.1038/s41436-018-0037-1

PubMed Abstract | CrossRef Full Text | Google Scholar