Miaofei Xu

Efficient typing of copy number variations in a segmental duplication-mediated rearrangement hotspot using multiplex competitive amplification

Local genomic architecture, such as segmental duplications (SDs), can induce copy number variations (CNVs) hotspots in the human genome, many of which manifest as genomic disorders. Significant technological advances have been achieved for genome-wide CNV investigations, but these costly methods are not suitable for genotyping certain disease-associated CNVs or other loci of interest in populations. Recently, two independent studies showed that the murine meiosis expressed gene 1 (Meig1) was critical to spermatogenesis. We found that the human orthologue MEIG1 is flanked by an SD pair, between which non-allelic homologous recombination (NAHR) can cause recurrent CNVs. To study this potential CNV hotspot and its role in spermatogenesis, we developed a new CNV genotyping method, AccuCopy, based on multiplex competitive amplification to investigate 320 patients with spermatogenic impairment and 93 healthy controls. Three MEIG1 duplications (two in patients and one in controls) were identified, whereas no deletion was found. As NAHR results in more recurrent deletions than duplications at a locus, the over representation of recurrent MEIG1 duplications suggests a potential purifying selection operating on this hotspot, possibly via fecundity. We also showed that AccuCopy is an efficient and reliable method for multiplex CNV genotyping.

Additional genomic duplications in AZFc underlie the b2/b3 deletion-associated risk of spermatogenic impairment in Han Chinese population

The azoospermia factor c (AZFc) region on the Y chromosome is a genetically dynamic locus in the human genome. Numerous genomic rearrangements, including deletion, duplication and inversion, have been identified in AZFc. The complete deletion of AZFc can cause spermatogenic impairment. However, the roles of partial AZFc deletions (e.g. b2/b3 deletion) in spermatogenesis are controversial and variable among human populations. Secondary duplication has been hypothesized to be a compensatory factor for partial AZFc deletions. To further study genomic duplications in AZFc as a potential genetic modifier underlying the phenotypic variations of partial AZFc deletions in spermatogenesis, we conducted comprehensive molecular analyses in 711 idiopathic infertile men and 390 healthy controls. Unexpectedly, we found that additional AZFc duplications accompanying the b2/b3 deletion, instead of the b2/b3 deletion alone, led to the b2/b3 deletion-associated risk of spermatogenic impairment previously reported in Han Chinese population. In addition, partial AZFc duplication also rendered a risk factor in the non-deletion patients. DAZ is a multi-copy AZFc gene (DAZ1-DAZ4) implicated in spermatogenesis. Genetic variations do exist between DAZ copies. Intriguingly, we found that the DAZ1/2 cluster was the main duplicated copies in the partial AZFc duplications associated with spermatogenic impairment, suggesting a potential different role of spermatogenesis between DAZ copies. Our findings demonstrated that additional AZFc duplications did not compensate but convey the susceptibility of the b2/b3 deletion to spermatogenic impairment in the tested population. Notably, genomic duplications and deletions in AZFc deserve comprehensive investigations to uncover spermatogenic roles of the AZFc region.

Mitochondria-related miR-151a-5p reduces cellular ATP production by targeting CYTB in asthenozoospermia

Mitochondria, acting as the energy metabolism factory, participate in many key biological processes, including the maintenance of sperm viability. Mitochondria-related microRNA (miRNA), encoded by nuclear genome or mitochondrial genome, may play an important regulatory role in the control of mitochondrial function. To investigate the potential role of mitochondria-related miRNAs in asthenozoospermia, we adopted a strategy consisting of initial screening by TaqMan Low Density Array (TLDA) and further validation with quantitative reverse transcriptase polymerase chain reaction (qRT-PCR). Validation of the profiling results was conducted in two independent phases. Eventually, two seminal plasma miRNAs (sp-miRs) (miR-101-3p, let-7b-5p) were found to be significantly decreased, while sp-miR-151a-5p was significantly increased in severe asthenozoospermia cases compared with healthy controls. To further study their potential roles in asthenozoospermia, we then evaluated mitochondrial function of GC-2 cells transfected with these potentially functional miRNAs. Our results demonstrated that transfection with miR-151a-5p mimics decreased the mitochondrial respiratory activity. Besides, Adenosine Triphosphate (ATP) level was decreased when transfected with miR-151a-5p mimics. In addition, Cytochrome b (Cytb) mRNA and protein levels were also decreased when miR-151a-5p was overexpressed. These results indicate that miR-151a-5p may participate in the regulation of cellular respiration and ATP production through targeting Cytb.

Evaluation of Five Candidate Genes from GWAS for Association with Oligozoospermia in a Han Chinese Population

Background Oligozoospermia is one of the severe forms of idiopathic male infertility. However, its pathology is largely unknown, and few genetic factors have been defined. Our previous genome-wide association study (GWAS) has identified four risk loci for non-obstructive azoospermia (NOA). Objective To investigate the potentially functional genetic variants (including not only common variants, but also less-common and rare variants) of these loci on spermatogenic impairment, especially oligozoospermia. Design, Setting, and Participants A total of 784 individuals with oligozoospermia and 592 healthy controls were recruited to this study from March 2004 and January 2011. Measurements We conducted a two-stage study to explore the association between oligozoospermia and new makers near NOA risk loci. In the first stage, we used next generation sequencing (NGS) in 96 oligozoospermia cases and 96 healthy controls to screen oligozoospermia-susceptible genetic variants. Next, we validated these variants in a large cohort containing 688 cases and 496 controls by SNPscan for high-throughput Single Nucleotide Polymorphism (SNP) genotyping. Results and Limitations Totally, we observed seven oligozoospermia associated variants (rs3791185 and rs2232015 in PRMT6, rs146039840 and rs11046992 in Sox5, rs1129332 in PEX10, rs3197744 in SIRPA, rs1048055 in SIRPG) in the first stage. In the validation stage, rs3197744 in SIRPA and rs11046992 in Sox5 were associated with increased risk of oligozoospermia with an odds ratio (OR) of 4.62 (P  =  0.005, 95%CI 1.58-13.4) and 1.82 (P  =  0.005, 95%CI 1.01-1.64), respectively. Further investigation in larger populations and functional characterizations are needed to validate our findings. Conclusions Our study provides evidence of independent oligozoospermia risk alleles driven by variants in the potentially functional regions of genes discovered by GWAS. Our findings suggest that integrating sequence data with large-scale genotyping will serve as an effective strategy for discovering risk alleles in the future.

Mitochondria-related miR-141-3p contributes to mitochondrial dysfunction in HFD-induced obesity by inhibiting PTEN

Mitochondria-related microRNAs (miRNAs) have recently emerged as key regulators of cell metabolism and can modulate mitochondrial fusion and division. In order to investigate the roles of mitochondria-related miRNAs played in obesity, we conducted comprehensive molecular analysis in vitro and in vivo. Based on high-fat-diet (HFD) induced obese mice, we found that hepatic mitochondrial function was markedly altered. Subsequently, we evaluated the expression levels of selected mitochondria-related miRNAs and found that miR-141-3p was up-regulated strikingly in HFD mice. To further verify the role of miR-141-3p in obesity, we carried out gain-and-loss-of-function study in human HepG2 cells. We found that miR-141-3p could modulate ATP production and induce oxidative stress. Through luciferase report gene assay, we identified that phosphatase and tensin homolog (PTEN) was a target of miR-141-3p. Inhibiting PTEN could alter the mitochondrial function, too. Our study suggested that mitochondria-related miR-141-3p induced mitochondrial dysfunction by inhibiting PTEN.

Gene copy number alterations in the azoospermia-associated AZFc region and their effect on spermatogenic impairment

The azoospermia factor c (AZFc) region in the long arm of human Y chromosome is characterized by massive palindromes. It harbors eight multi-copy gene families that are expressed exclusively or predominantly in testis. To assess systematically the role of the AZFc region and these eight gene families in spermatogenesis, we conducted a comprehensive molecular analysis (including Y chromosome haplogrouping, AZFc deletion typing and gene copy quantification) in 654 idiopathic infertile men and 781 healthy controls in a Han Chinese population. The b2/b3 partial deletion (including both deletion-only and deletion-duplication) was consistently associated with spermatogenic impairment. In the subjects without partial AZFc deletions, a notable finding was that the frequency of DAZ and/or BPY2 copy number alterations in the infertile group was significantly higher than in the controls. Combined patterns of DAZ and/or BPY2 copy number abnormality were associated with spermatogenic impairment when compared with the pattern of all AZFc genes with common level copies. In addition, in Y chromosome haplogroup O1 (Y-hg O1), the frequency of copy number alterations of all eight gene families was significantly higher in the case group than that in the control group. Our findings indicate that the DAZ, BPY2 genes may be prominent players in spermatogenesis, and genomic rearrangements may be enriched in individuals belonging to Y-hg O1. Our findings emphasize the necessity of routine molecular analysis of AZFc structural variation during the workup of azoospermia and/or oligozoospermia, which may diminish the genetic risk of assisted reproduction.

Pathogenic variants screening in five non-obstructive azoospermia-associated genes

Non-obstructive azoospermia (NOA) is one of the most severe forms of male infertility and a recent, genome-wide association study (GWAS) has identified four risk loci associated with NOA. However, a large portion of the heritability of NOA has not been well explained by GWAS. By hypothesizing that rare, low-frequency and common genetic variants might point toward a causal relation between candidate genes and NOA, we performed a two-stage study including deep exon sequencing in 96 NOA cases and 96 healthy controls and a replication study in a larger population containing 522 NOA cases and 484 healthy controls. In the solexa sequencing stage, a total of two rare mutations (chr20. 1902132 and chr20. 1902301 in SIRPA), four common mutations (rs1048055 and rs2281807 in SIRPG, rs11046992 and rs146039840 in SOX5) were identified by using next generation sequencing (NGS). In the validation stage, subjects in the NOA group had a significantly decreased frequency of the heterozygous GA genotype in SIRPA (4.23%, 22 out of 520) than that in the control group (8.60%, 41 out of 477) [odds ratios (OR) 0.47, 95% confidence intervals (CI) 0.28-0.80] (P = 6.00 × 10(-3)). The rs1048055 in SIRPG was associated with a significantly increased risk of spermatogenic impairment, compared with the CC genotype (OR 3.93, 95% CI 1.59-9.70) (P = 3.00 × 10(-3)). Our study provides evidence of independent NOA risk alleles driven by variants in the protein-coding sequence of two of the genes (SIRPA and SIRPG) discovered by GWAS. Further investigation in larger populations and functional characterizations are needed to validate our findings.

DAZ duplications confer the predisposition of Y chromosome haplogroup K* to non-obstructive azoospermia in Han Chinese populations

STUDY QUESTION What are the genetic causes for the predisposition of certain Y chromosome haplogroups (Y-hgs) to spermatogenic impairment? SUMMARY ANSWER The AZFc(azoospermia factor c)/DAZ (deleted in azoospermia) duplications might underlie the susceptibility of Y-hg K* to spermatogenic impairment. WHAT IS KNOWN ALREADY The roles of Y chromosomal genetic background in spermatogenesis are controversial and vary among human populations. Individuals in predisposed Y-hgs may carry some genetic factors, which might be a potential genetic modifier for the Y-hg-specific susceptibility to spermatogenic impairment. STUDY DESIGN, SIZE, DURATION A total of 2444 individuals with azoospermia or oligozoospermia and 2456 healthy controls were recruited to this study from March 2004 and January 2011. PARTICIPANTS/MATERIALS, SETTING, METHODS We performed a two-stage association study to investigate the risk and/or protective Y-hgs for spermatogenic impairment. In addition, the genetic causes for the predisposition of certain Y-hg to spermatogenic impairment were investigated. Deletion typing and DAZ gene copy number quantification were performed for individuals in predisposed Y-hgs. MAIN RESULTS AND THE ROLE OF CHANCE Y-hgs K* and O3e* showed significantly different distribution between cases and controls consistently in two-stage studies. Combined analyses identified significant predisposition to non-obstructive azoospermia in Y-hg K* [odds ratio (OR) 8.58; 95% confidence interval (CI) 3.31-22.28; P = 1.40 × 10⁻⁵], but a protecting effect in Y-hg O3e* (OR 0.64; 95% CI 0.53-0.78; P = 4.20 × 10⁻⁵). Based on the dynamic nature of the Y chromosome, we hypothesized that Y-hgs K* and O3e* may be accompanied by modifying genetic factors for their predisposing or protecting effects in spermatogenesis. Accordingly, we quantified the multi-copy DAZ gene, which has variable copy numbers between individuals and plays an important role in spermatogenesis. In combined analysis, we found that the over-dosage of DAZ was significantly more frequent in Y-hg K* than in O3e* (OR 4.79; 95% CI 1.67-13.70; P = 6 × 10⁻³). LIMITATIONS, REASONS FOR CAUTION Owing to the inconsistency of genetic background, it remains to be determined whether the results derived from Han Chinese populations are applicable to other ethnic groups. WIDER IMPLICATIONS OF THE FINDINGS The findings of this study can advance the etiology of spermatogenic impairment, and also shed new light on Y chromosome evolution in human populations. Y-hg-specific genetic factors of modifying spermatogenic phenotypes deserve further investigation in larger and diverse populations.

Genetic variants in meiotic program initiation pathway genes are associated with spermatogenic impairment in a Han Chinese population.

Background The meiotic program initiation pathway genes (CYP26B1, NANOS1 and STRA8) have been proposed to play key roles in spermatogenesis. Objective To elucidate the exact role of the genetic variants of the meiosis initiation genes in spermatogenesis, we genotyped the potential functional genetic variants of CYP26B1, NANOS1 and STRA8 genes, and evaluated their effects on spermatogenesis in our study population. Design, Setting, and Participants In this study, all subjects were volunteers from the affiliated hospitals of Nanjing Medical University between March 2004 and July 2009 (NJMU Infertile Study). Total 719 idiopathic infertile cases were recruited and divided into three groups according to WHO semen parameters: 201 azoospermia patients (no sperm in the ejaculate even after centrifugation), 155 oligozoospermia patients (sperm counts <20×106/ml) and 363 infertility/normozoospermia subjects (sperm counts >20×106/ml). The control group consisted of 383 subjects with normal semen parameters, all of which had fathered at least one child without assisted reproductive technologies. Measurements Eight single nucleotide polymorphisms (SNPs) in CYP26B1, NANOS1 and STRA8 genes were determined by TaqMan allelic discrimination assay in 719 idiopathic infertile men and 383 healthy controls. Results and Limitations The genetic variant rs10269148 of STRA8 gene showed higher risk of spermatogenic impairment in the groups of abnormospermia (including azoospermia subgroup and oligozoospermia subgroup) and azoospermia than the controls with odds ratios and 95% confidence intervals of 2.52 (1.29–4.94) and 2.92 (1.41–6.06), respectively (P = 0.006, 0.002 respective). Notably, larger sample size studies and in vivo or in vitro functional studies are needed to substantiate the biological roles of these variants. Conclusions Our results provided epidemiological evidence supporting the involvement of genetic polymorphisms of the meiotic program initiation genes in modifying the risk of azoospermia and oligozoospermia in a Han-Chinese population.

A genome-wide association study of mitochondrial DNA in Chinese men identifies two risk single nucleotide substitutions for idiopathic oligoasthenospermia.

Mitochondrial DNA (mtDNA) is believed to be both the source and target of reactive oxygen species (ROS), and mtDNA genetic alterations have been reported to be associated with molecular defects in the oxidative phosphorylation (OXPHOS) system. In order to investigate the potentially susceptible mtDNA genetic variants to oligoasthenospermia, we conducted a two-stage study in 921 idiopathic infertile men with oligoasthenospermia and 766 healthy controls using comprehensive molecular analysis. In the screen stage, we used next generation sequencing (NGS) in 233 cases and 233 controls to screen oligoasthenospermia susceptible mitochondrial genetic variants. In total, seven variants (C5601T, T12338C, A12361G, G13928C, A15235G, C16179T and G16291A) were screened to be potentially associated with idiopathic oligoasthenospermia. In the validation stage, we replicated these variants in 688 cases and 533 healthy controls using SNPscan. Our results demonstrated that the genetic alteration of C16179T was associated with idiopathic male infertility (odds ratio (OR) 3.10, 95% CI 1.41-6.79) (p=3.10×10(-3)). To elucidate the exact role of the genetic variants in spermatogenesis, two main sperm parameters (sperm count and motility) were taken into account. We found that C16179T was associated with both low sperm count and motility, with ORs of 4.18 (95% CI 1.86-9.40) (p=1.90×10(-4)) and 3.17 (95% CI 1.40-7.16) (p=3.50×10(-3)), respectively. Additionally, A12361G was found to be associated with low sperm count, with an OR of 3.30 (95% CI 1.36-8.04) (p=5.50×10(-3)). These results indicated that C16179T influenced both the process of spermatogenesis and sperm motility, while A12361G may just only participate in the process of spermatogenesis. Further investigation in larger populations and functional characterizations are needed to validate our findings.