Pathway analysis supports association of nonsyndromic cryptorchidism with genetic loci linked to cytoskeleton-dependent functions
Julia Spencer Barthold,Yanping Wang,Thomas F Kolon,Claude Kollin,Agneta Nordenskjöld,Alicia Olivant Fisher,T Ernesto Figueroa,Ahmad H BaniHani,Jennifer A Hagerty,Ricardo Gonzaléz,Paul H Noh,Rosetta M Chiavacci,Kisha R Harden,Debra J Abrams,Cecilia E Kim,Jin Li,Hakon Hakonarson,Marcella Devoto
DOI: https://doi.org/10.1093/humrep/dev180
Abstract:Study question: What are the genetic loci that increase susceptibility to nonsyndromic cryptorchidism, or undescended testis? Summary answer: A genome-wide association study (GWAS) suggests that susceptibility to cryptorchidism is heterogeneous, with a subset of suggestive signals linked to cytoskeleton-dependent functions and syndromic forms of the disease. What is known already: Population studies suggest moderate genetic risk of cryptorchidism and possible maternal and environmental contributions to risk. Previous candidate gene analyses have failed to identify a major associated locus, although variants in insulin-like 3 (INSL3), relaxin/insulin-like family peptide receptor 2 (RXFP2) and other hormonal pathway genes may increase risk in a small percentage of patients. Study design, size, duration: This is a case-control GWAS of 844 boys with nonsyndromic cryptorchidism and 2718 control subjects without syndromes or genital anomalies, all of European ancestry. Participants/materials, setting, methods: All boys with cryptorchidism were diagnosed and treated by a pediatric specialist. In the discovery phase, DNA was extracted from tissue or blood samples and genotyping performed using the Illumina HumanHap550 and Human610-Quad (Group 1) or OmniExpress (Group 2) platform. We imputed genotypes genome-wide, and combined single marker association results in meta-analyses for all cases and for secondary subphenotype analyses based on testis position, laterality and age, and defined genome-wide significance as P = 7 × 10(-9) to correct for multiple testing. Selected markers were genotyped in an independent replication group of European cases (n = 298) and controls (n = 324). We used several bioinformatics tools to analyze top (P < 10(-5)) and suggestive (P < 10(-3)) signals for significant enrichment of signaling pathways, cellular functions and custom gene lists after multiple testing correction. Main results and the role of chance: In the full analysis, we identified 20 top loci, none reaching genome-wide significance, but one passing this threshold in a subphenotype analysis of proximal testis position (rs55867206, near SH3PXD2B, odds ratio = 2.2 (95% confidence interval 1.7, 2.9), P = 2 × 10(-9)). An additional 127 top loci emerged in at least one secondary analysis, particularly of more severe phenotypes. Cytoskeleton-dependent molecular and cellular functions were prevalent in pathway analysis of suggestive signals, and may implicate loci encoding cytoskeletal proteins that participate in androgen receptor signaling. Genes linked to human syndromic cryptorchidism, including hypogonadotropic hypogonadism, and to hormone-responsive and/or differentially expressed genes in normal and cryptorchid rat gubernaculum, were also significantly overrepresented. No tested marker showed significant replication in an independent population. The results suggest heterogeneous, multilocus and potentially multifactorial susceptibility to nonsyndromic cryptorchidism. Limitations, reasons for caution: The present study failed to identify genome-wide significant markers associated with cryptorchidism that could be replicated in an independent population, so further studies are required to define true positive signals among suggestive loci. Wider implications of the findings: As the only GWAS to date of nonsyndromic cryptorchidism, these data will provide a basis for future efforts to understand genetic susceptibility to this common reproductive anomaly and the potential for additive risk from environmental exposures. Study funding/competing interests: This work was supported by R01HD060769 (the Eunice Kennedy Shriver National Institute for Child Health and Human Development (NICHD)), P20RR20173 (the National Center for Research Resources (NCRR), currently P20GM103464 from the National Institute of General Medical Sciences (NIGMS)), an Institute Development Fund to the Center for Applied Genomics at The Children's Hospital of Philadelphia, and Nemours Biomedical Research. The authors have no competing interests to declare.