A Systematic Review of the Impact of Mitochondrial Variations on Male Infertility

Abstract

According to current estimates, infertility affects one in four couples trying to conceive. Primary or secondary infertility can be due either to both partners or only to the man or the woman. Up to 15% of infertility cases in men can be attributed to genetic factors that can lead to irreversible partial or complete spermatogenic arrest. The increased use of assisted reproductive technology (ART) has provided not only insights into the causes of male infertility but also afforded a diagnostic tool to detect and manage this condition among couples. Genes control a variety of physiological attributes, such as the hypothalamic–pituitary–gonadal axis, development, and germ cell differentiation. In the era of ART, it is important to understand the genetic basis of infertility so as to provide the most tailored therapy and counseling to couples. Genetic factors involved in male infertility can be chromosome abnormalities or single-gene disorders, mitochondrial DNA (mtDNA) mutations, Y-chromosome deletions, multifactorial disorders, imprinting disorders, or endocrine disorders of genetic origin. In this review, we discuss the role of mitochondria and the mitochondrial genome as an indicator of sperm quality and fertility.

1. Mitochondrial DNA (mtDNA)

The mitochondrion was first defined by Von Kölliker in 1856, during his study of muscle tissue. About two billion years ago, mitochondria originated when a precursor to the modern eukaryotic cell engulfed an α-proteobacterium through an opportunistic relationship [1]. Since then, mitochondria have kept the double-membrane structure of their ancestors, but their overall shape and composition have changed consistently. Mitochondria have retained their genomes (mtDNA), reflecting their evolutionary origin as bacteria. While nuclear genes encode most mitochondrial proteins, some respiratory proteins and mitochondrial tRNAs are still encoded by the mitochondrial genome [2].

Another feature of mitochondria is the presence of specific ribosomes that allow localized protein synthesis [3].

mtDNA is present in stroma and is only transmitted through the female germline [4]. Sperm have only 100 mtDNA copies compared to eggs, which have 150,000 [5]. Mitochondria are located in the midsection of spermatozoa and have one copy of mtDNA per mitochondria [6,7,8].

Farge et al. (2019) reported that the gene structure and arrangement of mtDNA are highly conserved across different mammal species, including dolphins, bears, dogs, foxes, horses, llamas, mice, pandas, and many others. This includes approximately 16.6 kbs of close circular double-stranded DNA (heavy and light) strands [9].

Cristae, the main sites of mitochondrial energy conversion, are far into the matrix. ATP synthase in cristae membranes operates under a low pH gradient between the intermembrane space (pH 7.2–7.4) and the matrix (pH 7.9–8) [10].

There are four electron transport chain (ETC) complexes inside the inner membrane. Each ETC complex, except II, has genes encoded by the mitochondrial genome, while the rest of the ETC subunits are encoded by the nucleus. The mitochondrial respiratory chain consists of 13 mtDNA-encoded proteins. Complex I includes seven nicotinamide adenine dinucleotide hydride (NADH) dehydrogenase subunits: NADH dehydrogenase 1 (ND1), NADH dehydrogenase 2 (ND2), NADH dehydrogenase 3 (ND3), NADH dehydrogenase 4 (ND4), NADH dehydrogenase 4L (ND4L), NADH dehydrogenase 5 (ND5), and NADH dehydrogenase 6 (ND6); complex III contains cytochrome B (CYB); complex IV contains three subunits: cytochrome oxidase subunit I (COX I), cytochrome oxidase subunit II (COX II), and cytochrome oxidase subunit III (COX III), while complex V contains ATPase 6 and ATPase 8 [11].

Complex I is the main point of entry of electrons into the respiratory chain and is suggested as the rate-limiting step in general respiration. Therefore, it plays a critical role in energy metabolism [12]. For that reason, mitochondria change in location and number according to cell type. When they are present in large numbers in cells, this means that these cells require a lot of energy; for example, in oocytes, mitochondria can number up to 100,000 [10].

mtDNA is a naked molecule that lacks introns and histones, and both of its strands are transcribed to synthesize functional proteins. Due to asexual replication, very basic repair mechanisms, a lack of protective histones, and close proximity to free radical formation sites, mtDNA has a 10- to 20-fold higher mutation rate than nuclear DNA. Furthermore, mtDNA replicates rapidly in the absence of DNA repair machinery [7]. Thus, mtDNA is 100-fold more prone to mutation than nuclear DNA [13].

Scientists believe that human mtDNA is inherited only from the mother and that paternal mtDNA disappears after the cleavage stage [14].

In October 2018, Luo and his team provided evidence that paternal mtDNA can be passed on to offspring [15]. A recent study by Annis et al. (2019) dismissed this model of inheritance, which hinged on the notion of biparental mtDNA transmission to the offspring, as inaccurate [16].

2. Paternal Inheritance of Mitochondrial Genome (mtDNA)

Sperm-derived paternal mitochondria and their mtDNA reach the cytoplasm of an oocyte upon fertilization, typically disappear during early embryogenesis, and are never passed on to offspring.

However, the molecular mechanisms involved in this paternal mitochondrial clearance remain relatively unclear [17,18].

For decades, it has been agreed that human mtDNA is inherited entirely from the maternal line. The co-occurrence of mutant and wild-type variant alleles in the same individual (heterogeneity) and rapid changes in allele frequencies result in different disease severities in offspring [19].

Contrary to the long-held view that mtDNA is strictly inherited from the mother, a recent study challenges this view and provides evidence of additional paternal mtDNA transmission from father to offspring [15]. Paternal mtDNA is determined by quasi-Mendelian inheritance [16]. Furthermore, during intracytoplasmic sperm injection (ICSI), all sperm are injected into the cytoplasm of the oocyte and, interestingly, mtDNA is preserved, allowing offspring to share their father’s mtDNA [20].

A recent study provides evolutionary evidence in support of a role for paternal mtDNA in fertilization, with data collected from two neotropical primate families (Cebidae and Atelidae) suggesting that the midpiece-containing mtDNA has evolved to become larger and wider in the younger species (Atelidae) compared to the narrower and shorter midpiece in the ancestral species (Cebidae) [21].

A previous phenomenon supporting the role of sperm mtDNA in early embryonic development is the bi-uniparental inheritance of mitochondria in sea mussels, as females contain predominantly maternal mtDNA in their somatic cells, while males contain maternal mtDNA in their somatic cells and paternal mtDNA in gonads [22].

In a study performed by Cogswell et al. (2006), paternal mtDNA is seen to behave differently in eggs that will develop into males or females. In male-producing eggs, the paternal mtDNA tend to congregate together within the same cell, while in female-producing oocytes, the paternal mtDNA are scattered throughout the egg. During pre-cleavage stages of development, the paternal mitochondrion continues to have sex-specific roles and is seen to be localized in different ways, which supports its role in early development [23].

In a study on Caenorhabditis elegans (C. elegans), researchers found that paternal mtDNA was completely destroyed after being fertilized by autophagosomes [24]. Another study showed that the entire sperm mtDNA was destroyed in the eggs of pigs and monkeys following fertilization by the ubiquitin–proteasome system, thanks to a specific microtubule-associated protein called SQSTM1 [25].

3. Mitochondrial Genome Mutations/Variations in Humans

Recently, many mtDNA mutations that are associated with the development and progression of human disease have been identified [26]. As many as 2 million Americans are affected by mitochondrial diseases (Lemonick, 2006). According to the United Mitochondrial Disease Foundation (UMDF), “every 30 min a child is born with mitochondrial disease which develops at age 10 [sic]”. Mitochondrial disease occurs in about 1 in 4000 people in the United States. The Center for Mitochondrial and Metabolic Diseases at UC San Diego estimates that 1000–4000 newborns in the United States are born with mitochondrial disease each year [27].

Numerous epidemiological studies supporting this idea have also been conducted in Europe [19,28,29,30].

Genetic variants in mitochondrial genes are associated with many diseases (Table 1). Most of these diseases affect organs with high energy demands, such as the brain, skeletal muscle, eyes, and heart [31]. Because mtDNA is not protected by histones or other DNA-binding proteins, it is more susceptible to DNA damage caused by excess reactive oxygen species (ROS) and free radicals in the matrix [32]. Furthermore, the mtDNA repair machinery is less efficient compared to that used in nuclear DNA repair [33]. Together, these elements increase the mutation rate in mtDNA by a factor of 10 to 100 compared to nuclear DNA [13].

A previous study found that the incidence of the 4216 T>C variant was higher in diabetic patients than in controls and was statistically associated with type 2 diabetes [34]. Another study found an interesting male-specific association between the 4216 T>C variant and infection rates, leading to complex sepsis and death [35].

Likewise, the 13708 G>A variant is associated with multiple clinical manifestations. It has been shown to increase susceptibility to multiple sclerosis [36], and, in another study, it was found to increase hereditary Leber‘s amplified expression in optic neuropathy (LHON) disease [37]. Males who carry the 13708 G>A variant have increased chances of developing Alzheimer’s disease. In fact, it was found that this variant is more dangerous in male patients than in female patients [38].

Recent studies have reported significant associations between MTND3 polymorphisms and the risk of Parkinson’s disease, type 2 diabetes mellitus, and breast and esophageal cancer, but not gastric cancer [39,40,41]. Furthermore, among the identified MTND3 single-nucleotide polymorphisms (SNPs), rs2853826 (A10398G) (MT-ND3) has been reported to be associated with the increased production of ROS in mitochondria, leading to oxidative stress and mtDNA damage [41].

Some mutations in the MT-ND4L gene have been found to be connected to specific disease disorders, such as LHON [42]. In the MT-ND4 gene, variation is related to macular degeneration in the eye (AMD), mesial temporal lobe epilepsy (MTLE), cystic fibrosis, and even aging [43,44,45,46].

A significant association between rs28358280 (A10550G) (MT-ND4L) and body mass index (BMI) was identified for the first time with an increase in the G allele leading to higher BMI compared to when the allele alone was present [47]. Furthermore, rs2857285 is associated with aggressive ovarian cancer [48].

Four versions of the ND6 gene were found to be correlated with Leigh disease in patients from Italy, France, and Germany: 14459G>A, 14495A>G, 14482C>A, and 14568C>T [42]. The variant 14439G>A was associated with a child who had mitochondrial respiratory chain disorder [49]. Another variant, 14459G>A, was found to be related to Leigh disease [50].

Table 1. Summary of studies reporting different mitochondrial disorders.

Mitochondrial Genome Abnormalities		Disease	References
Complex 1 NADH dehydrogenase	4216 T>C in MT-ND 1 (missense variant)	Diabetes mellitus type 2 (T2D)	[51,52]
		Leber’s Hereditary Optic Neuropathy (LHON)	[53]
		Male-specific infection, leading to complicated sepsis and death	[35]
	5178 C>A in MT-ND2	T2D	[51]
	rs2853826 and rs414676521 in MT-ND3	Earlier age at onset in males, Machado–Joseph disease, breast cancer, T2D, Parkinson’s disease, esophageal cancer, gastric cancer, LHON	[39]
	120271 T>C and 12096 T>A in MT-ND4	Schizophrenia (SCZ), age-related muscular degeneration (AMD), mesial temporal lobe epilepsy (MTLE), cystic fibrosis	[43,44,45,46,54]
	rs28358280 in MT-ND4L	Body mass index	[47]
	rs2853495 in MT-ND4L	Ulcerative colitis and pancreatic cancer	[55,56]
	rs869096886 in MT-ND4L	SCZ	[57]
	rs2857285 in MT-ND4L	Ovarian cancer	[48]
	11777C>A in MT-ND4L	Late-onset encephalopathy	[58]
	13708 G>A in MT-ND 5 gene	SCZ, increase in the susceptibility to multiple sclerosis, enhanced expression of LHON, increase in the risk of Alzheimer’s disease specifically in the male patients, breast cancer	[36,38,53,59,60,61,62]
	14439G>A in MT-ND 6 gene (missense variant)	Mitochondrial respiratory chain disease	[49]
	14459 G>A in MT-ND 6 gene	Leigh syndrome	[63]
	14459G>A, 14495A>G; 14482C>A and 14568C>T in MT- ND 6	LHON disease among patients from Germany, France, and Italy	[42]
Complex III MT-CYB	rs2853506 (15218A>G)	Epileptogenesis	[64]
	rs2853508	Breast cancer	[65,66]
	rs41518645	LHON	[67]
Complex IV	Mutations in MT-CO III genes	Recurrent myoglobinuria, LHON, severe encephalopathy, isolated myopathy	[68]
the rRNA: RNR 1 (12 S), RNR2 (16 S)	1709G>A, 15851A>G	Parkinson’s disease	[60,69]
tRNAs	Variant at position, 15928	Alzheimer’ disease	[70,71]
	8344 A>G in tRNA Lys gene	Myoclonus epilepsy and ragged–red fiber (MERRE) diseases	[72]
	8363 G>A in tRNA Lys gene	Correlated with autism spectrum disorders (ASD)	[73]
	8326 A>G in tRNA Lys, 15995 G>A in tRNA pro	Cystic fibrosis	[74]
ATP Synthetase 6 gene (ATPase 6)	9176 T>C	Mild myopathic change	[75]
	8839G>C	Retinitis pigmentosa syndrome (NARP)	[76]
	8914C>T	Mitochondrial encephalomyopathies	[77]
	8593 A>G	Leigh syndrome with a deficiency in mitochondrial energy production	[78]
ATP Synthetase 8 gene (ATP 8)	Mutations in MT-ATP8	LHON, MELAS, Leigh syndrome, NARP	[76,79]
COII gene	7750 C>A	SCZ	[80]

MT-: mitochondrial, NADH: nicotinamide adenine dinucleotide hydride dehydrogenase, ND1: NADH dehydrogenase 1, ND2: NADH dehydrogenase 2, ND3: NADH dehydrogenase 3, ND4: NADH dehydrogenase 4, ND4L: NADH dehydrogenase 4L, ND5: NADH dehydrogenase 5, and ND6: NADH dehydrogenase 6; CYB: Cytochrome b, COX II: cytochrome oxidase subunit II, and COX III: cytochrome oxidase subunit III.

Table 1. Summary of studies reporting different mitochondrial disorders.

Mitochondrial Genome Abnormalities		Disease	References
Complex 1 NADH dehydrogenase	4216 T>C in MT-ND 1 (missense variant)	Diabetes mellitus type 2 (T2D)	[51,52]
		Leber’s Hereditary Optic Neuropathy (LHON)	[53]
		Male-specific infection, leading to complicated sepsis and death	[35]
	5178 C>A in MT-ND2	T2D	[51]
	rs2853826 and rs414676521 in MT-ND3	Earlier age at onset in males, Machado–Joseph disease, breast cancer, T2D, Parkinson’s disease, esophageal cancer, gastric cancer, LHON	[39]
	120271 T>C and 12096 T>A in MT-ND4	Schizophrenia (SCZ), age-related muscular degeneration (AMD), mesial temporal lobe epilepsy (MTLE), cystic fibrosis	[43,44,45,46,54]
	rs28358280 in MT-ND4L	Body mass index	[47]
	rs2853495 in MT-ND4L	Ulcerative colitis and pancreatic cancer	[55,56]
	rs869096886 in MT-ND4L	SCZ	[57]
	rs2857285 in MT-ND4L	Ovarian cancer	[48]
	11777C>A in MT-ND4L	Late-onset encephalopathy	[58]
	13708 G>A in MT-ND 5 gene	SCZ, increase in the susceptibility to multiple sclerosis, enhanced expression of LHON, increase in the risk of Alzheimer’s disease specifically in the male patients, breast cancer	[36,38,53,59,60,61,62]
	14439G>A in MT-ND 6 gene (missense variant)	Mitochondrial respiratory chain disease	[49]
	14459 G>A in MT-ND 6 gene	Leigh syndrome	[63]
	14459G>A, 14495A>G; 14482C>A and 14568C>T in MT- ND 6	LHON disease among patients from Germany, France, and Italy	[42]
Complex III MT-CYB	rs2853506 (15218A>G)	Epileptogenesis	[64]
	rs2853508	Breast cancer	[65,66]
	rs41518645	LHON	[67]
Complex IV	Mutations in MT-CO III genes	Recurrent myoglobinuria, LHON, severe encephalopathy, isolated myopathy	[68]
the rRNA: RNR 1 (12 S), RNR2 (16 S)	1709G>A, 15851A>G	Parkinson’s disease	[60,69]
tRNAs	Variant at position, 15928	Alzheimer’ disease	[70,71]
	8344 A>G in tRNA Lys gene	Myoclonus epilepsy and ragged–red fiber (MERRE) diseases	[72]
	8363 G>A in tRNA Lys gene	Correlated with autism spectrum disorders (ASD)	[73]
	8326 A>G in tRNA Lys, 15995 G>A in tRNA pro	Cystic fibrosis	[74]
ATP Synthetase 6 gene (ATPase 6)	9176 T>C	Mild myopathic change	[75]
	8839G>C	Retinitis pigmentosa syndrome (NARP)	[76]
	8914C>T	Mitochondrial encephalomyopathies	[77]
	8593 A>G	Leigh syndrome with a deficiency in mitochondrial energy production	[78]
ATP Synthetase 8 gene (ATP 8)	Mutations in MT-ATP8	LHON, MELAS, Leigh syndrome, NARP	[76,79]
COII gene	7750 C>A	SCZ	[80]

MT-: mitochondrial, NADH: nicotinamide adenine dinucleotide hydride dehydrogenase, ND1: NADH dehydrogenase 1, ND2: NADH dehydrogenase 2, ND3: NADH dehydrogenase 3, ND4: NADH dehydrogenase 4, ND4L: NADH dehydrogenase 4L, ND5: NADH dehydrogenase 5, and ND6: NADH dehydrogenase 6; CYB: Cytochrome b, COX II: cytochrome oxidase subunit II, and COX III: cytochrome oxidase subunit III.

4. Mitochondrial Genome Mutations/Variations and Male Infertility

It is estimated that approximately 15% to 30% of cases of male infertility are caused by genetic defects [81]. Mitochondria have their own genome (more than 93,000 genes). Some genes related to fertility and longevity are also known to be linked to mitochondria (Table 2).

Sperm have only 100 mtDNA copies compared to oocytes, which have 150,000 [82]. Spermatozoa rely on mitochondrial oxidative phosphorylation (OxPhos) machinery to generate the energy required for their motility [83].

Single-nucleotide polymorphisms or large deletions are types of mutations that affect sperm mtDNA. Thus, mutant mtDNA in sperm can lead to respiratory dysfunction that impairs energy production and results in decreased motility, which affects the normal activity of sperm [84,85,86,87].

Complex 1 plays a key role in OxPhos by receiving electrons from NADH. The captured energy of these electrons is used to release protons into the intermembrane space, which is subsequently used to generate ATP [88]. In addition, the mitochondrial genes ATPase 6, ATPase 8, COX 3, COX 2, CYB, ND3, ND4, ND5, and ND6 play important roles in the construction of mature sperm and in progressive flagellar motility after ejaculation [87].

Therefore, pathogenic variants in the ND gene are expected to affect complex 1 activity, resulting in insufficient energy production, which will negatively affect sperm motility [84].

Large mtDNA deletions, including 4977 bp, 7345 bp, 7436 bp, and 7599 bp, have been associated with asthenozoospermia [89,90,91,92]. These deficiencies include deletions of several mitochondrial genes, namely ATPase8 (lost at 7599 bp only), ATPase 6, COX III, CYB, and NADH dehydrogenase (ND) 3, 4, 4L, 5, and 6. These deletions also include the loss of eight tRNA genes [91]. Gene deficiencies from mtDNA deletions play critical roles in OxPhos in mitochondria; thus, their removal from mtDNA reduces energy conservation, which, in turn, negatively affects sperm-flagellar motility and contributes to asthenozoospermia [92].

The deletion of 4977 bp is considered the most common among mtDNA deficiencies [93]. This deletion involves the removal of seven genes and five transfer RNAs located between 8483 bp and 13459 bp in mtDNA, and the deleted site is located between two 13 bp repeats (5′-ACCTCCCTCACCA-3′) [94].

A previous study reported that sperm motility is inversely correlated with mtDNA deletion. According to the study’s results, mtDNA deletions were found in immobilized spermatozoa, whereas spermatozoa with normal motility did not have these deletions [95]. Moreover, it has been noted that the percentage of this deletion in mtDNA is higher in normal-motility sperm than in low-motility sperm [90].

Chari et al. (2015) found sperm with an abnormal motility which have a 4866 bp deletion in their mtDNA [96]. Another study found that 7436 bp sperm mtDNA deletions occurred more frequently in asthenozoospermic men than in normozoospermic men. Therefore, the use of this deletion as an indicator of reduced sperm motility has been proposed [89].

The change in mtDNA is likely to affect the speed and quality of sperm, according to studies performed by La Vignera et al. (2019). Their study showed that mtDNA alterations could have a significant impact on the quality of sperm as well as on its ability to swim [97].

In contrast, another study found no significant difference in mtDNA deletions in sperm from asthenozoospermic and normozoospermic men. Therefore, the researchers chose to ignore the role of these deletions in male infertility [14].

A previous study also showed that mtDNA deletions did not affect sperm motility, as the deletions did not differ significantly in the occurrence of mtDNA deficiencies in poor-quality and high-quality sperm [98].

A study in New Zealand comparing the incidence of mitochondrial variants in infertile and fertile men found 11 nucleotide substitutions in a group of infertile patients, compared with only 7 nucleotide substitutions in men with normal sperm. In addition, they identified two SNPs in ATPase6, and ND4 mitochondrial genes were found to be associated with asthenozoospermia. Furthermore, 9055G>A occurred at a frequency of 10.7% in the ATP synthase 6 gene, and this missense variant resulted in the substitution of alanine by threonine [99].

In the Güney study, three different gene variants were identified in the ND1 gene in patients who were infertile (C4159, G4153A, and T4114G). These gene variants were only found in patients who were unable to generate sperm, and not in patients who were able to produce sperm [100].

A previous study reported that the presence of massive mtDNA deletion is associated with asthenozoospermia because it results in the loss of some mitochondrial genes responsible for mitochondrial respiration, which provides sperm with the energy they need to move, thereby affecting male fertility [101].

An 11696G>A variant in the MT-ND4 gene caused by a missense change was also found to be related to reduced sperm motility. This variant caused a change from valine at position 313 to isoleucine, which altered the structure of the protein [102].

Barbhuiya et al. (2016) found that the genes with the highest number of SNPs were the ATPase 6 gene (21 SNPs), followed by the ND2 gene (12 SNPs), and the ATPase 8 gene (9 SNPs). The ND4 gene showed the lowest number of SNP point mutations [103].

A genetic mutation causing a change in SNP T4216C in both fertile and infertile males was found by Khan et al. (2016) [104]. Zhang et al. also noted a decrease in the risk of asthenozoospermia for SNP C3398T, as it has a low frequency and small sample size [105]. Mughal et al. discovered that there is a strong link between the 15 bp deletion of COX III (at location 9390 to 9413) and infertility in men, with a P value of 0.033 [106].

Three different missense genetic variants were found in Complex I by Alsmadi et al. (2021) by screening mitochondrial genes ND6, ND5, ND2, and ND1 in the sperm of male partners. The variants were 13708 G>A, 4216 T>C, and 12506T>A. These three variants were correlated negatively with sperm motility and ICSI outcomes. They were significantly different from wild types in terms of sperm motility, fertilization rate, embryo quality score, and median embryo cleavage score [107].

Recently, we scanned sterile and fertile males for polymorphisms by the direct sequencing of MTND3, MTND4L, and MTND4 genes in our laboratory [108].

An SNP called rs2853495 in the MTND4 gene was found to have an association with male infertility in the genotype frequency test. Moreover, G11719A (rs2853495) and A11251G (rs869096886) in gene MTND4 were also tied to male infertility. Therefore, even if the genotype of the allele is different, the presence of the allele itself could still be tied to male infertility [109].

In addition, a total of 49 SNPs were identified and genotyped: 13 SNPs in MT-CYB, 14 SNPs in MT-ATP6, and 10 SNPs in MT-ATP8 [110].

For three variants in MT-CYB, genotype frequencies were significantly different between fertile and sterile groups, and two SNPs showed a significant association between male infertility and allele frequencies: rs41504845 (C15833T) (p = 0.0147) and rs527236194 (T15784C) (p = 0.0014)). Furthermore, for MT-CO3 and MT-ATP6, only rs7520428 showed statistically significant differences between the subfertile and fertile groups in both genotype and allele frequency tests (p < 0.0001 for both) [110].

Table 2. Summary of studies reporting different mitochondrial genome abnormalities associated with male infertility.

Mitochondrial Genome Abnormalities	Description	Effect on Male Infertility	References
4977 bp deletion	Most common deletion, located between 8483 bp and 13459 bp and characterized by the presence of two 13-bp repeated sequences (5′-ACCTCCCTCA CCA-3′)	Removal of seven genes and five tRNAs in mitochondrial DNA (mtDNA) associated with asthenozoospermia.	[90,93,111,112,113]
7599 bp deletion	Located between 8642 and 16243-bp and characterized by the presence of 7 nucleotides’ direct repeat (5′CATCAAC-3′) on both sides	-Removal of several mitochondrial genes: ATP8 (lost with 7599 bp only), ATP6, cytochrome oxidase (COX) III, cytochrome b (CYB), NADH dehydrogenase (ND) 3,4, 4L, 5, and 6. -Reduction in the obtained energy, which in turn has a negative effect on sperm flagellum movement and leads to asthenozoospermia.	[91,92]
7345 bp deletion	Located between 9009 and 1654-bp
4216 T>G	Located on MT-ND1 gene	Negative correlation with sperm motility.	[107]
3243A>G		Positive correlation with the mt DNA copy number in embryo after ICSI as an adaptation for inefficient ATP production via oxidative phosphorylation due to mutated mtDNA.	[114]
A point mutation in the ND1 gene at locus 4216		Association with recurrent pregnancy loss.	[115]
Copy-number variations (CNV)		mtDNA copy number affects implantation rate after ICSI.	[116]
11719G>A	Located on MT-ND4 gene	Association with male infertility.	[109]
11251A>G
9055 G>A		Association with poor sperm quality.	[99]
11696G>A (missense variant)		Association with reduced sperm motility.	[102]
11719 G>A		Association with poor semen quality.	[99]
11994 C>T (missense variant)		Negative correlation with sperm motility.	[117]
12506T>A	Located on MT-ND5 gene	Negative correlation with sperm motility.	[107]
13708 G>A (missense variant)
14172 T>C	Located on MT-ND6 gene	Significant difference between the total fertilization failure group and control.	[118]
14368 C>T
G15301A	Located on MT-CYB	These SNPs showed a statistically significant link to male infertility.	[110]
A 15326G
A 15487 T
15 bp deletion of cytochrome c oxidase III	Location 9390 to 9413	This deletion linked to human male infertility.	[106]

MT-: mitochondrial, NADH: nicotinamide adenine dinucleotide hydride dehydrogenase, ND1: NADH dehydrogenase 1, ND3: NADH dehydrogenase 3, ND4: NADH dehydrogenase 4, ND4L: NADH dehydrogenase 4L, ND5: NADH dehydrogenase 5, and ND6: NADH dehydrogenase 6; CYB: Cytochrome b, and COX III: cytochrome oxidase subunit III.

Table 2. Summary of studies reporting different mitochondrial genome abnormalities associated with male infertility.

Mitochondrial Genome Abnormalities	Description	Effect on Male Infertility	References
4977 bp deletion	Most common deletion, located between 8483 bp and 13459 bp and characterized by the presence of two 13-bp repeated sequences (5′-ACCTCCCTCA CCA-3′)	Removal of seven genes and five tRNAs in mitochondrial DNA (mtDNA) associated with asthenozoospermia.	[90,93,111,112,113]
7599 bp deletion	Located between 8642 and 16243-bp and characterized by the presence of 7 nucleotides’ direct repeat (5′CATCAAC-3′) on both sides	-Removal of several mitochondrial genes: ATP8 (lost with 7599 bp only), ATP6, cytochrome oxidase (COX) III, cytochrome b (CYB), NADH dehydrogenase (ND) 3,4, 4L, 5, and 6. -Reduction in the obtained energy, which in turn has a negative effect on sperm flagellum movement and leads to asthenozoospermia.	[91,92]
7345 bp deletion	Located between 9009 and 1654-bp
4216 T>G	Located on MT-ND1 gene	Negative correlation with sperm motility.	[107]
3243A>G		Positive correlation with the mt DNA copy number in embryo after ICSI as an adaptation for inefficient ATP production via oxidative phosphorylation due to mutated mtDNA.	[114]
A point mutation in the ND1 gene at locus 4216		Association with recurrent pregnancy loss.	[115]
Copy-number variations (CNV)		mtDNA copy number affects implantation rate after ICSI.	[116]
11719G>A	Located on MT-ND4 gene	Association with male infertility.	[109]
11251A>G
9055 G>A		Association with poor sperm quality.	[99]
11696G>A (missense variant)		Association with reduced sperm motility.	[102]
11719 G>A		Association with poor semen quality.	[99]
11994 C>T (missense variant)		Negative correlation with sperm motility.	[117]
12506T>A	Located on MT-ND5 gene	Negative correlation with sperm motility.	[107]
13708 G>A (missense variant)
14172 T>C	Located on MT-ND6 gene	Significant difference between the total fertilization failure group and control.	[118]
14368 C>T
G15301A	Located on MT-CYB	These SNPs showed a statistically significant link to male infertility.	[110]
A 15326G
A 15487 T
15 bp deletion of cytochrome c oxidase III	Location 9390 to 9413	This deletion linked to human male infertility.	[106]

MT-: mitochondrial, NADH: nicotinamide adenine dinucleotide hydride dehydrogenase, ND1: NADH dehydrogenase 1, ND3: NADH dehydrogenase 3, ND4: NADH dehydrogenase 4, ND4L: NADH dehydrogenase 4L, ND5: NADH dehydrogenase 5, and ND6: NADH dehydrogenase 6; CYB: Cytochrome b, and COX III: cytochrome oxidase subunit III.

5. Mitochondrial Genome Mutations/Variations and Fertilization/Pregnancy Outcomes

Until a few years ago, maternally restricted mtDNA inheritance was the only accepted idea because paternal mtDNA disappears after the embryonic cleavage stage [14]. In 2018, Luo et al. presented strong proof of biparental mtDNA inheritance, showing evidence of parental mtDNA transmission from father to offspring according to mitochondrial disease inheritance patterns in three independent multigenerational families [15]. Ecker found that sons with ICSI had the same SNPs in mitochondrial genes (COX1, ND1, ND4, and ND5) as their fathers [20]. Paternal mtDNA similarity is sometimes as high as 99% [20].

Another study observed that mtDNA myopathy could be passed from father to son by ICSI and found that sperm mtDNA mutations were retained in the embryo [119]. This suggests that, at least in some cases, paternal mtDNA can be transmitted to offspring. Paternal mtDNA mutations can be diluted due to mitochondrial heterogeneity, as maternal mtDNA copies are much more numerous than paternal mtDNA [120].

Diez-Juan et al. demonstrated a strong correlation between implantation rate and mtDNA copy number in euploid embryos. They found that successfully implanted embryos had lower mtDNA content; this result applies to embryo transfer at the division or blastocyst stage. On the other hand, poorly implanted embryos are estimated to have high mtDNA copy numbers. Therefore, researchers believe that mtDNA copy number can be reliably used to predict successful engraftment [116].

A study of a Chinese population found that men with haplogroup Z were more prone to IVF failure, but fertilization was observed between fertile and infertile cohorts with haplogroup D or haplogroup G. There was no significant difference in failure. Therefore, it was concluded that the SNPs of the MT-ND2 gene were not associated with overall fertility [121].

In a recent study, variations in mtDNA correlated with lower levels of embryo development compared to quality at the blastocyst stage. The quality of the embryo at the blastocyst stage correlated with better sperm motility [122].

Forty-five ND6 gene variants were found, two of which had a significant difference in total fertilization failure between the groups studied [118].

Alsmadi et al. (2021) found that pregnancy rates were negatively affected by the presence of 4216T>C in ND1, and 13708G>A and 12506T>A variants in ND5. The lowest pregnancy rates were due to men with severe asthenozoospermia and mtDNA SNPs, while the highest pregnancy rates were due to normozoospermic men without mtDNA SNPs [107].

These results are consistent with a previous study that found a point mutation in the T>C variant at locus 4216 in the ND1 gene to be associated with recurrent miscarriage [115].

One variant in the ND4 gene was investigated in two separate studies. In one of them, a missense variant, 11994C>T, was found to be strongly associated with oligoasthenozoospermia in India [117]. In the other study, no association was found between this variant and oligoasthenozoospermia [123].

6. Conclusions

Human male infertility is a complex phenotype. Research into the field of genetics associated with this phenotype is increasing, and many genes, including their mutations and/or deletions, have been identified. In addition, mitochondria and their genomes have received attention due to their importance in sperm function and successful fertilization. To date, various mitochondrial variants associated with male infertility have been observed. However, more research is needed to understand how these changes lead to adverse outcomes in men.