Search tips
Search criteria 


Logo of oncotargetLink to Publisher's site
Oncotarget. 2016 December 13; 7(50): 82594–82601.
Published online 2016 October 31. doi:  10.18632/oncotarget.12993
PMCID: PMC5347716

Idiopathic male infertility in the Han population in China is affected by polymorphism in the VDAC2 gene



It has been proved that human voltage-dependent anion channel 2 (VDAC2) plays a significant role in sperm function and male fertility. This study was primarily aimed at exploring whether VDAC2 is a risk factor for idiopathic male infertility.


We determined a significantly increased risk of idiopathic infertility with abnormal semen parameters in association with the variant rs2804535 and a decreased risk of idiopathic infertility with abnormal semen parameters in association with the variant rs11001334. However, among subjects with normal semen parameters, no significant differences could be found in these genotypes. Moreover, we could not find any differences in the variants rs7896741 and rs1259503, which showed no risk of male infertility, whether normal or abnormal.

Materials and methods

All of the experimental subjects, including 523 men who cannot conceive children and 277 fertile controls, underwent complete historical and physical examinations. Each participant donated an ejaculate for semen analysis and 5 ml of peripheral blood for genomic DNA extraction. A computer-assisted semen analysis system was used for the semen analysis. Four single-nucleotide polymorphisms were identified and analyzed using TaqMan SNP Genotyping Assays.


The result shows that the relationships between different variants in the VDAC2 gene and male fertility differ, and the individuals who carry those variants may have a decreased or increased risk of abnormal semen parameters associated with male infertility.

Keywords: VDAC2, genetic diversity, semen parameters, idiopathic male infertility


Studies have been conducted on voltage-dependent anion channel (VDAC) proteins, i.e., channel-forming proteins in the outer mitochondrial membrane (OMM) because of their protective action in modulating metabolite exchange and the transfer of ATP as well as ions between the cytoplasm and mitochondria [1]. There are three isoforms of VDAC, VDAC1, VDAC2 and VDAC3, which show >70% similarity in their sequence. Even though all VDAC isoforms can appear everywhere, the major form is expressed in sperm outer dense fibers (ODFs), an important separate part of the cytoskeleton without a membrane that is involved in flagellar motility [2]; ODFs are mainly focused on the top of the membrane and principal piece of the sperm [35]. These proteins also play roles as transportation proteins in sperm functions by mediating transmembrane calcium transport [4]. A recent study showed that VDAC2 mRNA tends to be more enriched in less motile sperm [6]. A proteomic study revealed that VDAC2 has a strong relationship with male fertility [7]. On the other hand, it is still unknown that whether differences in VDAC2 increase the risk of male fertility or not.

The VDAC2 gene encodes human VDAC2, which is also named outer mitochondrial membrane protein porin 2 (UniProtKB ID: P45880; Genetic ID: 7417). This gene is approximately 20.6 kb in length, is located on chromosome 10q22.2 and includes 13 exons [8]. Using HapMap genotype data, which were obtained from Han Chinese subjects who live in Beijing, four tag single-nucleotide polymorphisms (tagSNPs) (rs2804535, rs7896741, rs11001334 and rs1259503) were selected. Then, we began the analysis of those four genotyping-identified tagSNPs in the VDAC2 gene and measured how those genetic variants affect semen parameters with respect to male infertility; however, all steps were completed with the help of the hospital, which recruited 523 men with infertility and 277 fertile men from the Han-Chinese population.


The mean age of the 523 males in the Han population in this study who could not have babies was 32.0 ± 5.6 years of age (range 20–55 years). The time between experimentation and semen collection was 5.0 ± 2.7 days (range 2–20 days). In addition, the sperm parameters are presented in Table Table1.1. Briefly, there was nothing uncommon or special in the stratification for smoking and drinking for all the features above. However, the sperm number per ejaculate, sperm concentration, and sperm motility exhibited significant differences among the three age groups (P = 0.001, 0.001 and 0.034, respectively). The sperm number per ejaculate and sperm concentration were significantly greater in the 25 kg/m2 group than the < 25 kg/m2 group (P = 0.036 and 0.028, respectively). In regard to the duration of sexual abstinence, the ≥7 days group had a higher semen volume and sperm number per ejaculate than the < 7 days group (P = 0.000 and 0.001, respectively). The remaining parameters (sperm concentration, sperm number per ejaculate, sperm motility, and semen volume) were similar after stratification for age, BMI and the duration of sexual abstinence.

Table 1
Association between the selected individual characteristics and sperm parameters in 523 males with definite idiopathic infertility

Table Table22 illustrates the sperm parameters and the genotype frequency for each variant. The relationship between the sperm number per ejaculate and the subjects was not clear. The sperm concentration, semen volume, and sperm motility were largely similar between the genotypes of the two variants (rs2804535 and rs7896741). For rs11001334, the sperm concentration in subjects with the TT genotype was significantly lower than that in those harboring the CC genotype (P = 0.045) and significantly higher in subjects with the CT and TT genotypes than those harboring the CC genotype (P = 0.037). For rs1259503, the semen volume in subjects with the GC and CC genotypes was significantly lower than that in the subjects harboring the GG genotype (P = 0.039) (Table (Table2).2). As shown in Table Table3,3, we initially divided the cases into the case Group I and case Group II by semen parameters and found that no significant relationships were found for the frequency of all four genotype variants.

Table 2
Sperm parameters according to the genetic variants of VADC2 gene in 523 males with definite idiopathic infertility
Table 3
Genotype frequencies of VADC2 genetic variants among the cases and controls and their association with male infertility

As shown in Table Table4,4, the rs2804535 genotype frequencies were 67.86% (TT), 32.14% (TC) and 0% (CC) in the cases of Subgroup 2 and 85.06% (TT), 14.11% (TC) and 0.83% (CC) in the controls, and they showed statistically significant differences (P = 0.023). A co-dominant allele effect might exist when we put the rs2804535 TC and the rs2804535 CC together, and we can also see a greater number of the rs2804535 TC + CC variant genotypes in the group (P = 0.025). Logistic regression analysis showed us that the frequency of the heterozygous rs2804535 TC genotype (OR = 2.86, 95% CI = 1.19–6.83) and the rs2804535 TC + CC variant genotypes (OR = 2.69, 95% CI = 1.13–6.43) was higher in the subjects in Subgroup 2. However, the frequency of rs2804535 was not higher in the other subgroups.

Table 4
Genotype frequencies of VADC2 genetic variants among the cases and controls and their association with male infertility

Furthermore, the rs11001334 genotype frequencies were 88.10% (CC), 11.61% (CT) and 0.28% (TT) in the cases in Subgroup 1 and 81.74% (CC), 18.26% (CT) and 0% (TT) in the controls, with an obvious statistical difference (P = 0.030). When we combined rs11001334 CT and rs11001334 TT, there was a co-dominant allele effect, and the combination of the rs11001334 CT + TT variant was lower in the experimental group (P = 0.032). Logistic regression analysis indicated that the heterozygous rs11001334 CT genotype (OR = 0.59, 95% CI = 0.37–0.96) and the rs11001334 CT + TT variant genotype (OR=0.61, 95% CI=0.38–0.96) were present at remarkably lower frequency in the subjects in Subgroup 1. However, there was no evidence of rs7896741 and rs1259503 among the study and the experimental subjects.


Although it has been suggested that human VDAC2 plays an important role in sperm function and male fertility [5, 7, 9, 10], there have been no related studies of the hidden role that VDAC2 plays in male infertility, particularly idiopathic male infertility. This research was primarily aimed at analyzing 4 different tagSNPs (rs2804535, rs7896741, rs11001334 and rs1259503) of the VDAC2 gene to determine their relationship with male infertility. For the variant rs11001334, we determined that there was a negative relationship between this SNP and sperm parameters. In addition, there was a positive relationship between the same parameters and variant rs2804535 (P = 0.030 and 0.023). However, the men with normal semen parameters shared the same genotypes. Furthermore, among the rs7896741 and rs1259503 variants, there were no signs of a relationship with the risk of male infertility. These results illustrated that the different polymorphisms in the VDAC2 gene are not equal because they have extremely different relationships (positive or negative) with unique semen parameters in infertile males.

The SNP location helps predict disease etiology [11]. Despite the fact that the most common polymorphisms are situated in exonic sequences, an increasing number are located in introns [1216]. Recently, mutations deep in intronic sequences were reported to alter splicing efficiency by the means of disrupting or producing a splicing motif, affecting the binding of the splicing machinery, which could lead to a disorder in impacted individuals [12]. By analyzing the locations of the four tagSNPs (rs2804535, rs7896741, rs11001334 and rs1259503), we found that the SNP rs2804535 (T > C) was located within the fifth intron of the human VDAC2 gene. Thus, we hypothesized that the intronic SNP rs2804535 (T > C) may modulate the splicing efficiency of the VDAC2 gene. Specific motifs in the VDAC2 protein are required for mitochondrial Bak import [17, 18]. Therefore, this SNP may affect VDAC2 protein function by modulating the splicing efficiency of the VDAC2 gene and may ultimately result in idiopathic male infertility.

Interestingly, a study reported that one SNP in exon 5 of a non-coding RNA increased the transcriptional level of the novel gene and may play some role in making people sick [19]. In our study, the SNP rs11001334 (C > T) is a non-coding transcript exon variant located in exon 12 and may affect the transcriptional level of a non-coding RNA, which then may alter the transcription of the VDAC2 gene and may ultimately result in idiopathic male infertility. Furthermore, the rs11001334 (C > T) genotype had a strong relationship with sperm concentration and motility. Through our data, we confirmed that different polymorphisms in VDAC2 are responsible for the inability to conceive a baby along with abnormal semen parameters. The physiological impact of these SNPs on idiopathic male infertility requires further study.



This research was started with the approval of the Medical Ethics Committee of Nanjing Maternity and Child Health Care Hospital. From March 2014 to June 2015, 523 male patients who were proven to be infertile (while their wives were not) and 277 fertile men from the Nanjing Maternity and Child Health Care Hospital were continuously recruited to participate in the experiment. All the patients were of Han-Chinese ethnicity, had been informed of the aim of this study, and had been provided time to consider enrolling. Patients with idiopathic male infertility were asked to not to have sex or conceive a child for at least one year. Individuals who suffered from orchitis, congenital bilateral absence of the vas deferens, obstruction, cryptorchidism, Y chromosome microdeletions, and cytogenetic abnormalities were not ideal candidates for the study, as previously described [20]. The experimental subjects had already fathered at least one child without using any assistive techniques. All the patients provided a medical history and underwent physical examinations, which also required that they provide a volume of blood of 5 ml with the purpose of genomic DNA extraction as well as an ejaculate for performing semen analysis following a questionnaire.

Semen analysis

The analysis of semen was carried out in a manner similar to what has been described previously [21].

There are four indexes that were used to evaluate the association with VDAC2 gene polymorphisms, including sperm concentration, semen volume, sperm motility and sperm number per ejaculate.

VDAC2 SNP selection and genotyping

Using the HapMap genotype data that were collected from non-related Han Chinese individuals from the city of Beijing (HapMap Data Rel 24/phase II nov08, on NCBI B36 assembly, dbSNP b126;, SNPs within the 20641 bp human VDAC2 gene as well as the 1500 bp upstream and 1500 bp downstream that had a low frequency (>0.05%) in the Han Chinese population in Beijing and were located precisely on chromosome 10 from 76640569 to 76661210 were selected. Nineteen SNPs were identified in the region. Linkage disequilibrium (LD) in the zone was evaluated with the r2 values, which show the likelihood a certain SNP is capable of predicting the presence of other SNPs [22], using Haploview 4.0 software. By applying Tagger and a tagging threshold of r2 higher than 0.80, four tagSNPs (rs2804535, rs7896741, rs11001334 and rs1259503) that could predict the other SNPs from 76640569 to 76661210 were selected, with a mean r2 value of 1.0 (Figure (Figure11).

Figure 1
The overlook of the chosen tagging single nucleotide polymorphisms (tagSNPs) and the individual characteristics

Genomic DNA was isolated from peripheral blood leukocytes of the 523 Han Chinese men who suffered from the inability to conceive a baby in accordance with the prescribed protocols (Genomic DNA kit; Tiangen, Beijing, China). Using the Taq amplification method, TaqMan SNP Genotyping Assays were used for genotyping in an ABI Prism 7900 HT Fast Real-Time PCR system (Applied Biosystems, USA). The PCR amplification protocol included an initial step at 95°C for 10 minutes, then 40 cycles at the temperature of 95°C for approximately 15 seconds, 56°C for 10 seconds, and 60°C for 1 minute, and an additional step at 60°C for 10 minutes. Approximately 10 percent of the people who volunteered for the study were randomly chosen. In our study, the SNP genotyping coverage was perfect, as was the concordance of duplicates.

Statistical analysis

The statistical analysis was successfully completed in the Stata statistical package (Version 7.0, StataCorp, LP). The online software SHEsis is very effective in determining Hardy-Weinberg equilibrium (HWE) ( [23]. Logarithmic transformation was performed for the sperm concentration to confirm a normal distribution, and the outlying data were excluded. Cut-offs for the semen parameters were established according to the WHO reference values (World Health Organization, 1999 World Health Organization; Cambridge, UK: Cambridge University Press, 1999; WHO Laboratory Manual for the Examination of Human Semen and Sperm-Cervical Mucus Interaction, 4th ed.) for sperm concentration (< 20 × 106/ml), semen volume (< 2 ml), sperm motility (< 50% motile sperm), and sperm number per ejaculate (< 40 × 106/ml). The experimental subjects who previously conceived a child were split into 2 groups: case Group I, defined as male-factor infertility that met all four parameters, while case Group II failed to meet one or more of the semen parameters and was further broken into four subgroups: Subgroup I (sperm motility < 50% motile sperm), Subgroup II (semen volume < 2 ml), Subgroup III (sperm concentration < 20 × 106/ml) and Subgroup IV (sperm number per ejaculate < 40 × 106/ml). In the case of Group II, each of the subjects may have been classified into more than one subgroup. The method of comparing to the mean age, body mass index (BMI) and duration of sexual abstinence was used to analyze the difference of variance in the case of Group II. The Mann-Whitney U test was used to attempt to determine the differences in the subjects' daily life behaviors, such as smoking and the frequency of the VDAC2 allele in case groups. Multivariate logistic regression analysis was conducted to find the odds ratios (ORs) in the proven infertile men and the 95 percent confidence intervals (95% CI), adjusting for age, smoking and drinking habits, and body mass index (BMI) when appropriate. The χ2 test was used to precisely identify the difference in selected VDAC2 alleles between the HapMap study project and the experimental subjects. Differences were considered significant at P < 0.05.


This study was supported by grants from Medical Science and technology development Foundation Nanjing Department of Health (YKK14129,YKK15154), Nanjing science and technology project (201402024, 201503047), Natural Science Foundation of Jiangsu Province of China (BK20140083), and Science and Technology Development Fund of Nanjing Medical University (2014NJMUZD052).



The authors declare no conflicts of interest.

Contributed by

Authors' contributions

DZ and JHM conceived of and designed the experiments. DZ, JHM, DXQ, JYL, PX, JL and LJP performed the experiments. JHM, JYL and DXQ analyzed the data. LJP, JHM and DZ wrote and revised the manuscript.


1. Shoshan-Barmatz V, De Pinto V, Zweckstetter M, Raviv Z, Keinan N, Arbel N. VDAC, a multi-functional mitochondrial protein regulating cell life and death. Mol Aspects Med. 2010;31:227–85. [PubMed]
2. Hinsch KD, De Pinto V, Aires VA, Schneider X, Messina A, Hinsch E. Voltage-dependent anion-selective channels VDAC2 and VDAC3 are abundant proteins in bovine outer dense fibers, a cytoskeletal component of the sperm flagellum. J Biol Chem. 2004;279:15281–8. [PubMed]
3. Triphan X, Menzel VA, Petrunkina AM, Cassara MC, Wemheuer W, Hinsch KD, Hinsch E. Localisation and function of voltage-dependent anion channels (VDAC) in bovine spermatozoa. Pflugers Arch. 2008;455:677–86. [PubMed]
4. Liu B, Wang P, Wang Z, Zhang W. The use of anti-VDAC2 antibody for the combined assessment of human sperm acrosome integrity and ionophore A23187-induced acrosome reaction. PLoS One. 2011;6:e16985. [PMC free article] [PubMed]
5. Kwon WS, Park YJ, Mohamed el SA, Pang MG. Voltage-dependent anion channels are a key factor of male fertility. Fertil Steril. 2013;99:354–61. [PubMed]
6. Liu B, Wang P, Wang Z, Jia Y, Niu X, Wang W, Zhang W. Analysis and difference of voltage-dependent anion channel mRNA in ejaculated spermatozoa from normozoospermic fertile donors and infertile patients with idiopathic asthenozoospermia. J Assist Reprod Genet. 2010;27:719–24. [PMC free article] [PubMed]
7. Park YJ, Kwon WS, Oh SA, Pang MG. Fertility-related proteomic profiling bull spermatozoa separated by percoll. J Proteome Res. 2012;11:4162–8. [PubMed]
8. Blachly-Dyson E, Zambronicz EB, Yu WH, Adams V, McCabe ER, Adelman J, Colombini M, Forte M. Cloning and functional expression in yeast of two human isoforms of the outer mitochondrial membrane channel, the voltage-dependent anion channel. J Biol Chem. 1993;268:1835–41. [PubMed]
9. Hiyama G, Matsuzaki M, Mizushima S, Dohra H, Ikegami K, Yoshimura T, Shiba K, Inaba K, Sasanami T. Sperm activation by heat shock protein 70 supports the migration of sperm released from sperm storage tubules in Japanese quail (Coturnix japonica) Reproduction. 2014;147:167–78. [PubMed]
10. Zhang Y, Ou Y, Cheng M, Saadi HS, Thundathil JC, van der Hoorn FA. KLC3 is involved in sperm tail midpiece formation and sperm function. Dev Biol. 2012;366:101–10. [PubMed]
11. Lieben L. Disease genetics: SNP location helps predict disease aetiology. Nat Rev Genet. 2015 [PubMed]
12. Seo S, Takayama K, Uno K, Ohi K, Hashimoto R, Nishizawa D, Ikeda K, Ozaki N, Nabeshima T, Miyamoto Y, Nitta A. Functional analysis of deep intronic SNP rs13438494 in intron 24 of PCLO gene. PLoS One. 2013;8:e76960. [PMC free article] [PubMed]
13. He Q, Lu Y, Hu S, Huang Q, Li S, Huang Y, Hu Q, Wu L, Chen W. An intron SNP rs807185 in ATG4A decreases the risk of lung cancer in a southwest Chinese population. Eur J Cancer Prev. 2015 [PubMed]
14. Chaudhry AS, Prasad B, Shirasaka Y, Fohner A, Finkelstein D, Fan Y, Wang S, Wu G, Aklillu E, Sim SC, Thummel KE, Schuetz EG. The CYP2C19 Intron 2 Branch Point SNP is the Ancestral Polymorphism Contributing to the Poor Metabolizer Phenotype in Livers with CYP2C19*35 and CYP2C19*2 Alleles. Drug Metab Dispos. 2015;43:1226–35. [PubMed]
15. Abelleyro MM, Rossetti LC, de Los Angeles Curto M, Radic CP, Marchione VD, De Brasi CD. F8 intron 22 inversions and SNP rs73563631 in unrelated families with severe haemophilia A: clinical features and gene testing implications. Thromb Haemost. 2015:115. [PubMed]
16. Wang X, Zhong J, Gao Y, Ju Z, Huang J. A SNP in intron 8 of CD46 causes a novel transcript associated with mastitis in Holsteins. BMC Genomics. 2014;15:630. [PMC free article] [PubMed]
17. Ma SB, Nguyen TN, Tan I, Ninnis R, Iyer S, Stroud DA, Menard M, Kluck RM, Ryan MT, Dewson G. Bax targets mitochondria by distinct mechanisms before or during apoptotic cell death: a requirement for VDAC2 or Bak for efficient Bax apoptotic function. Cell Death Differ. 2014;21:1925–35. [PMC free article] [PubMed]
18. Naghdi S, Varnai P, Hajnoczky G. Motifs of VDAC2 required for mitochondrial Bak import and tBid-induced apoptosis. Proc Natl Acad Sci U S A. 2015;112:E5590–9. [PubMed]
19. Ishii N, Ozaki K, Sato H, Mizuno H, Saito S, Takahashi A, Miyamoto Y, Ikegawa S, Kamatani N, Hori M, Saito S, Nakamura Y, Tanaka T. Identification of a novel non-coding RNA, MIAT, that confers risk of myocardial infarction. J Hum Genet. 2006;51:1087–99. [PubMed]
20. Ding X, Zhang J, Bian Z, Xia Y, Lu C, Gu A, Li Y, Song L, Wang S, Wang X. Variants in the Eppin gene show association with semen quality in Han-Chinese population. Reprod Biomed Online. 2010;20:125–31. [PubMed]
21. Ding X, Zhang J, Fei J, Bian Z, Li Y, Xia Y, Lu C, Song L, Wang S, Wang X. Variants of the EPPIN gene affect the risk of idiopathic male infertility in the Han-Chinese population. Hum Reprod. 2010;25:1657–65. [PubMed]
22. Barrett JC, Fry B, Maller J, Daly MJ. Haploview: analysis and visualization of LD and haplotype maps. Bioinformatics. 2005;21:263–5. [PubMed]
23. Shi YY, He L. SHEsis, a powerful software platform for analyses of linkage disequilibrium, haplotype construction, and genetic association at polymorphism loci. Cell Res. 2005;15:97–8. [PubMed]

Articles from Oncotarget are provided here courtesy of Impact Journals, LLC