Rare copy number variants in over 100,000 European ancestry subjects reveal multiple disease associations.

Copy number variants (CNVs) are suggested to have a widespread impact on the human genome and phenotypes. To understand the role of CNVs across human diseases, we examine the CNV genomic landscape of 100,028 unrelated individuals of European ancestry, using SNP and CGH array datasets. We observe an average CNV burden of ~650 kb, identifying a total of 11,314 deletion, 5625 duplication, and 2746 homozygous deletion CNV regions (CNVRs). In all, 13.7% are unreported, 58.6% overlap with at least one gene, and 32.8% interrupt coding exons. These CNVRs are significantly more likely to overlap OMIM genes (2.94-fold), GWAS loci (1.52-fold), and non-coding RNAs (1.44-fold), compared with random distribution (P < 1 × 10-3). We uncover CNV associations with four major disease categories, including autoimmune, cardio-metabolic, oncologic, and neurological/psychiatric diseases, and identify several drug-repurposing opportunities. Our results demonstrate robust frequency definition for large-scale rare variant association studies, identify CNVs associated with major disease categories, and illustrate the pleiotropic impact of CNVs in human disease.

ARAF recurrent mutation causes central conducting lymphatic anomaly treatable with a MEK inhibitor.

The treatment of lymphatic anomaly, a rare devastating disease spectrum of mostly unknown etiologies, depends on the patient manifestations1. Identifying the causal genes will allow for developing affordable therapies in keeping with precision medicine implementation2. Here we identified a recurrent gain-of-function ARAF mutation (c.640T>C:p.S214P) in a 12-year-old boy with advanced anomalous lymphatic disease unresponsive to conventional sirolimus therapy and in another, unrelated, adult patient. The mutation led to loss of a conserved phosphorylation site. Cells transduced with ARAF-S214P showed elevated ERK1/2 activity, enhanced lymphangiogenic capacity, and disassembly of actin skeleton and VE-cadherin junctions, which were rescued using the MEK inhibitor trametinib. The functional relevance of the mutation was also validated by recreating a lymphatic phenotype in a zebrafish model, with rescue of the anomalous phenotype using a MEK inhibitor. Subsequent therapy of the lead proband with a MEK inhibitor led to dramatic clinical improvement, with remodeling of the patient's lymphatic system with resolution of the lymphatic edema, marked improvement in his pulmonary function tests, cessation of supplemental oxygen requirements and near normalization of daily activities. Our results provide a representative demonstration of how knowledge of genetic classification and mechanistic understanding guides biologically based medical treatments, which in our instance was life-saving.

Pathogenic variant in EPHB4 results in central conducting lymphatic anomaly.

Central conducting lymphatic anomaly (CCLA) is one of the complex lymphatic anomalies characterized by dilated lymphatic channels, lymphatic channel dysmotility and distal obstruction affecting lymphatic drainage. We performed whole exome sequencing (WES) of DNA from a four-generation pedigree and examined the consequences of the variant by transfection of mammalian cells and morpholino and rescue studies in zebrafish. WES revealed a heterozygous mutation in EPHB4 (RefSeq NM_004444.4; c.2334 + 1G>C) and RNA-Seq demonstrated that the EPHB4 mutation destroys the normal donor site, which leads to the use of a cryptic splice donor that results in retention of the intervening 12-bp intron sequence. Transient co-expression of the wild-type and mutant EPHB4 proteins showed reduced phosphorylation of tyrosine, consistent with a loss-of-function effect. Zebrafish ephb4a morpholino resulted in vessel misbranching and deformities in the lymphatic vessel development, indicative of possible differentiation defects in lymphatic vessels, mimicking the lymphatic presentations of the patients. Immunoblot analysis using zebrafish lysates demonstrated over-activation of mTORC1 as a consequence of reduced EPHB4 signaling. Strikingly, drugs that inhibit mTOR signaling or RAS-MAPK signaling effectively rescued the misbranching phenotype in a comparable manner. Moreover, knock-in of EPHB4 mutation in HEK293T cells also induced mTORC1 activity. Our data demonstrate the pathogenicity of the identified EPHB4 mutation as a novel cause of CCLA and suggesting that ERK inhibitors may have therapeutic benefits in such patients with complex lymphatic anomalies.

Fasoracetam in adolescents with ADHD and glutamatergic gene network variants disrupting mGluR neurotransmitter signaling.

The glutamatergic neurotransmitter system may play an important role in attention-deficit hyperactivity disorder (ADHD). This 5-week, open-label, single-blind, placebo-controlled study reports the safety, pharmacokinetics and responsiveness of the metabotropic glutamate receptor (mGluR) activator fasoracetam (NFC-1), in 30 adolescents, age 12-17 years with ADHD, harboring mutations in mGluR network genes. Mutation status was double-blinded. A single-dose pharmacokinetic profiling from 50-800 mg was followed by a single-blind placebo at week 1 and subsequent symptom-driven dose advancement up to 400 mg BID for 4 weeks. NFC-1 treatment resulted in significant improvement. Mean Clinical Global Impressions-Improvement (CGI-I) and Severity (CGI-S) scores were, respectively, 3.79 at baseline vs. 2.33 at week 5 (P < 0.001) and 4.83 at baseline vs. 3.86 at week 5 (P < 0.001). Parental Vanderbilt scores showed significant improvement for subjects with mGluR Tier 1 variants (P < 0.035). There were no differences in the incidence of adverse events between placebo week and weeks on active drug. The trial is registered at https://clinicaltrials.gov/ct2/show/study/NCT02286817 .

Erratum: A genome-wide association study of anorexia nervosa suggests a risk locus implicated in dysregulated leptin signaling.

A correction to this article has been published and is linked from the HTML version of this paper. The error has been fixed in the paper.

A genome-wide association study of anorexia nervosa suggests a risk locus implicated in dysregulated leptin signaling.

We conducted a genome-wide association study (GWAS) of anorexia nervosa (AN) using a stringently defined phenotype. Analysis of phenotypic variability led to the identification of a specific genetic risk factor that approached genome-wide significance (rs929626 in EBF1 (Early B-Cell Factor 1); P = 2.04 × 10-7; OR = 0.7; 95% confidence interval (CI) = 0.61-0.8) with independent replication (P = 0.04), suggesting a variant-mediated dysregulation of leptin signaling may play a role in AN. Multiple SNPs in LD with the variant support the nominal association. This demonstrates that although the clinical and etiologic heterogeneity of AN is universally recognized, further careful sub-typing of cases may provide more precise genomic signals. In this study, through a refinement of the phenotype spectrum of AN, we present a replicable GWAS signal that is nominally associated with AN, highlighting a potentially important candidate locus for further investigation.

Association of a rare NOTCH4 coding variant with systemic sclerosis: a family-based whole exome sequencing study.

BACKGROUND: Systemic sclerosis (SSc) is a rheumatologic disease with a multifactorial etiology. Genome-wide association studies imply a polygenic, complex mode of inheritance with contributions from variation at the human leukocyte antigen locus and non-coding variation at a locus on chromosome 6p21, among other modestly impactful loci. Here we describe an 8-year-old female proband presenting with diffuse cutaneous SSc/scleroderma and a family history of SSc in a grandfather and maternal aunt. METHODS: We employed whole exome sequencing (WES) of three members of this family. We examined rare missense, nonsense, splice-altering, and coding indels matching an autosomal dominant inheritance model. We selected one missense variant for Sanger sequencing confirmation based on its predicted impact on gene function and location in a known SSc genetic locus. RESULTS: Bioinformatic analysis found eight candidate variants meeting our criteria. We identified a very rare missense variant in the regulatory NODP domain of NOTCH4 located at the 6p21 locus, c.4245G > A:p.Met1415Ile, segregating with the phenotype. This allele has a frequency of 1.83 × 10-5 by the data of the Exome Aggregation Consortium. CONCLUSION: This family suggests a novel mechanism of SSc pathogenesis in which a rare and penetrant coding variation can substantially elevate disease risk in contrast to the more modest non-coding variation typically found at this locus. These results suggest that modulation of the NOTCH4 gene might be responsible for the association signal at chromosome 6p21 in SSc.

Genomic copy number variation association study in Caucasian patients with nonsyndromic cryptorchidism.

BACKGROUND: Copy number variation (CNV) is a potential contributing factor to many genetic diseases. Here we investigated the potential association of CNV with nonsyndromic cryptorchidism, the most common male congenital genitourinary defect, in a Caucasian population. METHODS: Genome wide genotyping were performed in 559 cases and 1772 controls (Group 1) using Illumina HumanHap550 v1, HumanHap550 v3 or Human610-Quad platforms and in 353 cases and 1149 controls (Group 2) using the Illumina Human OmniExpress 12v1 or Human OmniExpress 12v1-1. Signal intensity data including log R ratio (LRR) and B allele frequency (BAF) for each single nucleotide polymorphism (SNP) were used for CNV detection using PennCNV software. After sample quality control, gene- and CNV-based association tests were performed using cleaned data from Group 1 (493 cases and 1586 controls) and Group 2 (307 cases and 1102 controls) using ParseCNV software. Meta-analysis was performed using gene-based test results as input to identify significant genes, and CNVs in or around significant genes were identified in CNV-based association test results. Called CNVs passing quality control and signal intensity visualization examination were considered for validation using TaqMan CNV assays and QuantStudio® 3D Digital PCR System. RESULTS: The meta-analysis identified 373 genome wide significant (p < 5X10-4) genes/loci including 49 genes/loci with deletions and 324 with duplications. Among them, 17 genes with deletion and 1 gene with duplication were identified in CNV-based association results in both Group 1 and Group 2. Only 2 genes (NUCB2 and UPF2) containing deletions passed CNV quality control in both groups and signal intensity visualization examination, but laboratory validation failed to verify these deletions. CONCLUSIONS: Our data do not support that structural variation is a major cause of nonsyndromic cryptorchidism.

GRIN2D Recurrent De Novo Dominant Mutation Causes a Severe Epileptic Encephalopathy Treatable with NMDA Receptor Channel Blockers.

N-methyl-D-aspartate receptors (NMDARs) are ligand-gated cation channels that mediate excitatory synaptic transmission. Genetic mutations in multiple NMDAR subunits cause various childhood epilepsy syndromes. Here, we report a de novo recurrent heterozygous missense mutation-c.1999G>A (p.Val667Ile)-in a NMDAR gene previously unrecognized to harbor disease-causing mutations, GRIN2D, identified by exome and candidate panel sequencing in two unrelated children with epileptic encephalopathy. The resulting GluN2D p.Val667Ile exchange occurs in the M3 transmembrane domain involved in channel gating. This gain-of-function mutation increases glutamate and glycine potency by 2-fold, increases channel open probability by 6-fold, and reduces receptor sensitivity to endogenous negative modulators such as extracellular protons. Moreover, this mutation prolongs the deactivation time course after glutamate removal, which controls the synaptic time course. Transfection of cultured neurons with human GRIN2D cDNA harboring c.1999G>A leads to dendritic swelling and neuronal cell death, suggestive of excitotoxicity mediated by NMDAR over-activation. Because both individuals' seizures had proven refractory to conventional antiepileptic medications, the sensitivity of mutant NMDARs to FDA-approved NMDAR antagonists was evaluated. Based on these results, oral memantine was administered to both children, with resulting mild to moderate improvement in seizure burden and development. The older proband subsequently developed refractory status epilepticus, with dramatic electroclinical improvement upon treatment with ketamine and magnesium. Overall, these results suggest that NMDAR antagonists can be useful as adjuvant epilepsy therapy in individuals with GRIN2D gain-of-function mutations. This work further demonstrates the value of functionally evaluating a mutation, enabling mechanistic understanding and therapeutic modeling to realize precision medicine for epilepsy.

Early Infantile Epileptic Encephalopathy in an STXBP1 Patient with Lactic Acidemia and Normal Mitochondrial Respiratory Chain Function.

A wide range of clinical findings have been associated with mutations in Syntaxin Binding Protein 1 (STXBP1), including multiple forms of epilepsy, nonsyndromic intellectual disability, and movement disorders. STXBP1 mutations have recently been associated with mitochondrial pathology, although it remains unclear if this phenotype is a part of the core feature for this gene disorder. We report a 7-year-old boy who presented for diagnostic evaluation of intractable epilepsy, episodic ataxia, resting tremor, and speech regression following a period of apparently normal early development. Mild lactic acidemia was detected on one occasion at the time of an intercurrent illness. Due to the concern for mitochondrial disease, ophthalmologic evaluation was performed that revealed bilateral midperiphery pigmentary mottling. Optical coherence tomography (OCT) testing demonstrated a bilaterally thickened ganglion cell layer in the perifovea. Skeletal muscle biopsy analysis showed no mitochondrial abnormalities or respiratory chain dysfunction. Exome sequencing identified a de novo c.1651C>T (p.R551C) mutation in STXBP1. Although mitochondrial dysfunction has been reported in some individuals, our proband had only mild lactic acidemia and no skeletal muscle tissue evidence of mitochondrial disease pathology. Thus, mitochondrial dysfunction is not an obligate feature of STXBP1 disease.

Variants in CXCR4 associate with juvenile idiopathic arthritis susceptibility.

BACKGROUND: Juvenile idiopathic arthritis (JIA) is the most common chronic rheumatic disease among children, the etiology of which involves a strong genetic component, but much of the underlying genetic determinants still remain unknown. Our aim was to identify novel genetic variants that predispose to JIA. METHODS: We performed a genome-wide association study (GWAS) and replication in a total of 1166 JIA cases and 9500 unrelated controls of European ancestry. Correlation of SNP genotype and gene expression was investigated. Then we conducted targeted resequencing of a candidate locus, among a subset of 480 cases and 480 controls. SUM test was performed to evaluate the association of the identified rare functional variants. RESULTS: The CXCR4 locus on 2q22.1 was found to be significantly associated with JIA, peaking at SNP rs953387. However, this result is subjected to subpopulation stratification within the subjects of European ancestry. After adjusting for principal components, nominal significant association remained (p < 10(-4)). Because of its interesting known function in immune regulation, we carried out further analyses to assess its relationship with JIA. Expression of CXCR4 was correlated with CXCR4 rs953387 genotypes in lymphoblastoid cell lines (p = 0.014) and T-cells (p = 0.0054). In addition, rare non-synonymous and stop-gain sequence variants in CXCR4, putatively damaging for CXCR4 function, were significantly enriched in JIA cases (p = 0.015). CONCLUSION: Our results suggest the association of CXCR4 variants with JIA, implicating that this gene may be involved in the pathogenesis of autoimmune disease. However, because this locus is subjected to population stratification within the subjects of European ancestry, additional replication is still necessary for this locus to be considered a true risk locus for JIA. This cell-surface chemokine receptor has already been targeted in other diseases and may serve as a tractable therapeutic target for a specific subset of pediatric arthritis patients with additional replication and functional validation of the locus.

Systematic data-querying of large pediatric biorepository identifies novel Ehlers-Danlos Syndrome variant.

BACKGROUND: Ehlers Danlos Syndrome is a rare form of inherited connective tissue disorder, which primarily affects skin, joints, muscle, and blood cells. The current study aimed at finding the mutation that causing EDS type VII C also known as "Dermatosparaxis" in this family. METHODS: Through systematic data querying of the electronic medical records (EMRs) of over 80,000 individuals, we recently identified an EDS family that indicate an autosomal dominant inheritance. The family was consented for genomic analysis of their de-identified data. After a negative screen for known mutations, we performed whole genome sequencing on the male proband, his affected father, and unaffected mother. We filtered the list of non-synonymous variants that are common between the affected individuals. RESULTS: The analysis of non-synonymous variants lead to identifying a novel mutation in the ADAMTSL2 (p. Gly421Ser) gene in the affected individuals. Sanger sequencing confirmed the mutation. CONCLUSION: Our work is significant not only because it sheds new light on the pathophysiology of EDS for the affected family and the field at large, but also because it demonstrates the utility of unbiased large-scale clinical recruitment in deciphering the genetic etiology of rare mendelian diseases. With unbiased large-scale clinical recruitment we strive to sequence as many rare mendelian diseases as possible, and this work in EDS serves as a successful proof of concept to that effect.

Association of Arrhythmia-Related Genetic Variants With Phenotypes Documented in Electronic Medical Records.

IMPORTANCE: Large-scale DNA sequencing identifies incidental rare variants in established Mendelian disease genes, but the frequency of related clinical phenotypes in unselected patient populations is not well established. Phenotype data from electronic medical records (EMRs) may provide a resource to assess the clinical relevance of rare variants. OBJECTIVE: To determine the clinical phenotypes from EMRs for individuals with variants designated as pathogenic by expert review in arrhythmia susceptibility genes. DESIGN, SETTING, AND PARTICIPANTS: This prospective cohort study included 2022 individuals recruited for nonantiarrhythmic drug exposure phenotypes from October 5, 2012, to September 30, 2013, for the Electronic Medical Records and Genomics Network Pharmacogenomics project from 7 US academic medical centers. Variants in SCN5A and KCNH2, disease genes for long QT and Brugada syndromes, were assessed for potential pathogenicity by 3 laboratories with ion channel expertise and by comparison with the ClinVar database. Relevant phenotypes were determined from EMRs, with data available from 2002 (or earlier for some sites) through September 10, 2014. EXPOSURES: One or more variants designated as pathogenic in SCN5A or KCNH2. MAIN OUTCOMES AND MEASURES: Arrhythmia or electrocardiographic (ECG) phenotypes defined by International Classification of Diseases, Ninth Revision (ICD-9) codes, ECG data, and manual EMR review. RESULTS: Among 2022 study participants (median age, 61 years [interquartile range, 56-65 years]; 1118 [55%] female; 1491 [74%] white), a total of 122 rare (minor allele frequency <0.5%) nonsynonymous and splice-site variants in 2 arrhythmia susceptibility genes were identified in 223 individuals (11% of the study cohort). Forty-two variants in 63 participants were designated potentially pathogenic by at least 1 laboratory or ClinVar, with low concordance across laboratories (Cohen κ = 0.26). An ICD-9 code for arrhythmia was found in 11 of 63 (17%) variant carriers vs 264 of 1959 (13%) of those without variants (difference, +4%; 95% CI, -5% to +13%; P = .35). In the 1270 (63%) with ECGs, corrected QT intervals were not different in variant carriers vs those without (median, 429 vs 439 milliseconds; difference, -10 milliseconds; 95% CI, -16 to +3 milliseconds; P = .17). After manual review, 22 of 63 participants (35%) with designated variants had any ECG or arrhythmia phenotype, and only 2 had corrected QT interval longer than 500 milliseconds. CONCLUSIONS AND RELEVANCE: Among laboratories experienced in genetic testing for cardiac arrhythmia disorders, there was low concordance in designating SCN5A and KCNH2 variants as pathogenic. In an unselected population, the putatively pathogenic genetic variants were not associated with an abnormal phenotype. These findings raise questions about the implications of notifying patients of incidental genetic findings.

The Philadelphia Neurodevelopmental Cohort: A publicly available resource for the study of normal and abnormal brain development in youth.

The Philadelphia Neurodevelopmental Cohort (PNC) is a large-scale study of child development that combines neuroimaging, diverse clinical and cognitive phenotypes, and genomics. Data from this rich resource is now publicly available through the Database of Genotypes and Phenotypes (dbGaP). Here we focus on the data from the PNC that is available through dbGaP and describe how users can access this data, which is evolving to be a significant resource for the broader neuroscience community for studies of normal and abnormal neurodevelopment.

Genetic sharing and heritability of paediatric age of onset autoimmune diseases.

Autoimmune diseases (AIDs) are polygenic diseases affecting 7-10% of the population in the Western Hemisphere with few effective therapies. Here, we quantify the heritability of paediatric AIDs (pAIDs), including JIA, SLE, CEL, T1D, UC, CD, PS, SPA and CVID, attributable to common genomic variations (SNP-h(2)). SNP-h(2) estimates are most significant for T1D (0.863±s.e. 0.07) and JIA (0.727±s.e. 0.037), more modest for UC (0.386±s.e. 0.04) and CD (0.454±0.025), largely consistent with population estimates and are generally greater than that previously reported by adult GWAS. On pairwise analysis, we observed that the diseases UC-CD (0.69±s.e. 0.07) and JIA-CVID (0.343±s.e. 0.13) are the most strongly correlated. Variations across the MHC strongly contribute to SNP-h(2) in T1D and JIA, but does not significantly contribute to the pairwise rG. Together, our results partition contributions of shared versus disease-specific genomic variations to pAID heritability, identifying pAIDs with unexpected risk sharing, while recapitulating known associations between autoimmune diseases previously reported in adult cohorts.

Meta-analysis of shared genetic architecture across ten pediatric autoimmune diseases.

Genome-wide association studies (GWASs) have identified hundreds of susceptibility genes, including shared associations across clinically distinct autoimmune diseases. We performed an inverse χ(2) meta-analysis across ten pediatric-age-of-onset autoimmune diseases (pAIDs) in a case-control study including more than 6,035 cases and 10,718 shared population-based controls. We identified 27 genome-wide significant loci associated with one or more pAIDs, mapping to in silico-replicated autoimmune-associated genes (including IL2RA) and new candidate loci with established immunoregulatory functions such as ADGRL2, TENM3, ANKRD30A, ADCY7 and CD40LG. The pAID-associated single-nucleotide polymorphisms (SNPs) were functionally enriched for deoxyribonuclease (DNase)-hypersensitivity sites, expression quantitative trait loci (eQTLs), microRNA (miRNA)-binding sites and coding variants. We also identified biologically correlated, pAID-associated candidate gene sets on the basis of immune cell expression profiling and found evidence of genetic sharing. Network and protein-interaction analyses demonstrated converging roles for the signaling pathways of type 1, 2 and 17 helper T cells (TH1, TH2 and TH17), JAK-STAT, interferon and interleukin in multiple autoimmune diseases.

Copy Number Variations in CTNNA3 and RBFOX1 Associate with Pediatric Food Allergy.

Food allergy is a significant public health concern, especially among children. Previous candidate gene studies suggested a few susceptibility loci for food allergy, but no study investigated the contribution of copy number variations (CNVs) to food allergy on a genome-wide scale. To investigate the genetics of food allergy, we performed CNV assessment using high-resolution genome-wide single nucleotide polymorphism arrays. CNV calls from a total of 357 cases with confirmed food allergy and 3980 controls were analyzed within a discovery cohort, followed by a replication analysis composed of 167 cases and 1573 controls. We identified that CNVs in CTNNA3 were significantly associated with food allergy in both the discovery cohort and the replication cohort. Of particular interest, CTNNA3 CNVs hit exons or intron regions rich in histone marker H3K4Me1. CNVs in a second gene (RBFOX1) showed a significant association (p = 7.35 × 10(-5)) with food allergy at the genome-wide level in our meta-analysis of the European ancestry cohorts. The presence of these CNVs was confirmed by quantitative PCR. Furthermore, knockdown of CTNNA3 resulted in upregulation of CD63 and CD203c in mononuclear cells upon PMA stimulation, suggesting a role in sensitization to allergen. We uncovered at least two plausible genes harboring CNV loci that are enriched in pediatric patients with food allergies. The novel gene candidates discovered in this study by genome-wide CNV analysis are compelling drug and diagnostic targets for food allergy.

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.

The Philadelphia Neurodevelopmental Cohort: constructing a deep phenotyping collaborative.

BACKGROUND: An integrative multidisciplinary approach is required to elucidate the multiple factors that shape neurodevelopmental trajectories of mental disorders. The Philadelphia Neurodevelopmental Cohort (PNC), funded by the National Institute of Mental Health Grand Opportunity (GO) mechanism of the American Recovery and Reinvestment Act, was designed to characterize clinical and neurobehavioral phenotypes of genotyped youths. Data generated, which are recently available through the NIMH Database of Genotypes and Phenotypes (dbGaP), have garnered considerable interest. We provide an overview of PNC recruitment and clinical assessment methods to allow informed use and interpretation of the PNC resource by the scientific community. We also evaluate the structure of the assessment tools and their criterion validity. METHODS: Participants were recruited from a large pool of youths (n = 13,958) previously identified and genotyped at The Children's Hospital of Philadelphia. A comprehensive computerized tool for structured evaluation of psychopathology domains (GOASSESS) was constructed. We administered GOASSESS to all participants and used factor analysis to evaluate its structure. RESULTS: A total of 9,498 youths (aged 8-21; mean age = 14.2; European American = 55.8%; African American = 32.9%; Other = 11.4%) were enrolled. Factor analysis revealed a strong general psychopathology factor, and specific 'anxious-misery', 'fear', and 'behavior' factors. The 'behavior' factor had a small negative correlation (-0.21) with overall accuracy of neurocognitive performance, particularly in tests of executive and complex reasoning. Being female had a high association with the 'anxious-misery' and low association with the 'behavior' factors. The psychosis spectrum was also best characterized by a general factor and three specific factors: ideas about 'special abilities/persecution,' 'unusual thoughts/perceptions', and 'negative/disorganized' symptoms. CONCLUSIONS: The PNC assessment mechanism yielded psychopathology data with strong factorial validity in a large diverse community cohort of genotyped youths. Factor scores should be useful for dimensional integration with other modalities (neuroimaging, genomics). Thus, PNC public domain resources can advance understanding of complex inter-relationships among genes, cognition, brain, and behavior involved in neurodevelopment of common mental disorders.

Comorbidity of physical and mental disorders in the neurodevelopmental genomics cohort study.

OBJECTIVES: To examine patterns of associations between a broad range of mental and physical conditions by using a large, systematically obtained pediatric registry. METHODS: The sample included 9014 youth ages 8 to 21 years (4349 males and 4665 females; 3585 aged <13 years, 3678 aged 13 to 18 years, and 1751 aged 19 to 21 years) from the Philadelphia Neurodevelopmental Cohort identified through pediatric clinics at the Children's Hospital of Philadelphia health care network by the Center for Applied Genomics. Measures were as follows: physical condition based on electronic medical records and interview data on 42 physical conditions of 14 organ systems/specialties and mental disorders based on an abbreviated version of the structured Kiddie-Schedule for Affective Disorders and Schizophrenia psychiatric diagnostic interview. RESULTS: There was a direct association between the severity of the physical condition and most classes of mental disorders, as well as with functional impairment. Models adjusted for sociodemographic correlates, other physical and mental disorders, and false discovery and revealed broad patterns of associations between neurodevelopmental disorders with behavior disorders (odds ratio [OR]: 1.5; 95% confidence interval [CI]: 1.3-1.8; P < .004) and attention-deficit/hyperactivity disorder (OR: 3.1; 95% CI: 2.7-3.6; P < .0001), and neurologic/central nervous system conditions (OR: 1.3; 95% CI: 1.1-1.9; P < .05) with mood disorders and attention-deficit/hyperactivity disorder (OR: 1.3; 95% CI: 1.1-1.5; P < .001), and autoimmune/inflammatory conditions with mood disorders (OR: 1.4; 95% CI: 1.1-1.8, P < .05). CONCLUSIONS: Findings show the strong overlap between physical and mental conditions and their impact on severity and functional impairment in youth. Specific patterns of comorbidity have important implications for etiology. Prospective tracking of cross-disorder morbidity will be important to establish more effective mechanisms for prevention and intervention.