| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||




From the * Centre for Cellular and Molecular
Biology, Hyderabad, India; the
University of
Madras, Chennai, India; the
Institute of
Reproductive Medicine, Salt Lake, Kolkata, India; and the
Infertility Institute and Research Centre,
Hyderabad, India.
| Correspondence to: K. Thangaraj, Centre for Cellular and Molecular Biology, Uppal Road, Hyderabad-500 007, India (e-mail: thangs{at}ccmb.res.in). |
| Received for publication November 26, 2007; accepted for publication May 18, 2008. |
| Abstract |
|---|
|
|
|---|
Key words: Haplotype, SNPs, ubiquitin, ubiquitin-conjugating enzyme
Ubiquitin (UB) is a small and highly conserved protein, present apparently in all eukaryotic cells. It is covalently attached to abnormal and short-lived proteins that are destined for ATP-dependent proteolysis in eukaryotic cells. UB is first activated by the ubiquitin-activating enzyme E1. Ubiquitin-conjugating enzyme E2 (ubiquitin carrier protein) transfers the activated ubiquitin to the substrate bound to E3, a member of the ubiquitin ligase family (Ciechanover, 1996). This ubiquitin pathway, by which proteins are tagged for degradation and ubiquitin-mediated proteolysis, plays a major role in a broad array of basic cellular processes (Tongaonkar et al, 2000). Among these are regulation of the cell cycle, differentiation and development, biogenesis of organelles, DNA repair, modulation of cell surface receptors, ion channels, signaling pathways, stress responses, and removal of damaged or misfolded proteins (Hochstrasser, 2001). Several investigations have clearly indicated the correlation between the UB pathway of protein degradation and spermatogenesis (Roest et al, 1996; Bebington et al, 2001).
Ube2b knockout mice were found to be infertile because of abnormal sperm morphology and immotility (Roest et al, 1996). Derailment of spermatogenesis became overt during the postmeiotic condensation of chromatin in spermatids. Roest et al (1996) stated that the findings in the knockout mice provided a parallel between yeast sporulation and mammalian spermatogenesis and strongly implicated UBE2B-dependent ubiquitination in chromatin remodeling in the human.
In another study, Baarends et al (2003) determined that the primary spermatocytes of Ube2b knockout mice underwent increased apoptosis during meiotic prophase. In the absence of UBE2B, the structure and telomere localization of synaptonemal complexes were altered within the nuclei of pachytene and diplotene spermatocytes. The number of foci containing the mismatch DNA repair protein MLH1 was increased and reflected a consistent 20%–25% increase in crossover frequency in mutant spermatocytes. The authors concluded that the UB-conjugating activity of UBE2B is required within the synaptonemal complex and for meiotic recombination in spermatocytes. Thus, it is now evident that Ube2b knockout mice are infertile because of abnormal sperm morphology and motility. Therefore, we undertook this study to assess the role of the UBE2B gene in infertile Indian men because there has been no report on this gene before this study.
| Materials and Methods |
|---|
|
|
|---|
|
The cases were classified as having no sperm in the ejaculate (ie, azoospermia, n = 350; and low sperm count [oligozoospermia, sperm count < 20 million/mL], n = 180). Of the latter category, 75 patients presented with <50% motile sperm (oligoasthenozoospermic) and 105 individuals with low motility and abnormal sperm shape and size (<30% sperm with normal morphology), known as oligoasthenoteratozoospermia. Microdeletion analysis of the Y chromosome was carried out for all the cases (Thangaraj et al, 2003a, unpublished data). Patients who did not exhibit obstruction, endocrinological defect, pelvic injury, major illness, karyotype abnormality, or Y chromosome microdeletion were included in the study. The age of infertile men ranged between 26 and 42 years. Three hundred randomly selected and ethnically matched fertile men (sperm count > 20 x 106/mL) were included in this study as controls (Table 1). All the 300 control men had fathered at least one child each, and their fertility was proven by short tandem repeat (STR)–based (Profiler plus; Applied Biosystems, Foster City, California) DNA fingerprinting. Blood samples (5 mL) from each infertile and fertile control man were collected with their informed written consent. Genomic DNA was isolated using the protocol described in our earlier study (Thangaraj et al, 2002). This study was approved by The Institutional Ethical Committee.
Analysis of the UBE2B Gene![]()
The DNA sequence of the UBE2B gene (accession number
ENSG00000119048) was downloaded from the Ensembl database
(www.ensembl.org),
and 6 pairs of primers covering the entire coding region and the exon-intron
boundaries and 2 internal primers (forward) for sequencing were designed with
GeneTool software (Table 2).
The primers were synthesized with the use of an ABI392 oligosynthesizer
(Applied Biosystems, Foster City, California). PCRs were performed according
to the following amplification conditions: 94°C for 7 minutes; 35 cycles
of denaturation at 94°C for 45 seconds, annealing at 60.5°C or
65.5°C (Table 2) for 30
seconds, extension at 72°C for 1 minute, and a final extension at 72°C
for 7 minutes (Table 1). PCR
products were checked on 2% agarose gel, stained with ethidium bromide, and
visualized under an ultraviolet transilluminator. The amplicons were further
subjected to direct sequencing (Thangaraj
et al, 2002) with a BigDye Terminator cycle sequencing kit
(version 3.1, Applied Biosystems) and 3730 DNA analyzer (Applied
Biosystems).
|
Statistical Analysis![]()
PEPI (Programs for Epidemiologic Analysis) CASECONT
(Abramson and Gahlinger, 2001)
was used to perform statistical analysis to determine whether any
single-nucleotide polymorphism (SNP) identified in this study is associated
with male infertility. The odds ratio (OR) and 95% confidence interval (CI)
were calculated to measure the strength of association between genotype
frequencies and infertility. P values <.05 were considered
statistically significant.
Haplotype Analysis![]()
Haplotypes were inferred from unphased genotypes with the program PHASE
version 2.1 (Stephens et al,
2001; Stephens and Scheet,
2005), which was run with default settings (100 iterations, with
100 burn-in, and a thinning interval of 1). A case control association
analysis was performed by PHASE v2.1 to test the hypothesis that the
haplotypes were randomly distributed between the infertile and fertile men.
2 tests were used to test individual haplotypes for their
over- or underrepresentation in infertile vs fertile men.
In Silico Analysis![]()
To evaluate whether SNPs present exclusively in infertile men and are the
potential cause for the defect in recognition/binding sites for the proteins
responsible for proper splicing, we have analyzed those sites with ESEfinder,
release 3.0
(http://rulai.cshl.edu/cgi-bin/tools/ESE3/esefinder.cgi),
which evaluates the presence of exonic splicing enhancer (ESE) at the mutation
sites. The ASD Workbench wrapper
(http://www.ebi.ac.uk/asd-srv/wb.cgi)
tools PPT analysis, BP analysis, Splicing Rainbow, and Regulatory Sequences
were used to identify the presence of polypyrimidine tracks, branch point
sites, binding sites for splicing factors, and exonic splicing
enhancers/silencers (ESE/ESS) or intronic splicing enhancers/silencers
(ISE/ISS), respectively, at the mutation sites. ESE-finder 3.0 searches the
ESE motifs corresponding to 4 serine/arginine-rich (SR) proteins, wherein a
Splicing Rainbow tool searches for the SR proteins as well as hnRNP
motifs.
| Results |
|---|
|
|
|---|
|
|
The transversion g.-293:T>G (rs17167484) was observed in 37 infertile
(6.2%) and 6 fertile men (2%) (panel A in the
Figure;
Table 3). The difference in the
frequency distribution of this SNP between infertile and fertile men was
statistically significant (
2 = 9.676, P = .002; OR =
3.68, CI = 1.46–9.81). Another substitution, g.20016:A>G (rs3777373)
was observed in 137 infertile (23%) and 16 (5.3%) fertile men (panel D in the
Figure;
Table 3), with a statistically
significant difference between the 2 groups (
2 = 53.625,
P < .001; OR = 6.19, CI = 3.52–11.04).
Haplotype Analysis![]()
Analyses of the case and control genotypes with the use of the
g.-293:T>G and g.20016:A>G sites (PHASE v2.1) revealed 4 haplotypes.
Distribution of haplotypes was not uniform between infertile and fertile men.
Furthermore, in fertile men frequency of the TA haplotype was high (95.83%)
when compared with infertile men (82.24%), suggesting that this haplotype is
advantageous for fertility. Our results suggest that haplotype TG conferred
significantly increased risk for male infertility (OR = 5.07, 95% CI =
1.29–23.29, P = .007) if the main haplotype TA was used as the
baseline genotype (Table
4).
|
In Silico Analysis![]()
To evaluate the functional significance of the SNPs, which were exclusively
present in infertile men, these intronic and exonic sites, along with their
flanking sequence, were analyzed by various bioinformatics tools. Analysis of
g.5197:T>G with the Splicing Rainbow tool predicted a binding site for
SRp55 in both normal and mutated sequences, but the potential binding site for
SC35, which was present in the normal sequence, was absent in the mutated
sequence (Table 5). This
analysis also predicted the binding sites for SRp40 and hnRNP U in the mutated
sequence, but not in the normal sequence. A similar analysis performed for the
exon 4 SNP g.9157:A>G showed a potential binding site for the SPp40, SC35,
and ASF/SF2 in both normal and mutated sequences, but the binding site for
Tra2Beta was predicted in normal, but not the mutant, sequence. We also
analyzed intronic SNPs with the PPT and BP analysis tools, but none of the SNP
were present in the branching point site or caused any change in the predicted
polypyrimidine track.
|
| Discussion |
|---|
|
|
|---|
We identified 2 novel SNPs (g.5197:T>G, g.9157:A>G) observed only in infertile men and 3 known SNPs (g.-293:T>G [rs17167484], g.19765:T>G [rs11538104], and g.20016:A>G [rs3777373]), of which the g.19765:T>G substitution was found only in azoospermic men, whereas the other 2 were found in both infertile and fertile men. However, the differences in the frequency of these mutations (g.-293:T>G and g.20016:A>G) between cases and control groups were statistically significant (Table 3). Because the Ube2b knockout mouse has been shown to have maturation arrest, we tried to correlate the histological phenotype in 18 azoospermic men (for which the biopsy results were available) with the mutations observed; however, we did not find any correlation. One possible reason could be availability of fewer samples for histological studies. The genotype frequencies of the known SNPs (g.-293:T>G, g.19765:T>G, and g.20016:A>G) in various populations, including this study, are shown in Table 6. The g.19765:T>G transversion was found only in azoospermic samples (1.71%) in this study; hence, we strongly suggest that this SNP might be associated with spermatogenic failure in Indian populations. Considering the unique ethnic origin of Indian populations, it is not surprising. India is known for its human diversity, having 4635 anthropologically well defined populations. Each population follows strict endogamy, social and cultural traditions; therefore, we presume that these SNPs might have originated in situ in Indians. Similar findings are quite common between different ethnic populations. We have earlier demonstrated that the A386G in the DAZL gene is not associated with male infertility on the Indian subcontinent (Thangaraj et al, 2006), whereas Teng et al (2002) earlier reported that this mutation was associated with severe spermatogenic failure in the Taiwan population. Similarly, expansion of the CAG repeat length was reported to be associated with male infertility (Patrizio et al, 2001), whereas, we did not find any significant difference in the CAG repeat length between infertile and fertile men (Thangaraj et al, 2002). This trend is not only observed in the genes related to infertility, but it also true for other genes involved in sex determination (Singh et al, 2006), alcoholism (Bhaskar et al, 2007; Rao et al, 2007), and mitochondrial disorders (Rani et al, 2006; Vanniarajan et al, 2006). Interestingly, g.9157:A>G and g.19765: T>G polymorphisms were present only in azoospermic patients, which suggests that these SNPs are specifically associated with azoospermia, but not with other phenotypes of infertility (Table 3).
|
In this investigation, 4 haplotypes were observed based on 2 SNPs (Table 4). Distribution of haplotypes was not uniform between the infertile and fertile men. The high frequency of TA haplotype in fertile men when compared with infertile men suggests that this haplotype is advantageous for fertility, whereas the TG haplotype confers susceptibility to spermatogenic failure. The remaining 2 haplotypes (GA and GG) were observed at low frequencies in both fertile and infertile men. The absence of GG haplotype in oligoasthenozoospermic men could be attributed to their small sample size. These findings indicate that substitution of a single variant in a haplotype could significantly alter gene function. Although both of the SNPs used for the haplotype analysis were present in fertile men as well, their frequency was significantly low. We observed similar findings in other genes; for example, a total of 5 haplotypes were observed for the NPY gene, only 2 of which were found to occur with a high frequency in all populations (Bhaskar et al, 2007).
Gene expression profile depends on the proper transcription, which includes DNA to mRNA formation; pre-mRNA splicing and specific consensus sequences (5' splice site, branch site, and 3' donor site) especially are found in virtually all exon-intron boundaries (Sun and Chasin, 2000). Among the splicing factors involved in splice site choice, members of SR protein family have been studied extensively. These proteins are rich in serine and arginine residues and have an RNA recognition motif domain that recognizes several classes of specific RNA motifs, including ESEs and ISEs. The hnRNPs are another class of proteins that recognize negative regulatory elements known as ESSs and ISSs. These sequences have been demonstrated to play key roles in both alternate and constitutive splice site selection (Gabut et al, 2005). Our in silico analysis of intronic and exonic SNPs specifically present in infertile men predict that these SNPs affect the binding site of splicing factors (SR and hnRNP proteins), which will in turn affect the proper splicing and translation of protein. Therefore, we predict that these SNPs bring about either loss or appearance of potential binding sites for various splicing factors.
This is for the first time we report SNPs and haplotypes in the human UBE2B gene and its role in human male infertility. We detected SNPs in all categories of infertile men (azoospermic, oligoasthenoteratozoospermia, and oligoasthenozoospermic) in the Indian population. From this study, it is clear that disturbance in the UBE2B gene not only causes structural abnormalities of spermatozoa, which is seen in knockout mice (Roest et al, 1996), but also could be responsible for sperm production in humans. Therefore, our study predicted the additional genetic causes of male infertility.
| Acknowledgments |
|---|
| Footnotes |
|---|
| References |
|---|
|
|
|---|
Baarends WM, Wassenaar E, Hoogerbrugge JW, Van Cappellen G, Roest
HP, Vreeburg J, Ooms M, Hoeijmakers JH, Grootegoed JA. Loss of HR6B ubiquitin
conjugating activity results in damaged synaptonemal complex structure and
increased crossing over frequency during the male meiotic prophase.
Mol Cell Biol. 2003; 23: 1151
–1162.
Bebington C, Doherty FJ, Fleming SD. The possible biological and
reproductive functions of ubiquitin. Hum Reprod
Update. 2001;7: 102
–111.
Bhaskar LVKS, Thangaraj K, Shah AM, Pardhasaradhi G, Kumar KP, Reddy AG, Papa Rao A, Mulligan CJ, Singh L, Rao VR. Allelic variation in the NPY gene in 14 Indian populations. J Hum Genet. 2007; 52: 592 –598.[CrossRef][Medline]
Ciechanover A. Ubiquitin-mediated proteolysis and male sterility. Nat Med . 1996;2: 1188 –1190.[CrossRef][Medline]
Gabut M, Mine M, Marsac C, Brivet M, Tazi J, Soret J. The SR
protein SC35 is responsible for aberrant splicing of the E1alpha pyruvate
dehydrogenase mRNA in a case of mental retardation with lactic acidosis.
Mol Cell Biol. 2005; 25: 3286
–3294.
Grootegoed JA, Baarends WM, Roest HP, Hoeijmakers JH. Knockout mouse model and gametogenic failure. Mol Cell Endocrinol. 1998;145: 161 –166.[CrossRef][Medline]
Hackstein JHP, Hochstenbach R, Pearson PL. Towards on understanding of the genetics of human male infertility: lessons from flies. Trends Genet. 2000; 16: 565 –572.[CrossRef][Medline]
Hochstrasser M. Spring for sumo. New functions bloom for a ubiquitin-like protein. Cell. 2001; 107: 5 –8.[CrossRef][Medline]
International HapMap Project. http://www.hapmap.org.
Krausz C, Forti G. Clinical aspects of male infertility. Results Probl Cell Differ . 2000; 28: 1 –21.[Medline]
Patrizio P, Leonard DGB, Chen KL, Hernandez-Ayup S, Trounson AO. Larger trinucleotide repeat size in the androgen receptor gene of infertile men with extremely severe oligozoospermia. J Androl. 2001; 22: 444 –448.[Abstract]
Rani DS, Vanniarajan A, Gupta NJ, Chakravarty B, Singh L, Thangaraj K. A novel missense mutation C11994T in the mitochondrial ND4 gene as a cause of low sperm motility in the Indian sub continent. Fertil Steril. 2006;86: 1783 –1785.[CrossRef][Medline]
Rao VR, Bhaskar LV, Annapurna C, Reddy AG, Thangaraj K, Papa Rao A, Singh L. Single nucleotide polymorphisms in alcohol dehydrogenase genes among some Indian populations. Am J Hum Biol. 2007; 19: 338 –344.[CrossRef][Medline]
Reynolds N, Cooke HJ. Role of the DAZ genes in male fertility. Reprod Biomed Online. 2005; 10: 72 –80.[Medline]
Roest HP, van Klaveren J, de Wit J, van Gurp CG, Koken MH, Vermey M, van Roijen JH, Hoogerbrugge JW, Vreeburg JT, Baarends WM, Bootsma D, Grootegoed JA, Hoeijmakers JH. Inactivation of the HR6B ubiquitin-conjugating DNA repair enzyme in mice causes male sterility associated with chromatin modification. Cell. 1996; 86: 799 –810.[CrossRef][Medline]
Shah K, Sivapalan G, Gibbons N, Tempest H, Griffin DK. The genetic basis of infertility. Reproduction. 2003; 126: 13 –25.[Abstract]
Singh R, Thangaraj K, Gupta NJ, Leelavathy N, Rani DS, Nambiar RG,
Kalavathy V, Santhiya ST, Rajangam S, Gopinath PM, Chakravarty B, Singh L. A
novel human sex-determining gene linked to Xp11.21–11.23. J
Clin Endocrinol Metab. 2006; 91: 4028
–4036.
Stephens M, Scheet P. Accounting for decay of linkage disequilibrium in haplotype inference and missing-data imputation. Am Hum Genet. 2005; 76: 449 –462.[CrossRef]
Stephens M, Smith NJ, Donnelly P. A new statistical method for haplotype reconstruction from population data. Am J Hum Genet. 2001;68: 978 –989.[CrossRef][Medline]
Sun H, Chasin LA. Multiple splicing defects in an intronic false
exon. Mol Cell Biol. 2000; 20: 6414
–6425.
Teng YN, Lin YM, Lin YH, Tsao SY, Hsu CC, Lin SJ, Tsai WC, Kuo PL.
Association of a single-nucleotide polymorphism of the
deleted-in-azoospermia-like gene with susceptibility to spermatogenic failure.
J Clin Endocrinol Metab. 2002; 87: 5258
–5264.
Thangaraj K, Chaubey G, Kivisild T, Reddy AG, Singh VK, Rasalkar
AA, Singh L. Reconstructing the origin of Andaman Islanders.
Science. 2005;308: 996
.
Thangaraj K, Deepa SR, Pavani K, Gupta NJ, Reddy P, Reddy AG, Chakravarty BN, Singh L. A to G transitions at 260, 386 and 437 in DAZL gene are not associated with spermatogenic failure in Indian population. Int J Androl. 2006; 29: 510 –514.[CrossRef][Medline]
Thangaraj K, Gupta NJ, Pavani K, Reddy AG, Subramanian S, Rani DS,
Ghosh B, Chakravarty B, Singh L. Y chromosome deletion in azoospermic men in
India. J Androl. 2003a; 24: 588
–597.
Thangaraj K, Joshi MB, Reddy AG, Gupta NJ, Chakravarty B, Singh L.
CAG repeat expansion in the androgen receptor gene is not associated with male
infertility in Indian populations. J Androl. 2002; 23: 815
–818.
Thangaraj K, Joshi MB, Reddy AG, Rasalkar AA, Singh L. Sperm
mitochondrial mutations as a cause of low sperm motility. J
Androl. 2003b;24: 388
–392.
Tongaonkar P, Chen L, Lambertson D, Ko B, Madura K. Evidence for an
interaction between ubiquitin-conjugating enzymes and the 26S proteasome.
Mol Cell Biol. 2000; 20: 4691
–4698.
Vanniarajan A, Rajshekher GP, Joshi MB, Reddy AG, Singh L, Thangaraj K. Novel mitochondrial mutation in the ND4 gene associated with Leigh syndrome. Acta Neurol Scand. 2006; 114: 350 –353.[CrossRef][Medline]
Venables JP, Cook HJ. Lessons from knockout and transgenic mice for infertility in men. J Endocrinol Invest. 2000; 23: 584 –591.[Medline]
Yang XJ, Shinka T, Nozawa S, Yan HT, Yoshiike M, Umeno M, Sato Y,
Chen G, Iwamoto T, Nakahori Y. Survey of the 2 polymorphisms in DAZL,
an autosomal candidate for the azoospermic factor, in Japanese infertile men
and implications for male infertility. Mol Hum Reprod. 2005; 11: 513
–515.
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |