Transcriptome analysis of the dihydrotestosterone-exposed fetal rat gubernaculum identifies common androgen and insulin-like 3 targets

Hyperandrogenism and insulin resistance-induced fetal loss: evidence for placental mitochondrial abnormalities and elevated reactive oxygen species production in pregnant rats that mimic the clinical features of polycystic ovary syndrome

KEY POINTS

Women with polycystic ovary syndrome (PCOS) commonly suffer from miscarriage, but the underlying mechanisms remain unknown. Herein, pregnant rats chronically treated with 5α-dihydrotestosterone (DHT) and insulin exhibited hyperandrogenism and insulin resistance, as well as increased fetal loss, and these features are strikingly similar to those observed in pregnant PCOS patients. Fetal loss in our DHT+insulin-treated pregnant rats was associated with mitochondrial dysfunction, disturbed superoxide dismutase 1 and Keap1/Nrf2 antioxidant responses, over-production of reactive oxygen species (ROS) and impaired formation of the placenta. Chronic treatment of pregnant rats with DHT or insulin alone indicated that DHT triggered many of the molecular pathways leading to placental abnormalities and fetal loss, whereas insulin often exerted distinct effects on placental gene expression compared to co-treatment with DHT and insulin. Treatment of DHT+insulin-treated pregnant rats with the antioxidant N-acetylcysteine improved fetal survival but was deleterious in normal pregnant rats. Our results provide insight into the fetal loss associated with hyperandrogenism and insulin resistance in women and suggest that physiological levels of ROS are required for normal placental formation and fetal survival during pregnancy.

ABSTRACT

Women with polycystic ovary syndrome (PCOS) commonly suffer from miscarriage, but the underlying mechanism of PCOS-induced fetal loss during pregnancy remains obscure and specific therapies are lacking. We used pregnant rats treated with 5α-dihydrotestosterone (DHT) and insulin to investigate the impact of hyperandrogenism and insulin resistance on fetal survival and to determine the molecular link between PCOS conditions and placental dysfunction during pregnancy. Our study shows that pregnant rats chronically treated with a combination of DHT and insulin exhibited endocrine aberrations such as hyperandrogenism and insulin resistance that are strikingly similar to those in pregnant PCOS patients. Of pathophysiological significance, DHT+insulin-treated pregnant rats had greater fetal loss and subsequently decreased litter sizes compared to normal pregnant rats. This negative effect was accompanied by impaired trophoblast differentiation, increased glycogen accumulation, and decreased angiogenesis in the placenta. Mechanistically, we report that over-production of reactive oxygen species (ROS) in the placenta, mitochondrial dysfunction, and disturbed superoxide dismutase 1 (SOD1) and Keap1/Nrf2 antioxidant responses constitute important contributors to fetal loss in DHT+insulin-treated pregnant rats. Many of the molecular pathways leading to placental abnormalities and fetal loss in DHT+insulin treatment were also seen in pregnant rats treated with DHT alone, whereas pregnant rats treated with insulin alone often exerted distinct effects on placental gene expression compared to insulin treatment in combination with DHT. We also found that treatment with the antioxidant N-acetylcysteine (NAC) improved fetal survival in DHT+insulin-treated pregnant rats, an effect related to changes in Keap1/Nrf2 and nuclear factor-κB signalling. However, NAC administration resulted in fetal loss in normal pregnant rats, most likely due to PCOS-like endocrine abnormality induced by the treatment. Our results suggest that the deleterious effects of hyperandrogenism and insulin resistance on fetal survival are related to a constellation of mitochondria-ROS-SOD1/Nrf2 changes in the placenta. Our findings also suggest that physiological levels of ROS are required for normal placental formation and fetal survival during pregnancy.

Subphenotype meta-analysis of testicular cancer genome-wide association study data suggests a role for RBFOX family genes in cryptorchidism susceptibility

STUDY QUESTION

Can subphenotype analysis of genome-wide association study (GWAS) data from subjects with testicular germ cell tumor (TGCT) provide insight into cryptorchidism (undescended testis, UDT) susceptibility?

SUMMARY ANSWER

Suggestive intragenic GWAS signals common to UDT, TGCT case-case and TGCT case-control analyses occur in genes encoding RBFOX RNA-binding proteins (RBPs) and their neurodevelopmental targets.

WHAT IS KNOWN ALREADY

UDT is a strong risk factor for TGCT, but while genetic risk factors for TGCT are well-known, genetic susceptibility to UDT is poorly understood and appears to be more complex.

STUDY DESIGN, SIZE, DURATION

We performed a secondary subphenotype analysis of existing GWAS data from the Testicular Cancer Consortium (TECAC) and compared these results with our previously published UDT GWAS data, and with data previously acquired from studies of the fetal rat gubernaculum.

PARTICIPANTS/MATERIALS, SETTING, METHODS

Studies from the National Cancer Institute (NCI), United Kingdom (UK) and University of Pennsylvania (Penn) that enrolled white subjects were the source of the TGCT GWAS data. We completed UDT subphenotype case-case (TGCT/UDT vs TGCT/non-UDT) and case-control (TGCT/UDT vs control), collectively referred to as 'TECAC' analyses, followed by a meta-analysis comprising 129 TGCT/UDT cases, 1771 TGCT/non-UDT cases, and 3967 unaffected controls. We reanalyzed our UDT GWAS results comprising 844 cases and 2718 controls by mapping suggestive UDT and TECAC signals (defined as P < 0.001) to genes using Ingenuity Pathway Analysis (IPA®). We compared associated pathways and enriched gene categories common to all analyses after Benjamini-Hochberg multiple testing correction, and analyzed transcript levels and protein expression using qRT-PCR and rat fetal gubernaculum confocal imaging, respectively.

MAIN RESULTS AND THE ROLE OF CHANCE

We found suggestive signals within 19 genes common to all three analyses, including RBFOX1 and RBFOX3, neurodevelopmental paralogs that encode RBPs targeting (U)GCATG-containing transcripts. Ten of the 19 genes participate in neurodevelopment and/or contribute to risk of neurodevelopmental disorders. Experimentally predicted RBFOX gene targets were strongly overrepresented among suggestive intragenic signals for the UDT (117 of 628 (19%), P = 3.5 × 10-24), TECAC case-case (129 of 711 (18%), P = 2.5 × 10-27) and TECAC case-control (117 of 679 (17%), P = 2 × 10-21) analyses, and a majority of the genes common to all three analyses (12 of 19 (63%), P = 3 × 10-9) are predicted RBFOX targets. Rbfox1, Rbfox2 and their encoded proteins are expressed in the rat fetal gubernaculum. Predicted RBFOX targets are also enriched among transcripts differentially regulated in the fetal gubernaculum during normal development (P = 3 × 10-31), in response to in vitro hormonal stimulation (P = 5 × 10-45) and in the cryptorchid LE/orl rat (P = 2 × 10-42).

LARGE SCALE DATA

GWAS data included in this study are available in the database of Genotypes and Phenotypes (dbGaP accession numbers phs000986.v1.p1 and phs001349.v1p1).

LIMITATIONS, REASONS FOR CAUTION

These GWAS data did not reach genome-wide significance for any individual analysis. UDT appears to have a complex etiology that also includes environmental factors, and such complexity may require much larger sample sizes than are currently available. The current methodology may also introduce bias that favors false discovery of larger genes.

WIDER IMPLICATIONS OF THE FINDINGS

Common suggestive intragenic GWAS signals suggest that RBFOX paralogs and other neurodevelopmental genes are potential UDT risk candidates, and potential TGCT susceptibility modifiers. Enrichment of predicted RBFOX targets among differentially expressed transcripts in the fetal gubernaculum additionally suggests a role for this RBP family in regulation of testicular descent. As RBFOX proteins regulate alternative splicing of Calca to generate calcitonin gene-related peptide, a protein linked to development and function of the gubernaculum, additional studies that address the role of these proteins in UDT are warranted.

STUDY FUNDING/COMPETING INTEREST(S)

The Eunice Kennedy Shriver National Institute for Child Health and Human Development (R01HD060769); National Center for Research Resources (P20RR20173), National Institute of General Medical Sciences (P20GM103464), Nemours Biomedical Research, the Testicular Cancer Consortium (U01CA164947), the Intramural Research Program of the NCI, a support services contract HHSN26120130003C with IMS, Inc., the Abramson Cancer Center at Penn, National Cancer Institute (CA114478), the Institute of Cancer Research, UK and the Wellcome Trust Case-Control Consortium (WTCCC) 2. None of the authors reports a conflict of interest.

Insulin-Like Peptide 3 (INSL3) Serum Concentration During Human Male Fetal Life

Context: Insulin-like peptide 3 (INSL3), a protein hormone produced by Leydig cells, may play a crucial role in testicular descent as male INSL3 knockout mice have bilateral cryptorchidism. Previous studies have measured human fetal INSL3 levels in amniotic fluid only. Objective: To measure INSL3 serum levels and mRNA in fetal umbilical cord blood and fetal testes, respectively. Design: INSL3 concentrations were assayed on 50 μl of serum from male human fetal umbilical cord blood by a non-commercial highly sensitive and specific radioimmunoassay. For secondary confirmation, quantitative real-time PCR was used to measure INSL3 relative mRNA expression in 7 age-matched human fetal testes. Setting: UT Southwestern Medical Center, Dallas, TX and Medical University of South Carolina, Charleston, SC. Patients or other Participants: Twelve human male umbilical cord blood samples and 7 human male testes were obtained from fetuses 14-21 weeks gestation. Male sex was verified by leukocyte genomic DNA SRY PCR. Interventions: None. Main Outcome Measures: Human male fetal INSL3 cord blood serum concentrations and testicular relative mRNA expression. Results: INSL3 serum concentrations during human male gestational weeks 15-20 were 2-4 times higher than published prepubertal male levels and were 5-100 times higher than previous reports of INSL3 concentrations obtained from amniotic fluid. Testicular fetal INSL3 mRNA relative expression was low from weeks 14-16, rose significantly weeks 17 and 18, and returned to low levels at week 21. Conclusions: These findings further support the role of INSL3 in human testicular descent and could prove relevant in uncovering the pathophysiology of cryptorchidism.

Molecular genetics of hypospadias and cryptorchidism recent developments

During the last decade, a tremendous amount of work has been devoted to the study of the molecular genetics of isolated hypospadias and cryptorchidism, two minor forms of disorders of sex development (DSD). Beyond the genes involved in gonadal determination and sex differentiation, including those underlying androgen biosynthesis and signaling, new genes have been identified through genome-wide association study and familial clustering. Even if no single genetic defect can explain the whole spectrum of DSD, these recent studies reinforce the strong role of the genetic background in the occurrence of these defects. The timing of signaling disruption may explain the different phenotypes.

Regional Variation in Androgen Receptor Expression and Biomechanical Properties May Contribute to Cryptorchidism Susceptibility in the LE/orl Rat

Background: The process of testicular descent requires androgen and insulin-like 3, hormones secreted by fetal Leydig cells. Knowledge concerning distinct and common functions of these hormones in regulating development of the fetal gubernaculum remains limited and/or conflicting. The current studies were designed to better define characteristics of androgen receptor (AR) expression, function and regulation, as well as the biomechanical properties of normal and cryptorchid gubernaculum during fetal development. Methods: We studied fetal gubernacula from Long Evans outbred (LE/wt) rats and an inbred (LE/orl) strain with an inherited form of cryptorchidism associated with an AR signaling defect. Gubernacular cells or whole organs obtained from LE/wt and LE/orl fetal gubernacula underwent AR immunostaining and quantitative image analysis. The effects of dihydrotestosterone (DHT) on AR expression, muscle fiber morphology, hyaluronan (HA) levels and glycosaminoglycan (GAG) content were measured in LE/wt gubernacula. Finally, the spatial mechanics of freshly harvested LE/wt and LE/orl fetal gubernacula were compared using micropipette aspiration. Results: AR is expressed in the nucleus of mesenchymal core, tip and cord cells of the embryonic (E) day 17 and 21 fetal gubernaculum, and is enhanced by DHT in primary cultures of gubernacular mesenchymal cells. Enhanced AR expression at the tip was observed in LE/wt but not LE/orl gubernacula. In in vitro studies of whole mount fetal gubernaculum, DHT did not alter muscle fiber morphology, HA content or GAG production. Progressive swelling with reduced cellular density of the LE/wt gubernaculum at E19-21 was associated with increased central stiffness in LE/wt but not in LE/orl fetuses. Conclusions: These data confirm nuclear AR expression in gubernacular mesenchyme with distal enhancement at the tip/cord region in LE/wt but not LE/orl rat fetuses. DHT enhanced cellular AR expression but had no major effects on muscle morphology or matrix composition in the rat fetal gubernaculum in vitro. Regional increased stiffness and decreased cell density between E19 and E21 were observed in LE/wt but not LE/orl fetal gubernacula. Developmental differences in cell-specific AR expression in LE/orl fetal gubernacula may contribute to the dysmorphism and aberrant function that underlies cryptorchidism susceptibility in this strain.

Perspective: A Neuro-Hormonal Systems Approach to Understanding the Complexity of Cryptorchidism Susceptibility

Nonsyndromic cryptorchidism is a common multifactorial, condition with long-term risks of subfertility and testicular cancer. Revealing the causes of cryptorchidism will likely improve prediction and prevention of adverse outcomes. Herein we provide our current perspective of cryptorchidism complexity in a synthesis of cumulative clinical and translational data generated by ourselves and others. From our recent comparison of genome-wide association study (GWAS) data of cryptorchidism with or without testicular germ cell tumor, we identified RBFOX family genes as candidate susceptibility loci. Notably, RBFOX proteins regulate production of calcitonin gene-related peptide (CGRP), a sensory neuropeptide linked to testicular descent in animal models. We also re-analyzed existing fetal testis transcriptome data from a rat model of inherited cryptorchidism (the LE/orl strain) for enrichment of Leydig cell progenitor genes. The majority are coordinately downregulated, consistent with known reduced testicular testosterone levels in the LE/orl fetus, and similarly suppressed in the gubernaculum. Using qRT-PCR, we found dysregulation of dorsal root ganglia (DRG) sensory transcripts ipsilateral to undescended testes. These data suggest that LE/orl cryptorchidism is associated with altered signaling in possibly related cell types in the testis and gubernaculum as well as DRG. Complementary rat and human studies thus lead us to propose a multi-level, integrated neuro-hormonal model of testicular descent. Variants in genes encoding RBFOX family proteins and/or their transcriptional targets combined with environmental exposures may disrupt this complex pathway to enhance cryptorchidism susceptibility. We believe that a systems approach is necessary to provide further insight into the causes and consequences of cryptorchidism.

Identification of early gene expression changes in primary cultured neurons treated with topoisomerase I poisons

Topoisomerase 1 (TOP1) poisons like camptothecin (CPT) are currently used in cancer chemotherapy but these compounds can have damaging, off-target effects on neurons leading to cognitive, sensory and motor deficits. To understand the molecular basis for the enhanced sensitivity of neurons to CPT, we examined the effects of compounds that inhibit TOP1-CPT, actinomycin D (ActD) and β-lapachone (β-Lap)-on primary cultured rat motor (MN) and cortical (CN) neurons as well as fibroblasts. Neuronal cells expressed higher levels of Top1 mRNA than fibroblasts but transcript levels are reduced in all cell types after treatment with CPT. Microarray analysis was performed to identify differentially regulated transcripts in MNs in response to a brief exposure to CPT. Pathway analysis of the differentially expressed transcripts revealed activation of ERK and JNK signaling cascades in CPT-treated MNs. Immediate-early genes like Fos, Egr-1 and Gadd45b were upregulated in CPT-treated MNs. Fos mRNA levels were elevated in all cell types treated with CPT; Egr-1, Gadd45b and Dyrk3 transcript levels, however, increased in CPT-treated MNs and CNs but decreased in CPT-treated fibroblasts. These transcripts may represent new targets for the development of therapeutic agents that mitigate the off-target effects of chemotherapy on the nervous system.

Polygenic inheritance of cryptorchidism susceptibility in the LE/orl rat

STUDY HYPOTHESIS

Susceptibility to inherited cryptorchidism in the LE/orl rat may be associated with genetic loci that influence developmental patterning of the gubernaculum by the fetal testis.

STUDY FINDING

Cryptorchidism in the LE/orl rat is associated with a unique combination of homozygous minor alleles at multiple loci, and the encoded proteins are co-localized with androgen receptor (AR) and Leydig cells in fetal gubernaculum and testis, respectively.

WHAT IS KNOWN ALREADY

Prior studies have shown aberrant perinatal gubernacular migration, muscle patterning defects and reduced fetal testicular testosterone in the LE/orl strain. In addition, altered expression of androgen-responsive, cytoskeletal and muscle-related transcripts in the LE/orl fetal gubernaculum suggest a role for defective AR signaling in cryptorchidism susceptibility.

STUDY DESIGN, SAMPLES/MATERIALS, METHODS

The long-term LE/orl colony and short-term colonies of outbred Crl:LE and Crl:SD, and inbred WKY/Ncrl rats were maintained for studies. Animals were intercrossed (LE/orl X WKY/Ncrl), and obligate heterozygotes were reciprocally backcrossed to LE/orl rats to generate 54 F2 males used for genotyping and/or linkage analysis. At least five fetuses per gestational time point from two or more litters were used for quantitative real-time RT-PCR (qRT-PCR) and freshly harvested embryonic (E) day 17 gubernaculum was used to generate conditionally immortalized cell lines. We completed genotyping and gene expression analyses using genome-wide microsatellite markers and single nucleotide polymorphism (SNP) arrays, PCR amplification, direct sequencing, restriction enzyme digest with fragment analysis, whole genome sequencing (WGS), and qRT-PCR. Linkage analysis was performed in Haploview with multiple testing correction, and qRT-PCR data were analyzed using ANOVA after log transformation. Imaging was performed using custom and commercial antibodies directed at candidate proteins in gubernaculum and testis tissues, and gubernaculum cell lines.

MAIN RESULTS AND THE ROLE OF CHANCE

LE/orl rats showed reduced fertility and fecundity, and higher risk of perinatal death as compared with Crl:LE rats, but there were no differences in breeding outcomes between normal and unilaterally cryptorchid males. Linkage analysis identified multiple peaks, and with selective breeding of outbred Crl:LE and Crl:SD strains for alleles within two of the most significant (P < 0.003) peaks on chromosomes 6 and 16, we were able to generate a non-LE/orl cryptorchid rat. Associated loci contain potentially functional minor alleles (0.25-0.36 in tested rat strains) including an exonic deletion in Syne2, a large intronic insertion in Ncoa4 (an AR coactivator) and potentially deleterious variants in Solh/Capn15, Ankrd28, and Hsd17b2. Existing WGS data indicate that homozygosity for these combined alleles does not occur in any other sequenced rat strain. We observed a modifying effect of the Syne2(del) allele on expression of other candidate genes, particularly Ncoa4, and for muscle and hormone-responsive transcripts. The selected candidate genes/proteins are highly expressed, androgen-responsive and/or co-localized with developing muscle and AR in fetal gubernaculum, and co-localized with Leydig cells in fetal testis.

LIMITATIONS, REASONS FOR CAUTION

The present study identified multiple cryptorchidism-associated linkage peaks in the LE/orl rat, containing potentially causal alleles. These are strong candidate susceptibility loci, but further studies are needed to demonstrate functional relevance to the phenotype.

WIDER IMPLICATIONS OF THE FINDINGS

Association data from both human and rat models of spontaneous, nonsyndromic cryptorchidism support a polygenic etiology of the disease. Both the present study and a human genome-wide association study suggest that common variants with weak effects contribute to susceptibility, and may exist in genes encoding proteins that participate in AR signaling in the developing gubernaculum. These findings have potential implications for the gene-environment interaction in the etiology of cryptorchidism.

LARGE SCALE DATA

Sequences were deposited in the Rat Genome Database (RGD, http://rgd.mcw.edu/).

STUDY FUNDING AND COMPETING INTERESTS

This work was supported by: R01HD060769 from the Eunice Kennedy Shriver National Institute for Child Health and Human Development (NICHD), 2P20GM103446 and P20GM103464 from the National Institute of General Medical Sciences (NIGMS), and Nemours Biomedical Research. The authors have no competing interests to declare.

Fetal Rat Gubernaculum Mesenchymal Cells Adopt Myogenic and Myofibroblast-Like Phenotypes

PURPOSE

Gubernaculum-cremaster complex development is hormonally regulated and abnormal in a cryptorchid rat model. Using cell tracking techniques and imaging we studied myogenic phenotypes and fates in the fetal rat gubernaculum-cremaster complex.

MATERIALS AND METHODS

Embryonic day 17 gubernaculum-cremaster complexes were labeled with CellTracker™ or the DNA synthesis marker EdU (5-ethynyl-2'-deoxyuridine), or immobilized in Matrigel® and grown in culture. Embryonic day 17 to 21 gubernaculum-cremaster complex sections and cells were imaged using wide field and deconvolution immunofluorescence microscopy, and muscle and/or myofibroblast specific antibodies. Deconvolved image stacks were used to create a 3-dimensional model of embryonic day 21 gubernaculum-cremaster complex muscle.

RESULTS

PAX7 (paired box 7) positive and myogenin positive muscle precursors were visible in a desmin-rich myogenic zone between muscle layers that elongated and became thicker during development. Gubernaculum-cremaster complex inner mesenchymal cells expressed desmin and αSMA (α smooth muscle actin) at lower levels than in the myogenic zone. After pulse labeling with CellTracker or EdU mesenchymal cells became incorporated into differentiated muscle. Conversely, mesenchymal cells migrated beyond Matrigel immobilized gubernaculum-cremaster complexes, expressed PAX7 and fused to form striated myotubes. Mesenchymal gubernaculum-cremaster complex cell lines proliferated more than 40 passages and showed contractile behavior but did not form striated muscle. Our 3-dimensional gubernaculum-cremaster complex model had 2 orthogonal ventral layers and an arcing inner layer of muscle.

CONCLUSIONS

Our data suggest that mesenchymal cells in the peripheral myogenic zone of the fetal gubernaculum-cremaster complex contribute to formation of a distinctively patterned cremaster muscle. Nonmyogenic, desmin and αSMA positive gubernaculum-cremaster complex mesenchymal cells proliferate and have a myofibroblast-like phenotype in culture. Intrinsic mechanical properties of these divergent cell types may facilitate perinatal inversion of the gubernaculum-cremaster complex.

Pathway analysis supports association of nonsyndromic cryptorchidism with genetic loci linked to cytoskeleton-dependent functions

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.

Phenotype specific association of the TGFBR3 locus with nonsyndromic cryptorchidism

PURPOSE

Based on a genome-wide association study of testicular dysgenesis syndrome showing a possible association with TGFBR3, we analyzed data from a larger, phenotypically restricted cryptorchidism population for potential replication of this signal.

MATERIALS AND METHODS

We excluded samples based on strict quality control criteria, leaving 844 cases and 2,718 controls of European ancestry that were analyzed in 2 separate groups based on genotyping platform (ie Illumina® HumanHap550, version 1 or 3, or Human610-Quad, version 1 BeadChip in group 1 and Human OmniExpress 12, version 1 BeadChip platform in group 2). Analyses included genotype imputation at the TGFBR3 locus, association analysis of imputed data with correction for population substructure, subsequent meta-analysis of data for groups 1 and 2, and selective genotyping of independent cases (330) and controls (324) for replication. We also measured Tgfbr3 mRNA levels and performed TGFBR3/betaglycan immunostaining in rat fetal gubernaculum.

RESULTS

We identified suggestive (p ≤ 1× 10(-4)) association of markers in/near TGFBR3, including rs9661103 (OR 1.40; 95% CI 1.20, 1.64; p = 2.71 × 10(-5)) and rs10782968 (OR 1.58; 95% CI 1.26, 1.98; p = 9.36 × 10(-5)) in groups 1 and 2, respectively. In subgroup analyses we observed strongest association of rs17576372 (OR 1.42; 95% CI 1.24, 1.60; p = 1.67 × 10(-4)) with proximal and rs11165059 (OR 1.32; 95% CI 1.15, 1.38; p = 9.42 × 10(-4)) with distal testis position, signals in strong linkage disequilibrium with rs9661103 and rs10782968, respectively. Association of the prior genome-wide association study signal (rs12082710) was marginal (OR 1.13; 95% CI 0.99, 1.28; p = 0.09 for group 1), and we were unable to replicate signals in our independent cohort. Tgfbr3/betaglycan was differentially expressed in wild-type and cryptorchid rat fetal gubernaculum.

CONCLUSIONS

These data suggest complex or phenotype specific association of cryptorchidism with TGFBR3 and the gubernaculum as a potential target of TGFβ signaling.

Cryptorchidism in the orl rat is associated with muscle patterning defects in the fetal gubernaculum and altered hormonal signaling

Cryptorchidism, or undescended testis, is a common male genital anomaly of unclear etiology. Hormonal stimulation of the developing fetal gubernaculum by testicular androgens and insulin-like 3 (INSL3) is required for testicular descent. In studies of the orl fetal rat, one of several reported strains with inherited cryptorchidism, we studied hormone levels, gene expression in intact and hormone-stimulated gubernaculum, and imaging of the developing cremaster muscle facilitated by a tissue clearing protocol to further characterize development of the orl gubernaculum. Abnormal localization of the inverted gubernaculum was visible soon after birth. In the orl fetus, testicular testosterone, gubernacular androgen-responsive transcript levels, and muscle-specific gene expression were reduced. However, the in vitro transcriptional response of the orl gubernaculum to androgen was largely comparable to wild type (wt). In contrast, increases in serum INSL3, gubernacular INSL3-responsive transcript levels, expression of the INSL3 receptor, Rxfp2, and the response of the orl gubernaculum to INSL3 in vitro all suggest enhanced activation of INSL3/RXFP2 signaling in the orl rat. However, DNA sequence analysis did not identify functional variants in orl Insl3. Finally, combined analysis of the present and previous studies of the orl transcriptome confirmed altered expression of muscle and cellular motility genes, and whole mount imaging revealed aberrant muscle pattern formation in the orl fetal gubernaculum. The nature and prevalence of developmental muscle defects in the orl gubernaculum are consistent with the cryptorchid phenotype in this strain. These data suggest impaired androgen and enhanced INSL3 signaling in the orl fetus accompanied by defective cremaster muscle development.

Insulin-Like Factor 3 and the HPG Axis in the Male

The hypothalamic-pituitary-gonadal (HPG) axis comprises pulsatile GnRH from the hypothalamus impacting on the anterior pituitary to induce expression and release of both LH and FSH into the circulation. These in turn stimulate receptors on testicular Leydig and Sertoli cells, respectively, to promote steroidogenesis and spermatogenesis. Both Leydig and Sertoli cells exhibit negative feedback to the pituitary and/or hypothalamus via their products testosterone and inhibin B, respectively, thereby allowing tight regulation of the HPG axis. In particular, LH exerts both acute control on Leydig cells by influencing steroidogenic enzyme activity, as well as chronic control by impacting on Leydig cell differentiation and gene expression. Insulin-like peptide 3 (INSL3) represents an additional and different endpoint of the HPG axis. This Leydig cell hormone interacts with specific receptors, called RXFP2, on Leydig cells themselves to modulate steroidogenesis, and on male germ cells, probably to synergize with androgen-dependent Sertoli cell products to support spermatogenesis. Unlike testosterone, INSL3 is not acutely regulated by the HPG axis, but is a constitutive product of Leydig cells, which reflects their number and/or differentiation status and their ability therefore to produce various factors including steroids, together this is referred to as Leydig cell functional capacity. Because INSL3 is not subject to the acute episodic fluctuations inherent in the HPG axis itself, it serves as an excellent marker for Leydig cell differentiation and functional capacity, as in puberty, or in monitoring the treatment of hypogonadal patients, and at the same time buffering the HPG output.