1. Introduction
Water buffaloes (Bubalus bubalis) are vital dairy animals with significant economic and cultural value in several regions across the world, such as Asia, South America, and the Mediterranean countries, including Egypt. The water buffalo, specifically river buffalo subspecies, was domesticated and introduced to Egypt and Italy during the eighth century by the Arabs [1]. According to FAOSTAT, 2023 “
Genome-wide association study (GWAS) has recently gained significant progress, facilitating rapid advancements in high-throughput genotyping and molecular technologies. The widespread adoption of GWAS has resulted in its integration into breeding programs for various livestock species, including swine [4], chickens [5], buffalo [6,7], and dairy cattle [8]. GWAS plays a vital role in the identification of markers and genomic regions associated with several economically important quantitative traits. A significant advancement has been made by introducing a novel approach known as single-step genomic best linear unbiased prediction (ssGBLUP) [9]. This approach has a distinct benefit and enhanced precision relative to the ordinary GWAS in simultaneously estimating breeding values for genotyped and non-genotyped animals and integrating genomic, pedigree, and phenotypic data into a unified model. Consequently, this has the potential to improve the accuracy of predicting the genetic merit of the animal [10]. The GEBVs obtained from ssGBLUP are converted into marker effects and marker weights, which are subsequently utilized in an iterative process to update the ssGBLUP solutions.
In practice, the application of GWAS has been employed to determine potential genes that are linked to semen traits in cattle [11]. Nevertheless, there is a paucity of research on genomic studies assessed for semen traits in buffalo. Therefore, the objectives of the present study were (1) to implement a ssGWAS methodology for semen traits in the Egyptian buffalo bulls using the Axiom® Buffalo Genotyping 90 K array and (2) to identify the genomic regions and the potential candidate genes associated with the studied semen traits. This study contributes to understanding the genetic/genomic architecture underlying the reproductive traits in buffalo bulls. Moreover, this study holds practical implications for buffalo breeders aiming to enhance the reproductive performance of their animals through selective breeding strategies.
2. Materials and Methods
2.1. Ethical Statement
This study is part of a project entitled “A genomic approach to improve production and reproduction traits in Egyptian buffalo”, funded by the Science, Technology & Innovation Funding Authority (STDF), Egyptian Ministry for Scientific Research, Egypt. The Institutional Animal Care and Use Committee of the Agricultural Research Center (ARC-IACUC), Egypt, reviewed and granted approval for the project’s protocols and procedures. The study was conducted under the approval number ARC-APRI-50-23.
2.2. Animals and Phenotypes
A total of 8465 semen ejaculates were collected throughout 14 consecutive years from 2009 to 2022 from 115 buffalo bulls. These bulls were reared in two different herds; 29 buffalo bulls (with 7626 ejaculates) originated in the International Livestock Management Training Center (IMTC) at Sakha, Kafr El-Sheikh Governorate, and 86 buffalo bulls (with 839 ejaculates) originated in Mahalet Mousa farm (MMF), Kafr El-Sheikh Governorate. The two farms belong to the Animal Production Research Institute (APRI), Ministry of Agriculture, Dokki, Egypt. Out of these 115 buffalo bulls with phenotypes, a total of 96 bulls have been genotyped. The original pedigree file involved 10,748 animals, with 179 sires and 2507 dams. The bulls of IMTC weighing 350–400 kg live body weight were fed a ration of 4 kg of concentrate feed, 3 kg of clover hay, and 4 kg of rice straw per day, and semen was collected twice a week at 8 a.m. from each bull. However, the bulls of MMF weighed 400–500 kg live body weight, fed concentrate feed contained 16% crude protein, rice straw, maize silage (during the summer season), and Egyptian Berseem, Trifolium alexandrinum (during the winter season). In that farm, every bull was fed twice a day, at 7 a.m. and 3 p.m., and semen was collected from each bull twice weekly (on Sunday and Wednesday). All the bulls used for semen collection were free of pathogenic diseases and healthy, with a likely body condition score of 3.5. These differences in the management followed in the two farms have been considered by adding the herd fixed effect in the model of analysis. In the two farms, the semen was collected through an artificial vagina (IMV Technologies, L’Aigle, France), maintained at 42–45 °C, and immediately transferred to the laboratory in a water bath (37 °C) for further semen evaluation. The semen traits examined were ejaculate volume (VOL), sperm mass motility % (MM), living sperm % (LS), abnormal sperm % (AS), and sperm concentration (CONC). The ejaculate volume was precisely measured in milliliters to the nearest 0.1 mL using a graduated glass tube. Semen mass motility (%) was calculated using the proportion of spermatozoa wave motion in a drop of semen placed on a glass slide. The live sperms (%) and abnormal sperms (%) were measured according to the methods of Hackett et al. (1965) [12] and Blom et al. (1983) [13]. The sperm concentration (109 spermatozoa/mL) in each ejaculate was measured using a Neubauer hemocytometer.
2.3. Genotyping and Quality Control
Blood samples were collected from the jugular vein and transferred into 15 mL Falcon tubes containing 1 mL of 0.5 M EDTA as an anticoagulant. The tubes were immediately chilled on ice to preserve the samples. Genomic DNA was subsequently extracted from the whole blood samples using the QIAamp® DNA Blood Mini Kit (QIAGEN, Hilden, Germany). The Axiom® Buffalo Genotyping 90 K array was utilized for genotyping, which was conducted by the Laboratorio Genetica e Servizi LGS-AGROTIS, Cremona, Italy (Agrotis S.r.L.-LGS, Cremona, Italy “
The initial genotypes comprised 123,040 SNP markers, evenly distributed across 24 autosomes and the X-chromosome (indicated as chromosome no. 25). The overall genotyping rate was 0.98. In order to ensure the quality of the data, a series of quality control steps were performed using PLINK 1.9 [14]. These procedures were performed involving the removal of SNPs with unknown and/or duplicate positions, low minor allele frequency (MAF < 0.01), or excessive missing genotypes per SNP (>0.10). Additionally, individuals with a low call rate (<0.10) were excluded from the analysis. Following the application of these filtering restrictions, the final dataset consisted of 192 genotyped animals (96 genotyped bulls and 96 genotyped cows) and 67,282 SNPs, with an improved genotyping rate of 0.99. The main objective of this study was to evaluate the semen quality in buffalo bulls. As a result, an additional 96 females were included in the genotyped dataset. This inclusion aimed to improve the accuracy of prediction during the construction of the H-matrix for single-step genomic evaluation.
2.4. Statistical Analyses
The mixed model procedure implemented in the R program was utilized to evaluate the significance of all available fixed effects [15]. Following the test, only the factors that demonstrated statistical significance were chosen for inclusion in the genetic analysis. Furthermore, prior to commencing the genetic analysis, a comprehensive examination of the data distribution was conducted to ensure the suitability for the subsequent steps.
The ss-GWAS was employed to evaluate each trait since this analysis procedure involves integrating phenotypes, genotypes, and pedigree information in one step using the BLUPF90 software family [16]. The animals were initially renumbered sequentially using the RENUMF90 software, and the resulting input files were then employed in the BLUPF90, PREGSF90, and POSTGSF90 programs for the purpose of conducting a single-step analysis. The univariate animal model was utilized to analyze the vector of y, which represents the phenotypic observations for each semen trait. The following univariate animal model was applied:
In this model, the vector comprises all the fixed effects, including the effects of the herd and the combined effect of the year and month of collection of semen samples. Additionally, linear regressions of the age of the bull at semen collection were incorporated into the model. The vectors of and are representative of random additive genetic and permanent environmental effects, respectively, whereas denotes the vector of random residuals. In order to examine the association between the observations denoted by the vector y and the different factors, namely, X, Z, and W incidence matrices were utilized. These three matrices are responsible for linking the observations of y to the fixed effects, random animal effects, and random permanent environmental effects, respectively. The analysis was conducted with assumptions of , , and , where H is a combined matrix from A (pedigree relationship matrix) and G (genomic relationship matrix), I is an identity matrix; and σ2α, σ2p, and σ2ε are additive, permanent environmental and residual variances, respectively. The inverse of the H matrix (H−1) was constructed according to Aguilar et al. (2010) [10] as: , where G−1 is the inverse of the genomic relationship matrix and is the inverse of numerator relationship matrix. The G matrix was constructed based on [17] as , where Z represents the difference between M matrix (matrix of genotypes with columns indicating the markers and rows representing the animals) and P matrix (frequency matrix of the second allele pj, expressed as 2pj).
The SNP effects were estimated using the ssGBLUP framework proposed by Wang et al. (2012) [9]. The mathematical representation of the function is given in matrix form as follows: , where û is the vector of SNP effect; D is the diagonal matrix containing weighting factors for the SNP effects; Z is the genotype matrix, and âg is the vector of breeding values estimated for genotyped animals. To calculate the variance explained by each SNP, the formula was used, where û corresponds to the SNP effect as described earlier, and p denotes the allele frequency of the SNP [18]. The percentage of genetic variance explained by a window segment of 30 adjacent SNPs was calculated as follows: , where ai i represents the genetic value of the i-th region that consists of 30 adjacent SNP and is the total genetic variance [19].
2.5. Candidate Genes Identification
In the current study, a sliding-window approach implemented in the ss-GWAS was used to assess the proportional contribution of genomic variance attributed to each genomic segment for the aforementioned semen traits. The top ten window segments (each spanning 30 consecutive SNPs) were selected, i.e., accounting for the highest proportion of the genomic variance to identify candidate genes for each semen trait. Furthermore, in order to comprehend the chromosomal impact, the percentages of variance explained by the windows belonging to each chromosome were aggregated. This approach provides valuable insights into the genomic architecture underpinning male fertility traits and contributes to pinpointing potential genomic regions of interest. The gene information was extracted from each window segment proximal for every trait and from the comprehensive database and library provided by Thermo Fisher Scientific “
2.6. Enrichment Analyses and Functional Annotation
The functional enrichment and metabolic pathways analysis were conducted using version 6.8 of the “Database for Annotation, Visualization, and Integrated Discovery” (DAVID) v 6.8 [20]. The parameters suggested by the authors were used in the computation for the functional analysis. In addition, the functions for each annotated gene were obtained independently by utilizing the database containing all annotated functions from Ensembl Rapid Release and DAVID.
3. Results
3.1. Descriptive Statistics of Semen Traits in the Egyptian Buffalo Bulls
Table 1 shows the main statistical features of the five studied semen traits in the Egyptian buffalo bulls. The overall means for VOL, MM, LS, AS, and CONC were 3.78 mL, 62.15%, 61.07%, 4.47%, and 0.76 × 109 spermatozoa/mL, respectively. The coefficient of variation (CV%) for semen traits exhibited a moderate range, spanning from 25.53% to 58.39%, as depicted in (Table 1).
3.2. Genotype Filtering and Quality Control
Following the quality control steps, the dataset comprised 67,282 SNPs, reflecting a substantially improved genotyping rate of 0.99. For the raw and filtered number of SNPs, the initial number of SNPs ranged from 1779 to 8533, while the filtered number of SNPs ranged from 1152 to 5526 (Table 2). The raw and filtered minor allele frequencies were moderate and ranged from 0.24 to 0.32 for initial MAF with an average of 0.28 and ranged from 0.28 to 0.34 for filtered MFA with an average of 0.32. The gap size between adjacent SNPs was high and ranged from 35.59 to 51.51 kb with an average of 38.9 ± 28 kb (Table 2). The raw and filtered whole chromosomal lengths were 2617 and 2614 Mb, respectively, in the whole buffalo genome (Table 2).
3.3. ss-GWAS for Semen Traits in the Egyptian Buffalo Bulls
As shown in Table 3, the X-chromosome accounted for 4.18% (Figure 1), 4.59% (Figure 2), 5.16% (Figure 3), 5.19% (Figure 4), and 4.31% (Figure 5) of the genomic variance related to VOL, MM, LS, AS, and CONC traits, respectively. Collectively, across the five semen traits, the X-chromosome explained 23.43% of the genomic variance, contrasting with just 8.4% attributed to the joint effect of the remaining eight chromosomes (3, 4, 6, 7, 14, 16, 20, and 22). Accordingly, the X-chromosome was the most critical chromosome compared with other autosomes. Within the set of regions under investigation, it has been observed that two specific genomic regions, spanning a chromosome length of 69–77 Mb and 81–92 Mb, display associations with scrotal circumference. Additionally, another distinct region located at 40–55 Mb exhibits polymorphisms that are correlated with the percentage of normal sperm.
3.4. Functional Annotation and Identification of Putative Candidate Genes for Semen Traits
The genomic regions that have been identified are located within or close to a collective total of 50 genes. The annotated genes associated with the semen traits in the Egyptian buffalo are described in Table 4. Regarding these genes, 33 genes are situated in the X-chromosome, while the remaining 17 are placed in the other chromosomes (3, 4, 6, 7, 14, 16, 20, and 22). However, all the genes investigated here are considered protein-coding genes, except one pseudo-coding protein described as dynein light chain 1, cytoplasmic-like gene (LOC102401551, Table 4).
The novel annotated gene (LOC102398394) with the Ensembl gene ID ENSBBUG00015018582 within the genomic region length of 48.19 to 121.68 Mb associated with VOL and it is known as melanoma-associated antigen D2 (Table 3 and Table 4). The novel gene (LOC112582129) annotated by the Ensembl gene ID of ENSBBUG00015019729 is located on the X-chromosome within the genomic region length of 48.19 to 121.68 Mb, known as cancer/testis antigen 47A-like gene and associated with VOL trait (Table 3 and Table 4). The novel gene (LOC102391626), the Ensembl gene ID of ENSBBUG00015015043, is located in the X-chromosome within the genomic region length of 29.12 to 127.08 Mb and associated with MM, LS, AS, and CONC traits (Table 3). This gene is known as actin-related protein T1 (Table 4). The novel gene (LOC102404920) nominated as a sodium/hydrogen exchanger 2-like gene was identified on a genomic window of X-chromosome, and it is associated with LS, AS, and CONC traits in the Egyptian buffalo (Table 3 and Table 4). A novel gene (LOC102397107) with an Ensembl gene ID of ENSBBUG00015009483 has been observed (Table 4). This gene is located oon chromosome 16 at a position of 33.86 Mb. It has been found to be associated with MM, LS, and AS traits (Table 3).
4. Discussion
The overall mean values for VOL, MM, LS, AS, and CONC exhibited a degree of similarity to those previously reported by Salem et al. (2023) [21] based on a subset of the current data from the IMTC herd. The authors reported that the measured values for VOL, MM, LS, AS, and CONC were 3.89 mL, 62.37%, 60.64%, 3.94%, and 0.67 × 109 spermatozoa/mL, respectively. Additionally, the present results fall within the range found in earlier studies in dairy and buffalo bulls [22,23]. In the Egyptian buffalo bulls, Mahmoud et al. (2013) [24] revealed that the overall means for VOL, MM, and LS were 2.9 mL, 70.9%, and 65.8%, respectively. Khattab et al. (2015) [22] reported that the actual means for VOL, MM, and LS were 3.26 mL, 58.89%, and 65.96%, respectively. Similarly, the mean values for VOL, LS, AS, and CONC were reported by Gabr et al. (2018) [25] as 2.6 mL, 67.3%, 18.4%, and 526.28 × 106 sperm/mL, respectively. Additionally, Rushdi et al. (2017) [26] reported values of 66.20%, 70.58%, and 15.15% for MM, LS, and AS, respectively. In the Indian Murrah buffalo bulls, Kumar et al. (2023) [27] observed that the actual means of ejaculate volume and sperm concentration were 2.82 mL and 1040.12 million/mL, respectively. Variations in the mean values for semen traits observed in this study and those reported by different researchers working on different breeds of dairy or buffalo bulls may be explained by differences in genetic makeup, reproductive and health status of bulls, age of bulls, frequency of collection, teamwork in collection, nutrition, season and year of collection, and management [21]. These CV figures agreed with those published by Khattab et al. (2015) [22], who demonstrated that the CV% range for Egyptian buffalo bulls’ semen traits was 21.86% to 38.61%.
A comprehensive process of genotype filtering and quality control was conducted to ensure the robustness of the dataset. Approximately 45% of the raw genotypes were excluded from further analysis due to various reasons. The aforementioned factors include unknown or duplicated positions, instances of missing genotypes exceeding the threshold of 10% per SNP, and MAF falling below the established threshold of 0.01. This meticulous filtration process guarantees the reliability and integrity of the subsequent analyses conducted in the current study. The current MAF values were in accordance with those of El-Halawany et al. (2017) [6], who conducted GWAS analyses on the Egyptian buffalo, indicated divergent genotyped SNPs, stating an average gap spacing of 62.66 ± 67.10 and an MAF of 0.30 ± 0.13. The present filtered SNPs are consistent with those GWAS for milk production traits in the Egyptian buffalo reported by Abdel-Shafy et al. (2020) [7], who revealed an average distance between SNPs of 40.8 ± 32.0 kb and an average MAF of 0.29 ± 0.13 across chromosomes. This variance underscores the potential influence of diverse genetic backgrounds or experimental conditions on SNP distribution and MAF estimation.
The optimization of buffalo breeding practices can be significantly improved by selectively identifying buffalo that demonstrate superior reproductive performance. Consequently, the identification of genetic markers for male fertility using GWAS can offer substantial benefits in evaluating semen traits, including VOL, MM, LS, AS, and CONC, which would be particularly relevant with the increased application of artificial insemination. In the current study, a sliding-window approach implemented in the ss-GWAS was used to assess the proportional contribution of genomic variance attributed to each genomic segment for the aforementioned semen traits. We selected the top ten window segments (each spanning 30 consecutive SNPs) that account for the highest proportion of the genomic variance for each trait. Furthermore, in the pursuit of comprehending the chromosomal impact, the percentages of variance explained by the windows belonging to each chromosome were aggregated. This approach provides valuable insights into the genomic architecture underpinning male fertility traits and aids in pinpointing potential genomic regions of interest.
Research on the area of dairy cattle breeding has emphasized the significance of the X-chromosome in relation to bull conception rates and semen-quality traits, including sperm concentration and motility [28,29]. Moreover, the X-chromosome relevance in exploring genes associated with spermatogenesis and male fertility has been proposed by Mueller et al. (2013) [30], highlighting the influence of the X-chromosome on genome-wide recombination and fertility. The investigation reported by Liu et al. (2019) [31] revealed that the enrichment of genes related to sexual development and reproduction is based on the sex chromosomes. Fortes et al. (2020) [32] evidenced a significant influence of the X-chromosome on bull fertility within two bovine populations. A study by Fortes et al. (2013) [33] has documented that the X-chromosome in bovines exhibits multiple regions linked to sperm morphology and testis development. To further this comprehension, a study conducted by Suchocki et al. (2015) [34] examined Polish Holstein-Friesian bulls and investigated semen characteristics using a GWAS. The researchers discovered a highly significant correlation between the X-chromosome and various semen traits, including sperm concentration, semen volume, motility score, and spermatozoa count. This finding indicates that the X-chromosome plays a crucial role in determining the quality of semen. In a similar context, Pacheco et al. (2020) [35] identified three candidate genes (FAM9B, TBL1X, and PIH1D3) within the X-specific genomic region that exhibited strong correlations with sire conception rate, testosterone concentration, spermatogenesis, and sperm motility, respectively, in dairy bulls.
The alignment in the genomic loci (29.35 to 127.09 Mb in the X-chromosome) highlights the similarity between the influences of the two traits (MM and LS). This consistency provides evidence for the substantial genetic correlation previously reported by Salem et al. (2023) [21] in the Egyptian buffalo bulls (0.99 ± 0.01), based on a subset of current data from the IMTC herd. A study by Khattab et al. (2015) [22] revealed a robust genetic correlation of 1.0 ± 0.03 between MM and LS traits, indicating that the enhancement of one trait could potentially lead to improvements in the other. The findings reported by Khattab et al. (2022) [23] for Friesian bulls gave a consistent genetic correlation of 0.64 ± 0.01 between MM and LS traits. According to the existing literature, no prior studies have investigated GWAS concerning semen traits in buffalo bulls. However, numerous genomic regions associated with semen characteristics have been identified in dairy and beef cattle through GWAS. In this regard, these genomic regions were previously identified in the Holstein Friesian bulls [36], in the US Jersey dairy bulls [29], in the Chinese Holstein bulls [28], in the Italian Holstein bulls [37], and in crossbred beef bulls [38].
The E3 ubiquitin-protein ligase SIAH1-like gene identified in the X-chromosome (48.19 to 121.68 Mb) plays a vital role in the spermatogenesis processes [39]. In the same genomic region, the novel annotated gene (known as melanoma-associated antigen D2) is considered one of the melanoma antigen proteins family (MAGE). Most MAGE genes are found in the X-chromosome and are expressed during the early stages of spermatogenesis [40]. Based on the tissue expression pattern, the MAGE gene can be divided into types I and II [41]. Type II MAGEs, on the other hand, have widespread expression in numerous tissues. According to several studies, MAGE genes are crucial for embryogenesis and the development of germ cells [41]. Moreover, in the X-chromosome, TRO, commonly known as Trophinin, was identified, and this gene could improve sperm motility [42].
GRIA3 gene was identified in the X-chromosome (48.19 to 121.68 Mb), and the mutation of this gene led to an increased risk of 20–30% in decreasing the libido [43]. In the same genomic region, a novel gene was identified (cancer/testis antigen 47A-like gene). This gene exhibits distinctive and limited expression patterns during different stages of spermatogenesis, including cell cycle progression, meiosis, and spermiogenesis. AR gene (located in the X-chromosome) is a ligand-dependent transcription factor that modulates testosterone signaling and is essential for spermatogenesis and fertility [44]. The effects of androgens, which are necessary for the growth and maintenance of the male reproductive system, are mediated through the AR gene. Destructive mutations in AR can cause disorders ranging from minor deviations to complete failure of normal male phenotypic development [45]. The key cell regulators of androgen signaling androgen receptors are DNA-binding transcription factors primarily activated by testosterone and 5-dihydrotestosterone.
The AR candidate gene had significant effects on sperm quality traits (abnormal spermatozoa rate) in boars [46]. Sarakul et al. (2018) [47] found that the COL4A6 gene was associated with semen volume, number of sperm, and sperm motility in the Thai multi-breed dairy population. ZC4H2 gene located in the genomic region of 29.12 to 127.08 Mb in the X-chromosome was identified as an RNF220 interacting protein, a partner protein of RNF220 that stabilizes RNF220 [48]. In some cases, it serves as a mediator for the interaction between RNF220 and target proteins [49]. The ZC3H12B gene (located in the 29.12 to 127.08 Mb genomic region of the X-chromosome) is a negative regulator of macrophage activation and may play a role in inflammatory disorders and host immunity in the same genomic region. A novel annotated gene was identified, and this gene is known as actin-related protein T1. This gene encodes a protein related to the cytoskeletal protein beta-actin, and it presumably participates in the production of spermatids since it makes up a significant portion of the calyx in the perinuclear theca of the mammalian sperm heads “
The novel annotated sodium/hydrogen exchanger 2-like gene was identified in the X-chromosome. The inhibition of acrosome production decreases its expression, ultimately causing male infertility. Interestingly, a distinct PFKFB activity was detected in the testes and a human testes-specific, suggesting that the unique metabolic requirements of sperm may necessitate a divergent PFKFB family member [52]. The GPR174 gene nominated as G protein-coupled receptor 174 was previously reported to play an essential role in inducing the acrosome reaction of bull sperm [53]. The DACH2 gene is a marker of immature Sertoli cells and may have a role in reproduction and development [54]. Regarding the X-chromosome, the BCAP31 gene could be a significant marker for spermatogenesis and apoptosis and plays a pivotal role in the control of testis function in humans [55] and sheep [56].
KIF2B gene is located on the genomic region length of 58.52 to 137.06 Mb of chromosome number three; known as a meiosis-related gene, it has been discovered to be hypomethylated in the sperm but hypermethylated in the somatic cells. The hypomethylated KIF2B has a positive correlation with the fertility of bulls [57]. The KIF2B gene affects spindle bipolarity by causing microtubule depolymerization to alter mitotic dynamics during spermatogenesis [58]. In the same genomic region of chromosome number three, the PTK2B was identified. This gene belongs to the family of focal adhesion kinases and can be activated in response to outside stimuli that work through receptor tyrosine kinases [59]. This gene is crucial for modifying actin at sperm-binding sites and encourages sperm integration into the oocyte [60]. Yang et al. (2019) [61] reported that the PTK2B gene plays a vital role in regulating bovine mastitis. In the genomic region length of 60.92 to 62.73 Mb of chromosome number three, the CA10 gene was identified, and it could have a crucial role in the interconversion of carbon dioxide, bicarbonate, and other physiological processes in many tissues [62]. CA10 gene might regulate buffaloes’ fat metabolism and reproductive development [63]. Moreover, Carvalho et al. (2019) [64] reported that the CA10 gene was jointly associated with other genes controlling hot carcass weight in Nelore cattle.
In chromosome number six, within the genomic region length of 51.58 to 52.47 Mb, the BARHL2 gene (BarH-like homeobox 2) was identified. It is a transcriptional regulator and a member of the BarH family of homeodomain proteins and has been shown to affect cell fate specification, differentiation, survival, and migration [65]. The KCTD8 gene (located in chromosome number seven within the genomic region length of 53.23 to 54.27 Mb) is associated with the milk fat percentage, as reported in GWAS in Indian buffaloes [63]. In the same genomic region of chromosome number seven, YIPF was identified. This gene could be suggested to play a significant role in the pathophysiology of varicocele. The potential correlation between the expression of YIPF and the testicular inflammatory response, as well as the fluctuation in TNF expression level induced by varicocele, has been suggested [66]. In the genomic region 34.04 to 35.39 Mb of chromosome 14, the BMP2 gene was identified, and it is physiologically expressed in the testis and acts as proliferative signals upon Sertoli cells and spermatogonia [67].
The novel annotated UBQLN3 gene (ubiquilin-3) was identified within the genomic region of 33.86 to 36.95 Mb. All ubiquitin-like proteins or ubiquilin-3 have a carboxy-terminal ubiquitin-associated (UBA) domain and an amino-terminal UBL domain. The testis is where this gene primarily expresses itself [68]. However, it has been proposed that this particular gene may have an impact on the progression of the cell cycle during the process of spermatogenesis and has been linked to male fertility [69,70].
5. Conclusions
The present study focuses on the implementation of single-step genome-wide association studies (ss-GWAS) to analyze semen traits in Egyptian buffalo bulls, representing the inaugural application of ss-GWAS in this particular population. The X-chromosome mainly explained 23.43% of the total genomic variance across the five semen traits. In contrast, the remaining eight chromosomes (3, 4, 6, 7, 14, 16, 20, and 22) accounted collectively for 8.40% of the total genomic variance. This finding suggests that semen traits exhibit polygenic inheritance patterns rather than simple ones. Several candidate genes associated with male fertility and spermatogenesis were identified in the X-chromosome. In practice, five novel candidate genes were detected across the whole genomic regions studied for semen traits in the Egyptian buffalo bulls (four novel genes in the X-chromosome and one novel gene in chromosome 16). A more precise buffalo genomic map is required to identify the candidate genes affecting semen traits, considering a larger sample size that may help to improve the power of candidate gene detection and to obtain more accuracy in the genomic breeding values.
A.G.E.N.: Designed the experiment, handled project administration, acquired funding, wrote the first draft, and revised the manuscript. M.H.K.: Designed the experiment and revised the manuscript. M.M.I.S.: Contributed to data analysis and manuscript writing. A.M.S.A.: Collected data from the farm, handled blood sample collection, facilitated other farm work, and organized DNA extraction. A.F.A.: Collected data from the farm, organized blood sample collection, and facilitated other farm work. M.M.H.: Assisted in collecting data from the farm and handling blood samples. H.A.-S.: Organized DNA extraction, analyzed the data, contributed to manuscript writing, and revised the manuscript. All authors have read and agreed to the published version of the manuscript.
The project protocols and procedures were checked and approved by the Institutional Animal Care and Use Committee of the Agricultural Research Center (ARC-IACUC), Egypt. The study was conducted under the approval number: ARC-APRI-50-23.
Not applicable.
The datasets of this study are confidential.
All the authors of the manuscript thank and acknowledge their respective Universities and Institutes. The authors thank the Science, Technology & Innovation Funding Authority (STDF); Egyptian Ministry for Scientific Research for supporting the present study.
The authors declare that they have no competing interest.
Footnotes
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.
Enlarge this image.Figure 1. ssGWAS results of ejaculate volume (VOL) in the Egyptian buffalo bulls. Single nucleotide polymorphisms on different chromosomes are denoted by different colors. Each dot represents one SNP. The X-axis represents 25 chromosomes. The Y-axis represents the percentage of the explained variance by SNP.
Enlarge this image.Figure 2. ssGWAS results of mass motility (MM) in the Egyptian buffalo bulls. Single nucleotide polymorphisms on different chromosomes are denoted by different colors. Each dot represents one SNP. The X-axis represents 25 chromosomes. The Y-axis represents the percentage of the explained variance by SNP.
Enlarge this image.Figure 3. ssGWAS results of sperm livability (LS) in the Egyptian buffalo bulls. Single nucleotide polymorphisms on different chromosomes are denoted by different colors. Each dot represents one SNP. The X-axis represents 25 chromosomes. The Y-axis represents the percentage of the explained variance by SNP.
Enlarge this image.Figure 4. ssGWAS results of percentage of abnormal sperms (AS) in the Egyptian buffalo bulls. Single nucleotide polymorphisms on different chromosomes are denoted by different colors. Each dot represents one SNP. The X-axis represents 25 chromosomes. The Y-axis represents the percentage of the explained variance by SNP.
Enlarge this image.Figure 5. ssGWAS results of sperm concentration (CONC) in the Egyptian buffalo bulls. Single nucleotide polymorphisms on different chromosomes are denoted by different colors. Each dot represents one SNP. The X-axis represents 25 chromosomes. The Y-axis represents the percentage of the explained variance by SNP.
Descriptive statistics for semen traits in the Egyptian buffalo bulls.
| Trait 1 | Mean | SD | Min | Max | CV (%) |
|---|---|---|---|---|---|
| Ejaculate volume, mL | 3.78 | 1.79 | 0.50 | 11.50 | 47.35 |
| Mass motility, % | 62.15 | 16.30 | 5.00 | 95.00 | 26.23 |
| Livability, % | 61.07 | 15.59 | 3.00 | 88.00 | 25.53 |
| Abnormality, % | 4.47 | 2.61 | 1.00 | 44.00 | 58.39 |
| Concentration 2 | 0.76 | 0.35 | 0.20 | 3.75 | 46.05 |
1: Total number of records is 8465 ejaculates; 2: Sperm concentration (109 sperm/mL).
Summary statistics of genotyped SNPs in the Egyptian buffalo using Axiom Buffalo Genotyping Array.
| Chr. | Number of SNPs | MAF | Chr. Length [Mb] | Gap Size [kb] # | |||||
|---|---|---|---|---|---|---|---|---|---|
| Initial | Filtered * | Initial | Filtered * | Initial | Filtered * | Mean | SD | Max | |
| 1 | 8533 | 5526 | 0.25 | 0.29 | 201.95 | 201.95 | 36.55 | 22.36 | 716.57 |
| 2 | 7786 | 4993 | 0.24 | 0.29 | 188.83 | 188.83 | 37.83 | 27.05 | 699.06 |
| 3 | 7192 | 4564 | 0.27 | 0.31 | 175.41 | 175.39 | 38.44 | 27.10 | 762.66 |
| 4 | 6499 | 4192 | 0.27 | 0.29 | 165.14 | 165.14 | 39.40 | 29.32 | 719.43 |
| 5 | 5027 | 3229 | 0.29 | 0.30 | 127.53 | 127.50 | 39.50 | 40.53 | 1721.08 |
| 6 | 5064 | 3264 | 0.29 | 0.29 | 120.32 | 120.24 | 36.85 | 23.83 | 477.34 |
| 7 | 4722 | 3027 | 0.29 | 0.30 | 117.10 | 117.10 | 38.70 | 23.74 | 418.96 |
| 8 | 5031 | 3244 | 0.28 | 0.29 | 119.61 | 119.50 | 36.85 | 22.03 | 512.54 |
| 9 | 4341 | 2825 | 0.27 | 0.30 | 110.10 | 109.97 | 38.94 | 28.94 | 689.09 |
| 10 | 3929 | 2539 | 0.28 | 0.29 | 104.07 | 104.07 | 41.01 | 24.94 | 414.49 |
| 11 | 4012 | 2572 | 0.27 | 0.29 | 102.20 | 102.05 | 39.69 | 27.67 | 609.86 |
| 12 | 4461 | 2880 | 0.27 | 0.28 | 106.40 | 106.37 | 36.95 | 19.29 | 180.88 |
| 13 | 3385 | 2154 | 0.28 | 0.29 | 90.40 | 90.22 | 41.90 | 45.88 | 1229.15 |
| 14 | 3342 | 2146 | 0.27 | 0.29 | 83.09 | 82.88 | 38.64 | 23.93 | 501.28 |
| 15 | 3428 | 2192 | 0.28 | 0.28 | 82.03 | 82.03 | 37.44 | 21.27 | 197.54 |
| 16 | 3185 | 2040 | 0.29 | 0.30 | 84.41 | 84.36 | 41.37 | 43.17 | 1139.89 |
| 17 | 2937 | 1873 | 0.28 | 0.29 | 73.12 | 72.68 | 38.82 | 24.00 | 398.63 |
| 18 | 2737 | 1728 | 0.29 | 0.30 | 65.76 | 65.69 | 38.04 | 32.20 | 560.14 |
| 19 | 2968 | 1940 | 0.28 | 0.30 | 71.58 | 71.54 | 36.90 | 18.94 | 235.90 |
| 20 | 2637 | 1691 | 0.30 | 0.31 | 68.47 | 68.45 | 40.50 | 28.75 | 386.39 |
| 21 | 2557 | 1703 | 0.28 | 0.29 | 60.72 | 60.58 | 35.59 | 18.07 | 231.64 |
| 22 | 2529 | 1633 | 0.28 | 0.30 | 61.97 | 61.69 | 37.80 | 20.60 | 173.44 |
| 23 | 2172 | 1410 | 0.32 | 0.34 | 51.44 | 51.18 | 36.32 | 21.94 | 388.82 |
| 24 | 1779 | 1152 | 0.30 | 0.34 | 42.13 | 42.00 | 36.49 | 22.36 | 264.13 |
| X-chromosome | 4891 | 2765 | 0.31 | 0.31 | 143.40 | 142.37 | 51.51 | 62.41 | 1427.38 |
| UnPos | 17,896 | --- | 0.25 | --- | --- | --- | --- | --- | --- |
| Overall | 123,040 | 67,282 | 0.28 | 0.30 | 2617 | 2614 | 38.88 | 28.01 | 1721.08 |
Chr: chromosome; Mb: mega base; kb: kilo base; UnPos: un-positioned SNPs; MAF: minor allele frequency; * Number of SNPs after applying filtering restrictions (see materials and methods for further details) and duplicates SNP position; # Distance between adjacent SNPs after cleaning the data.
The percentage of genomic variance explained by relevant genomic region 1 associated with semen traits in the Egyptian buffalo bulls.
| Chromsome No. | Genomic Map Position (Mb) | Trait (Genomic Explained Variance %) | Gene ID |
|---|---|---|---|
| X-chromosome | 48.19–121.68 | VOL (4.18%) | LOC102414874, LOC102401551, LOC102398394, TRO, CHRDL1, GRIA3, LOC112582129, AR, OPHN1, COL4A6 |
| X-chromosome | 29.35–127.08 | MM (4.59%) | ZC4H2, ZC3H12B, LOC102391626, SMARCA1, DMD, GNL3L, ITIH6, P2RY10, LOC102415058, RTL9, AMMECR1 |
| X-chromosome | 29.12–127.08 | LS (5.16%) | ZC4H2, ZC3H12B, LOC102391626, SMARCA1, DMD, GNL3L, ITIH6, P2RY10, LOC102415058, RTL9, AMMECR1, POU3F4, LOC102404920 |
| X-chromosome | 29.12–127.08 | AS (5.19%) | PFKFB1, TBX22, GPR174, DACH2, LOC112582031, LOC102391626, SMARCA1, POU3F4, LOC102404920, RTL9, AMMECR1, DMD |
| X-chromosome | 47.85–142.21 | CONC (4.31%) | BCAP31, P2RY10, LOC102415058, LOC102414874, LOC102401551, LOC102391626, SMARCA1, POU3F4, LOC102404920, LOC102394334, ZIC3, MAGEH1, FOXR2 |
| 3 | 58.52–137.06 | VOL (0.99%) | KIF2B, LOC102390200, PTK2B |
| 3 | 60.92–62.73 | CONC (1.14%) | CA10 |
| 4 | 43.42–44.85 | VOL (0.42%) | CSF2RB |
| 6 | 51.58–52.47 | CONC (0.76%) | ZNF326, BARHL2 |
| 7 | 53.23–54.25 | MM (0.45%) | KCTD8, YIPF7 |
| 7 | 53.26–54.27 | AS (0.57%) | KCTD8, YIPF7 |
| 14 | 34.03–35.39 | VOL (0.40%) | FERMT1, BMP2 |
| 16 | 33.86–36.95 | MM (0.51%) | LOC102402176, LOC102397107 |
| 16 | 33.86–36.95 | LS (0.47%) | LOC102402176, LOC102397107 |
| 16 | 33.86–36.95 | AS (0.49%) | LOC102402176, LOC102397107 |
| 20 | 46.22–47.48 | MM (0.51%) | BLM, CRTC3 |
| 20 | 46.22–47.48 | LS (0.57%) | BLM, CRTC3 |
| 20 | 46.22–47.48 | CONC (0.61%) | BLM, CRTC3 |
| 22 | 49.91–50.83 | AS (0.51%) | LOC102389067, CDH7 |
1: The top 10 SNP-windows explained the highest genomic variance; VOL: ejaculate volume; MM: mass motility; LS: sperm livability; AS: percentage of abnormal sperms; CONC: sperm concentration.
The function of candidate genes located in relevant genomic region 1 associated with semen traits in the Egyptian buffalo bulls.
| Candidate Gene 2 | Chromosome No. | Gene Function |
|---|---|---|
| LOC102414874 | X-chromosome | E3 ubiquitin-protein ligase SIAH1-like |
| LOC102401551 | X-chromosome | Dynein light chain 1, cytoplasmic-like |
| LOC102398394 | X-chromosome | Melanoma-associated antigen D2 |
| TRO | X-chromosome | Trophinin |
| CHRDL1 | X-chromosome | Chordin like 1 |
| GRIA3 | X-chromosome | Glutamate ionotropic receptor AMPA type subunit 3 |
| LOC112582129 | X-chromosome | Cancer/testis antigen 47A-like |
| AR | X-chromosome | Androgen receptor |
| OPHN1 | X-chromosome | Oligophrenin 1 |
| COL4A6 | X-chromosome | Collagen type IV alpha 6 chain |
| ZC4H2 | X-chromosome | Zinc finger C4H2-type containing |
| ZC3H12B | X-chromosome | Zinc finger CCCH-type containing 12B |
| LOC102391626 | X-chromosome | Actin-related protein T1 |
| SMARCA1 | X-chromosome | SWI/SNF related, matrix associated, actin dependent regulator of chromatin, subfamily a, member 1 |
| DMD | X-chromosome | Dystrophin |
| GNL3L | X-chromosome | G protein nucleolar 3, like |
| ITIH6 | X-chromosome | Inter-alpha-trypsin inhibitor heavy chain family member 6 |
| P2RY10 | X-chromosome | Putative P2Y purinoceptor 10 |
| LOC102415058 | X-chromosome | P2Y receptor family member 10 |
| RTL9 | X-chromosome | Retrotransposon Gag like 9 |
| AMMECR1 | X-chromosome | AMMECR nuclear protein 1 |
| POU3F4 | X-chromosome | POU class 3 homeobox 4 |
| LOC102404920 | X-chromosome | Sodium/hydrogen exchanger 2-like |
| PFKFB1 | X-chromosome | 6-phosphofructo-2-kinase/fructose-2,6-biphosphatase 1 |
| TBX22 | X-chromosome | T-box transcription factor 22 |
| GPR174 | X-chromosome | G protein-coupled receptor 174 |
| DACH2 | X-chromosome | Dachshund family transcription factor 2 |
| LOC112582031 | X-chromosome | Protein FAM209B-like |
| BCAP31 | X-chromosome | B cell receptor associated protein 31 |
| LOC102394334 | X-chromosome | Putative MAGE domain-containing protein MAGEA13P |
| ZIC3 | X-chromosome | Zic family member 3 |
| MAGEH1 | X-chromosome | MAGE family member H1 |
| FOXR2 | X-chromosome | Forkhead box R2 |
| KIF2B | 3 | Kinesin family member 2B |
| LOC102390200 | 3 | Peptidyl-prolyl cis-trans isomerase A-like |
| PTK2B | 3 | Protein tyrosine kinase 2 beta |
| CA10 | 3 | Carbonic anhydrase 10 |
| CSF2RB | 4 | Colony stimulating factor 2 receptor subunit beta |
| ZNF326 | 6 | Zinc finger protein 326 |
| BARHL2 | 6 | BarH like homeobox 2 |
| KCTD8 | 7 | Potassium channel tetramerization domain containing 8 |
| YIPF7 | 7 | Yip1 domain family member 7 |
| FERMT1 | 14 | FERM domain containing kindlin 1 |
| BMP2 | 14 | Bone morphogenetic protein 2 |
| LOC102402176 | 16 | Ubiquilin-1 or Ubiquitin-like proteins |
| LOC102397107 | 16 | Ubiquilin-3 |
| BLM | 20 | BLM RecQ-like helicase |
| CRTC3 | 20 | CREB-regulated transcription coactivator 3 |
| LOC102389067 | 22 | Serpin B8 |
| CDH7 | 22 | Cadherin 7 |
1: The top 10 SNP-windows explained the highest genomic variance; 2: All types of the candidate genes are known as protein-coding genes except one gene described as pseudo protein; dynein light chain 1, cytoplasmic-like; Source: NCBI “
References
1. Humberto, A.; Minervino, H.; Zava, M.; Vecchio, D.; Borghese, A. Bubalus Bubalis: A Short Story. Front. Vet. Sci.; 2020; 7, 570413. [DOI: https://dx.doi.org/10.3389/fvets.2020.570413]
2. Vale, W.G. Effects of Environment on Buffalo Reproduction. Ital. J. Anim. Sci.; 2007; 6, pp. 130-142. [DOI: https://dx.doi.org/10.4081/ijas.2007.s2.130]
3. Abdel-Shafy, H.; Awad, M.A.A.; El-Regalaty, H.; Ismael, A.; El-Assal, S.E.D.; Abou-Bakr, S. A Single-Step Genomic Evaluation for Milk Production in Egyptian Buffalo. Livest. Sci.; 2020; 234, [DOI: https://dx.doi.org/10.1016/j.livsci.2020.103977]
4. Do, D.N.; Janss, L.L.G.; Jensen, J.; Kadarmideen, H.N. SNP Annotation-Based Whole Genomic Prediction and Selection: An Application to Feed Efficiency and Its Component Traits in Pigs. J. Anim. Sci.; 2015; 93, pp. 2056-2063. [DOI: https://dx.doi.org/10.2527/jas.2014-8640] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/26020301]
5. Fulton, J.E. Genomic Selection for Poultry Breeding. Anim. Front.; 2012; 2, pp. 30-36. [DOI: https://dx.doi.org/10.2527/af.2011-0028]
6. El-Halawany, N.; Abdel-Shafy, H.; Shawky, A.E.M.A.; Abdel-Latif, M.A.; Al-Tohamy, A.F.M.; Abd El-Moneim, O.M. Genome-Wide Association Study for Milk Production in Egyptian Buffalo. Livest. Sci.; 2017; 198, pp. 10-16. [DOI: https://dx.doi.org/10.1016/j.livsci.2017.01.019]
7. Abdel-Shafy, H.; Awad, M.A.A.; El-Regalaty, H.; El-Assal, S.E.D.; Abou-Bakr, S. Prospecting Genomic Regions Associated with Milk Production Traits in Egyptian Buffalo. J. Dairy Res.; 2020; 87, pp. 389-396. [DOI: https://dx.doi.org/10.1017/S0022029920000953]
8. Hayes, B.J.; Bowman, P.J.; Chamberlain, A.J.; Goddard, M.E. Invited Review: Genomic Selection in Dairy Cattle: Progress and Challenges. J. Dairy Sci.; 2009; 92, pp. 433-443. [DOI: https://dx.doi.org/10.3168/jds.2008-1646]
9. Wang, H.; Misztal, I.; Aguilar, I.; Legarra, A.; Muir, W.M. Genome-Wide Association Mapping Including Phenotypes from Relatives without Genotypes. Genet. Res.; 2012; 94, pp. 73-83. [DOI: https://dx.doi.org/10.1017/S0016672312000274]
10. Aguilar, I.; Misztal, I.; Johnson, D.L.; Legarra, A.; Tsuruta, S.; Lawlor, T.J. Hot Topic: A Unified Approach to Utilize Phenotypic, Full Pedigree, and Genomic Information for Genetic Evaluation of Holstein Final Score. J. Dairy Sci.; 2010; 93, pp. 743-752. [DOI: https://dx.doi.org/10.3168/jds.2009-2730]
11. Abril-Parreño, L.; Carthy, T.R.; Keogh, K.; Štiavnická, M.; O’Meara, C.; Lonergan, P.; Kenny, D.A.; Fair, S. Genome-Wide Association Study Reveals Candidate Markers Related to Field Fertility and Semen Quality Traits in Holstein-Friesian Bulls. Animal; 2023; 17, 100841. [DOI: https://dx.doi.org/10.1016/j.animal.2023.100841] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/37224615]
12. Hackett, A.J.; Macpherson, J.W. A Method for Differential Staining of Bovine Spermatozoa after Extension in Sterile Milk. Can. Vet. J.; 1965; 6, 117. [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/14337429]
13. Blom, E. Pathological Conditions in Genital Organs and Sperm as a Cause for the Rejection of Breeding Bulls for Import into and Export from Denmark (an Andrologic Retrospective, 1958–1982). Nord. Vet. Med.; 1983; 35, pp. 105-130. [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/6878032]
14. Chang, C.C.; Chow, C.C.; Tellier, L.C.A.M.; Vattikuti, S.; Purcell, S.M.; Lee, J.J. Second-Generation PLINK: Rising to the Challenge of Larger and Richer Datasets. Gigascience; 2015; 4, pp. 1-16. [DOI: https://dx.doi.org/10.1186/s13742-015-0047-8] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/25722852]
15. Luke, S.G. Evaluating Significance in Linear Mixed-Effects Models in R. Behav. Res. Methods; 2017; 49, pp. 1494-1502. [DOI: https://dx.doi.org/10.3758/s13428-016-0809-y] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/27620283]
16. Misztal, I.; Tsuruta, S.; Strabel, T.; Auvray, B.; Druet, T.; Lee, D.H. Blupf90 and Related Programs (Bgf90). Proceedings of the 7th World Congress on Genetics Applied to Livestock Production; Montpellier, France, 19–23 August 2002; pp. 2001-2002.
17. VanRaden, P.M. Efficient Methods to Compute Genomic Predictions. J. Dairy Sci.; 2008; 91, pp. 4414-4423. [DOI: https://dx.doi.org/10.3168/jds.2007-0980] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/18946147]
18. Zhang, Z.; Liu, J.; Ding, X.; Bijma, P.; de Koning, D.J.; Zhang, Q. Best Linear Unbiased Prediction of Genomic Breeding Values Using a Trait-Specific Marker-Derived Relationship Matrix. PLoS ONE; 2010; 5, e12648. [DOI: https://dx.doi.org/10.1371/journal.pone.0012648]
19. Wang, H.; Misztal, I.; Aguilar, I.; Legarra, A.; Fernando, R.L.; Vitezica, Z.; Okimoto, R.; Wing, T.; Hawken, R.; Muir, W.M. Genome-Wide Association Mapping Including Phenotypes from Relatives without Genotypes in a Single-Step (SsGWAS) for 6-Week Body Weight in Broiler Chickens. Front. Genet.; 2014; 5, 134. [DOI: https://dx.doi.org/10.3389/fgene.2014.00134]
20. Jiao, X.; Sherman, B.T.; Huang, D.W.; Stephens, R.; Baseler, M.W.; Lane, H.C.; Lempicki, R.A. DAVID-WS: A Stateful Web Service to Facilitate Gene/Protein List Analysis. Bioinformatics; 2012; 28, pp. 1805-1806. [DOI: https://dx.doi.org/10.1093/bioinformatics/bts251]
21. Salem, M.M.I.; Amin, A.M.S.; Ashour, A.F.; El Nagar, A.G. Estimation of Genetic Parameters for Semen Traits in Egyptian Buffalo Bulls. Trop. Anim. Health Prod.; 2023; 55, 264. [DOI: https://dx.doi.org/10.1007/s11250-023-03680-7]
22. Khattab, A.S.; Zeweil, H.S.; Ahmed, M.H.; Asmaa, Z.M.; Mourd, K.A. A Comparison of Different Selection Indexes for Some Semen Traits on Egyptian Buffalo Bulls. Middle East J. Appl. Sci.; 2015; 5, pp. 368-373.
23. Khattab, A.S.; Adenaike, A.S.; Shymma, M.E.K.; El-Gendy, M.M.; Peters, S.O.; Hussein, A.M. Genetic Evaluation of Semen Traits in Friesian Bulls Raised in Egypt. Theriogenology; 2022; 179, pp. 39-44. [DOI: https://dx.doi.org/10.1016/j.theriogenology.2021.11.011] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/34823060]
24. Mahmoud, K.G.M.; El-Sokary, A.A.E.; El-Roos, M.E.A.A.; Ghaffar, A.D.A.; Nawito, M. Sperm Characteristics in Cryopreserved Buffalo Bull Semen and Field Fertility. Iran. J. Appl. Anim. Sci.; 2013; 3, pp. 777-783.
25. Gabr, A.A.W.; El Basuini, M.F. Effect of Tonophosphan, Zinc Oxide, and Ascorbic Acid on Semen, Sexual Desire, and the Fertility Rate of Egyptian Buffalo Bulls. Ann. Agric. Sci.; 2018; 63, pp. 215-221. [DOI: https://dx.doi.org/10.1016/j.aoas.2018.12.001]
26. Rushdi, H.E.; Saad, M.F.; Saeed, A.M. Association of Chromosomal Aberrations and Semen Quality in Egyptian Buffalo Bulls. Egypt. J. Anim. Prod.; 2017; 54, pp. 95-103. [DOI: https://dx.doi.org/10.21608/ejap.2017.93272]
27. Kumar, D.; Mehta, J.S.; Jerome, A.; Kumar, P.; Kumar, D.; Shivani, B.; Patil, C.S.; Bala, R.; Verma, N.; Nain, S. et al. Genetic Analysis for Semen Quality Traits in Buffalo Bulls. Trop. Anim. Health Prod.; 2023; 55, 313. [DOI: https://dx.doi.org/10.1007/s11250-023-03756-4]
28. Qin, C.; Yin, H.; Zhang, X.; Sun, D.; Zhang, Q.; Liu, J.; Ding, X.; Zhang, Y.; Zhang, S. Genome-Wide Association Study for Semen Traits of the Bulls in Chinese Holstein. Anim. Genet.; 2017; 48, pp. 80-84. [DOI: https://dx.doi.org/10.1111/age.12433]
29. Rezende, F.M.; Dietsch, G.O.; Peñagaricano, F. Genetic Dissection of Bull Fertility in US Jersey Dairy Cattle. Anim. Genet.; 2018; 49, pp. 393-402. [DOI: https://dx.doi.org/10.1111/age.12710]
30. Mueller, J.L.; Skaletsky, H.; Brown, L.G.; Zaghlul, S.; Rock, S.; Graves, T.; Auger, K.; Warren, W.C.; Wilson, R.K.; Page, D.C. Independent Specialization of the Human and Mouse X Chromosomes for the Male Germ Line. Nat. Genet.; 2013; 45, pp. 1083-1087. [DOI: https://dx.doi.org/10.1038/ng.2705]
31. Liu, W.S. Mammalian Sex Chromosome Structure, Gene Content, and Function in Male Fertility. Annu. Rev. Anim. Biosci.; 2019; 7, pp. 103-124. [DOI: https://dx.doi.org/10.1146/annurev-animal-020518-115332]
32. Fortes, M.R.S.; Porto-Neto, L.R.; Satake, N.; Nguyen, L.T.; Freitas, A.C.; Melo, T.P.; Scalez, D.C.B.; Hayes, B.; Raidan, F.S.S.; Reverter, A. et al. X Chromosome Variants Are Associated with Male Fertility Traits in Two Bovine Populations. Genet. Sel. Evol.; 2020; 52, 46. [DOI: https://dx.doi.org/10.1186/s12711-020-00563-5]
33. Fortes, M.R.S.; Reverter, A.; Kelly, M.; Mcculloch, R.; Lehnert, S.A. Genome-Wide Association Study for Inhibin, Luteinizing Hormone, Insulin-like Growth Factor 1, Testicular Size and Semen Traits in Bovine Species. Andrology; 2013; 1, pp. 644-650. [DOI: https://dx.doi.org/10.1111/j.2047-2927.2013.00101.x]
34. Suchocki, T.; Szyda, J. Genome-Wide Association Study for Semen Production Traits in Holstein-Friesian Bulls. J. Dairy Sci.; 2015; 98, pp. 5774-5780. [DOI: https://dx.doi.org/10.3168/jds.2014-8951] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/26051317]
35. Pacheco, H.A.; Rezende, F.M.; Peñagaricano, F. Gene Mapping and Genomic Prediction of Bull Fertility Using Sex Chromosome Markers. J. Dairy Sci.; 2020; 103, pp. 3304-3311. [DOI: https://dx.doi.org/10.3168/jds.2019-17767] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/32063375]
36. Hering, D.M.; Oleński, K.; Ruść, A.; Kaminski, S. Genome-Wide Association Study for Semen Volume and Total Number of Sperm in Holstein-Friesian Bulls. Anim. Reprod. Sci.; 2014; 151, pp. 126-130. [DOI: https://dx.doi.org/10.1016/j.anireprosci.2014.10.022] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/25465359]
37. Ramirez-Diaz, J.; Cenadelli, S.; Bornaghi, V.; Bongioni, G.; Montedoro, S.M.; Achilli, A.; Capelli, C.; Rincon, J.C.; Milanesi, M.; Passamonti, M.M. et al. Identification of Genomic Regions Associated with Total and Progressive Sperm Motility in Italian Holstein Bulls. J. Dairy Sci.; 2023; 106, pp. 407-420. [DOI: https://dx.doi.org/10.3168/jds.2021-21700] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/36400619]
38. Sweett, H.; Fonseca, P.A.S.; Suárez-Vega, A.; Livernois, A.; Miglior, F.; Cánovas, A. Genome-Wide Association Study to Identify Genomic Regions and Positional Candidate Genes Associated with Male Fertility in Beef Cattle. Sci. Rep.; 2020; 10, 20102. [DOI: https://dx.doi.org/10.1038/s41598-020-75758-3]
39. Dickins, R.A.; Frew, I.J.; House, C.M.; O’Bryan, M.K.; Holloway, A.J.; Haviv, I.; Traficante, N.; de Kretser, D.M.; Bowtell, D.D.L. The Ubiquitin Ligase Component Siah1a Is Required for Completion of Meiosis I in Male Mice. Mol. Cell. Biol.; 2002; 22, pp. 2294-2303. [DOI: https://dx.doi.org/10.1128/MCB.22.7.2294-2303.2002]
40. Weon, J.L.; Potts, P.R. The MAGE Protein Family and Cancer. Curr. Opin. Cell Biol.; 2015; 37, pp. 1-8. [DOI: https://dx.doi.org/10.1016/j.ceb.2015.08.002]
41. Lee, A.K.; Potts, P.R. A Comprehensive Guide to the MAGE Family of Ubiquitin Ligases. J. Mol. Biol.; 2017; 429, pp. 1114-1142. [DOI: https://dx.doi.org/10.1016/j.jmb.2017.03.005]
42. Park, S.K.; Yoon, J.; Wang, L.; Shibata, T.K.; Motamedchaboki, K.; Shim, K.J.; Chang, M.S.; Lee, S.H.; Tamura, N.; Hatakeyama, S. et al. Enhancement of Mouse Sperm Motility by Trophinin-Binding Peptide. Reprod. Biol. Endocrinol.; 2012; 10, 101. [DOI: https://dx.doi.org/10.1186/1477-7827-10-101] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/23194061]
43. Lipshultz, L.I.; Pastuszak, A.W.; Goldstein, A.T.; Giraldi, A.; Perelman, M.A. Management of Sexual Dysfunction in Men and Women: An Interdisciplinary Approach; Springer: Berlin/Heidelberg, Germany, 2016; ISBN 9781493931002
44. Wang, R.S.; Yeh, S.; Tzeng, C.R.; Chang, C. Androgen Receptor Roles in Spermatogenesis and Fertility: Lessons from Testicular Cell-Specific Androgen Receptor Knockout Mice. Endocr. Rev.; 2009; 30, pp. 119-132. [DOI: https://dx.doi.org/10.1210/er.2008-0025]
45. Ito, H.; Langenhorst, T.; Ogden, R.; Inoue-Murayama, M. Androgen Receptor Gene Polymorphism in Zebra Species. Meta Gene; 2015; 5, pp. 120-123. [DOI: https://dx.doi.org/10.1016/j.mgene.2015.06.006]
46. Lin, C.; Tholen, E.; Jennen, D.; Ponsuksili, S.; Schellander, K.; Wimmers, K. Evidence for Effects of Testis and Epididymis Expressed Genes on Sperm Quality and Boar Fertility Traits. Reprod. Domest. Anim.; 2006; 41, pp. 538-543. [DOI: https://dx.doi.org/10.1111/j.1439-0531.2006.00710.x] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/17107514]
47. Sarakul, M.; Elzo, M.A.; Koonawootrittriron, S.; Suwanasopee, T.; Jattawa, D.; Laodim, T. Characterization of Biological Pathways Associated with Semen Traits in the Thai Multibreed Dairy Population. Anim. Reprod. Sci.; 2018; 197, pp. 324-334. [DOI: https://dx.doi.org/10.1016/j.anireprosci.2018.09.002] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/30213568]
48. Ma, P.; Song, N.N.; Cheng, X.; Zhu, L.; Zhang, Q.; Zhang, L.L.; Yang, X.; Wang, H.; Kong, Q.; Shi, D. et al. ZC4H2 Stabilizes RNF220 to Pattern Ventral Spinal Cord through Modulating Shh/Gli Signaling. J. Mol. Cell Biol.; 2020; 12, pp. 337-344. [DOI: https://dx.doi.org/10.1093/jmcb/mjz087] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/31336385]
49. Ma, P.; Mao, B. The Many Faces of the E3 Ubiquitin Ligase, RNF220, in Neural Development and Beyond. Dev. Growth Differ.; 2022; 64, pp. 98-105. [DOI: https://dx.doi.org/10.1111/dgd.12756]
50. Fu, J.; Wang, Y.; Fok, K.L.; Yang, D.; Qiu, Y.; Chan, H.C.; Koide, S.S.; Miao, S.; Wang, L. Anti-ACTL7a Antibodies: A Cause of Infertility. Fertil. Steril.; 2012; 97, pp. 1226-1233.e8. [DOI: https://dx.doi.org/10.1016/j.fertnstert.2012.02.023]
51. Fortes, M.R.S.; Reverter, A.; Hawken, R.J.; Bolormaa, S.; Lehnert, S.A. Candidate Genes Associated with Testicular Development, Sperm Quality, and Hormone Levels of Inhibin, Luteinizing Hormone, and Insulin-like Growth Factor 1 in Brahman Bulls. Biol. Reprod.; 2012; 87, 58. [DOI: https://dx.doi.org/10.1095/biolreprod.112.101089]
52. El-Maghrabi, M.R.; Correiat, J.J.; Heil, P.J.; Pate, T.M.; Cobb, C.E.; Pilkis, S.J. Tissue distribution, immunoreactivity, and physical properties of 6-phosphofructo-2-kinase/fructose-2, 6-bisphosphatase. Proc. Natl. Acad. Sci. USA; 1986; 83, pp. 5005-5009. [DOI: https://dx.doi.org/10.1073/pnas.83.14.5005]
53. Umehara, T.; Tsujita, N.; Shimada, M.; ActiUmehara, T.; Tsujita, N.; Shimada, M. Activation of Toll-like Receptor 7/8 Encoded by the X Chromosome Alters Sperm Motility and Provides a Novel Simple Technology for Sexing Sperm. PLoS Biol.; 2019; 17, e3000398. [DOI: https://dx.doi.org/10.1371/journal.pbio.3000398] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/31408454]
54. Winge, S.B.; Dalgaard, M.D.; Belling, K.G.; Jensen, J.M.; Nielsen, J.E.; Aksglaede, L.; Schierup, M.H.; Brunak, S.; Skakkebæk, N.E.; Juul, A. et al. Transcriptome Analysis of the Adult Human Klinefelter Testis and Cellularity-Matched Controls Reveals Disturbed Differentiation of Sertoli- and Leydig Cells Article. Cell Death Dis.; 2018; 9, 586. [DOI: https://dx.doi.org/10.1038/s41419-018-0671-1] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/29789566]
55. Guan, Y.; Liang, G.; Martin, G.B.; Guan, L.L. Functional Changes in mRNA Expression and Alternative Pre-mRNA Splicing Associated with the Effects of Nutrition on Apoptosis and Spermatogenesis in the Adult Testis. BMC Genomics; 2017; 18, 64. [DOI: https://dx.doi.org/10.1186/s12864-016-3385-8] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/28068922]
56. Sun, Y.; Sun, P.; Hu, Y.; Shan, L.; Geng, Q.; Gong, Y.; Fan, H.; Zhang, T.; Zhou, Y. Elevated Testicular Apoptosis Is Associated with Elevated Sphingosine Driven by Gut Microbiota in Prediabetic Sheep. BMC Biol.; 2022; 20, 121. [DOI: https://dx.doi.org/10.1186/s12915-022-01326-y] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/35606800]
57. Özbek, M.; Hitit, M.; Kaya, A.; Jousan, F.D.; Memili, E. Sperm Functional Genome Associated With Bull Fertility. Front. Vet. Sci.; 2021; 8, 610888. [DOI: https://dx.doi.org/10.3389/fvets.2021.610888]
58. Ma, D.D.; Wang, D.H.; Yang, W.X. Kinesins in Spermatogenesis. Biol. Reprod.; 2017; 96, pp. 267-276. [DOI: https://dx.doi.org/10.1095/biolreprod.116.144113] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/28203733]
59. Shi, C.S.; Kehrl, J.H. Pyk2 Amplifies Epidermal Growth Factor and C-Src-Induced Stat3 Activation. J. Biol. Chem.; 2004; 279, pp. 17224-17231. [DOI: https://dx.doi.org/10.1074/jbc.M311875200]
60. Wang, H.; Hong, X.; Kinsey, W.H. Sperm-Oocyte Signaling: The Role of IZUMO1R and CD9 in PTK2B Activation and Actin Remodeling at the Sperm Binding Site. Biol. Reprod.; 2021; 104, pp. 1292-1301. [DOI: https://dx.doi.org/10.1093/biolre/ioab048]
61. Yang, F.; Chen, F.; Li, L.; Yan, L.; Badri, T.; Lv, C.; Yu, D.; Zhang, M.; Jang, X.; Li, J. et al. Three Novel Players: PTK2B, SYK, and TNFRSF21 Were Identified to Be Involved in the Regulation of Bovine Mastitis Susceptibility via GWAS and Post-Transcriptional Analysis. Front. Immunol.; 2019; 10, 1579. [DOI: https://dx.doi.org/10.3389/fimmu.2019.01579]
62. Pastorekova, S.; Parkkila, S.; Pastorek, J.; Supuran, C.T. Carbonic Anhydrases: Current State of the Art, Therapeutic Applications and Future Prospects. J. Enzyme Inhib. Med. Chem.; 2004; 19, pp. 199-229. [DOI: https://dx.doi.org/10.1080/14756360410001689540]
63. de Camargo, G.M.F.; Aspilcueta-Borquis, R.R.; Fortes, M.R.S.; Porto-Neto, R.; Cardoso, D.F.; Santos, D.J.A.; Lehnert, S.A.; Reverter, A.; Moore, S.S.; Tonhati, H. Prospecting Major Genes in Dairy Buffaloes. BMC Genomics; 2015; 16, 872. [DOI: https://dx.doi.org/10.1186/s12864-015-1986-2] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/26510479]
64. Carvalho, M.E.; Baldi, F.S.; Alexandre, P.A.; Santana, M.H.A.; Ventura, R.V.; Bueno, R.S.; Bonin, M.N.; Rezende, F.M.; Coutinho, L.L.; Eler, J.P. et al. Genomic Regions and Genes Associated with Carcass Quality in Nelore Cattle. Genet. Mol. Res.; 2019; 18, pp. 1-15. [DOI: https://dx.doi.org/10.4238/gmr18226]
65. Ding, Q.; Joshi, P.S.; Xie, Z.H.; Xiang, M.; Gan, L. BARHL2 Transcription Factor Regulates the Ipsilateral/ Contralateral Subtype Divergence in Postmitotic DI1 Neurons of the Developing Spinal Cord. Proc. Natl. Acad. Sci. USA; 2012; 109, pp. 1566-1571. [DOI: https://dx.doi.org/10.1073/pnas.1112392109] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/22307612]
66. Zhang, L.; Zhao, X.; Wang, W. LncRNA and MRNA Sequencing of the Left Testis in Experimental Varicocele Rats Treated with Morinda Officinalis Polysaccharide. Exp. Ther. Med.; 2021; 22, 1136. [DOI: https://dx.doi.org/10.3892/etm.2021.10570] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/34466146]
67. Puglisi, R.; Montanari, M.; Chiarella, P.; Stefanini, M.; Boitani, C. Regulatory Role of BMP2 and BMP7 in Spermatogonia and Sertoli Cell Proliferation in the Immature Mouse. Eur. J. Endocrinol.; 2004; 151, pp. 511-520. [DOI: https://dx.doi.org/10.1530/eje.0.1510511]
68. Conklin, D.; Holderman, S.; Whitmore, T.E.; Maurer, M.; Feldhaus, A.L. Molecular Cloning, Chromosome Mapping and Characterization of UBQLN3 a Testis-Specific Gene That Contains an Ubiquitin-like Domain. Gene; 2000; 249, pp. 91-98. [DOI: https://dx.doi.org/10.1016/S0378-1119(00)00122-0]
69. Zhao, W.; Li, Z.; Ping, P.; Wang, G.; Yuan, X.; Sun, F. Outer Dense Fibers Stabilize the Axoneme to Maintain Sperm Motility. J. Cell. Mol. Med.; 2018; 22, pp. 1755-1768. [DOI: https://dx.doi.org/10.1111/jcmm.13457]
70. Abu-Halima, M.; Belkacemi, A.; Ayesh, B.M.; Simone Becker, L.; Sindiani, A.M.; Fischer, U.; Hammadeh, M.; Keller, A.; Meese, E. MicroRNA-Targeting in Spermatogenesis: Over-Expressions of MicroRNA-23a/b-3p and Its Affected Targeting of the Genes ODF2 and UBQLN3 in Spermatozoa of Patients with Oligoasthenozoospermia. Andrology; 2021; 9, pp. 1137-1144. [DOI: https://dx.doi.org/10.1111/andr.13004]
© 2023 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https://creativecommons.org/licenses/by/4.0/). Notwithstanding the ProQuest Terms and Conditions, you may use this content in accordance with the terms of the License.
Details
; Salem, Mohamed M I 2
; Amin, Amin M S 3 ; Khalil, Maher H 1 ; Ashour, Ayman F 3 ; Hegazy, Mohammed M 3
; Abdel-Shafy, Hamdy 4
1 Department of Animal Production, Faculty of Agriculture at Moshtohor, Benha University, Benha 13736, Egypt;
2 Department of Animal and Fish Production, Faculty of Agriculture (El-Shatby), Alexandria University, Alexandria 21545, Egypt;
3 Animal Production Research Institute, Agricultural Research Center, Dokki, Giza 12619, Egypt;
4 Department of Animal Production, Faculty of Agriculture, Cairo University, El-Gamma Street, Giza 12613, Egypt;