|Ahead of print publication
Association between diacylglycerol kinase kappa variants and hypospadias susceptibility in a Han Chinese population
Hua Xie1, Xiao-Ling Lin2, Song Zhang3, Ling Yu1, Xiao-Xi Li1, Yi-Chen Huang1, Yi-Qing Lyu1, Hai-Tao Chen4, Jianfeng Xu2,4,5, Fang Chen1
1 Department of Urology, Shanghai Children's Hospital, Shanghai Jiao Tong University, Shanghai 200062, China
2 Fudan Institute of Urology, Huashan Hospital, Fudan University, Shanghai 200040, China
3 State Key Laboratory of Genetic Engineering, School of Life Sciences, Fudan University, Shanghai 200433, China
4 Center for Genomic Translational Medicine and Prevention, School of Public Health, Fudan University, Shanghai 200433, China
5 Program for Personalized Cancer Care, NorthShore University HealthSystem, Evanston, IL 60201, USA
|Date of Submission||11-Dec-2016|
|Date of Acceptance||14-Mar-2017|
|Date of Web Publication||09-Jun-2017|
Department of Urology, Shanghai Children's Hospital, Shanghai Jiao Tong University, Shanghai 200062, China
Fudan Institute of Urology, Huashan Hospital, Fudan University, Shanghai 200040, China; Center for Genomic Translational Medicine and Prevention, School of Public Health, Fudan University, Shanghai 200433, China; Program for Personalized Cancer Care, NorthShore University HealthSystem, Evanston, IL 60201, USA
Source of Support: None, Conflict of Interest: None
Previous genome-wide association studies have identified variants in the diacylglycerol kinase kappa (DGKK) gene associated with hypospadias in populations of European descent. However, no variants of DGKK were confirmed to be associated with hypospadias in a recent Han Chinese study population, likely due to the limited number of single-nucleotide polymorphisms (SNPs) included in the analysis. In this study, we aimed to address the inconsistent results and evaluate the association between DGKK and hypospadias in the Han Chinese population through a more comprehensive analysis of DGKK variants. We conducted association analyses for 17 SNPs in or downstream of DGKK with hypospadias among 322 cases (58 mild, 113 moderate, 128 severe, and 23 unknown) and 1008 controls. Five SNPs (rs2211122, rs4554617, rs7058226, rs7063116, and rs5915254) in DGKK were significantly associated with hypospadias (P < 0.05), with odds ratios (ORs) of 1.64-1.76. When only mild and moderate cases were compared to controls, 10 SNPs in DGKK were significant (P < 0.05), with ORs of 1.56-2.13. No significant SNP was observed when only severe cases were compared to controls. This study successfully implicated DGKK variants in hypospadias risk among a Han Chinese population, especially for mild/moderate cases. Severe forms of hypospadias are likely due to other genetic factors.
Keywords: association study; diacylglycerol kinase kappa; hypospadias; polymorphisms
Article in PDF
|How to cite this URL:|
Xie H, Lin XL, Zhang S, Yu L, Li XX, Huang YC, Lyu YQ, Chen HT, Xu J, Chen F. Association between diacylglycerol kinase kappa variants and hypospadias susceptibility in a Han Chinese population. Asian J Androl [Epub ahead of print] [cited 2017 Dec 14]. Available from: http://www.ajandrology.com/preprintarticle.asp?id=207709
Hua Xie, Xiao-Ling Lin
These authors contributed equally to the work.
| Introduction|| |
Hypospadias, affecting approximately 1 out of every 750 births in Europe, is a common congenital disease characterized by urogenital malformation, in which the position of urethral orifice is abnormal. Hypospadias are caused by incomplete urethral fusion during gestational weeks 8 to 16. The phenotype of hypospadias is divided into mild (glandular), moderate (midpenile), and severe (in scrotum or perineum) depending on the abnormal location of the urethral opening.,, The prevalence of hypospadias has increased in developed nations since the 1960s, and its incidence has plateaued in recent years.,, However, in China, the prevalence of hypospadias has shown an increasing trend, particularly in well-developed areas., Although the etiology of hypospadias is largely unknown, genetic factors have been demonstrated to play an important role in the development of hypospadias.,
Hypospadias is a complex disease affected by both genetic and environmental factors. The first genome-wide association study (GWAS) of hypospadias, among a European population, was reported in 2011, which included 436 Dutch cases and 494 controls. In the study, two single-nucleotide polymorphisms (SNPs) (rs1934179 and rs7063116) were identified in the diacylglycerol kinase kappa (DGKK) gene, which is located on the X chromosome and is strongly associated with hypospadias. Subsequently, a fine-mapping study was conducted in an American population to determine whether the association of DGKK with hypospadias could be replicated in a more racially diverse population in 2013. Results from this second study confirmed the relationship between hypospadias and the above-mentioned SNPs. Fifteen significant SNPs were also found associated with mild and moderate cases. In 2014, a study of the DGKK gene was performed in a Chinese population for the first time. Fourteen tag SNPs in the DGKK gene were tested, but all failed to show a statistically significant association with hypospadias. Considering the small number of tag SNPs and limited sample size in the previous Chinese study, we aimed to conduct more comprehensive research of the DGKK gene variants with all 17 previously identified SNPs in a larger case-control study. Our hypothesis was that DGKK genetic variants would be associated with hypospadias risk in a Han Chinese population.
| Materials and Methods|| |
This study included 322 unrelated cases and 1008 controls, all of which were Han Chinese. Patients with hypospadias were recruited from the Department of Urology at Shanghai Children's Hospital and were diagnosed by the Department of Urology from January 2013 to January 2014. Only patients with hypospadias without other system abnormalities were included. The severity of hypospadias cases was classified into one of three categories: mild (glandular), moderate (penile), or severe (in the scrotum or perineum) according to the position of the urethral opening. Of the 322 patients ultimately enrolled, there were 58 mild cases, 113 moderate cases, 128 severe cases, and 23 cases with unknown classification. The 1008 healthy controls were recruited from the Chinese Consortium for Prostate Cancer Genetics (ChinaPCa).
Each patient was informed of the purpose of this study, and written consent was obtained from all participants or their parent/legal guardian. Ethical approval was obtained from the Shanghai Children's Hospital in China.
DNA extraction and genotyping
Genomic DNA was extracted from peripheral blood samples using the Gentra Puregene Blood Kit (Qiagen, Dusseldorf, Germany). Polymerase chain reaction (PCR) and extension reactions were performed according to the manufacturers' protocol. The SNP genotypes were obtained using the MassARRAY iPLEX platform (Sequenom, San Diego, USA) and the genotyping call rates of all SNPs were >97%. Duplicated and water samples were included in each 96-well plate as PCR negative controls.
Selection of tagging SNPs
A total of 17 tagging SNPs were selected based on the Han Chinese population (CHB) in HapMap data and Genome Variation Server (http://gvs.gs.washington.edu/GVS/) with the criteria of r2 > 0.8 and minor allele frequency (MAF) >0.05. The SNP rs5915330 was located in CCNB3, and the rest were in DGKK. Two SNPs (rs7063116 and rs1934179) were obtained from the first GWAS study of a European population.
PLINK software was used to test the association of each SNP with risk of hypospadias. Association analyses were performed on all cases grouped together as well as separate subgroups based on phenotype severity. Mild and moderate cases were grouped together due to their phenotype similarities. The allelic OR and 95% confidence interval (CI) were calculated using logistic regression models to estimate relative risks. We performed clump and haplotype association analyses based on PLINK software. All P values were two-sided tests, and P < 0.05 was considered statistically significant. P = 0.003 (0.05/17) was the significance threshold through the strict Bonferroni correction. The pairwise linkage disequilibrium (LD) structure (r2) value was calculated using Haploview 4.2 software (https://www.broadinstitute.org/haploview/haploview).
| Results|| |
To validate whether the genetic variants of DGKK gene were associated with hypospadias risk in a Han Chinese population, we genotyped 17 tagging SNPs. Genotype information is detailed in [Table 1]. Sixteen SNPs are from the DGKK gene and one SNP (rs5915330) is located on the CCNB3 gene, downstream of DGKK. With all patients grouped together, 5 of the 16 DGKK SNPs (rs2211122, rs4554617, rs7058226, rs7063116, and rs5915254) were significantly associated with hypospadias, which had P < 0.05 and were associated with increased risk, with ORs ranging from 1.64 to 1.76. Another 5 SNPs (rs5915330, rs4143304, rs1934188, rs12171755, and rs1934179) in DGKK possessed marginal P values that were slightly more than 0.05. The other seven SNPs, including the variant in CCNB3, did not reach statistical significance with P < 0.05. Because of the similar results for mild and moderate cases, we merged them into one group for analysis (data not shown). Among mild/moderate cases, ten SNPs (rs4599945, rs4143304, rs1934188, rs12171755, rs1934179, rs2211122, rs4554617, rs7058226, rs7063116, and rs5915254) had P < 0.05 and their ORs were higher compared to the analysis including all cases. For severe cases, no P values reached statistical significance and ORs tended to be closer to 1.00.
|Table 1: Association of diacylglycerol kinase kappa single-nucleotide polymorphisms with different phenotype severities of hypospadias in a Han Chinese population |
Click here to view
To explore whether there was a high linkage disequilibrium that existed in these SNPs, we conducted a clump analysis ([Table 2]). The results showed that rs4143304, rs1934188, rs12171755, rs1934179, rs2211122, rs4554617, rs7058226, rs7063116, and rs5915254 were in the same block. All of them were found to be associated with hypospadias when only mild and moderate cases were compared with controls. Four out of these nine SNPs were estimated in the same haplotype block using Haploview 4.2 software ([Figure 1]).
|Figure 1: The pairwise linkage disequilibrium (LD) structure of 17 diacylglycerol kinase kappa (DGKK) SNPs and estimated haplotype blocks. Bottom, pairwise LD (r2) values were calculated based on data from the study samples and the color intensity of each SNP corresponded to pairwise D'/log of the odds (LOD). White squares: D' <1 and LOD <2; squares in shades of pink/red: D' <1 and LOD ≥2; bright red squares: D' of 1 and LOD ≥2. SNP: single-nucleotide polymorphism; D': the value of D prime between the two loci.|
Click here to view
|Table 2: Clump analysis based on linkage disequilibrium between diacylglycerol kinase kappa single-nucleotide polymorphisms |
Click here to view
We performed haplotype analysis in a block that included five SNPS (rs2211122, rs4554617, rs7058226, rs7063116, and rs5915254), which were in high linkage disequilibrium and identified three haplotypes ([Table 3]). The haplotype omnibus test revealed overall significant associations between these SNPs and hypospadias (P = 0.004). Meanwhile, individual haplotype analyses were consistent with the omnibus test and both ATGGA and GGAAG haplotypes reached significant P values (P = 0.002 and 0.001, respectively). ATGGA reflected the major allele and was the most common for each SNP, while GGAAG reflected the minor allele and was the next most common for each SNP.
|Table 3: Association of diacylglycerol kinase kappa haplotypes with hypospadias among all cases |
Click here to view
| Discussion|| |
In the present study, we found five SNPs (rs2211122, rs4554617, rs7058226, rs7063116, and rs5915254) that were significantly associated with hypospadias in DGKK (P < 0.05). When only mild and moderate cases were compared with controls, ten SNPs in DGKK were significantly associated with hypospadias (P < 0.05). No risk SNP was found to be associated with severe cases. Heritability of hypospadias is approximately 65%-75%, and the risk was estimated to be increased 12- to 20-fold among first-degree relatives.,, A previous study also demonstrated that genetic factors have a more important role in causing familial hypospadias than intrauterine environmental factors. GWAS analysis is a useful method in elucidating the genetic contributions of common variants. The first GWAS of hypospadias was conducted in a European population and identified two SNPs (rs1934179 and rs7063116) of DGKK, which had compelling evidence for association with hypospadias. These two susceptibility genetic loci and multiple other genetic loci of DGKK were further confirmed in moderate and mild cases in an independent American study population. However, the first study in a Han Chinese population indicated that the role of DGKK genetic variants was not likely to have major influences on hypospadias for the failure of identifying significant susceptibility genetic variants in DGKK. In our study, we confirmed the previous findings in the European population, and identified several novel variants in DGKK in a Han Chinese population, supporting the evidence that DGKK variants do contribute to urethral development as a major susceptibility gene for hypospadias.
According to previous studies, the severity of hypospadias might affect the association with DGKK SNPs; thus, we divided the 322 cases into mild, moderate, and severe subgroups. For the two previously reported SNPs,, SNP rs7063116 is outside of the DGKK coding sequence, located upstream of the coding region, and rs1934179 is located in an intron. We validated that SNP rs7063116 increased risk of hypospadias, except in severe cases, while SNP rs1934179 possessed only a marginal P value (0.05) in all cases, and did not show a significant difference after the strict Bonferroni correction. We also identified additional, neighboring SNPs (rs2211122, rs4554617, rs7058226, and rs5915254) that significantly increased risk for hypospadias. Two SNPs (rs4143304 and rs1934188) previously reported in an American population also showed significant association with hypospadias in our Han Chinese population. SNP rs12171755 was a novel genetic risk variant, which was significantly associated with mild to moderate hypospadias among the Han Chinese population in our study. Despite that none of these SNPs were found to have a relationship with hypospadias in severe cases in a Han Chinese population, our results indicated that genetic variants of DGKK play a major role in moderate and mild cases.
The DGKK gene, located on chromosome Xp11.22, encodes diacylglycerol kinase, which plays an important role in signal transduction by modulating the balance between diacylglycerol and phosphatidic acid. DGKK mRNA is found to be most abundant in the testis and second in the placenta. In addition, real-time quantitative PCR analyses showed that DGKK was expressed in the preputial skin of 10 healthy boys and 14 hypospadias cases. However, little is known about the functions of DGKK, and the biological mechanisms underlying the association between DGKK genetic variants and hypospadias risk remain elusive.
There are still several limitations in our study. First, the limited sample size may not have enough power to validate the genetic variants that exert effects on hypospadias in a Han Chinese population. Furthermore, the recruitment methods of cases and controls were different, which may have caused selection bias; thus, a larger sample size is warranted to confirm the association in a Han Chinese population. Second, little is known about the biological mechanisms of DGKK, and functional studies are needed in future projects.
| Conclusions|| |
In summary, we validated the finding of the first GWAS study of hypospadias and identified several novel SNPs in a Han Chinese population. These results from our fine-mapping study indicated that genetic variants of DGKK were associated with hypospadias in mild and moderate Han Chinese cases. Our findings support that further investigations and more comprehensive studies are warranted to address the functions and biological mechanisms of DGKK.
| Author Contributions|| |
HX, XLL, JX, and FC participated in conceiving and designing the study. HX, XLL, and SZ, participated in drafting the manuscript; XLL, SZ, and HTC participated in statistical analyses; HX, LY, XXL, YCH, and YQL participated in administrative, technical, and material support. All authors read and approved the final manuscript.
| Competing Interests|| |
All authors declare no competing interests.
| Acknowledgments|| |
This study was sponsored by the Shanghai Municipal Commission of Health and Family Planning (No. 20134283) and the Shanghai Municipal Commission of Science and Technology (No. 14411950403).
| References|| |
van der Zanden LF, van Rooij IA, Feitz WF, Knight J, Donders AR, et al.
Common variants in DGKK are strongly associated with risk of hypospadias. Nat Genet
2011; 43: 48-50.
Fredell L, Kockum I, Hansson E, Holmner S, Lundquist L, et al.
Heredity of hypospadias and the significance of low birth weight. J Urol
2002; 167: 1423-7.
Carmichael SL, Mohammed N, Ma C, Iovannisci D, Choudhry S, et al.
Diacylglycerol kinase K variants impact hypospadias in a California study population. J Urol
2013; 189: 305-11.
Ma Q, Tang Y, Lin H, Xu M, Xu G, et al.
Diacylglycerol kinase kappa (DGKK) variants and hypospadias in Han Chinese: association and meta-analysis. BJU Int
2015; 116: 634-40.
Paulozzi LJ. International trends in rates of hypospadias and cryptorchidism. Environ Health Perspect
1999; 107: 297-302.
Bergman JE, Loane M, Vrijheid M, Pierini A, Nijman RJ, et al.
Epidemiology of hypospadias in Europe: a registry-based study. World J Urol
2015; 33: 2159-67.
Caione P. Prevalence of hypospadias in European countries: is it increasing? Eur Urol
2009; 55: 1027-9.
Fisch H, Lambert SM, Hensle TW, Hyun G. Hypospadias rates in new york state are not increasing. J Urol
2009; 181: 2291-4.
Li Y, Mao M, Dai L, Li K, Li X, et al.
Time trends and geographic variations in the prevalence of hypospadias in China. Birth Defects Res A Clin Mol Teratol
2012; 94: 36-41.
Sun G, Tang D, Liang J, Wu M. Increasing prevalence of hypospadias associated with various perinatal risk factors in Chinese newborns. Urology
2009; 73: 1241-5.
George M, Schneuer FJ, Jamieson SE, Holland AJ. Genetic and environmental factors in the aetiology of hypospadias. Pediatr Surg Int
2015; 31: 519-27.
Soderhall C, Korberg IB, Thai HT, Cao J, Chen Y, et al.
Fine mapping analysis confirms and strengthens linkage of four chromosomal regions in familial hypospadias. Eur J Hum Genet
2015; 23: 516-22.
Purcell S, Neale B, Todd-Brown K, Thomas L, Ferreira MA, et al.
PLINK: a tool set for whole-genome association and population-based linkage analyses. Am J Hum Genet
2007; 81: 559-75.
Barrett JC, Fry B, Maller J, Daly MJ. Haploview: analysis and visualization of LD and haplotype maps. Bioinformatics
2005; 21: 263-5.
Carmichael SL, Shaw GM, Lammer EJ. Environmental and genetic contributors to hypospadias: a review of the epidemiologic evidence. Birth Defects Res A Clin Mol Teratol
2012; 94: 499-510.
Harris EL. Genetic epidemiology of hypospadias. Epidemiol Rev
1990; 12: 29-40.
Schnack TH, Zdravkovic S, Myrup C, Westergaard T, Christensen K, et al.
Familial aggregation of hypospadias: a cohort study. Am J Epidemiol
2008; 167: 251-6.
Imai S, Kai M, Yasuda S, Kanoh H, Sakane F. Identification and characterization of a novel human type II diacylglycerol kinase, DGK kappa. J Biol Chem
2005; 280: 39870-81.
[Table 1], [Table 2], [Table 3]