|Ahead of print publication
The association of stromal antigen 3 (STAG3) sequence variations with spermatogenic impairment in the male Korean population
Yeojung Nam1, Kyung Min Kang2, Se Ra Sung2, Ji Eun Park2, Yun-Jeong Shin2, Seung Hun Song3, Ju Tae Seo4, Tae Ki Yoon5, Sung Han Shim1,2
1 Department of Biomedical Science, College of Life Science, CHA University, Seongnam 13488, Korea
2 Genetics Laboratory, Fertility Center of CHA Gangnam Medical Center, CHA University, Seoul 06135, Korea
3 Department of Urology, CHA Gangnam Medical Center, CHA University, Seoul 06135, Korea
4 Department of Urology, Cheil General Hospital, Seoul 04619, Korea
5 Department of Obstetrics and Gynecology, CHA Gangnam Medical Center, Seoul 04637, Korea
|Date of Submission||12-Oct-2018|
|Date of Acceptance||18-Jan-2019|
|Date of Web Publication||17-May-2019|
Sung Han Shim,
Department of Biomedical Science, College of Life Science, CHA University, Seongnam 13488, Korea; Genetics Laboratory, Fertility Center of CHA Gangnam Medical Center, CHA University, Seoul 06135, Korea
Source of Support: None, Conflict of Interest: None
The stromal antigen 3 (STAG3) gene, encoding a meiosis-specific cohesin component, is a strong candidate for causing male infertility, but little is known about this gene so far. We identified STAG3 in patients with nonobstructive azoospermia (NOA) and normozoospermia in the Korean population. The coding regions and their intron boundaries of STAG3 were identified in 120 Korean men with spermatogenic impairments and 245 normal controls by using direct sequencing and haplotype analysis. A total of 30 sequence variations were identified in this study. Of the total, seven were exonic variants, 18 were intronic variants, one was in the 5'-UTR, and four were in the 3'-UTR. Pathogenic variations that directly caused NOA were not identified. However, two variants, c.3669+35C>G (rs1727130) and +198A>T (rs1052482), showed significant differences in the frequency between the patient and control groups (P = 0.021, odds ratio [OR]: 1.79, 95% confidence interval [CI]: 1.098–2.918) and were tightly linked in the linkage disequilibrium (LD) block. When pmir-rs1052482A was cotransfected with miR-3162-5p, there was a substantial decrease in luciferase activity, compared with pmir-rs1052482T. This result suggests that rs1052482 was located within a binding site of miR-3162-5p in the STAG3 3'-UTR, and the minor allele, the rs1052482T polymorphism, might offset inhibition by miR-3162-5p. We are the first to identify a total of 30 single-nucleotide variations (SNVs) of STAG3 gene in the Korean population. We found that two SNVs (rs1727130 and rs1052482) located in the 3'-UTR region may be associated with the NOA phenotype. Our findings contribute to understanding male infertility with spermatogenic impairment.
Keywords: linkage disequilibrium; meiotic-specific gene; single-nucleotide variations; spermatogenic impairment; stromal antigen 3 gene
Article in PDF
|How to cite this URL:|
Nam Y, Kang KM, Sung SR, Park JE, Shin YJ, Song SH, Seo JT, Yoon TK, Shim SH. The association of stromal antigen 3 (STAG3) sequence variations with spermatogenic impairment in the male Korean population. Asian J Androl [Epub ahead of print] [cited 2019 Oct 23]. Available from: http://www.ajandrology.com/preprintarticle.asp?id=258492
Yeojung Nam, Kyung Min Kang
These authors contributed equally to this work.
| Introduction|| |
Azoospermia affects approximately 1% of the male population, accounts for over 15% of all male infertility,,, and includes genital tract obstruction (obstructive azoospermia) and spermatogenic impairment (nonobstructive azoospermia, NOA)., NOA is mainly caused by severely impaired spermatogenesis and is reported to account for more than 70% of azoospermia in Korean patients. Chromosomal abnormalities, such as Klinefelter syndrome, balanced chromosomal rearrangements, and Yq microdeletions, are well known genetic causes of NOA. In many cases, the genetic etiology remains unknown. Matzuk and Lamb reviewed many genes involved in spermatogenesis and mutations in some of those genes were identified in patients with NOA.,,
The stromal antigen 3 (STAG3) gene was mapped to chromosome 7 and consists of 34 exons encoding a protein involved in the meiotic cohesion complex. Human STAG3 is highly expressed in the testis and several other organs including the ovary., During meiosis, STAG3 forms a cohesion core with three other proteins, structural maintenance of chromosome 3 (SMC3), structural maintenance of chromosomes 1β (SMC1β), and Rad21 cohesin complex component like 1 (Rad21L1)., In mice, defective Stag3 protein causes aberrant meiotic chromosomal features and infertility., In humans, a homozygous 1-bp deletion in STAG3 has been found in a consanguineous family with premature ovarian failure (POF), and a homozygous donor splice-site mutation has been found in two sisters with premature ovarian insufficiency (POI). Therefore, STAG3 has been suggested as a strong candidate gene target for causing male infertility.,, To date, no homozygous or compound heterozygotic mutations of STAG3 have been identified in patients with spermatogenic impairment, and studies of STAG3 mutations have not been examined in infertile male populations.
In this study, we investigated whether STAG3 variations may be a genetic cause of spermatogenic impairment in Korean men.
| Patients and Methods|| |
A total of 120 Korean men with spermatogenic impairment (43 oligozoospermic and 77 azoospermic) and 245 normal controls were obtained from the Cha Gangnam Medical Center at Cha University, Seoul, Korea, between January 2010 and December 2012. General and clinical characteristics of patients and controls are presented in [Table 1]. Patients with tubule obstruction, chromosomal abnormalities, or microdeletion of the Y chromosome AZF region were excluded. Normal controls had a normal sperm concentration and no history of infertility. Testicular size was measured by a Prader orchidometer (Pro-Health Product Ltd., Guangzhou, China). Serum testosterone, follicle-stimulating hormone (FSH), and luteinizing hormone (LH) were measured on a Cobas e601 analyzer (Roche Dignostics, Penzberg, Germany) using electrochemiluminescence immunoassay (ECLIA) method. Semen analysis was performed according to the World Health Organization criteria (WHO, 2010). The study was approved by the Institutional Review Board of Cha Gangnam Medical Center, Seoul, Korea, and written informed consent was obtained from all participants.
Genomic DNA was isolated from peripheral blood of the patient and control samples with the QuickGene DNA blood kit (KURABO industries, Neyagawa, Japan), according to the manufacturer's instructions. DNA yield was quantified with the NanoDrop™ spectrophotometer (Thermo Scientific, Wilmington, DE, USA). The extracted DNA was stored at −80°C until further analysis.
Polymerase chain reaction
Coding regions of STAG3 (NM_001282716.1) were amplified for genetic screening by polymerase chain reaction (PCR), using primers for the 34 exons and their intron boundaries that were designed by Primer 3 (http://primer3.ut.ee). As this gene has multiple pseudogenes, we performed long-range PCR that produced eleven fragments and designed eleven primer pairs to cover the 34 exons and avoid pseudogene sequences. The locations and sequences of primer sets are presented in [Supplementary Table 1 [Additional file 1]].
All eleven PCR products were purified with ExoSAP-IT (USB Corporation, Cleveland, OH, USA). To sequence the eleven fragments, we designed 30 sequencing primers to cover 34 exons [Supplementary Table 1]. All the samples were amplified by PCR and sequenced bidirectionally using the BigDye® Terminator v3.1 Cycle Sequencing Kit (Applied Biosystems, Austin, TX, USA) and BigDye® X-Terminator™ solutions (Applied Biosystems, Bedford, MA, USA) with standard conditions. The sample supernatant was loaded on the ABI 3130XL Genetic Analyzer (Applied Biosystems, Foster City, CA, USA) and processed with a BigDye® X-Terminator run module. All assays were repeated once for confirmation and the results matched over 99.0%.
Statistical analyses and database search
For each sequence variation, the data were statistically analyzed with Statistical Package for the Social Sciences software (SPSS version 22, IBM Software Group, Chicago, IL, USA). To evaluate the association between patient and control groups, odds ratios (OR), 95% confidence intervals (95% CI), and the applied P values were calculated from both the Chi-squared test and Fisher's exact test (two-tailed). Applied P < 0.05 were considered statistically significant. Three databases, Polyphen-2, SIFT, and Mutation Tester, were used to predict potentially damaging effects due to amino acid changes.
Multiple hypothesis testing was performed with the Benjamini–Hochberg method to control false discovery rate (FDR) in the logistic regression analysis. Calculating the FDR is a way to address problems associated with multiple comparisons, and FDR provides a measure of the expected proportion of false-positives in the data.
Haplotype block structure was established by HaploView 4.1 software (https://www.broadinstitute.org) using the method of block definition of Gabriel et al. Haplotype association tests were also conducted with this software.
The entire 3'-untranslated region (3'-UTR) of STAG3 was amplified from genomic DNA, which contained the rs1052482 A or T alleles, using primers that included Sac I and Xba Irestriction sites. Primer sequences of the STAG3 3'-UTR were: forward: 5'-GAGCTCccgttgctgtgtcctgtgta, reverse: 5'-TCTAGAgaccaagaacctgacctcca (for a predicted 476 bp product). PCR products were cloned into the pmirGLO vector (Promega, Madison, WI, USA) via the Sac I and Xba I sites and all constructs were verified by DNA sequencing.
Dual luciferase assay
HEK293T cells were seeded in 48-well plates (3 × 10 per well). After 24 h, 200 ng pmir-rs1052482A, pmir-rs1052482T, or pmir-empty vector was transiently transfected with Lipofectamine 2000 (Invitrogen, Carlsbad, CA, USA). The Renilla vector was used as an internal control for transfection efficiency. After 48 h, the transfected cells were harvested and lysed with a dual-luciferase reporter assay system (Promega), and the activities of Firefly and Renilla were measured in a Luminometer, Centro XS LB960 (Berthold Technologies, Bad Wildbad, Germany), and MikroWin2000 software (https://mikrowin-20001.software.informer.com). Transfection experiments were performed in triplicate, and activity measurements were done for three times. Relative luciferase activity was determined by normalizing firefly luciferase activity against Renilla luciferase activity. An average value of firefly/Renilla was calculated and then normalized to the average value of the empty vector to yield the vector-normalized ratio.
Computational prediction of putative targets for STAG3 mRNA was performed by searching mirmap.ezlab.org, www.targetscan.org, www.microrna.org, and www.mirdb.org for the target prediction algorithms. From the potential miRNAs interacting with STAG3 mRNA 3'-UTR, we selected seven candidate miRNAs (miR-148a, miR-2909, miR-3162-5p, miR-33a-5p, miR-33b-5p, miR-4739, and miR-6508-3p) that allowed rs1052482A or T alleles to be included in the seed sequence. Seven miRNAs and their inhibitors were constructed by Bioneer (Daejeon, Korea). The prediction score of the seven miRNAs and the 3'-UTR of STAG3 and details for their in silico interactions are presented in [Supplementary Figure 1 [Additional file 2]] and [Supplementary Table 2 [Additional file 3]]. Mimics (5–10 pmol) were transiently cotransfected with pmir-vectors with Lipofectamine 2000 according to the manufacturer's instructions.
| Results|| |
Thirty single-nucleotide variations (SNVs) in STAG3 were identified in this study. The locations, types, frequencies, and P values of the variations are presented in [Supplementary Table 3 [Additional file 4]]. The distributions of genotypes of all the SNVs followed the Hardy–Weinberg equilibrium in patients and controls. Seven were exonic, 18 were intronic, one was in the 5'-UTR, and four were in the 3'-UTR of the 30 variations. Three exonic variants were nonsynonymous and the other four variants were synonymous. Three variants were found in a patient but not in controls (c.1269C>T p. Asp423, +112G>A, +315C>T); two variants, c.3669+35C>G and +198A>T, showed significant differences in the frequency between patient and control groups (P = 0.021, OR: 1.79, 95% CI: 1.098–2.918). Haplotype analysis by HaploView 4.1 showed that nineteen variants were separated into five linkage disequilibrium (LD) blocks [Figure 1]. As shown in [Table 2], in particular, the frequencies of G-C-A (case: control = 0.504: 0.606, P = 0.009), G-G-T (case: control = 0.162: 0.075, P < 0.001), and C-C-A (case: control = 0.075: 0.018, P < 0.001) haplotypes in block 5 were significantly different between patients and controls.
|Table 2: Haplotype analysis between the SNPs of STAG3 and nonobstructive azoospermia |
Click here to view
[Table 3] shows a genetic model of the 3 SNVs (rs2246713, rs1727130, and rs1052482) between the cases and the controls. We found that the individuals with the CG genotype of rs1727130 and AT genotype of rs1052482 had an increased risk susceptibility to NOA in the codominant model, and those with the minor allele G of rs1727130 and T of rs1052482 had an increased risk susceptibility to NOA in the dominant model (P = 0.039, OR: 1.64, 95% CI: 1.030–2.608).
|Table 3: The genotype distributions of STAG3 rs2246713, rs1727130, and rs1052482 in the cases and the controls |
Click here to view
|Figure 1: LD pattern in the locus of STAG3 gene. Nineteen variants were separated into five LD blocks. Numbers in the squares indicate D' index (level of LD) between the corresponding SNPs. STAG3: stromal antigen 3; SNPs: single-nucleotide polymorphism; LD: linkage disequilibrium.|
Click here to view
To determine whether variations in the 3'-UTR region affected miRNA-mediated gene expression regulation, miRNAs predicted to interact with STAG3 were examined by in silico analysis, and seven miRNAs, miR-148a, miR-2909, miR-3162-5p, miR-33a-5p, miR-33b-5p, miR-4739, and miR-6508-3p, were found potentially to interact with the 3'-UTR of STAG3 mRNA. The effect of these seven miRNAs on +198A>T variation (rs1052482) was examined by a luciferase assay. There was no significant difference in relative luciferase activities between rs1052482A and rs1052482T [Figure 2]a, but a substantial decrease in luciferase activity was observed in cells cotransfected with pmir-rs1052482A and miR-3162-5p (mean ± standard deviation [s.d.]: 51.6% ± 2.5%, P = 0.002), compared with control vector pmir-GLO or pmir-rs1052482T (mean ± s.d.: 85.9% ± 3.6%, P = 0.061) [Figure 2]b. A sequential decrease was observed in cells cotransfected with pmir-rs1052482A and miR-3162-5p in proportion to the amount of the mimic (rs1052482A [mean ± s.d.], 63.7% ± 1.7%: 45.0% ± 2.7%; rs1052482T [mean ± s.d.], 89.2 ± 5.0%: 78.3% ± 1.4%, respectively), and this reduction was inhibited when the mir-3162-5p inhibitor was cotransfected with pmir-rs1052482A and miR-3162-5p [Figure 2]c. The data indicate that miR-3162-5p targets the rs1052482A sequence more efficiently than that of rs1052482T.
|Figure 2: Relative luciferase activities of pmir-rs1052482A and pmir-rs1052482T in HEK293T cells. a - c represent the mean RLU values ± s.d. of triplicates. ( a ) No significant differences in relative luciferase activity among three vectors. ( b ) Luciferase activity was substantially decreased in cells cotransfected with pmir-rs1052482A and miR-3162-5p compared to the cells cotransfected with pmir-control and miR-3162-5p (P < 0.01). ( c ) A sequential decrease was observed in cells cotransfected with pmir-rs1052482A and miR-3162-5p in proportion to the amount of the mimic and this reduction was rescued when the mir-3162-5p inhibitor was cotransfected with pmir-rs1052482A and miR-3162-5p. -: absence of the indicated one; +: presence of the indicated one; *P < 0.05 (pmir-rs1052482A + mir3162-5p [5 pmol] compared to pmir-rs1052482A); **P < 0.01 (pmir-rs1052482A + mir3162-5p [10 pmol] compared to pmir-rs1052482A). CTL: control; RLU: relative light unit; s.d.: standard deviation.|
Click here to view
| Discussion|| |
Many autosomal genes such as ring finger protein 212(Rnf212), testis expressed 15(Tex15), syntaxin 2(Stx2), and siah E3 ubiquitin protein ligase 1A(Siah1a)are reported to be crucial factors for the meiotic process and spermatogenesis in mouse studies.,,, Human homologs of these genes also play a role in meiosis, and variations of these genes are thought to induce spermatogenic impairment.,,, In recent studies, the homologous deletion of Stag 3 has been shown to induce sterility associated with the premature arrest of meiotic prophase I in both male and female mice., Therefore, we have investigated the association between STAG3 gene variations and male infertility.
We identified 30 variations in the coding regions and intron boundariesof STAG3 in patients with NOA and in control samples. Of the 30 variations, seven were exonic and three were found only in different infertile patients. For these variations, we evaluated the potential pathogenic effects by three prediction methods, Polyphen-2, SIFT, and Mutation Taster [Table 4]. Most variations (six of nine) were considered benign and three variations did not show consistent results in the three predictive programs. Minor allele frequencies (MAFs) of these variations were not significantly different between patients and controls, and these MAFs were similar to those of other Asian populations on the NCBI SNP database from the 1000 Genomes Project. Considering these data, the above-mentioned variations are unlikely to be related to the spermatogenic impairment in our Korean male population. Interestingly, two variations, c.3669+35C>G and +198A>T (rs1727130 and rs1052482) located in 3'-UTR, had a significantly different frequency between the patient and control groups. However, there is a discrepancy in this result. According to Yu et al., there is no significant difference in the frequencies of allele and genotype at SNP rs1052482 between patients with NOA and controls and they suggested that this SNP is not associated with azoospermia. Two variants of rs1727130 and rs1052482 are close together (372 bp apart) and tightly linked, as shown in the LD block analysis. The allele distribution of 2 SNVs between the patient and control groups is more evident in the LD block analysis. On the basis of the data, we propose that multiple SNVs linked to a block can interact with each other to regulate gene function, rather than allowing each SNV to function independently.
|Table 4: In silico analysis for exonic variations and three variations found only in patients |
Click here to view
The presence of SNVs in the 3'-UTR of genes may interfere with mRNA stability and translation through effects on polyadenylation and regulatory protein–mRNA and miRNA–mRNA interactions, or may locally alter secondary structures of mRNAs, affecting the accessibility of binding sites for interacting transelements.,,
We investigated whether rs1052482, the SNV in the 3'-UTRof STAG3, could affect interaction with miRNAs and thus affect posttranscriptional repression of STAG3. When pmir-rs1052482A was cotransfected with miR-3162-5p, a substantial decrease was observed in luciferase activity compared with pmir-rs1052482T. This result suggests that rs1052482 is located within a binding site for miR-3162-5p in the STAG3 3'-UTR, and the minor rs1052482T allele may offset the inhibition by miR-3162-5p.
According to Fukuda et al.,, STAG3 interacts with the three different α-kleisin subunits present in mammalian meiotic cells, depending on the temporal and spatial distribution. STAG3 combined with meiotic recombination protein (REC8), one of the α-kleisin subunits, and promoted synapsis between homologous chromosomes, while the same complexes inhibited synaptonemal complex assembly between sister chromatids. Therefore, we hypothesize that STAG3 with the rs1052482T variation reduces the normal inhibitory function of mir-3162-5p, thereby increasing the amount of STAG3 protein and ultimately disturbing synapses between homologous chromosomes or sister chromatids. However, it is unclear how elevated STAG3 may affect meiotic chromosome dynamics. In previous studies, mir-3162-5p was identified as a regulator of prostate or cervical antigen in cell carcinomas.,, However, little is known about mir-3162-5p's regulation of human meiosis, because of the difficulties in uncovering the spatiotemporal and sequential expression of miRNAs in human germ cells, and identifying which miRNAs are the actual operators for the onset of meiosis or spermatogenesis.
| Conclusion|| |
We have identified 30 SNVs of STAG3 in the Korean population. Pathogenic variations that directly cause NOA were not identified. However, we found that two SNVs, rs1727130 and rs1052482, located in the 3'-UTR region may be associated with the NOA phenotype through the regulation of miRNA. Further studies are needed to determine whether variations in the 3'-UTR region of STAG3 actually affect gene expression through miRNAs, including mir-3162-5p in germ cells. While there is still much to learn about the exact mechanisms regulating human meiosis or spermatogenesis, our findings contribute to the understanding of spermatogenic impairment, as well as the identification of predictive susceptibility biomarkers.
| Author Contributions|| |
SH Shim conceived and designed the article; SH Shim, YN, and KMK drafted the manuscript; YN, KMK, and SRS designed and performed the experiments; SRS, JEP, and YJS analyzed data and interpreted findings. SH Song collected the samples and performed andrology workup. JTS and TKY prepared the publication. All authors read, edited, and approved the final manuscript.
| Competing Interests|| |
All authors declared no competing interests.
| Acknowledgments|| |
This research was supported by a grant from the Korea Health Technology R&D Project through the Korea Health Industry Development Institute (KHIDI), funded by the Ministry of Health & Welfare, Republic of Korea (Grant No. H14C0106020014). Funding for this study was obtained by SH Shim.
Supplementary Information is linked to the online version of the paper on the Asian Journal of Andrology website.
| References|| |
Cooper TG, Hellenkemper B, Jonckheere J, Callewaert N, Grootenhuis AJ, et al.
Azoospermia: virtual reality or possible to quantify? J Androl
2006; 27: 483–90.
Esteves SC, Miyaoka R, Agarwal A. An update on the clinical assessment of the infertile male. [corrected]. Clinics (Sao Paulo)
2011; 66: 691–700.
Gudeloglu A, Parekattil SJ. Update in the evaluation of the azoospermic male. Clinics (Sao Paulo)
2013; 68 Suppl 1: 27–34.
Kolettis PN. The evaluation and management of the azoospermic patient. J Androl
2002; 23: 293–305.
Esteves SC, Prudencio C, Seol B, Verza S, Knoedler C, et al.
Comparison of sperm retrieval and reproductive outcome in azoospermic men with testicular failure and obstructive azoospermia treated for infertility. Asian J Androl
2014; 16: 602–6.
Lee HD, Lee HS, Park SH, Jo DG, Choe JH, et al.
Causes and classification of male infertility in Korea. Clin Exp Reprod Med
2012; 39: 172–5.
Esteves SC. Clinical management of infertile men with nonobstructive azoospermia. Asian J Androl
2015; 17: 459–70.
Matzuk MM, Lamb DJ. The biology of infertility: research advances and clinical challenges. Nat Med
2008; 14: 1197–213.
Sun F, Turek P, Greene C, Ko E, Rademaker A, et al.
Abnormal progression through meiosis in men with nonobstructive azoospermia. Fertil Steril
2007; 87: 565–71.
Sun F, Trpkov K, Rademaker A, Ko E, Martin RH. Variation in meiotic recombination frequencies among human males. Hum Genet
2005; 116: 172–8.
Sun F, Kozak G, Scott S, Trpkov K, Ko E, et al.
Meiotic defects in a man with non-obstructive azoospermia: case report. Hum Reprod
2004; 19: 1770–3.
Pezzi N, Prieto I, Kremer L, Perez Jurado LA, Valero C, et al.
STAG3, a novel gene encoding a protein involved in meiotic chromosome pairing and location of STAG3-related genes flanking the Williams-Beuren syndrome deletion. FASEB J
2000; 14: 581–92.
Houmard B, Small C, Yang L, Naluai-Cecchini T, Cheng E, et al.
Global gene expression in the human fetal testis and ovary. Biol Reprod
2009; 81: 438–43.
Nogues C, Fernandez C, Rajmil O, Templado C. Baseline expression profile of meiotic-specific genes in healthy fertile males. Fertil Steril
2009; 92: 578–82.
Rankin S. Complex elaboration: making sense of meiotic cohesin dynamics. FEBS J
2015; 282: 2426–43.
Garcia-Cruz R, Brieno MA, Roig I, Grossmann M, Velilla E, et al.
Dynamics of cohesin proteins REC8, STAG3, SMC1β and SMC3 are consistent with a role in sister chromatid cohesion during meiosis in human oocytes. Hum Reprod
2010; 25: 2316–27.
Fukuda T, Fukuda N, Agostinho A, Hernandez-Hernandez A, Kouznetsova A, et al.
STAG3-mediated stabilization of REC8 cohesin complexes promotes chromosome synapsis during meiosis. EMBO J
2014; 33: 1243–55.
Caburet S, Arboleda VA, Llano E, Overbeek PA, Barbero JL, et al.
Mutant cohesin in premature ovarian failure. N Engl J Med
2014; 370: 943–9.
He WB, Banerjee S, Meng LL, Du J, Gong F, et al.
Whole-exome sequencing identifies a homozygous donor splice-site mutation inSTAG3 that causes primary ovarian insufficiency. Clin Genet
2018; 93: 340–4.
Llano E, Gomez HL, Garcia-Tunon I, Sanchez-Martin M, Caburet S, et al.
STAG3 is a strong candidate gene for male infertility. Hum Mol Genet
2014; 23: 3421–31.
Bayes M, Prieto I, Noguchi J, Barbero JL, Perez Jurado LA. Evaluation of the Stage3 gene and the synaptonemal complex in a rat model (as/as) for male infertility. Mol Reprod Dev
2001; 60: 414–7.
Murray KS, James A, McGeady JB, Reed ML, Kuang WW, et al.
The effect of the new 2010 World Health Organization criteria for semen analyses on male infertility. Fertil Steril
2012; 98: 1428–31.
Hsueh HM, Chen JJ, Kodell RL. Comparison of methods for estimating the number of true null hypotheses in multiplicity testing. J Biopharm Stat
2003; 13: 675–89.
Gabriel SB, Schaffner SF, Nguyen H, Moore JM, Roy J, et al.
The structure of haplotype blocks in the human genome. Science
2002; 296: 2225–9.
Reynolds A, Qiao H, Yang Y, Chen JK, Jackson N, et al.
RNF212 is a dosage-sensitive regulator of crossing-over during mammalian meiosis. Nat Genet
2013; 45: 269–78.
Yang F, Eckardt S, Leu NA, McLaughlin KJ, Wang PJ. Mouse TEX15 is essential for DNA double-strand break repair and chromosomal synapsis during male meiosis. J Cell Biol
2008; 180: 673–9.
Fujiwara Y, Ogonuki N, Inoue K, Ogura A, Handel MA, et al.
t-SNARE Syntaxin2 (STX2) is implicated in intracellular transport of sulfoglycolipids during meiotic prophase in mouse spermatogenesis. Biol Reprod
2013; 88: 141.
Dickins RA, Frew IJ, House CM, O'Bryan MK, Holloway AJ, et al.
The ubiquitin ligase component Siahla is required for completion of meiosis I in male mice. Mol Cell Biol
2002; 22: 2294–303.
Kong A, Thorleifsson G, Stefansson H, Masson G, Helgason A, et al.
Sequence variants in the RNF212 gene associate with genome-wide recombination rate. Science
2008; 319: 1398–401.
Okutman O, Muller J, Baert Y, Serdarogullari M, Gultomruk M, et al.
Exome sequencing reveals a nonsense mutation in TEX15 causing spermatogenic failure in a Turkish family. Hum Mol Genet
2015; 24: 5581–8.
Choi Y, Jeon S, Choi M, Lee MH, Park M, et al.
Mutations in SOHLH1 gene associate with nonobstructive azoospermia. Hum Mutat
2010; 31: 788–93.
Hopkins J, Hwang G, Jacob J, Sapp N, Bedigian R, et al.
Meiosis-specific cohesin component, Stage3 is essential for maintaining centromere chromatid cohesion, and required for DNA repair and synapsis between homologous chromosomes. PLoS Genet
2014; 10: e1004413.
Yu CH, Xie T, Zhang RP, A ZC. Association of the common SNPs in RNF212, STAG3 and RFX2 gene with male infertility with azoospermia in Chinese population. Eur J Obstet Gynecol Reprod Biol
2018; 221: 109–12.
Bartel DP. MicroRNAs: genomics, biogenesis, mechanism, and function. Cell
2004; 116: 281–97.
Moreno-Moya JM, Vilella F, Simon C. MicroRNA: key gene expression regulators. Fertil Steril
2014; 101: 1516–23.
Kotaja N. MicroRNAs and spermatogenesis. Fertil Steril
2014; 101: 1552–62.
Fukuda T, Pratto F, Schimenti JC, Turner JM, Camerini-Otero RD, et al.
Phosphorylation of chromosome core components may serve as axis marks for the status of chromosomal events during mammalian meiosis. PLoS Genet
2012; 8: e1002485.
Stegeman S, Amankwah E, Klein K, O'Mara TA, Kim D, et al.
A large-scale analysis of genetic variants within putative miRNA binding sites in prostate cancer. Cancer Discov
2015; 5: 368–79.
Yousef GM. miRSNP-based approach identifies a miRNA that regulates prostate-specific antigen in an allele-specific manner. Cancer Discov
2015; 5: 351–2.
Chen J, Yao D, Li Y, Chen H, He C, et al.
Serum microRNA expression levels can predict lymph node metastasis in patients with early-stage cervical squamous cell carcinoma. Int J Mol Med
2013; 32: 557–67.
[Figure 1], [Figure 2]
[Table 1], [Table 2], [Table 3], [Table 4]