Novel evidence of CNV deletion in KCTD13 related to the severity of isolated hypospadias in Chinese population

Background: CNV in KCTD13 has been identified to influence androgen receptor function via its changes in gene dosage, which might contribute to hypospadias. However, there is lack of population-level evidence to assess the contribution of KCTD13 CNV to hypospadias.

Methods: 349 isolated hypospadias patients were recruited and their genotyping was performed using real-time qPCR. We use Database of Genomic Variants (DGV) and CNV calls from SNP-array intensity data in 1,008 Chinese healthy men as reference.

Results: 11.17% of patients were identified to have KCTD13 CNV deletion, significantly higher than 0.05% in DGV (P < 0.001), but no cases found to have CNV duplication. Meanwhile, no CNV calls encompassing KCTD13 region were detected in Chinese healthy men. Incidence of KCTD13 CNV deletion was significantly increased with the severity of hypospadias, P_{_trend} = 9.00 × 10⁻⁶. Compared to distal hypospadias, ORs for the proximal and midshaft were 10.07 (2.91–34.84) and 6.08 (1.69–21.84) respectively. In addition, the association between genital characteristics (stretched penile length and glans width) and KCTD13 CNV showed no significance in hypospadias children (P > 0.05).

Conclusions: We demonstrate KCTD13 CNV deletion is strongly associated with hypospadias and its severity, but duplication is not, characterizing KCTD13 genetic variation in more detail than previously described.

Introduction

Hypospadias is one of genitourinary (GU) tract anomalies, characterized by abnormal urethral development during embryonic development, resulting in the opening of the urethra being located anywhere along the underside of the penis (1). Its prevalence is approximately 20.9 per 10,000 births and has been increasing over time (2). Based on the abnormal location of urethral meatus, hypospadias is classified into distal, midshaft and proximal types (1). It may occur as an isolated disorder or accompany with other malformations (1, 3). Genetic factors have been demonstrated to play a principal role in the etiology of hypospadias, accounting for 77% of the heritability (4). Despite monogenic and chromosomal causes of hypospadias account for about 30% of all cases, the majority of genetic factors remain unknown (5–7). Previous studies are mainly based on sequencing and genotyping, and research on structural variations in hypospadias is limited (8, 9). Although structural variations are known to contribute to the genetic etiology of birth defects (10), few validated genetic markers have been associated with the risk of hypospadias in the Chinese population (11).

Copy number variation (CNV), a form of structural variations, is often observed as submicroscopic duplications or deletions of large genomic fragments exceeding 1 kb in length (12). They contribute to the development of human diseases by disrupting coding sequences, interfering with long-term gene regulation, or altering gene dosage (12). In a cohort of disorders of sex sexual development (DSD), 17.6% (3/17) of isolated hypospadias patients were found to carry a certain proportion of submicroscopic CNV, with duplication and deletion events noted on 16p11.2 for hypospadias (13). However, the contribution of CNV to hypospadias is still unknown. Investigation of CNV in the development of isolated hypospadias may provide valuable supplementary information for genetic research on this condition.

A recent study identified CNV events in potassium channel tetramerization domain containing 13 (KCTD13) on 16p11.2, which has been suggested as a potential candidate gene for GU anomalies in American population (14). The study observed a significantly higher prevalence of KCTD13 CNV among individuals with congenital anomalies of GU tract (2.58%) compared to the general population (0.10%) (14). However, due to the limited sample size of hypospadias cases in this study, the association between KCTD13 CNV and isolated hypospadias remains unclear. Therefore, we hypothesize that KCTD13 CNV may be identified in hypospadias in the Chinese population and could be associated with severity of hypospadias.

Patients and methods

Study cohort and design

In this study, we employed a case design (Table 1). A total of 349 patients with isolated hypospadias (no cryptorchidism and no micropenis) were recruited from Shanghai Children's Hospital of China between 2014 and 2021, ranging in age from 0.50 to 14.75 years (mean age: 2.98 ± 2.84 years; median age: 1.83 years). Clinical diagnosis was established by pediatric urologists and/or endocrinologists through direct clinical examination. Patients were excluded if they: (i) carried known karyotype abnormalities or abnormalities in the sex-determining region Y (SRY) gene abnormalities; (ii) were diagnosed with hypospadias as a result of disorders of sex development (DSD); (iii) showed endocrine disorder; (iv) were diagnosed with additional clinical features such as cryptorchidism or micropenis. According to the abnormal location of the urethral meatus, patients were classified into three categories: distal (glandular, coronal, and subcoronal: n = 129), midshaft (distal penile, midshaft, and proximal penile: n = 113), and proximal (penoscrotal, scrotal and perineal: n = 107). Stretched penile length (SPL) was standardized was measured by standardizing the dorsal length from the pubic bone to the apical end of the penis using a steel ruler after maximum traction (15), and a vernier caliper was used to measure the maximum width of the glans (16).

Table 1

Table 1. General characteristics of the study population.

All demographic and clinical information of participants were collected prior to surgery. This study was approved by the Ethics Committee of the Shanghai Children's Hospital in China (2022R045-E01), and written informed consents were obtained from all participants.

DNA extraction

Peripheral blood was collected from all subjects in vacuum blood tubes containing EDTA-K2. Genomic DNA was extracted from peripheral blood samples using the Gentra Puregene Blood Kit (Qiagen, Dusseldorf, Germany). The concentration and purity of DNA were determined using a Nanodrop ND2000 UV–vis spectrophotometer (Thermo Fisher Scientific, USA), and genomic DNA was stored at −80°C prior to CNV detection.

Determination of KCTD13 CNV

All samples were genotyped using the same genotyping system at the Central Laboratory of Shanghai Children's Hospital. The Applied Biosystems protocols for the TaqMan quantitative real-time polymerase chain reaction (qPCR) method were used to assess CNV (KCTD13: Assay Hs01712568_cn) in each sample. TaqMan CNV reactions were performed in triplicate, and all the above reagents were obtained from Thermo Fisher Scientific (Thermo Fisher Scientific, Waltham, MA, USA). The overall concordance rate for the CNV of KCTD13 among 10% of samples (39 duplicate samples) was 100%. The CNV assay was performed using an ABI QuantStudio6 Real-time PCR system. Same genotyping was analyzed using CopyCaller V2.1 software from Applied Biosystems.

Generation of CNV calls from intensity data of the genome-wide genotyping

We generated CNV calls from a genome-wide dataset of healthy controls in Chinese Consortium for Prostate Cancer Genetics (ChinaPCa), including 1,008 male Han Chinese from the community population. Details on study characteristics were reported in our previous study (17). A total of 731,458 SNPs have been genotyped using Illumina Human OmniExpress BeadChips (Illumina, San Diego, California).

A standard quality control procedure was used for generating CNV calls. SNPs with the following conditions were excluded: (i) genotype call rate <90%; (ii) minor allele frequency (MAF) < 0.01; or (iii) p < 0.001 for the Hardy–Weinberg equilibrium (HWE) test. Log R ratio (LRR) and B allele frequency (BAF) were generated using GenomeStudio (San Diego, CA, USA) and CNVs were called using cnvPartition (San Diego, CA, USA). Samples with the following conditions were removed: (i) standard deviation (SD) of LRR > 0.3; (ii) BAF_drift value >0.01; (iii) waviness factor (WF) > 0.05 or WF < −0.05; (iv) samples with more than 50 CNV calls. After individual quality control, a total of 982 controls remained. Minimum number of probes of 5 and confidence score of 35 were set as threshold for CNV calling when using cnvPartition analysis.

Statistical analysis

Mann-Whitney U test and t-test were utilized to investigate the statistical differences of quantitative variables with CNV. The Cochran-Armitage trend test was used to evaluate the trends of KCTD13 CNV detection rates in patients with different hypospadias ranges. Logistic regression analysis was performed to evaluate the independent influencing factors of KCTD13 CNV. A p value of <0.05 (two-sided) was considered to be statistical significant. All statistical analyses were performed using the Statistical Package for the Social Sciences (SPSS v.26.0, SPSS Inc.) and GraphPad Prism (v10.0.3, GraphPad Software). The Database of Genomic Variants (DGV), a database of human genomic structural variants (http://dgv.tcag.ca/), was used to as a reference for general population.

Results

The study cohort comprised 349 hypospadias patients, with 129 (36.96%) having distal, 113 (32.38%) for midshaft, and 107 (30.66%) proximal hypospadias. Table 1 showed the demographic information, including age, body mass index (BMI), stretched penile length (SPL) and glans width in the study population. The median age of patients with hypospadias was 1.83 (1.21–3.29) years old, while BMI was 15.97 (14.79–17.31) kg/m². Additionally, the median SPL and glans width were 5.00 (4.50–5.50) cm and 1.00 (0.90–1.20) cm, respectively.

Table 2 showed the carried rates of CNV detection of KCTD13 in the cohort. Among all patients, 11.17% (39 out of 349) were found to carry KCTD13 CNV deletion (≤1), and none with CNV duplication (>2) in KCTD13. Chi-squared test analysis showed a statistically significant difference in the frequency of copy number deletion between the three groups (P = 4.90 × 10⁻⁵, Table 3), and there was an increasing trend in the rate of carrying copy number deletion of the gene as the severity of hypospadias increased (P_{_trend} = 9.00 × 10⁻⁶, Figure 1).

Table 2

Table 2. Frequency distributions of KCTD13 CNV among different severity of hypospadias and its association with severity levels.

Table 3

Table 3. Comparison of general characteristics and clinical data of KCTD13 CNV status among hypospadias patients.

Figure 1

Figure 1. Patients with CNV in KCTD13 gene at 16p11.2 exhibit genitourinary tract anomalies. (A) Comparison of the incidences of KCTD13 CNV among our special disease cohort with isolated hypospadias, the general population from the Database of Genetic Variants (DGV) and Chinese healthy men according to existing SNP-array intensity data. (B) Distribution of KCTD13 copy number deletion in patients with different severity of hypospadias. P_{_trend} obtained using chi-squared tests of linear-by-linear association.

Combining general population samples from a public database (Database of Genomic Variants, DGV), we found that the frequency of KCTD13 CNV in the general population was low (0.05%, Figure 1). However, the samples in our study carried a much higher frequency of KCTD13 CNV, up to 11.17%, which is significantly higher than the former (P < 0.001). Furthermore, cnvPartity-based analyses of our existing SNP-array intensity data didn't generate any CNV calls spanning KCTD13 gene loci in the cohort of 1,008 healthy Chinese men.

Next, we performed the association of this CNV deletion with several clinical parameters. As shown in Table 3, the deletion was significantly correlated with lower levels of BMI (P = 0.031) and the more severe hypospadias (P < 0.001). However, the deletion was not related to SPL and Glans width, which are penis-related parameters. Further analysis showed that KCTD13 CNV deletion was still significantly associated with more severe hypospadias after adjustment for BMI by logistic regression model (Table 2). It indicated that patients with KCTD13 CNV deletion have a significantly increased risk for proximal hypospadias (OR = 10.07, 95% CI = 2.91–34.84, P = 2.70 × 10⁻⁴) and midshaft hypospadias (OR = 6.08, 95% CI = 1.69–21.84, P = 0.006) compared to the distal hypospadias.

Given the age, postnatal penile growth occurs concurrently with mini-puberty (from birth to 3 months of age) and then enters an extremely slow growth state before prepuberty (or during 1–3 years old). In our cohort, 62.18% (217/349) patients were aged between 1 and 3 years and 72.35% (157) of these had full record for SPL and glans width. Therefore, we selected these patients to reassess the correlation of KCTD13 CNV with penile parameters. As the result, no correlation was found between SPL, glans width and KCTD13 CNV in 1–3 years old children with isolated hypospadias (P > 0.05, Table 4).

Table 4

Table 4. Penile parameters in 1–3 years old children with hypospadias.

Discussion

In a previous study, KCTD13 CNV with deletion and duplication was observed to be associated with clinical cases of upper and lower tract genitourinary (GU) anomalies, including more than hypospadias, in the American population, and its changes in gene dosage could result in penile and testicular anomalies via diminished androgen receptor (AR) function (14). However, the precise role of KCTD13 CNV in hypospadias remains to be clarified. Our study determined the presence of KCTD13 CNV in cases of hypospadias and contributed to the evidence, indicating that the deletion of CNV in KCTD13 was significantly associated with increased risk to isolated hypospadias in a Chinese population. Although we did not find any association between KCTD13 CNV and penile parameters in children aged 1–3 years with isolated hypospadias, the deletion of this CNV was positively correlated with the severity of hypospadias. To the best of our knowledge, this is the first reported study to demonstrate a relationship between KCTD13 CNV and susceptibility to isolated hypospadias in the Chinese population.

In our study, the overall frequency of KCTD13 CNV in hypospadias was 11.17%, which is much higher than that observed in patients with GU anomalies in the American population (2.58%) according to Seth et al. (14). Meanwhile, we utilized the existing SNP-array intensity data in our previous study (17) and the general population within the DGV database as the reference and found that the frequency of KCTD13 CNV deletion in our cohort was also significantly higher than that in the general population (DGV: 0.05%), indicating its importance as a genetic marker for hypospadias susceptibility. Different from the previous research of Seth et al. (14), we strictly included isolated hypospadias as the study participants to rule out some unknown confounding variables, such as endocrine abnormality and symptom complex. All the detection of CNV type in KCTD13 was deletion (single gene copy) and its carriage was more common among Chinese (11.17%) compared to the American population (14).

With our research, it added to the evidence different copy number of KCTD13 reflected different contributions in the clinical phenotype. Previous research showed that a patient presenting both hypospadias and cryptorchidism exhibited KCTD13 copy number duplications, whereas only one case with isolated hypospadias displayed KCTD13 copy number deletion. In our study, all the determinations of KCTD13 CNV type in isolated hypospadias patients were deletion. Our findings further validated these previous results and suggested that children with isolated hypospadias are more likely to exhibit KCTD13 copy number deletion. The KCTD13 gene is located at the 16p11.2 locus, and prior studies have demonstrated that 16p11.2 copy number deletions are more prevalent than duplications (18). These results reinforce the findings of the current study. On the other hand, deletion of copy number of KCTD13 was also linked to other conditions like epilepsy and autism (19, 20). But limited research has been dedicated to investigate the association between CNV and the severity of diseases. In this study, we observed an association between KCTD13 CNV and the severity of disease. Patients with distal, midshaft, and proximal isolated hypospadias all exhibited a certain proportion of KCTD13 CNV, and the prevalence of KCTD13 copy number deletion increased as the severity of hypospadias worsened (P < 0.001).

Some studies have demonstrated that the severity of hypospadias seemed to correlate indicators of masculinization, such as SPL and glans width. In an observational study, the research found severe hypospadias patients had significantly shorter SPL compared to mild one when evaluating actual post-pubertal penile size and factors affecting it in hypospadias patients (21). Similarly, Bush et al. observed a smaller glans width among patients with hypospadias aged 0–24 months than the boys presenting for newborn circumcision (16). Building upon findings from previous mouse models (14), KCTD13-deficient mice exhibit testicular and penile abnormalities at maturity. Consequently, our study delved further into the analysis of the correlation between penile parameters and KCTD13 CNV. Considering the phenomenon of “small puberty” during infancy (from birth to 3 months), when the penis undergoes rapid growth (22), and the subsequent period of slow growth in children aged 1–3 years old (22–24), we selected patients within this age range for correlation analysis between penile parameters and KCTD13 CNV. However, our study did not identify any significant associations between SPL or glans width and KCTD13 CNV among isolated cases of hypospadias in children aged 1–3 years old. The limited number of samples available for comparison may have contributed to the absence of significant trends. The expansion of the sample in future studies may provide more clarity on the association between the aforementioned variables. Additionally, the previous study selected sexually mature Kctd13^−/−mice for comparison of penile development (14). While our study excluded interference caused by age development, it is important to note that the majority of children have not yet reached sexual maturity. Therefore, assessing whether their penile development is affected by KCTD13 CNV may not be appropriate, and long-term follow-up may be necessary to address this limitation.

KCTD13 functions as a substrate-specific adapter of the BCR (BTB-CUL3-RBX1) E3 ubiquitin-protein ligase complex (14, 25) and plays integral roles in regulating cytoskeleton, migration, proliferation, and neural development (25–27). A study has found that knockdown of intracellular KCTD13 reduces the translocation of AR into the nucleus upon androgen stimulation, which in turn affects the function of AR in transmitting androgen signals (14). Previous studies have indicated that the AR plays a crucial role in determining neonatal penile growth by controlling the proliferation of genital tubercle mesenchymal cells and promoting urethral fusion through the regulation of apoptosis. If these processes are compromised, it may result in the development of penile malformations like micropenis and hypospadias, and even lead to varying degrees of severity in hypospadias (28). However, abnormal expression, subcellular localization or function of AR can also contribute to hypospadias occurrence (28–30). Previous researches have demonstrated that E3 ubiquitin ligases can facilitate the degradation of AR (31, 32). The reduced expression of AR in testicular tissues of Kctd13^−/− mice (14) raises the hypothesis that the deletion of KCTD13 may lead to increased AR degradation, ultimately contributing to insufficient masculinization phenotypes. However, further studies are required to clarify the relationship between the two. Additionally, further research is still required to elucidate the molecular mechanisms underlying the association between KCTD13 and hypospadias.

Up to date, there were few CNV studies particularly related to hypospadias. These studies predominantly included cases of hypospadias accompanied with other symptom and could only identified large piece of CNV (>1.5 Mb) (Supplementary Table) (3, 13), which may not fully represent the genetic etiology of isolated cases. Compared to large piece of CNV, smaller ones could specify candidate genes predisposing to the diseases and help bridge the gap in understanding complex genetic disorders like hypospadias. Therefore, the contribution of CNV to the etiology of isolated hypospadias remains largely unexplored. Hypospadias is a clinically and genetically heterogeneous disorder (7), it not only presents as a local deformity but also serves as an indication of systemic endocrine dysfunction, particularly proximal hypospadias, which in the majority of cases are ultimately due to problems with the child's sex hormone levels (33). A consensus thought that patients with isolated hypospadias were more likely to have normal endocrine testicular function by comparing testicular function between patients with isolated hypospadias and those with hypospadias with micropenis, cryptorchidism, or ambiguous genitalia (34). That's the reasons we chose the isolated hypospadias as a typical object of study for this project, which could exclude the potential interference of endocrine factors and phenotypic heterogeneity. Moreover, unlike previous studies, the sample size in our study was considerably large. And for the first time, we have demonstrated the high prevalence of KCTD13 CNV among Chinese children with isolated hypospadias.

The limitation of this study lies in conducting only single gene copy number verification, without performing whole-genome CNV deletion analysis. Additionally, due to the majority of participants being young children with immature sexual development, there is a lack of data regarding penile development to comprehensively assess the degree of masculinization and predict its progression.

In the future, enhancing the long-term postoperative follow-up of children after surgery, particularly during puberty, will be beneficial for assessing and predicting the degree of masculinization. Furthermore, the prospective clinical utilization of KCTD13 copy number deletion may aid in evaluating and predicting the long-term development of masculinization in children with varying degrees of hypospadias, especially those with proximal hypospadias.

Conclusions

In conclusion, our study provides substantial evidence of an association between KCTD13 CNV deletion and susceptibility to isolated hypospadias. This approach allows us to identify the genetic markers for isolated hypospadias in Chinese boys and to analyze differences in the genetic determinants of hypospadias between the Chinese population and other populations.

Data availability statement

The data presented in the study are deposited in the Figshare repository. The data is available here https://doi.org/10.6084/m9.figshare.26934262.v1

Ethics statement

The studies involving humans were approved by Ethics Committee of the Shanghai Children's Hospital in China. The studies were conducted in accordance with the local legislation and institutional requirements. Written informed consent for participation in this study was provided by the participants’ legal guardians/next of kin.

Author contributions

YC: Formal Analysis, Methodology, Writing – original draft. LZ: Methodology, Visualization, Writing – review & editing. FC: Funding acquisition, Investigation, Writing – review & editing. ZC: Software, Writing – review & editing. YH: Investigation, Writing – review & editing. YL: Investigation, Writing – review & editing. MW: Investigation, Writing – review & editing. XL: Formal Analysis, Funding acquisition, Investigation, Writing – review & editing, Writing – original draft. HX: Investigation, Project administration, Supervision, Writing – review & editing, Writing – original draft.

Funding

The author(s) declare financial support was received for the research, authorship, and/or publication of this article.

This work was funded by grants from Science and Technology Commission of Shanghai Municipality to FC (No.23Y21900100), and Shanghai Children's Hospital to XL (No.2021YGZM03).

Acknowledgments

We thank all of the subjects included in this study.

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher's note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Supplementary material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fped.2024.1409264/full#supplementary-material

References

1. Bouty A, Ayers KL, Pask A, Heloury Y, Sinclair AH. The genetic and environmental factors underlying hypospadias. Sex Dev. (2015) 9(5):239–59. doi: 10.1159/000441988

PubMed Abstract | Crossref Full Text | Google Scholar

2. Yu X, Nassar N, Mastroiacovo P, Canfield M, Groisman B, Bermejo-Sánchez E, et al. Hypospadias prevalence and trends in international birth defect surveillance systems, 1980–2010. Eur Urol. (2019) 76(4):482–90. doi: 10.1016/j.eururo.2019.06.027

PubMed Abstract | Crossref Full Text | Google Scholar

3. Kon M, Suzuki E, Dung VC, Hasegawa Y, Mitsui T, Muroya K, et al. Molecular basis of non-syndromic hypospadias: systematic mutation screening and genome-wide copy-number analysis of 62 patients. Hum Reprod. (2015) 30(3):499–506. doi: 10.1093/humrep/deu364

PubMed Abstract | Crossref Full Text | Google Scholar

4. Schnack TH, Zdravkovic S, Myrup C, Westergaard T, Christensen K, Wohlfahrt J, et al. Familial aggregation of hypospadias: a cohort study. Am J Epidemiol. (2008) 167(3):251–6. doi: 10.1093/aje/kwm317

PubMed Abstract | Crossref Full Text | Google Scholar

5. Geller F, Feenstra B, Carstensen L, Pers TH, van Rooij IA, Körberg IB, et al. Genome-wide association analyses identify variants in developmental genes associated with hypospadias. Nat Genet. (2014) 46(9):957–63. doi: 10.1038/ng.3063

PubMed Abstract | Crossref Full Text | Google Scholar

6. Ságodi L, Kiss A, Kiss-Tóth E, Barkai L. Prevalence and possible causes of hypospadias. Orv Hetil. (2014) 155(25):978–85. doi: 10.1556/oh.2014.29858

Crossref Full Text | Google Scholar

7. van der Zanden LF, van Rooij IA, Feitz WF, Franke B, Knoers NV, Roeleveld N. Aetiology of hypospadias: a systematic review of genes and environment. Hum Reprod Update. (2012) 18(3):260–83. doi: 10.1093/humupd/dms002

PubMed Abstract | Crossref Full Text | Google Scholar

8. Scott CH, Amarillo IE. Identification of small regions of overlap from copy number variable regions in patients with hypospadias. Int J Mol Sci. (2022) 23(8):4246. doi: 10.3390/ijms23084246

PubMed Abstract | Crossref Full Text | Google Scholar

9. Singh N, Gupta DK, Sharma S, Sahu DK, Mishra A, Yadav DK, et al. Single-nucleotide and copy-number variance related to severity of hypospadias. Pediatr Surg Int. (2018) 34(9):991–1008. doi: 10.1007/s00383-018-4330-5

PubMed Abstract | Crossref Full Text | Google Scholar

10. Chen Z, Finnell RH, Lei Y, Wang H. Progress and clinical prospect of genomic structural variants investigation. Sci Bull (Beijing). (2024) 69(6):705–8. doi: 10.1016/j.scib.2024.01.035

PubMed Abstract | Crossref Full Text | Google Scholar

11. Lan F, Chen Z, Lin X. Systematic analysis and evaluation of chromosome aberrations in major birth defects associated with infertility. Intractable Rare Dis Res. (2024) 13(1):29–35. doi: 10.5582/irdr.2023.01100

PubMed Abstract | Crossref Full Text | Google Scholar

12. Stranger BE, Forrest MS, Dunning M, Ingle CE, Beazley C, Thorne N, et al. Relative impact of nucleotide and copy number variation on gene expression phenotypes. Science. (2007) 315(5813):848–53. doi: 10.1126/science.1136678

PubMed Abstract | Crossref Full Text | Google Scholar

13. Tannour-Louet M, Han S, Corbett ST, Louet JF, Yatsenko S, Meyers L, et al. Identification of de novo copy number variants associated with human disorders of sexual development. PLoS One. (2010) 5(10):e15392. doi: 10.1371/journal.pone.0015392

PubMed Abstract | Crossref Full Text | Google Scholar

14. Seth A, Rivera A, Chahdi A, Choi IS, Medina-Martinez O, Lewis S, et al. Gene dosage changes in KCTD13 result in penile and testicular anomalies via diminished androgen receptor function. Faseb J. (2022) 36(11):e22567. doi: 10.1096/fj.202200558R

PubMed Abstract | Crossref Full Text | Google Scholar

15. Cinaz P, Yeşilkaya E, Onganlar YH, Boyraz M, Bideci A, Camurdan O, et al. Penile anthropometry of normal prepubertal boys in Turkey. Acta Paediatr. (2012) 101(1):e33–36. doi: 10.1111/j.1651-2227.2011.02386.x

PubMed Abstract | Crossref Full Text | Google Scholar

16. Bush NC, DaJusta D, Snodgrass WT. Glans penis width in patients with hypospadias compared to healthy controls. J Pediatr Urol. (2013) 9(6 Pt B):1188–91. doi: 10.1016/j.jpurol.2013.05.004

PubMed Abstract | Crossref Full Text | Google Scholar

17. Xu J, Mo Z, Ye D, Wang M, Liu F, Jin G, et al. Genome-wide association study in Chinese men identifies two new prostate cancer risk loci at 9q31.2 and 19q13.4. Nat Genet. (2012) 44(11):1231–5. doi: 10.1038/ng.2424

PubMed Abstract | Crossref Full Text | Google Scholar

18. Tucker T, Giroux S, Clément V, Langlois S, Friedman JM, Rousseau F. Prevalence of selected genomic deletions and duplications in a French-Canadian population-based sample of newborns. Mol Genet Genomic Med. (2013) 1(2):87–97. doi: 10.1002/mgg3.12

PubMed Abstract | Crossref Full Text | Google Scholar

19. Gu J, Ke P, Guo H, Liu J, Liu Y, Tian X, et al. KCTD13-mediated ubiquitination and degradation of GluN1 regulates excitatory synaptic transmission and seizure susceptibility. Cell Death Differ. (2023) 30(7):1726–41. doi: 10.1038/s41418-023-01174-5

PubMed Abstract | Crossref Full Text | Google Scholar

20. Madison JM, Duong K, Vieux EF, Udeshi ND, Iqbal S, Requadt E, et al. Regulation of purine metabolism connects KCTD13 to a metabolic disorder with autistic features. iScience. (2021) 24(1):101935. doi: 10.1016/j.isci.2020.101935

PubMed Abstract | Crossref Full Text | Google Scholar

21. Moriya K, Nakamura M, Nishimura Y, Kitta T, Kanno Y, Chiba H, et al. Factors affecting post-pubertal penile size in patients with hypospadias. World J Urol. (2016) 34(9):1317–21. doi: 10.1007/s00345-016-1763-0

PubMed Abstract | Crossref Full Text | Google Scholar

22. Boas M, Boisen KA, Virtanen HE, Kaleva M, Suomi AM, Schmidt IM, et al. Postnatal penile length and growth rate correlate to serum testosterone levels: a longitudinal study of 1962 normal boys. Eur J Endocrinol. (2006) 154(1):125–9. doi: 10.1530/eje.1.02066

PubMed Abstract | Crossref Full Text | Google Scholar

23. Park SK, Ergashev K, Chung JM, Lee SD. Penile circumference and stretched penile length in prepubertal children: a retrospective, single-center pilot study. Investig Clin Urol. (2021) 62(3):324–30. doi: 10.4111/icu.20200495

PubMed Abstract | Crossref Full Text | Google Scholar

24. Wang YN, Zeng Q, Xiong F, Zeng Y. Male external genitalia growth curves and charts for children and adolescents aged 0 to 17 years in Chongqing, China. Asian J Androl. (2018) 20(6):567–71. doi: 10.4103/aja.aja_51_18

PubMed Abstract | Crossref Full Text | Google Scholar

25. Chen Y, Yang Z, Meng M, Zhao Y, Dong N, Yan H, et al. Cullin mediates degradation of RhoA through evolutionarily conserved BTB adaptors to control actin cytoskeleton structure and cell movement. Mol Cell. (2009) 35(6):841–55. doi: 10.1016/j.molcel.2009.09.004

PubMed Abstract | Crossref Full Text | Google Scholar

26. Escamilla CO, Filonova I, Walker AK, Xuan ZX, Holehonnur R, Espinosa F, et al. Kctd13 deletion reduces synaptic transmission via increased RhoA. Nature. (2017) 551(7679):227–31. doi: 10.1038/nature24470

PubMed Abstract | Crossref Full Text | Google Scholar

27. Srivastava D, Chakrabarti O. Mahogunin-mediated α-tubulin ubiquitination via noncanonical K6 linkage regulates microtubule stability and mitotic spindle orientation. Cell Death Dis. (2014) 5(2):e1064. doi: 10.1038/cddis.2014.1

PubMed Abstract | Crossref Full Text | Google Scholar

28. Zheng Z, Armfield BA, Cohn MJ. Timing of androgen receptor disruption and estrogen exposure underlies a spectrum of congenital penile anomalies. Proc Natl Acad Sci U S A. (2015) 112(52):E7194–7203. doi: 10.1073/pnas.1515981112

PubMed Abstract | Crossref Full Text | Google Scholar

29. Chao C, Jian YY, Zhao XY, Liu YL, Xie QK. The involvement of hsa_circ_0000417 in the development of hypospadias by regulating AR. Differentiation. (2020) 116:9–15. doi: 10.1016/j.diff.2020.09.003

PubMed Abstract | Crossref Full Text | Google Scholar

30. Chen L, Wang J, Lu W, Xiao Y, Ni J, Wang W, et al. Characterization with gene mutations in Han Chinese patients with hypospadias and function analysis of a novel AR genevariant. Front Genet. (2021) 12:673732. doi: 10.3389/fgene.2021.673732

PubMed Abstract | Crossref Full Text | Google Scholar

31. Qi J, Tripathi M, Mishra R, Sahgal N, Fazli L, Ettinger S, et al. The E3 ubiquitin ligase Siah2 contributes to castration-resistant prostate cancer by regulation of androgen receptor transcriptional activity. Cancer Cell. (2013) 23(3):332–46. doi: 10.1016/j.ccr.2013.02.016

PubMed Abstract | Crossref Full Text | Google Scholar

32. Rees I, Lee S, Kim H, Tsai FT. The E3 ubiquitin ligase CHIP binds the androgen receptor in a phosphorylation-dependent manner. Biochim Biophys Acta. (2006) 1764(6):1073–9. doi: 10.1016/j.bbapap.2006.03.013

PubMed Abstract | Crossref Full Text | Google Scholar

33. van der Horst HJ, de Wall LL. Hypospadias, all there is to know. Eur J Pediatr. (2017) 176(4):435–41. doi: 10.1007/s00431-017-2864-5

PubMed Abstract | Crossref Full Text | Google Scholar

34. Rey RA, Codner E, Iñíguez G, Bedecarrás P, Trigo R, Okuma C, et al. Low risk of impaired testicular sertoli and leydig cell functions in boys with isolated hypospadias. J Clin Endocrinol Metab. (2005) 90(11):6035–40. doi: 10.1210/jc.2005-1306

PubMed Abstract | Crossref Full Text | Google Scholar