Next Article in Journal
Clinical Mass Spectrometry Discovered Human IgG Sialylation as a Potential Biosignature for Kidney Function
Previous Article in Journal
Factors and Prognostic Significance of Impaired Exercise Tolerance in Women over 40 with Arterial Hypertension
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
JPM
Volume 11
Issue 8
10.3390/jpm11080760
jpm-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Genetic Association of the Functional WDR4 Gene in Male Fertility

by
Yu-Jia Wang
1,†,
Eko Mugiyanto
2,†,
Yun-Ting Peng
1,3,†,
Wan-Chen Huang
4,
Wan-Hsuan Chou
1,5,
Chi-Chiu Lee
4,
Yu-Shiuan Wang
2,
Lalu Muhammad Irham
1,6,
Dyah Aryani Perwitasari
6,
Ming-I Hsu
7,* and
Wei-Chiao Chang
1,3,5,8,*
1
Department of Clinical Pharmacy, School of Pharmacy, College of Pharmacy, Taipei Medical University, Taipei 11031, Taiwan
2
PhD Program in Clinical Drug Development of Chinese Herbal Medicine, College of Pharmacy, Taipei Medical University, Taipei 11031, Taiwan
3
Department of Pharmacy, Wan Fang Hospital, Taipei Medical University, Taipei 11696, Taiwan
4
Institute of Cellular and Organismic Biology, Academia Sinica, Taipei 115, Taiwan
5
Master Program for Clinical Pharmacogenomics and Pharmacoproteomics, School of Pharmacy, Taipei Medical University, Taipei 11031, Taiwan
6
Faculty of Pharmacy, University of Ahmad Dahlan, Yogyakarta 55164, Indonesia
7
Department of Obstetrics and Gynaecology, Wan Fang Hospital, Taipei Medical University, Taipei 11696, Taiwan
8
Integrative Research Center for Critical Care, Wan Fang Hospital, Taipei Medical University, Taipei 11696, Taiwan
*
Authors to whom correspondence should be addressed.
†
These authors equally contributed to this work.
J. Pers. Med. 2021, 11(8), 760; https://0-doi-org.brum.beds.ac.uk/10.3390/jpm11080760
Submission received: 18 May 2021 / Revised: 23 July 2021 / Accepted: 27 July 2021 / Published: 30 July 2021
(This article belongs to the Section Pharmacogenetics)
Browse Figures
Review Reports Versions Notes

Abstract

:
Infertility is one of the important problems in the modern world. Male infertility is characterized by several clinical manifestations, including low sperm production (oligozoospermia), reduced sperm motility (asthenozoospermia), and abnormal sperm morphology (teratozoospermia). WDR4, known as Wuho, controls fertility in Drosophila. However, it is unclear whether WDR4 is associated with clinical manifestations of male fertility in human. Here, we attempted to determine the physiological functions of WDR4 gene. Two cohorts were applied to address this question. The first cohort was the general population from Taiwan Biobank. Genomic profiles from 68,948 individuals and 87 common physiological traits were applied for phenome-wide association studies (PheWAS). The second cohort comprised patients with male infertility from Wan Fang Hospital, Taipei Medical University. In total, 81 male participants were recruited for the genetic association study. Clinical records including gender, age, total testosterone, follicle-stimulating hormone (FSH), luteinizing hormone (LH), total sperm number, sperm motility, and sperm morphology were collected. In the first cohort, results from PheWAS exhibited no associations between WDR4 genetic variants and 87 common physiological traits. In the second cohort, a total of four tagging single-nucleotide polymorphisms (tSNPs) from WDR4 gene (rs2298666, rs465663, rs2248490, and rs3746939) were selected for genotyping. We found that SNP rs465663 solely associated with asthenozoospermia. Functional annotations through the GTEx portal revealed the correlation between TT or TC genotype and low expression of WDR4. Furthermore, we used mouse embryonic fibroblasts cells from mwdr4 heterozygous (+/‒) mice for functional validation by western blotting. Indeed, low expression of WDR4 contributed to ROS-induced DNA fragmentation. In conclusion, our results suggest a critical role of WDR4 gene variant as well as protein expression in asthenozoospermia.

1. Introduction

Infertility is defined as a disease characterized by the inability of couples to conceive a pregnancy after 12 months of regular unprotected intercourse [1]. Infertility is a global health issue affecting 8~12% of couples worldwide, and male infertility factors are responsible for about 20~30% of infertility cases. Approximately 7% of the male population is affected by infertility [2]. Evaluating male infertility is mainly based on semen analyses. The most significant abnormalities of semen quality include a low sperm number (oligozoospermia), poor sperm motility (asthenozoospermia), and abnormal sperm morphology (teratozoospermia) [3,4]. Oxidative stress has been considered as a critical factor to infertility that possibly associated with the DNA fragmentation [5]. Moreover, the DNA fragmentation levels were reported to correlate with sperm motility [6].
Genetic variations are considered as one of the important factors in male infertility. For example, the azoospermia factor (AZF) region in the long arm of the Y chromosome contains three sub-regions referred to as AZFa, AZFb, and AZFc. Microdeletions in these sub-regions cause impaired sperm production [7]. In addition, HSFY, which is related to the heat-shock transcription factor (HSF) family, has a long open reading frame carrying an HSF-type DNA binding domain involved in azoospermia. Some point mutations in the testis-specific HSFY gene family, such as microdeletions, may be a cause of unexplained cases of idiopathic male infertility [1]. Furthermore, Ozdemir et al. evaluated the genetic variants of SRY and AZF in male infertility [8]. Moreover, three SNPs (rs4920566, rs11763979, and rs3741843) on taste receptor genes (TASR) are associated with male infertility [9]. Several genes were further identified to correlate with infertility-related azoospermia in humans, such as RNU7-6P, ZFP64 [10], SOX5, PLCH2 [11], and ART3 genes [12].
Wu et al. reported that the wh gene (Wuho gene with WD40 repeats, also named WDR4 in homo sapiens) is important in Drosophila spermatogenesis [13]. However, the physiological functions as well as the clinical manifestation for genetic variations of WDR4 gene are still unclear. Here, two cohorts were applied to address this question. The first cohort was the general population from Taiwan Biobank used for phenome-wide association study (PheWAS). The association between WDR4 variants and 87 common traits was comprehensively screened. The second cohort was male participants recruited from Wan Fang Hospital. The effects of WDR4 variants on human male infertility (oligozoospermia, asthenozoospermia, teratozoospermia) were investigated.

2. Materials and Methods

2.1. Participant Recruitment and Sample Collection

Men aged ≥ 20 years old were recruited from the infertility clinic of Wan Fang Hospital, Taipei Medical University between November 2013 and June 2015. Participants with a history of vasectomy or cancer were excluded from the study. Semen samples were collected from every participant by masturbation after 3~5 days of sexual abstinence. These samples were kept at room temperature and delivered to the reproductive laboratory for semen analysis within 1 h. Peripheral blood samples were collected for serum hormone analysis and genomic DNA extraction. The study protocol was approved by the Taipei Medical University-Joint Institutional Review Board (TMU-JIRB) (No: 201302040), and written informed consent was received from all participants.

2.2. Semen Analysis

Semen samples were examined after liquefaction. Semen analysis was conducted separately by two experienced technicians in the reproductive laboratory. Semen parameters including total sperm number, sperm motility, and sperm morphology were evaluated. The World Health Organization (WHO) 2010 recommendations for semen analysis were followed [14]. A semen sample (10 µL) was first transferred to a MAKLER counting chamber (Irvine Scientific, Santa Ana, CA, USA) to measure total number of spermatozoa. Oligozoospermia was defined as a number of sperm < 39 × 106 spermatozoa. The motility of each spermatozoon was graded into three categories: progressive motility (PR), non-progressive motility (NP), and immotile spermatozoa. Asthenozoospermia was defined as PR < 32%. For sperm morphological evaluation, each spermatozoon was classified into eight groups as follows: normal sperm (NS), amorphous head (AH), duplicated head (DH), large head (LH), small head (SH), normal head with other defect (NH), tapering head (TH), and other defect (OD). Teratozoospermia was defined as <4% of morphologically normal spermatozoa. In addition, combined oligozoospermia was defined as participants with sperm number less than 39 × 106 spermatozoa and other defects (astheno- and/or teratozoospermia), a rationale also applied to two other sperm characteristics (combined asthenozoospermia and combined teratozoospermia).

2.3. Genomic DNA Extraction

Peripheral blood samples were first centrifuged at 3000 rpm and 4 °C for 10 min to separate serum and blood cells. The buffy coat layer was extracted and washed with red blood cell (RBC) lysis buffer to isolate peripheral blood mononuclear cells (PBMCs). PBMCs were then lysed using a cell lysis buffer. Proteins were precipitated using a protein precipitation solution followed by 95% isopropanol and 80% alcohol to isolate total genomic DNA. Finally, the DNA purity and concentration were measured using NanoDrop (Thermo Fisher Scientific, Wilmington, DE, USA).

2.4. Genotyping of SNPs in the WDR4 Gene

Tagging (t)SNPs with a minimum allele frequency (MAF) of >10% in a Beijing Han Chinese (CHB) population were selected through UCSC (http://genome.ucsc.edu) and HapMap vers. 2010-08_phase II + III (http://hapmap.ncbi.nlm.hig.gov/). Haploview 4.2 was applied for tSNP selection. A total of four SNPs of the WDR4 gene were selected and genotyped (Figure 1). Characteristics of these SNPs are shown in Table S1. Genotyping was performed using the TaqMan Allelic Discrimination Assay (Applied Biosystems, Foster City, CA, USA). A polymerase chain reaction (PCR) was carried out with an ABI StepOnePlus Thermal Cycler (Applied Biosystems). The fluorescence from different probes was measured and analyzed by System SDS software vers. 2.2.2 (Applied Biosystems).

2.5. Functional Annotation Data Query

Tissue-specific cis-expression quantitative trait loci (cis-eQTL) were queried from the GTEx Portal (http://www.gtexportal.org/home/) to evaluate the effects of the SNPs on gene expressions in different human tissues. Data were obtained from the Genotype-Tissue Expression (GTEx) Portal on 4 May 2020.

2.6. Phenome-Wide Association Study (PheWAS)

A phenome-wide association study (PheWAS) was conducted to screen the association between WDR4 genetic variation and common traits. Imputed Axiom Genome-Wide TWB 2.0 Array data from 68,948 individuals were obtained from the Taiwan Biobank. The data were subjected to quality control based on a variant call rate (>98%), sex check, sample call rate (>98%), Hardy–Weinberg equilibrium (HWE; p > 10−6), heterozygosity (within the mean ± 3 standard deviations (SDs)), and identical by descent (IBD) check (<0.1875). After quality control and extracting common variants (with minor allele frequencies of >0.05), 3,982,815 variants and 59,448 individuals remained. Furthermore, 105 variants of the WDR4 gene were extracted to conduct the PheWAS with 87 phenotypes using Plink v1.9 (www.cog-genomics.org/plink/1.9/) [15,16].

2.7. Cell Culture

The mouse embryonic fibroblasts (MEF) used in this study were purchased from mWh heterozygous mice (+/−) [17] and maintained in DMEM (GIBCO, 11995) contained with 10% FBS, 1% Pen/Strep (GIBCO, 15140122). The 3 × 105 cells were seeded in 6 cm culture dish for 24 h before treated with H2O2.

2.8. Western Blotting

To detect the reactive oxygen species (ROS)-caused DNA damage, cells were treated with 100 μM H2O2 for 24 h. Cells were then lysed, and the total proteins were analyzed by 10% SDS–polyacrylamide gel. After running the gel, proteins were transferred to a PVDF membrane and blocked with 5% non-fat dry milk for 1 h at RT. Membranes were then washed with PBST (PBS contained with 0.1% Tween 20) and incubated with primary antibodies in 1:1000 dilutions overnight at 4 °C; second antibodies were incubated in 1:5000 dilutions for 1 h at RT. The peptides LKKKRQRSPFPGSPEQTK were synthesized to detect mouse Wdr4 protein (mWdr4). The antibody was purified by affinity chromatography with peptide antigens before use [17]. The protein band intensities were detected by an ECL-plus detection system (GE Healthcare, RPN2235).

2.9. Statistical Analysis

For the PheWAS, plink v1.9 was used for analysis (www.cog-genomics.org/plink/1.9/) [15,16]. Furthermore, the simpleM method was used for variants (18 independent variants after correction), and Pearson’s correlation was used for phenotypes (76 independent phenotypes after correction with r2 > 0.7) to adjust the number of independent tests for multiple testing correction in PheWAS [18]. In this study, tests with p < 3.65 × 10−5 (0.05/18 × 76) were regarded as significant for phenome-wide associations. For the WDR4 genetic association study with male infertility, R 3.2.0 (http://www.r-project.org) was used for the statistical analysis. HWE of SNPs was evaluated using a Chi-squared test. Magnitudes of associations between SNPs and semen quality parameters were examined through a linear regression analysis under a recessive model. Age, total testosterone, follicle-stimulating hormone (FSH), and luteinizing hormone (LH) were included as covariates in the regression model. p < 0.05 was considered statistically significant.

3. Results

3.1. PheWAS for WDR4 Variants in A Taiwanese Population

To further explore the roles of WDR4 in human traits, a PheWAS was conducted to investigate associations between WDR4 variants and common traits in a Taiwanese population. The associations between 105 variants of WDR4 gene and 87 phenotypes were tested through the PheWAS. A variety of associations between the WDR4 genetic variants and the common phenotypes were identified at the significance level of p < 0.05. However, none of them remained significant after multiple testing correction (Figure 2).

3.2. Participant Characteristics

Eighty-one males were recruited from the infertility clinic. Participants were aged 28~51 years, with a mean age of 36 years. Men with oligozoospermia, asthenozoospermia, and teratozoospermia accounted for 15%, 27%, and 7%, respectively. Their clinical characteristics are summarized in Table 1. Furthermore, correlation coefficient of age, serum hormones, and semen parameters were shown by Spearman’s rank correlation. In agreement with the previous study, our data demonstrated morphology decline with age, whereas follicle-stimulating hormone (FSH) levels rose (Figure S1) [19].

3.3. Associations between WDR4 Variants and Oligozoospermia

The effects of four SNPs on total sperm number were first evaluated. There were no statistically significant differences in genotypic distributions between oligozoospermia cases and the non-oligozoospermia group (Table 2).

3.4. Associations between WDR4 Variants and Asthenozoospermia

Next, the effects of four SNPs on total sperm motility were analyzed. As shown in Table 3, the genotypic distribution of rs465663 differed between asthenozoospermia and non-asthenozoospermia groups. The proportions of TT and TC genotypes were higher in asthenozoospermia compared to non-asthenozoospermia (p = 0.025).

3.5. Associations between WDR4 Variants and Teratozoospermia

The effects of four SNPs on normal sperm morphology were also examined. As shown in Table 4, no statistically significant differences in genotypic distributions were observed between teratozoospermia cases and the non-teratozoospermia group.

3.6. Functional Annotation from GTEx Portal for rs465663

The profile of WDR4 gene expression in various tissues was determined by the GTEx portal database (Figure 3). In the genetic association analysis, rs465663, an intronic variant of WDR4, was associated with asthenozoospermia. To further elucidate the possible functions of rs465663, cis-eQTL results were retrieved from the GTEx portal. rs465663 could affect the expression level of several genes in different tissue types. Subjects carrying the T allele showed a lower expression level of WDR4 in testes, whole blood, and esophageal mucosal tissue (Table S2).

3.7. Effects of WDR4 in DNA Fragmentation through γH2AX Expression

Previous studies showed that the sperm DNA fragmentation level was increased by ROS stress and was correlated with sperm motility [5,6]. The higher expression level of γH2AX (biomarker of DNA fragmentation level) was revealed to correlate with male infertility [20]. To investigate the role of WDR4 in male infertility, the MEF (mouse embryonic fibroblasts) cells were isolated from mwdr4 heterozygous (+/−) mice, and the H2O2 was used to induce DNA fragmentation for mimicking oxidative stress. The mwdr4 and the γH2AX protein expression levels were confirmed with western blotting (Figure 4). The results showed that the mWdr4 protein expression level in the mWdr4 heterozygous MEF group was low. Importantly, the H2O2-induced γH2AX protein level was highly increased in the mwdr4 heterozygous MEF group compared to the wild type MEF group. These results highlighted an important role of WDR4 in male infertility.

4. Discussion

Wuho is a member of the evolutionarily conserved WD repeat protein family that is expressed by the genes wuho in Drosophila, TRM82 in yeast, and WDR4 in humans [17]. The WDR4 domains usually contain four to eight repeating sequences, which are separated by approximately 40 amino acids. Each repeat consists of two sites, a poorly conserved site and a well-conserved site [21]. Wu at al. reported that the lack of WDR4 function is associated with dramatic germline-specific phenotypes by arresting the spermatogenesis at the spermatid elongating stage [13]. The study indicated that approximately 20% of the ovarioles in WDR4 mutant female have apparent defects in oogenesis with an over-proliferation of cystocytes. Additionally, WDR4 is associated with germline cell development through cytosolic tRNA modifications [22].
Male infertility is characterized (manifested) by low sperm production (oligozoospermia), reduced sperm motility (asthenozoospermia), and abnormal sperm morphology (teratozoospermia). In this study, we wanted to seek out whether any genetic variants associate with the clinical outcomes. Importantly, the analysis using continuous traits showed consistent results with dichotomous models (Tables S3–S5). Asthenozoospermia and oligozoospermia are the most common factors responsible for male infertility [23]. In line with this, our study indicated that asthenozoospermia comprised the most cases followed consecutively by oligozoospermia and teratozoospermia. Here, we identified variants of the WDR4 gene that related to male infertility. rs465663 significantly associated with the susceptibility to asthenozoospermia. The variation in rs465663 was located in the intronic region of chromosome 21, while proportions of the TT and the TC genotypes of rs465663 were higher than the CC genotype in cases with asthenozoospermia. Interestingly, the C allele is the minor allele of rs465663. The distribution of the minor allele of rs465663 in current study was lower than those reported in other populations, including Asians, Europeans, and Americans (Table S1). Trends of infertility globally showed that two populations of Europeans, represented by Central and Eastern Europe (8~12%), and Americans, represented by North America (4.5~6%), were higher compared to Sub-Saharan Africans (2.5~4.8%); unfortunately, those studies did not provide the percentage of Asians, possibly due to underreporting [24].
Our findings also emphasized that functional annotations through a bioinformatic approach using the GTEx portal revealed that subjects had lower expression of the TT genotype of WDR4 than the CC and the TC genotypes in several human tissues, including testes and whole blood. Meanwhile, the major allelic frequency of rs465663 of WDR4 had lower expression in testes. This result implied that rs465663 might influence expression of the WDR4 gene, especially in testes. Disease-related male infertility can be caused by testicular deficiencies and spermatogenesis [4]. This evidence provided additional clues that variations of the WDR4 gene may affect the testes, including their production of sperm (Figure 3). Furthermore, we noticed that ROS-induced DNA fragmentation level was significantly increased in the low WDR4 group (Figure 4). The results further support the possibility that lower expression of WDR4 in males resulted in infertility.
Several phenotypes are associated with WDR4 variants, such as rs370189685, which was correlated with the fasting plasma glucose level [25]. A study from Drosophila germ cells suggested that WDR4 is a regulator of Mei-p26, and it interacts with TRIM32 to control tissue homeostasis in other stem cell systems [26]. Other genes related to male infertility were identified to associate with azoospermia (e.g., SOHLH1 [27], SYCP3 [28], and TEX11 [29]). WDR4 gene with a missense mutation was reported to be related to primordial dwarfism through m7G46 methylation, which impaired transfer (t)RNAs [30]. Sperm carries thousands of different RNAs [31]. Interestingly, according to the REACTOME database (R-HSA-6782315) [32], variations in rs465663 of the WDR4 gene might involve in RNA metabolism and tRNA modifications in nuclei and cytosol. tRNAs play pivotal roles in protein synthetases. Mutations in tRNAs which modify enzymes are associated with human diseases, including cancer, type 2 diabetes (T2D), neurological disorders, and mitochondrion-linked disorders [33]. Furthermore, as reported by a previous study, tRNA modifications can affect proteostasis in humans [34]. A similar study supported that tRNA modifications in mice resulted in increased apoptosis in male germ cells and male infertility [35]. In addition, epigenomic assays are also essential for relating noncoding genetic variations to regulatory mechanisms underlying phenotypic changes, including genomic variations of the WDR4 gene. The Encyclopedia of DNA Elements (ENCODE) database revealed that rs465663 is located in an intronic region with known histone modifications [32]. Subsequent evidence revealed that an epigenomic histone modification was involved in male fertility [36]. Some genes, such as FAM50B and GNAS, were reported to involve in the quality of sperm in asthenozoospermia through histone modification-type methylation [37]. Taken together, previous studies provided comprehensive integration between the identified genetic variants and male infertility.
Our research is focusing on the genetic variants of WDR4 in male fertility s. However, some limitations still exist in this study. First, although the subject number (68,948 individuals) for PheWAS was good, no significant physiological traits were found. Regarding the second cohort (male participants) from the hospital, the small sample size limits the statistical power. Thus, larger sample sizes with different populations are necessary to confirm our findings. Second, as the candidate gene approach was performed, the influence of variants in other unexamined genes cannot be ruled out. Since the pathophysiology of infertility in humans is complexly regulated by many signaling pathways, using whole genome sequencing technology may yield further insights into the genomic variations in infertility.

5. Conclusions

This study offers important information related to genetic variants and expression level of the WDR4 which might affect infertility in the cases of asthenozoospermia. However, further functional studies and larger sample sizes are required to validate the variants.

Supplementary Materials

The following are available online at https://0-www-mdpi-com.brum.beds.ac.uk/article/10.3390/jpm11080760/s1, Figure S1: Spearman’s rank correlation coefficient between age, serum hormones, and semen parameters, Table S1: Characteristics of tagging single-nucleotide polymorphisms (tSNPs) of the WDR4 gene, Table S2: Cis-expression quantitative trait loci (cis-eQTL) of rs465663, Table S3: Associations between genetic variants of the WDR4 gene and total sperm number, Table S4: Associations between genetic variants of the WDR4 gene and sperm progressive motility, Table S5: Associations between genetic variants of the WDR4 gene and sperm morphology.

Author Contributions

Conceptualization: Y.-J.W., Y.-T.P., W.-H.C., M.-I.H. and W.-C.C.; sample collection: Y.-J.W. and Y.-T.P.; data curation: Y.-J.W., Y.-T.P., M.-I.H. and W.-C.C.; formal analyses: Y.-J.W., Y.-T.P., W.-H.C., M.-I.H. and W.-C.C.; experimental preparation: C.-C.L., Y.-S.W., W.-H.C.; data interpretation and discussion: Y.-J.W., Y.-T.P., E.M., W.-C.H., W.-H.C., C.-C.L., Y.-S.W., L.M.I., D.A.P., M.-I.H. and W.-C.C.; writing review & editing: Y.-J.W., Y.-T.P., E.M., W.-C.H., W.-H.C., C.-C.L., Y.-S.W., L.M.I., D.A.P., M.-I.H. and W.-C.C.; supervision: W.-C.H., D.A.P., M.-I.H. and W.-C.C. All authors have made significant contributions to this study. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by grants from the Health and welfare surcharge of tobacco products grant (MOHW110-TDU-B-212-144014; MOHW110-TDU-B-212-144020), Ministry of Science and Technology, Taiwan (MOST109-2314-B-038-131 and MOST110-2628-B-038-020), and Taipei Medical University, Taiwan (12310-106079; Yusuke Nakamura Chair Professorship).

Institutional Review Board Statement

This study was approved by the Taipei Medical University-Joint Institutional Review Board (TMU-JIRB no: 201302040). Samples were collected from participants after they had provided written informed consent.

Informed Consent Statement

Informed consent was obtained from all subjects involved in the study.

Acknowledgments

The paper is in memory of academician, Tao-shih Hsieh, Academia Sinica, Taiwan. Thanks to his guidance and initiation of this project.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Zegers-Hochschild, F.; Adamson, G.; de Mouzon, J.; Ishihara, O.; Mansour, R.; Nygren, K.; Sullivan, E.; Vanderpoel, S. International Committee for Monitoring Assisted Reproductive Technology (ICMART) and the World Health Organization (WHO) revised glossary of ART terminology, 2009. Fertil. Steril. 2009, 92, 1520–1524. [Google Scholar] [CrossRef]
  2. Vander Borght, M.; Wyns, C. Fertility and infertility: Definition and epidemiology. Clin. Biochem. 2018, 62, 2–10. [Google Scholar] [CrossRef] [PubMed]
  3. Esteves, S.C. Clinical relevance of routine semen analysis and controversies surrounding the 2010 World Health Organization criteria for semen examination. Int. Braz. J. Urol. 2014, 40, 433–453. [Google Scholar] [CrossRef] [PubMed]
  4. Jungwirth, A.; Giwercman, A.; Tournaye, H.; Diemer, T.; Kopa, Z.; Dohle, G.; Krausz, C.G. European Association of Urology Guidelines on Male Infertility: The 2012 Update. Eur. Urol. 2012, 62, 324–332. [Google Scholar] [CrossRef] [PubMed]
  5. Moustafa, M.H.; Sharma, R.; Thornton, J.; Mascha, E.; Abdel-Hafez, M.A.; Thomas, A.J.; Agarwal, A. Relationship between ROS production, apoptosis and DNA denaturation in spermatozoa from patients examined for infertility. Hum. Reprod. 2004, 19, 129–138. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  6. Elbashir, S.; Magdi, Y.; Rashed, A.; Ibrahim, M.A.; Edris, Y.; Abdelaziz, A.M. Relationship between sperm progressive motility and DNA integrity in fertile and infertile men. Middle East Fertil. Soc. J. 2018, 23, 195–198. [Google Scholar] [CrossRef]
  7. Hopps, C.V.; Mielnik, A.; Goldstein, M.; Palermo, G.D.; Rosenwaks, Z.; Schlegel, P.N. Detection of sperm in men with Y chromosome microdeletions of the AZFa, AZFb and AZFc regions. Hum. Reprod. 2003, 18, 1660–1665. [Google Scholar] [CrossRef] [Green Version]
  8. Ozdemir, O.; Gül, E.; Kilicarslan, H.; Gokce, G.; Beyaztas, F.Y.; Ayan, S.; Sezgin, I. SRY and AZF gene variation in male infertility: A cytogenetic and molecular approach. Int. Urol. Nephrol. 2007, 39, 1183–1189. [Google Scholar] [CrossRef]
  9. Gentiluomo, M.; Crifasi, L.; Luddi, A.; Locci, D.; Barale, R.; Piomboni, P.; Campa, D. Taste receptor polymorphisms and male infertility. Hum. Reprod. 2017, 32, 2324–2331. [Google Scholar] [CrossRef]
  10. Aston, K.I.; Carrell, D.T. Genome-wide study of single-nucleotide polymorphisms associated with azoospermia and severe oligozoospermia. J. Androl. 2009, 30, 711–725. [Google Scholar] [CrossRef]
  11. Hu, Z.; Xia, Y.; Guo, X.; Dai, J.; Li, H.; Hu, H.; Jiang, Y.; Lu, F.; Wu, Y.; Yang, X.; et al. A genome-wide association study in Chinese men identifies three risk loci for non-obstructive azoospermia. Nat. Genet. 2011, 44, 183–186. [Google Scholar] [CrossRef] [PubMed]
  12. Okada, H.; Tajima, A.; Shichiri, K.; Tanaka, A.; Tanaka, K.; Inoue, I. Genome-Wide Expression of Azoospermia Testes Demonstrates a Specific Profile and Implicates ART3 in Genetic Susceptibility. PLoS Genet. 2008, 4, e26. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  13. Wu, J.; Hou, J.H.; Hsieh, T.-S. A new Drosophila gene wh (wuho) with WD40 repeats is essential for spermatogenesis and has maximal expression in hub cells. Dev. Biol. 2006, 296, 219–230. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  14. WHO. WHO Laboratory Manual for the Examination and Processing of Human Semen, 5th ed.; World Health Organization: Geneva, Switzerland, 2010. [Google Scholar]
  15. Purcell, S.; Chang, C. Plink v1.9. Available online: https://www.cog-genomics.org/plink/ (accessed on 18 May 2021).
  16. 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, 7. [Google Scholar] [CrossRef]
  17. Cheng, I.-C.; Chen, B.C.; Shuai, H.-H.; Chien, F.-C.; Chen, P.; Hsieh, T.-S. Wuho Is a New Member in Maintaining Genome Stability through its Interaction with Flap Endonuclease 1. PLoS Biol. 2016, 14, e1002349. [Google Scholar] [CrossRef]
  18. Gao, X.; Starmer, J.; Martin, E.R. A multiple testing correction method for genetic association studies using correlated single nucleotide polymorphisms. Genet. Epidemiol. 2008, 32, 361–369. [Google Scholar] [CrossRef]
  19. Pasqualotto, F.F.; Sobreiro, B.P.; Hallak, J.; Pasqualotto, E.B.; Lucon, A.M. Sperm concentration and normal sperm morphology decrease and follicle-stimulating hormone level increases with age. BJU Int. 2005, 96, 1087–1091. [Google Scholar] [CrossRef]
  20. Zhong, H.-Z.; Lv, F.-T.; Deng, X.-L.; Hu, Y.; Xie, D.-N.; Lin, B.; Mo, Z.-N.; Lin, F.-Q. Evaluating γH2AX in spermatozoa from male infertility patients. Fertil. Steril. 2015, 104, 574–581. [Google Scholar] [CrossRef]
  21. Riedl, S.J.; Salvesen, G.S. The apoptosome: Signalling platform of cell death. Nat. Rev. Mol. Cell Biol. 2007, 8, 405–413. [Google Scholar] [CrossRef]
  22. Tahmasebi, S.; Khoutorsky, A.; Mathews, M.B.; Sonenberg, N. Translation deregulation in human disease. Nat. Rev. Mol. Cell Biol. 2018, 19, 791–807. [Google Scholar] [CrossRef]
  23. Ikechebelu, J.; Adinma, B.; Orie, E.; Ikegwuonu, S. High prevalence of male infertility in southeastern Nigeria. J. Obstet. Gynaecol. 2003, 23, 657–659. [Google Scholar] [CrossRef] [PubMed]
  24. Agarwal, A.; Mulgund, A.; Hamada, A.; Chyatte, M.R. A unique view on male infertility around the globe. Reprod. Biol. Endocrinol. 2015, 13, 37. [Google Scholar] [CrossRef] [Green Version]
  25. Nagy, R.; Boutin, T.S.; Marten, J.; Huffman, J.E.; Kerr, S.M.; Campbell, A.; Evenden, L.; Gibson, J.; Amador, C.; Howard, D.M.; et al. Exploration of haplotype research consortium imputation for genome-wide association studies in 20,032 Generation Scotland participants. Genome Med. 2017, 9, 23. [Google Scholar] [CrossRef] [Green Version]
  26. Rastegari, E.; Kajal, K.; Tan, B.-S.; Huang, F.; Chen, R.-H.; Hsieh, T.-S.; Hsu, H.-J. WD40 protein Wuho controls germline homeostasis via TRIM-NHL tumor suppressor Mei-p26 in Drosophila. Development 2020, 147, dev182063. [Google Scholar] [CrossRef] [PubMed]
  27. Choi, Y.; Jeon, S.; Choi, M.; Lee, M.-H.; Park, M.; Lee, D.R.; Jun, K.-Y.; Kwon, Y.; Lee, O.-H.; Song, S.-H.; et al. Mutations in SOHLH1 gene associate with nonobstructive Azoospermia. Hum. Mutat. 2010, 31, 788–793. [Google Scholar] [CrossRef] [PubMed]
  28. Miyamoto, T.; Hasuike, S.; Yogev, L.; Maduro, M.R.; Ishikawa, M.; Westphal, H.; Lamb, D.J. Azoospermia in patients heterozygous for a mutation in SYCP3. Lancet 2003, 362, 1714–1719. [Google Scholar] [CrossRef]
  29. Yatsenko, A.N.; Georgiadis, A.P.; Röpke, A.; Berman, A.J.; Jaffe, T.; Olszewska, M.; Westernströer, B.; Sanfilippo, J.; Kurpisz, M.; Rajkovic, A.; et al. X-Linked TEX11 Mutations, Meiotic Arrest, and Azoospermia in Infertile Men. N. Engl. J. Med. 2015, 372, 2097–2107. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  30. Shaheen, R.; Abdel-Salam, G.; Guy, M.P.; AlOmar, R.; Abdel-Hamid, M.S.; Afifi, H.H.; Ismail, S.; Emam, B.A.; Phizicky, E.M.; Alkuraya, F.S. Mutation in WDR4 impairs tRNA m7G46 methylation and causes a distinct form of microcephalic primordial dwarfism. Genome Biol. 2015, 16, 210. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  31. Dadoune, J.-P. Spermatozoal RNAs: What about their functions? Microsc. Res. Tech. 2009, 72, 536–551. [Google Scholar] [CrossRef]
  32. Oscanoa, J.; Sivapalan, L.; Gadaleta, E.; Ullah, A.Z.D.; Lemoine, N.R.; Chelala, C. SNPnexus: A web server for functional annotation of human genome sequence variation (2020 update). Nucleic Acids Res. 2020, 48, W185–W192. [Google Scholar] [CrossRef]
  33. Torres, A.G.; Batlle, E.; de Pouplana, L.R. Role of tRNA modifications in human diseases. Trends Mol. Med. 2014, 20, 306–314. [Google Scholar] [CrossRef] [PubMed]
  34. Pereira, M.; Francisco, S.; Varanda, A.S.; Santos, M.; Santos, M.A.S.; Soares, A.R. Impact of tRNA Modifications and tRNA-Modifying Enzymes on Proteostasis and Human Disease. Int. J. Mol. Sci. 2018, 19, 3738. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  35. Lin, F.-J.; Shen, L.; Jang, C.-W.; Falnes, P.Ø.; Zhang, Y. Ikbkap/Elp1 Deficiency Causes Male Infertility by Disrupting Meiotic Progression. PLoS Genet. 2013, 9, e1003516. [Google Scholar] [CrossRef] [PubMed]
  36. Giacone, F.; Cannarella, R.; Mongioì, L.M.; Alamo, A.; Condorelli, R.A.; Calogero, A.E.; La Vignera, S. Epigenetics of Male Fertility: Effects on Assisted Reproductive Techniques. World J. Men’s Health 2019, 37, 148–156. [Google Scholar] [CrossRef] [PubMed]
  37. Xu, J.; Zhang, A.; Zhang, Z.; Wang, P.; Qian, Y.; He, L.; Shi, H.; Xing, Q.; Du, J. DNA methylation levels of imprinted and nonimprinted genes DMRs associated with defective human spermatozoa. Andrologia 2016, 48, 1027–1035. [Google Scholar] [CrossRef] [PubMed]
Jpm 11 00760 g001 550
Figure 1. Graphic view of the genotyped WDR4 gene.
Figure 1. Graphic view of the genotyped WDR4 gene.
Jpm 11 00760 g001
Jpm 11 00760 g002 550
Figure 2. Phenome-wide scan for associations with WDR4 variants. The blue line represents p = 0.05, and the red line represents the significance level for phenome-wide associations after multiple testing corrections (p = 3.65 × 10−5).
Figure 2. Phenome-wide scan for associations with WDR4 variants. The blue line represents p = 0.05, and the red line represents the significance level for phenome-wide associations after multiple testing corrections (p = 3.65 × 10−5).
Jpm 11 00760 g002
Jpm 11 00760 g003 550
Figure 3. Profile of WDR4 gene expression in various tissues. Blue represents the male gender and red the female gender.
Figure 3. Profile of WDR4 gene expression in various tissues. Blue represents the male gender and red the female gender.
Jpm 11 00760 g003
Jpm 11 00760 g004 550
Figure 4. H2O2-induced γH2AX expression increased in wdr4 heterozygous (+/−) MEF cells. The MEFs cells were subcultured into 6 cm culture dish 24 h before being treated with 100 μM H2O2. The mwdr4, the γH2AX, and the GAPDH protein expression levels were determined by western blot after H2O2 treatment for 24 h.
Figure 4. H2O2-induced γH2AX expression increased in wdr4 heterozygous (+/−) MEF cells. The MEFs cells were subcultured into 6 cm culture dish 24 h before being treated with 100 μM H2O2. The mwdr4, the γH2AX, and the GAPDH protein expression levels were determined by western blot after H2O2 treatment for 24 h.
Jpm 11 00760 g004
Table
Table 1. Characteristics of the 81 male participants.
Table 1. Characteristics of the 81 male participants.
CharacteristicParticipantsNormal Range
Age (years) a36.21 ± 4.58
 Range (years)27~51
Semen analysis, no. (%)
 Oligozoospermia (total sperm number < 39 × 106)
  Oligozoospermia 12 (15)
    Isolated9
    Combined 3
  Non-oligozoospermia69 (85)
 Asthenozoospermia (PR < 32%)
  Asthenozoospermia22 (27)
    Isolated10
    Combined12
  Non-asthenozoospermia59 (73)
 Teratozoospermia (sperm with normal morphology <4%)
  Teratozoospermia 6 (7)
    Isolated0
    Combined6
  Non-teratozoospermia75 (93)
Serum hormone analysis a
 Total testosterone (ng/mL) 4.87 ± 1.592.51~10.60
 FSH (mIU/mL) 5.42 ± 4.241~14
 LH (mIU/mL) 2.40 ± 1.591.5~9.2
 SHBG (nmol/L)28.88 ± 14.8214.5~48.4
 Inhibin B (pg/mL)252.40 ± 162.44—
 TSH (μIU/mL)1.63 ± 0.820.34~5.60
 T3 (ng/dL)112.20 ± 18.0687~178
 T4 (μg/dL)8.15 ± 1.286.09~12.23
 Free T4 (ng/dL)0.90 ± 0.100.61~1.12
 AMH (ng/mL)12.51 ± 7.84—
 Zinc (μg/L)1076.10 ± 217.26800~1200
a Mean ± standard deviation. PR, progressive motility; FSH, follicle-stimulating hormone; LH, luteinizing hormone; SHBG, sex hormone-binding globulin; TSH, thyroid-stimulating hormone; T3, triiodothyronine; T4, thyroxine; AMH, anti-müllerian hormone. — No optimal recommendation available for male.
Table
Table 2. Associations between genetic variants of the WDR4 gene and oligozoospermia.
Table 2. Associations between genetic variants of the WDR4 gene and oligozoospermia.
SNPGenotypeCasesNon-OligozoospermiaRecessive Model
p-Value
No.%No.%
rs2298666GG1083.34870.60.665
GA216.71826.5
AA0022.9
rs465663TT975.03349.30.157
TC325.02435.8
CC001014.9
rs2248490CC866.73043.50.095
CG433.32942.0
GG001014.5
rs3746939AA866.73753.60.165
AC433.3 2739.1
CC0057.2
The p-value was adjusted for age, total testosterone, follicle-stimulating hormone, and luteinizing hormone. SNP, single-nucleotide polymorphism.
Table
Table 3. Associations between genetic variants of the WDR4 gene and asthenozoospermia.
Table 3. Associations between genetic variants of the WDR4 gene and asthenozoospermia.
SNPGenotypeCasesNon-AsthenozoospermiaRecessive Model
p-Value
No.%No.%
rs2298666GG17814169.50.399
GA4191627.1
AA0023.4
rs465663TT1571.42746.60.025 *
TC628.62136.2
CC001017.2
rs2248490CC1568.223390.677
CG522.72847.5
GG29.1813.6
rs3746939AA1568.23050.80.513
AC627.32542.4
CC14.546.8
The p-value was adjusted for age, total testosterone, follicle-stimulating hormone, and luteinizing hormone. * p < 0.05. SNP, single-nucleotide polymorphism.
Table
Table 4. Associations between genetic variants of the WDR4 gene and teratozoospermia.
Table 4. Associations between genetic variants of the WDR4 gene and teratozoospermia.
SNPGenotypeCasesNon-TeratozoospermiaRecessive Model
p-Value
No.%No.%
rs2298666GG583.35371.60.759
GA116.71925.7
AA0022.7
rs465663TT583.33750.70.395
TC116.72635.6
CC001013.7
rs2248490CC583.333440.154
CG116.73242.7
GG001013.3
rs3746939AA466.74154.70.186
AC233.32938.7
CC0056.7
The p-value was adjusted for age, total testosterone, follicle-stimulating hormone, and luteinizing hormone. SNP, single-nucleotide polymorphism.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Wang, Y.-J.; Mugiyanto, E.; Peng, Y.-T.; Huang, W.-C.; Chou, W.-H.; Lee, C.-C.; Wang, Y.-S.; Irham, L.M.; Perwitasari, D.A.; Hsu, M.-I.; et al. Genetic Association of the Functional WDR4 Gene in Male Fertility. J. Pers. Med. 2021, 11, 760. https://0-doi-org.brum.beds.ac.uk/10.3390/jpm11080760

AMA Style

Wang Y-J, Mugiyanto E, Peng Y-T, Huang W-C, Chou W-H, Lee C-C, Wang Y-S, Irham LM, Perwitasari DA, Hsu M-I, et al. Genetic Association of the Functional WDR4 Gene in Male Fertility. Journal of Personalized Medicine. 2021; 11(8):760. https://0-doi-org.brum.beds.ac.uk/10.3390/jpm11080760

Chicago/Turabian Style

Wang, Yu-Jia, Eko Mugiyanto, Yun-Ting Peng, Wan-Chen Huang, Wan-Hsuan Chou, Chi-Chiu Lee, Yu-Shiuan Wang, Lalu Muhammad Irham, Dyah Aryani Perwitasari, Ming-I Hsu, and et al. 2021. "Genetic Association of the Functional WDR4 Gene in Male Fertility" Journal of Personalized Medicine 11, no. 8: 760. https://0-doi-org.brum.beds.ac.uk/10.3390/jpm11080760

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop