👤 Nicholas J Schork

Endometriosis is a common condition associated with debilitating pelvic pain and infertility. A genome-wide association study meta-analysis, including 60,674 cases and 701,926 controls of European and East Asian descent, identified 42 genome-wide significant loci comprising 49 distinct association signals. Effect sizes were largest for stage 3/4 disease, driven by ovarian endometriosis. Identified signals explained up to 5.01% of disease variance and regulated expression or methylation of genes in endometrium and blood, many of which were associated with pain perception/maintenance (SRP14/BMF, GDAP1, MLLT10, BSN and NGF). We observed significant genetic correlations between endometriosis and 11 pain conditions, including migraine, back and multisite chronic pain (MCP), as well as inflammatory conditions, including asthma and osteoarthritis. Multitrait genetic analyses identified substantial sharing of variants associated with endometriosis and MCP/migraine. Targeted investigations of genetically regulated mechanisms shared between endometriosis and other pain conditions are needed to aid the development of new treatments and facilitate early symptomatic intervention. Show less

Recent genome-wide association studies (GWAS) and pathway analyses supported long-standing observations of an association between immune-mediated diseases and Parkinson disease (PD). The post-GWAS era provides an opportunity for cross-phenotype analyses between different complex phenotypes. To test the hypothesis that there are common genetic risk variants conveying risk of both PD and autoimmune diseases (ie, pleiotropy) and to identify new shared genetic variants and their pathways by applying a novel statistical framework in a genome-wide approach. Using the conjunction false discovery rate method, this study analyzed GWAS data from a selection of archetypal autoimmune diseases among 138 511 individuals of European ancestry and systemically investigated pleiotropy between PD and type 1 diabetes, Crohn disease, ulcerative colitis, rheumatoid arthritis, celiac disease, psoriasis, and multiple sclerosis. NeuroX data (6927 PD cases and 6108 controls) were used for replication. The study investigated the biological correlation between the top loci through protein-protein interaction and changes in the gene expression and methylation levels. The dates of the analysis were June 10, 2015, to March 4, 2017. The primary outcome was a list of novel loci and their pathways involved in PD and autoimmune diseases. Genome-wide conjunctional analysis identified 17 novel loci at false discovery rate less than 0.05 with overlap between PD and autoimmune diseases, including known PD loci adjacent to GAK, HLA-DRB5, LRRK2, and MAPT for rheumatoid arthritis, ulcerative colitis and Crohn disease. Replication confirmed the involvement of HLA, LRRK2, MAPT, TRIM10, and SETD1A in PD. Among the novel genes discovered, WNT3, KANSL1, CRHR1, BOLA2, and GUCY1A3 are within a protein-protein interaction network with known PD genes. A subset of novel loci was significantly associated with changes in methylation or expression levels of adjacent genes. The study findings provide novel mechanistic insights into PD and autoimmune diseases and identify a common genetic pathway between these phenotypes. The results may have implications for future therapeutic trials involving anti-inflammatory agents. Show less

Women are under-represented in research on suicidality to date. Although women have a lower rate of suicide completion than men, due in part to the less-violent methods used, they have a higher rate of suicide attempts. Our group has previously identified genomic (blood gene expression biomarkers) and clinical information (apps) predictors for suicidality in men. We now describe pilot studies in women. We used a powerful within-participant discovery approach to identify genes that change in expression between no suicidal ideation (no SI) and high suicidal ideation (high SI) states (n=12 participants out of a cohort of 51 women psychiatric participants followed longitudinally, with diagnoses of bipolar disorder, depression, schizoaffective disorder and schizophrenia). We then used a Convergent Functional Genomics (CFG) approach to prioritize the candidate biomarkers identified in the discovery step by using all the prior evidence in the field. Next, we validated for suicidal behavior the top-ranked biomarkers for SI, in a demographically matched cohort of women suicide completers from the coroner's office (n=6), by assessing which markers were stepwise changed from no SI to high SI to suicide completers. We then tested the 50 biomarkers that survived Bonferroni correction in the validation step, as well as top increased and decreased biomarkers from the discovery and prioritization steps, in a completely independent test cohort of women psychiatric disorder participants for prediction of SI (n=33) and in a future follow-up cohort of psychiatric disorder participants for prediction of psychiatric hospitalizations due to suicidality (n=24). Additionally, we examined how two clinical instruments in the form of apps, Convergent Functional Information for Suicidality (CFI-S) and Simplified Affective State Scale (SASS), previously tested in men, perform in women. The top CFI-S item distinguishing high SI from no SI states was the chronic stress of social isolation. We then showed how the clinical information apps combined with the 50 validated biomarkers into a broad predictor (UP-Suicide), our apriori primary end point, predicts suicidality in women. UP-Suicide had a receiver-operating characteristic (ROC) area under the curve (AUC) of 82% for predicting SI and an AUC of 78% for predicting future hospitalizations for suicidality. Some of the individual components of the UP-Suicide showed even better results. SASS had an AUC of 81% for predicting SI, CFI-S had an AUC of 84% and the combination of the two apps had an AUC of 87%. The top biomarker from our sequential discovery, prioritization and validation steps, BCL2, predicted future hospitalizations due to suicidality with an AUC of 89%, and the panel of 50 validated biomarkers (BioM-50) predicted future hospitalizations due to suicidality with an AUC of 94%. The best overall single blood biomarker for predictions was PIK3C3 with an AUC of 65% for SI and an AUC of 90% for future hospitalizations. Finally, we sought to understand the biology of the biomarkers. BCL2 and GSK3B, the top CFG scoring validated biomarkers, as well as PIK3C3, have anti-apoptotic and neurotrophic effects, are decreased in expression in suicidality and are known targets of the anti-suicidal mood stabilizer drug lithium, which increases their expression and/or activity. Circadian clock genes were overrepresented among the top markers. Notably, PER1, increased in expression in suicidality, had an AUC of 84% for predicting future hospitalizations, and CSNK1A1, decreased in expression, had an AUC of 96% for predicting future hospitalizations. Circadian clock abnormalities are related to mood disorder, and sleep abnormalities have been implicated in suicide. Docosahexaenoic acid signaling was one of the top biological pathways overrepresented in validated biomarkers, which is of interest given the potential therapeutic and prophylactic benefits of omega-3 fatty acids. Some of the top biomarkers from the current work in women showed co-directionality of change in expression with our previous work in men, whereas others had changes in opposite directions, underlying the issue of biological context and differences in suicidality between the two genders. With this study, we begin to shed much needed light in the area of female suicidality, identify useful objective predictors and help understand gender commonalities and differences. During the conduct of the study, one participant committed suicide. In retrospect, when the analyses were completed, her UP-Suicide risk prediction score was at the 100 percentile of all participants tested. Show less

Genome-wide association studies (GWAS) have identified multiple common variants associated with body mass index (BMI). In this study, we tested 23 genotyped GWAS-significant SNPs (p-value<5*10-8) for longitudinal associations with BMI during childhood (3-17 years) and adulthood (18-45 years) for 658 subjects. We also proposed a heuristic forward search for the best joint effect model to explain the longitudinal BMI variation. After using false discovery rate (FDR) to adjust for multiple tests, childhood and adulthood BMI were found to be significantly associated with six SNPs each (q-value<0.05), with one SNP associated with both BMI measurements: KCTD15 rs29941 (q-value<7.6*10-4). These 12 SNPs are located at or near genes either expressed in the brain (BDNF, KCTD15, TMEM18, MTCH2, and FTO) or implicated in cell apoptosis and proliferation (FAIM2, MAP2K5, and TFAP2B). The longitudinal effects of FAIM2 rs7138803 on childhood BMI and MAP2K5 rs2241423 on adulthood BMI decreased as age increased (q-value<0.05). The FTO candidate SNPs, rs6499640 at the 5 '-end and rs1121980 and rs8050136 downstream, were associated with childhood and adulthood BMI, respectively, and the risk effects of rs6499640 and rs1121980 increased as birth weight decreased. The best joint effect model for childhood and adulthood BMI contained 14 and 15 SNPs each, with 11 in common, and the percentage of explained variance increased from 0.17% and 9.0*10(-6)% to 2.22% and 2.71%, respectively. In summary, this study evidenced the presence of long-term major effects of genes on obesity development, implicated in pathways related to neural development and cell metabolism, and different sets of genes associated with childhood and adulthood BMI, respectively. The gene effects can vary with age and be modified by prenatal development. The best joint effect model indicated that multiple variants with effects that are weak or absent alone can nevertheless jointly exert a large longitudinal effect on BMI. Show less

Anorexia nervosa (AN) is a mental illness with high mortality that most commonly afflicts adolescent female individuals. Clinical symptoms include chronic food refusal, weight loss and body image distortions. We carried out a genome-wide association study on 1033 AN cases and 3733 pediatric control subjects, all of whom were of European ancestry and were genotyped on the Illumina HumanHap610 platform (Illumina, San Diego, CA, USA). We confirmed that common single-nucleotide polymorphisms (SNPs) within OPRD1 (rs533123, P=0.0015) confer risk for AN, and obtained suggestive evidence that common SNPs near HTR1D (rs7532266, P=0.04) confer risk for restricting-type AN specifically. However, no SNPs reached genome-wide significance in our data, whereas top association signals were detected near ZNF804B, CSRP2BP, NTNG1, AKAP6 and CDH9. In parallel, we performed genome-wide analysis on copy number variations (CNVs) using the signal intensity data from the SNP arrays. We did not find evidence that AN cases have more CNVs than control subjects, nor do they have over-representation of rare or large CNVs. However, we identified several regions with rare CNVs that were only observed in AN cases, including a recurrent 13q12 deletion (1.5 Mb) disrupting SCAS in two cases, and CNVs disrupting the CNTN6/CNTN4 region in several AN cases. In conclusion, our study suggests that both common SNPs and rare CNVs may confer genetic risk to AN. These results point to intriguing genes that await further validation in independent cohorts for confirmatory roles in AN. Show less

To identify loci for age at menarche, we performed a meta-analysis of 32 genome-wide association studies in 87,802 women of European descent, with replication in up to 14,731 women. In addition to the known loci at LIN28B (P = 5.4 × 10⁻⁶⁰) and 9q31.2 (P = 2.2 × 10⁻³³), we identified 30 new menarche loci (all P < 5 × 10⁻⁸) and found suggestive evidence for a further 10 loci (P < 1.9 × 10⁻⁶). The new loci included four previously associated with body mass index (in or near FTO, SEC16B, TRA2B and TMEM18), three in or near other genes implicated in energy homeostasis (BSX, CRTC1 and MCHR2) and three in or near genes implicated in hormonal regulation (INHBA, PCSK2 and RXRG). Ingenuity and gene-set enrichment pathway analyses identified coenzyme A and fatty acid biosynthesis as biological processes related to menarche timing. Show less