Past Issue

Volume 10, Number 3, Oct-Dec 2016, Pages: 303-309

Original Article(s)

Mitochondrial Genetic Variation in Iranian Infertile Men with Varicocele


Mohammad Mehdi Heidari, Ph.D., 1,*Mehri Khatami, Ph.D., 1Amirhossein Danafar, M.Sc., 2Tahere Dianat, M.Sc., 1Ghazaleh Farahmand, M.Sc., 3Ali Reza Talebi, Ph.D., 4
Department of Biology, Faculty of Science, Yazd University, Yazd, Iran
Department of Biology, Ashkezar Islamic Azad University, Ashkezar, Yazd, Iran
Department of Biology, Faculty of Science, Islamic Azad University Shahrekord, Shahrekord, Iran
Research and Clinical Center for Infertility and Department of Anatomy, Shahid Sadughi University of Medical Sciences, Yazd, Iran
*Corresponding Address: P.O.Box: 89195 -741 Department of Biology Faculty of Science Yazd University Yazd Iran Email: Heidarimm@yazd.ac.ir

Abstract

Background:

Several recent studies have shown that mitochondrial DNA mutations lead to major disabilities and premature death in carriers. More than 150 mutations in human mitochondrial DNA (mtDNA) genes have been associated with a wide spectrum of disorders. Varicocele, one of the causes of infertility in men wherein abnormal inflexion and distension of veins of the pampiniform plexus is observed within spermatic cord, can increase reactive oxygen species (ROS) production in semen and cause oxidative stress and sperm dysfunction in patients. Given that mitochondria are the source of ROS production in cells, the aim of this study was to scan nine mitochondrial genes (MT-COX2, MT-tRNALys , MT-ATP8, MT-ATP6, MT-COX3, MT-tRNAGly , MT-ND3, MT-tRNAArg and MT-ND4L) for mutations in infertile patients with varicocele.

Materials and Methods:

In this cross-sectional study, polymerase chain reaction-single strand conformation polymorphism (PCR-SSCP) and DNA sequencing were used to detect and identify point mutations respectively in 9 mitochondrial genes in 72 infertile men with varicocele and 159 fertile men. In brief, the samples showing altered electrophoretic patterns of DNA in the SSCP gel were sent for DNA sequencing to identify the exact nucleotide variation.

Results:

Ten type nucleotide variants were detected exclusively in mitochondrial DNA of infertile men. These include six novel nucleotide changes and four variants previously reported for other disorders.

Conclusion:

Mutations in mitochondrial genes may affect respiratory complexes in combination with environmental risk factors. Therefore these nucleotide variants probably lead to impaired ATP synthesis and mitochondrial function ultimately interfering with sperm motility and infertility.

Introduction

The major concern among married couples when they are unsuccessful to conceive after one year of regular unprotected intercourse is that they may be infertile. Male factors can be attributed to half of these cases (1,2). The most common surgically reversible cause of infertility is varicocele. Its prevalence is about 4.4-22.6% in the general population, 21-41% in men with primary infertility and 75-80% in men with secondary infertility (3,4). Despite the advances in molecular medicine, the pathophysiology of varicocele induced infertility remains unknown. Several proposed mechanisms include venous pressure changes and increased testicular temperature due to dilation and tortuosity of the pampiniform plexus of veins, oxidative stress, retrograde flow of renal or adrenal products, Leydig cell dysfunction and hyperthermia (5,6). In addition, a number of patients with varicocele have genetic abnormalities like Yq-microdeletions (7). Among them, oxidative stress-induced DNA damage appears to be a more likely cause which may severely affect sperm quality leading to infertility (8). This damage is one of the potential etiological factors in varicocele. A major source of partially reduced derivatives of molecular oxygen (O2) is mitochondria (9). The variety of reactive oxygen species (ROS) that mitochondria produce principally include hydrogen peroxide (H2O2), superoxide anion (O2.-) and the hydroxyl radical (OH.) (10,11). In normal physiology, ROS perform several roles in regulating cellular functions by interacting with cellular components (12). In fertile men, physiological levels of ROS play important roles in sperm function, acrosome reaction, capacitation, hyper-activation and the penetration of oocyte by spermatozoa. However, in varicocele patients ROS generation is abnormally enhanced (13,14).

Specific point mutations and deletions of mitochondrial DNA (mtDNA) have been associated with poor sperm motility and semen quality in several studies. Sperm mtDNA is highly sensitive to mutations due to increased ROS by-products generated during oxidative respiration (15). When large amounts of mutant mtDNA accumulate in the testes, reduction in ATP production, mitochondrial respiratory dysfunction and meiotic arrest are induced in spermatogonia cells (16). Each mitochondrion has 2-10 mitochondrial genomes responsible for coding the subunits of the OXPHOS complex. The OXPHOS machinery is made up of over 80 different polypeptides, of which the mtDNA encodes 13 polypeptides including complex I, III, IV and V subunits (17).

In this study, for first time, we further analyzed nine genes (MT-COX2, MT-tRNALys, MT-ATP8, MT-ATP6, MT-COX3, MT-tRNAGly, MT-ND3, MT-tRNAArg and MT-ND4L) in the mitochondrial genome by polymerase chain reaction-single strand conformation polymorphism (PCR-SSCP) assay and direct sequencing techniques to identify the possible association between mtDNA variation with varicocele in the Iranian population.

Materials and Methods

Patients

This study was a cross sectional study. Seventy two Iranian infertile men with clinical varicocele were recruited in the study. The varicocele diagnosis was made by the urologists for the patients by physical examination in standing position and via scrotal palpation in a temperature controlled room (23°C). Semen analysis was performed according to the WHO laboratory manual (18). Patients with varicocele were in 3 grades: i. Grade I (n=12), ii. Grade II (n=27) and iii. Grade III (n=33). The control group (healthy volunteers) consisted of 159 fertile and normospermic men from the Yazd Infertility Center who fathered at least one child. The ethnic and geographical origin of all patients and controls was the same. All participants were fully informed of the objectives of the study and those that signed the consent form were assigned to the study. All infertile men in the age group ranging from 22 to 36 years (mean, 29 years) were referred for evaluation of their infertility (1 year of unprotected intercourse and not leading to conception). The Yazd University Ethics Committee approved recruitment of patients and laboratory protocols in this study.

DNA extraction and mutation analysis

Peripheral blood samples were obtained from varicocele patients and the DNA was extracted using a standard salting-out procedure. Purified DNA samples from leukocytes were used for the PCR reactions. To amplify MT-COX2, MT-tRNALys , MT-ATP8, MT-ATP6, MT-COX3, MT-tRNAGly , MT-ND3, MT-tRNAArg and MT-ND4L mitochondrial genes, four pairs of PCR primers were designed, which were located in the flanking regions of each gene (Table 1,). Primer Design was based on the human mitochondrial sequence by primer design software (Primer Premier 5.0; Premier Biosoft Inc., Canada), and their secondary structure was examined with Gene Runner version 3.05 (Hastings Software Inc. Hastings, NY, USA, http://www.generunner.com). Each reaction was prepared to a final volume of 25 µl containing 1XMasterMix PCR (Yekta Tajhiz Azma Co., Iran), 0.2 mM of each primer and 0.5 µg DNA template. The PCR conditions were an initial denaturation of 95°C for 5 minutes followed by 35 cycles of denaturation at 95°C for 30 seconds, the annealing temperature (Table 1,) for 30 seconds and extension at 72°C for 30 seconds, which was extended for 5 minutes in the final cycle. The PCR products were electrophoresed on an ethidium bromide-stained 2% agarose gel.

Table 1- Primers used for mitochondrial genes

For the SSCP assay, PCR products were heatdenatured at 95 ºC for 5 minutes and chilled on ice for 5 minutes, and then loaded onto an 8% nondenaturing polyacrylamide/TBE 0.5x gel. Gels were stained with silver nitrate to reveal the bands of single strand DNA. Various band patterns of the amplified PCR products were marked and scored. The typical gene variants got sequenced using a commercial company (Macrogen, South Korea). All the data obtained from automated sequencing was checked with Sequencher. The online multiple sequence alignment software ClustalW2 (http:// www.ebi.ac.uk/tools/msa/clustalw2/) and BLAST analysis were used to determine the nature of mutations and percent homology of the sequences that have been obtained in the study with all other sequences of five other species (chimpanzee, monkey, cattle, zebrafish and drosophila).

Software and databases

We used the tool PolyPhen-2 (http://genetics.bwh.harvard.edu/pph2/) for prediction of the functional consequences of mutations and damaging effect of missense mutations on protein structure. The sequence alignment was performed using the blastp program available at the National Center for Biotechnology Information (NCBI) web site (http://www.ncbi.nlm.nih.gov/Blastp) and the ClustalW program (http://bioinfo.hku.hk/services/ analyseq/cgi-bin/clustalw_in.pl). For detection of structural features of mammalian mitochondrial tRNAs and human diseases linked to point mutations in mitochondrial tRNA genes, we used Mamit-tRNA (http://mamit-trna.u-strasbg.fr).

Statistical analysis

The GraphPad Prism software (GraphPad Software, Inc. USA) was used for statistical analysis. Distributions of continuous variables in groups were expressed as mean ± SD, and compared with unpaired Student’s t test. P<0.05 were regarded as statistically significant.

Results

The age difference between the 72 Iranian infertile men with varicocele (mean age of 30.76 ± 6.47) and 159 normal controls (mean age: 28.8 ± 6.01) was not significant (P=0.785). Mutation analyses for the mitochondrial MT-COX2, MTtRNALys, MT-ATP8, MT-ATP6, MT-COX3, MTtRNAGly, MT-ND3, MT-tRNAArg and MT-ND4L genes were carried out in all of patients and healthy controls by PCR-SSCP. Mobility shift of single strand DNA on polyacrylamide gel electrophoresis was the criterion for sequencing and the identification of DNA variation (Fig .1,). We found ten different nucleotide substitutions of which 4 caused an amino acid change, of which one occurred in tRNAArg. None of the ten mutations were found in healthy controls. All the mutations identified are summarized in Table 2. In addition, 6 were novel mutations of which four were silent mutations. Four reported polymorphisms, including m.8258T>C, m.9911C>A, m.9932G>A and m.10463T>C were found in six patients. The m.9911C>A variant in MT-COX3 was heteroplasmic. The novel 9 bp heteroplasmic insertion was found in the non-coding MT-NC7 locus in one patient.

- Silver staining SSCP analysis of fragment 3. A. Polyacrylamide gel electrophoresis. Lanes 1, 3 and 5 show 3 patients who did not have mutations, B. Lane 6 shows a patient with the m.9929 C>A mutation, C. Lanes 2, 4 and 8 show 2 patients with the m.10141A>C mutation, and, D. Lane 7 shows a patient with m.9911C>A and lanes 9, 10 and 11 are men without varicocele.

The m.9911C>A mutation, an aromatic amino acid phenylalanine codon (TTC), changes to leucine codon (TTA), a hydrophobic amino acid at position 235 (designated F235L) in 1 patient (Fig .2,) and the novel 9929C>A mutation changes a polar tyrosine to threonine.

Two novel mutations were detected in 5 patients with one (9929C>A) being a nonsense mutation and changes tyrosine to stop codon (Y241X) and the other (10141C>G) being a missense mutation that changes Asparagine to Lysine (N27K). Also, three synonymous polymorphisms were found that were not reported previously (Table 2,). The other variation was the m.10463T>C substitution (homoplasmic state) in the tRNAArggene that was found in 3 patients.

- Protein alignment of m.9911C>A missense mutation MT-COX3 and the arrow indicate the site of the F235L mutation.
Table 2-Mitochondrial variation found in infertile men with varicocele

Discussion

One of the most frequent causes of male infertility is varicocele, however, the pathogenic mechanisms by which it leads to changes in spermatogenesis are not clear (25). Some of these mechanisms may be related to mutations in mitochondrial complexes that affect flagellar movement and cause sperm dysmotility.

DNA alterations including point mutations and deletions of mtDNA have been reported in infertile patients with low sperm motility (26). The effect of mtDNA mutations on male infertility has also been studied. Shamsi et al. (27) reported that generation of ROS and mtDNA mutations are associated with pathogenic molecular mechanisms. Agarwal et al. (28) showed an increased oxidative stress in varicocele patients. Thangaraj et al. (29) demonstrated that sperm mitochondrial mutations is one of the causes of low sperm motility which is strongly dependent on ATP biosynthesis which is carried out by the mitochondrial OXPHOS system. Furthermore, it has been demonstrated that cells with some base substitutions in mtDNA can greatly influence semen quality (9,30,31).

It has been established that mitochondrial dysfunction caused by mtDNA mutations and oxidative damages is one of the important reasons for most types of infertility such as Varicocele (32). The mtDNAs alterations may accumulate in the spermatids or during gametogenesis and thereby impair the respiratory function and motility of spermatozoa (33).

We observed three heteroplasmic variations in 4 patients. A nine base pair heteroplasmic insertion in the non-coding MT-NC7 locus were found in 1 patient. Although this insertion (5ˊ-CCCCCTCTA-3ˊ) has been found in a noncoding region, it may cause mitochondrial rearrangements and DNA strand break affected by topoisomerases or DNA recombinase (34).

The heteroplasmic m.9911C>A and m.9929C>A transversions in MT-COX3 alter two conserved codons. Given that these variants change highly conserved amino acids and were not identified in normal controls, they may be considered as pathogenic mutations for the following reasons. First, these missense mutations are found in several patients. Second, these mutations are not reported as polymorphisms in the general population and are not detected in the control individuals from the same ethnic background. Third, the mutations are heteroplasmic in the lymphocyte cells. Fourth, we propose that these mutations may affect the polarity of the protein due to the replacement of a natural amino acid with a polar amino acid. Using PolyPhen2, we found that these mutations are expected to change protein function.

Here, we describe seven homoplasmic variants in 25 patients: two missense mutations in MT-COX2 (8258T>C) and MT-ND3 (10141C>A), 4 synonymous polymorphisms in MT-COX3, MTATPase6, MT-ND3 and MT-ND4L and one mutation (10463T>C) in MT-tRNAArg. This mutation is located at a moderate conserved region of the acceptor stem of tRNA arginine. This mutation was not observed in healthy control subjects but was previously reported as a polymorphism in mitochondrial encephalomyopathy (35) and may be one of the several predisposing factors for varicocele.

Conclusion

Because sperms require an optimal energy to reach the oviduct during fertilization, the appropriate bioenergetic function of mitochondria is critical for male infertility. Therefore, any changes in mitochondrial genome can cause improper functioning of respiratory chain that in combination with environmental risk factors lead to infertility in men. This first Iranian study revealed that some Iranian infertile men carry variants in the nine mitochondrial genes and suggests that variants in these genes may be associated with varicocele.

Acknowledgements

This work was financially supported by Yazd University. We thank all the patients for providing blood samples for this study. The authors hereby declare that there is no conflict of interests.

References

Raman JD, Walmsley K, Goldstein M. Inheritance of varicoceles. Urology. 2005; 65(6): 1186-1189.
Evans CJ, Fowkes FG, Hajivassiliou CA, Harper DR, Ruckley CV. Epidemiology of varicose veins.A review. Int Angiol. 1994; 13(3): 263-720.
Canales BK, Zapzalka DM, Ercole CJ, Carey P, Haus E, Aeppli D. Prevalence and effect of varicoceles in an elderly population. Urology. 2005; 66(3): 627-631.
Niederberger C. Inheritance of varicoceles. J Urol. 2005; 174(5): 1942-1942.
Benoff S, Gilbert BR. Varicocele and male infertility: part I.Preface. Hum Reprod Update. 2001; 7(1): 47-54.
Talebi AR, Moein MR, Tabibnejad N, Ghasemzadeh J. Effect of varicocele on chromatin condensation and DNA integrity of ejaculated spermatozoa using cytochemical tests. Andrologia. 2008; 40(4): 245-251.
Güney AI, Javadova D, Kirac D, Ulucan K, Koc G, Ergec D. Detection of Y chromosome microdeletions and mitochondrial DNA mutations in male infertility patients. Genet Mol Res. 2012; 11(2): 1039-1048.
Lü YQ, Chen B. Progress in researches on the mechanism of varicocele-induced male infertility. Zhonghua Nan Ke Xue. 2008; 14(5): 454-458.
Forkink M, Smeitink JA, Brock R, Willems PH, Koopman WJ. Detection and manipulation of mitochondrial reactive oxygen species in mammalian cells. Biochim Biophys Acta. 2010; 1797(6-7): 1034-1044.
Starkov AA. The role of mitochondria in reactive oxygen species metabolism and signaling. Ann N Y Acad Sci. 2008; 1147: 37-52.
Romeo C, Santoro G. Free radicals in adolescent varicocele testis. Oxid Med Cell Longev. 2014; 2014: 912878-912878.
Shamsi MB, Venkatesh S, Tanwar M, Talwar P, Sharma RK, Dhawan A. DNA integrity and semen quality in men with low seminal antioxidant levels. Mutat Res. 2009; 665(1-2): 29-36.
Moustafa MH, Sharma RK, Thornton J, Mascha E, Abdel-Hafez MA, Thomas AJ Jr. Relationship between ROS production, apoptosis and DNA denaturation in spermatozoa from patients examined for infertility. Hum Reprod. 2004; 19(1): 129-138.
Benoff SH, Millan C, Hurley IR, Napolitano B, Marmar JL. Bilateral increased apoptosis and bilateral accumulation of cadmium in infertile men with left varicocele. Hum Reprod. 2004; 19(3): 616-627.
Agarwal A, Makker K, Sharma R. Clinical relevance of oxidative stress in male factor infertility: an update. Am J Reprod Immunol. 2008; 59(1): 2-11.
Baklouti-Gargouri S, Ghorbel M, Ben Mahmoud A, Mkaouar-Rebai E, Cherif M, Chakroun N. Mitochondrial DNA mutations and polymorphisms in asthenospermic infertile men. Mol Biol Rep. 2013; 40(8): 4705-4712.
Tuppen HA, Blakely EL, Turnbull DM, Taylor RW. Mitochondrial DNA mutations and human disease. Biochim Biophys Acta. 2010; 1797(2): 113-128.
. Laboratory manual of the WHO for the examination of human semen and spermcervical mucus interaction. Ann Ist Super Sanita. 2001; 37(1): I-123.
Achilli A, Rengo C, Magri C, Battaglia V, Olivieri A, Scozzari R. The molecular dissection of mtDNA haplogroup H confirms that the Franco-Cantabrian glacial refuge was a major source for the European gene pool. Am J Hum Genet. 2004; 75(5): 910-918.
Tang S, Wang J, Zhang VW, Li FY, Landsverk M, Cui H. Transition to next generation analysis of the whole mitochondrial genome: a summary of molecular defects. Hum Mutat. 2013; 34(6): 882-893.
Ozawa T, Tanaka M, Ino H, Ohno K, Sano T, Wada Y, et al. Distinct clustering of point mutations in mitochondrial DNA among patients with mitochondrial encephalomyopathies and with Parkinson's disease. Biochem Biophys Res Commun. 1991; 176(2): 938-946.
Yeh JJ, Lunetta KL, van Orsouw NJ, Moore FD Jr, Mutter GL, Vijg J. Somatic mitochondrial DNA (mtDNA) mutations in papillary thyroid carcinomas and differential mtDNA sequence variants in cases with thyroid tumours. Oncogene. 2000; 19(16): 2060-2066.
Houshmand M, Larsson NG, Holme E, Oldfors A, Tulinius MH, Andersen O. Automatic sequencing of mitochondrial tRNA genes in patients with mitochondrial encephalomyopathy. Biochim Biophys Acta. 1994; 1226(1): 49-55.
Howell N, Kubacka I, Halvorson S, Howell B, McCullough DA, Mackey D. Phylogenetic analysis of the mitochondrial genomes from Leber hereditary optic neuropathy pedigrees. Genetics. 1995; 140(1): 285-302.
Dada R, Shamsi MB, Venkatesh S, Gupta NP, Kumar R. Attenuation of oxidative stress & DNA damage in varicocelectomy: implications in infertility management. Indian J Med Res. 2010; 132: 728-730.
Benoff S, Marmar JL, Hurley IR. Molecular and other predictors for infertility in patients with varicoceles. Front Biosci (Landmark Ed). 2009; 14: 3641-3672.
Shamsi MB, Kumar R, Bhatt A, Bamezai RN, Kumar R, Gupta NP. Mitochondrial DNA mutations in etiopathogenesis of male infertility. Indian J Urol. 2008; 24(2): 150-154.
Agarwal A, Sharma RK, Nallella KP, Thomas AJ Jr, Alvarez JG, Sikka SC. Reactive oxygen species as an independent marker of male factor infertility. Fertil Steril. 2006; 86(4): 878-885.
Thangaraj K, Joshi MB, Reddy AG, Rasalkar AA, Singh L. Sperm mitochondrial mutations as a cause of low sperm motility. J Androl. 2003; 24(3): 388-392.
Agarwal A, Said TM. Role of sperm chromatin abnormalities and DNA damage in male infertility. Hum Reprod Update. 2003; 9(4): 331-345.
Schoor RA, Elhanbly SM, Niederberger C. The pathophysiology of varicocele-associated male infertility. Curr Urol Rep. 2001; 2(6): 432-436.
Palanichamy MG, Zhang YP. Identifying potential pitfalls in interpreting mitochondrial DNA mutations of male infertility cases. Indian J Med Res. 2011; 134: 447-451.
Venkatesh S, Deecaraman M, Kumar R, Shamsi MB, Dada R. Role of reactive oxygen species in the pathogenesis of mitochondrial DNA (mtDNA) mutations in male infertility. Indian J Med Res. 2009; 129(2): 127-137.
Fritsch ES, Chabbert CD, Klaus B, Steinmetz LM. A genome-wide map of mitochondrial DNA recombination in yeast. Genetics. 2014; 198(2): 755-771.
Uusimaa J, Finnilä S, Remes AM, Rantala H, Vainionpää L, Hassinen IE. Molecular epidemiology of childhood mitochondrial encephalomyopathies in a Finnish population: sequence analysis of entire mtDNA of 17 children reveals heteroplasmic mutations in tRNAArg, tRNAGlu, and tRNALeu(UUR) genes. Pediatrics. 2004; 114(2): 443-450.