Exome sequencing of 1190 non-syndromic clubfoot cases reveals HOXD12 as a novel disease gene.

BACKGROUND: Clubfoot, presenting as a rigid inward and downward turning of the foot, is one of the most common congenital musculoskeletal anomalies. The aetiology of clubfoot is poorly understood and variants in known clubfoot disease genes account for only a small portion of the heritability. METHODS: Exome sequence data were generated from 1190 non-syndromic clubfoot cases and their family members from multiple ethnicities. Ultra-rare variant burden analysis was performed comparing 857 unrelated clubfoot cases with European ancestry with two independent ethnicity-matched control groups (1043 in-house and 56 885 gnomAD controls). Additional variants in prioritised genes were identified in a larger cohort, including probands with non-European ancestry. Segregation analysis was performed in multiplex families when available. RESULTS: Rare variants in 29 genes were enriched in clubfoot cases, including PITX1 (a known clubfoot disease gene), HOXD12, COL12A1, COL9A3 and LMX1B. In addition, rare variants in posterior HOX genes (HOX9-13) were enriched overall in clubfoot cases. In total, variants in these genes were present in 8.4% (100/1190) of clubfoot cases with both European and non-European ancestry. Among these, 3 are de novo and 22 show variable penetrance, including 4 HOXD12 variants that segregate with clubfoot. CONCLUSION: We report HOXD12 as a novel clubfoot disease gene and demonstrate a phenotypic expansion of known disease genes (myopathy gene COL12A1, Ehlers-Danlos syndrome gene COL9A3 and nail-patella syndrome gene LMX1B) to include isolated clubfoot.

PheWAS and Beyond: The Landscape of Associations with Medical Diagnoses and Clinical Measures across 38,662 Individuals from Geisinger.

Most phenome-wide association studies (PheWASs) to date have used a small to moderate number of SNPs for association with phenotypic data. We performed a large-scale single-cohort PheWAS, using electronic health record (EHR)-derived case-control status for 541 diagnoses using International Classification of Disease version 9 (ICD-9) codes and 25 median clinical laboratory measures. We calculated associations between these diagnoses and traits with ∼630,000 common frequency SNPs with minor allele frequency > 0.01 for 38,662 individuals. In this landscape PheWAS, we explored results within diseases and traits, comparing results to those previously reported in genome-wide association studies (GWASs), as well as previously published PheWASs. We further leveraged the context of functional impact from protein-coding to regulatory regions, providing a deeper interpretation of these associations. The comprehensive nature of this PheWAS allows for novel hypothesis generation, the identification of phenotypes for further study for future phenotypic algorithm development, and identification of cross-phenotype associations.

A genome-wide association study of polycystic ovary syndrome identified from electronic health records.

BACKGROUND: Polycystic ovary syndrome is the most common endocrine disorder affecting women of reproductive age. A number of criteria have been developed for clinical diagnosis of polycystic ovary syndrome, with the Rotterdam criteria being the most inclusive. Evidence suggests that polycystic ovary syndrome is significantly heritable, and previous studies have identified genetic variants associated with polycystic ovary syndrome diagnosed using different criteria. The widely adopted electronic health record system provides an opportunity to identify patients with polycystic ovary syndrome using the Rotterdam criteria for genetic studies. OBJECTIVE: To identify novel associated genetic variants under the same phenotype definition, we extracted polycystic ovary syndrome cases and unaffected controls based on the Rotterdam criteria from the electronic health records and performed a discovery-validation genome-wide association study. STUDY DESIGN: We developed a polycystic ovary syndrome phenotyping algorithm on the basis of the Rotterdam criteria and applied it to 3 electronic health record-linked biobanks to identify cases and controls for genetic study. In the discovery phase, we performed an individual genome-wide association study using the Geisinger MyCode and the Electronic Medical Records and Genomics cohorts, which were then meta-analyzed. We attempted validation of the significant association loci (P<1×10-6) in the BioVU cohort. All association analyses used logistic regression, assuming an additive genetic model, and adjusted for principal components to control for population stratification. An inverse-variance fixed-effect model was adopted for meta-analysis. In addition, we examined the top variants to evaluate their associations with each criterion in the phenotyping algorithm. We used the STRING database to characterize protein-protein interaction network. RESULTS: Using the same algorithm based on the Rotterdam criteria, we identified 2995 patients with polycystic ovary syndrome and 53,599 population controls in total (2742 cases and 51,438 controls from the discovery phase; 253 cases and 2161 controls in the validation phase). We identified 1 novel genome-wide significant variant rs17186366 (odds ratio [OR]=1.37 [1.23, 1.54], P=2.8×10-8) located near SOD2. In addition, 2 loci with suggestive association were also identified: rs113168128 (OR=1.72 [1.42, 2.10], P=5.2×10-8), an intronic variant of ERBB4 that is independent from the previously published variants, and rs144248326 (OR=2.13 [1.52, 2.86], P=8.45×10-7), a novel intronic variant in WWTR1. In the further association tests of the top 3 single-nucleotide polymorphisms with each criterion in the polycystic ovary syndrome algorithm, we found that rs17186366 (SOD2) was associated with polycystic ovaries and hyperandrogenism, whereas rs11316812 (ERBB4) and rs144248326 (WWTR1) were mainly associated with oligomenorrhea or infertility. We also validated the previously reported association with DENND1A1. Using the STRING database to characterize protein-protein interactions, we found both ERBB4 and WWTR1 can interact with YAP1, which has been previously associated with polycystic ovary syndrome. CONCLUSION: Through a discovery-validation genome-wide association study on polycystic ovary syndrome identified from electronic health records using an algorithm based on Rotterdam criteria, we identified and validated a novel genome-wide significant association with a variant near SOD2. We also identified a novel independent variant within ERBB4 and a suggestive association with WWTR1. With previously identified polycystic ovary syndrome gene YAP1, the ERBB4-YAP1-WWTR1 network suggests involvement of the epidermal growth factor receptor and the Hippo pathway in the multifactorial etiology of polycystic ovary syndrome.

Gene co-regulation and co-expression in the aryl hydrocarbon receptor-mediated transcriptional regulatory network in the mouse liver.

Four decades after its discovery, the aryl hydrocarbon receptor (AHR), a ligand-inducible transcription factor (TF) activated by the persistent environmental contaminant 2,3,7,8-Tetrachlorodibenzo-p-dioxin (TCDD), remains an enigmatic molecule with a controversial endogenous role. Here, we have assembled a global map of the AHR gene regulatory network in female C57BL/6 mice orally gavaged with 30 µg/kg of TCDD from a combination of previously published gene expression and genome-wide TF-binding data sets. Using Kohonen self-organizing maps and subspace clustering, we show that genes co-regulated by common upstream TFs in the AHR network exhibit a pattern of co-expression. Directly bound, indirectly bound, and non-genomic AHR target genes exhibit distinct expression patterns, with the directly bound targets associated with highest median expression. Interestingly, among the directly bound AHR target genes, the expression level increases with the number of AHR-binding sites in the proximal promoter regions. Finally, we show that co-regulated genes in the AHR network activate distinct groups of downstream biological processes. Although the specific findings described here are restricted to hepatic effects under short-term TCDD exposure, this work describes a generalizable approach to the reconstruction and analysis of transcriptional regulatory cascades underlying cellular stress response, revealing network hierarchy and the nature of information flow from the initial signaling events to phenotypic outcomes. Such reconstructed networks can form the basis of a new generation of quantitative adverse outcome pathways.

WHOLE GENOME SEQUENCING ANALYSIS OF BODY MASS INDEX IDENTIFIES NOVEL AFRICAN ANCESTRY-SPECIFIC RISK ALLELE.

Obesity is a major public health crisis associated with high mortality rates. Previous genome-wide association studies (GWAS) investigating body mass index (BMI) have largely relied on imputed data from European individuals. This study leveraged whole-genome sequencing (WGS) data from 88,873 participants from the Trans-Omics for Precision Medicine (TOPMed) Program, of which 51% were of non-European population groups. We discovered 18 BMI-associated signals (P < 5 × 10-9). Notably, we identified and replicated a novel low frequency single nucleotide polymorphism (SNP) in MTMR3 that was common in individuals of African descent. Using a diverse study population, we further identified two novel secondary signals in known BMI loci and pinpointed two likely causal variants in the POC5 and DMD loci. Our work demonstrates the benefits of combining WGS and diverse cohorts in expanding current catalog of variants and genes confer risk for obesity, bringing us one step closer to personalized medicine.

Exome-wide evaluation of rare coding variants using electronic health records identifies new gene-phenotype associations.

The clinical impact of rare loss-of-function variants has yet to be determined for most genes. Integration of DNA sequencing data with electronic health records (EHRs) could enhance our understanding of the contribution of rare genetic variation to human disease1. By leveraging 10,900 whole-exome sequences linked to EHR data in the Penn Medicine Biobank, we addressed the association of the cumulative effects of rare predicted loss-of-function variants for each individual gene on human disease on an exome-wide scale, as assessed using a set of diverse EHR phenotypes. After discovering 97 genes with exome-by-phenome-wide significant phenotype associations (P < 10-6), we replicated 26 of these in the Penn Medicine Biobank, as well as in three other medical biobanks and the population-based UK Biobank. Of these 26 genes, five had associations that have been previously reported and represented positive controls, whereas 21 had phenotype associations not previously reported, among which were genes implicated in glaucoma, aortic ectasia, diabetes mellitus, muscular dystrophy and hearing loss. These findings show the value of aggregating rare predicted loss-of-function variants into 'gene burdens' for identifying new gene-disease associations using EHR phenotypes in a medical biobank. We suggest that application of this approach to even larger numbers of individuals will provide the statistical power required to uncover unexplored relationships between rare genetic variation and disease phenotypes.

Genome-wide meta-analysis identifies 127 open-angle glaucoma loci with consistent effect across ancestries.

Primary open-angle glaucoma (POAG), is a heritable common cause of blindness world-wide. To identify risk loci, we conduct a large multi-ethnic meta-analysis of genome-wide association studies on a total of 34,179 cases and 349,321 controls, identifying 44 previously unreported risk loci and confirming 83 loci that were previously known. The majority of loci have broadly consistent effects across European, Asian and African ancestries. Cross-ancestry data improve fine-mapping of causal variants for several loci. Integration of multiple lines of genetic evidence support the functional relevance of the identified POAG risk loci and highlight potential contributions of several genes to POAG pathogenesis, including SVEP1, RERE, VCAM1, ZNF638, CLIC5, SLC2A12, YAP1, MXRA5, and SMAD6. Several drug compounds targeting POAG risk genes may be potential glaucoma therapeutic candidates.

Characterization of genetic and phenotypic heterogeneity of obstructive sleep apnea using electronic health records.

BACKGROUND: Obstructive sleep apnea (OSA) is defined by frequent episodes of reduced or complete cessation of airflow during sleep and is linked to negative health outcomes. Understanding the genetic factors influencing expression of OSA may lead to new treatment strategies. Electronic health records (EHRs) can be leveraged to both validate previously reported OSA-associated genomic variation and detect novel relationships between these variants and comorbidities. METHODS: We identified candidate single nucleotide polymorphisms (SNPs) via systematic literature review of existing research. Using datasets available at Geisinger (n = 39,407) and Vanderbilt University Medical Center (n = 24,084), we evaluated associations between 40 previously implicated SNPs and OSA diagnosis, defined using clinical codes. We also evaluated associations between these SNPs and OSA severity measures obtained from sleep reports at Geisinger (n = 6571). Finally, we used a phenome-wide association study approach to help reveal pleiotropic genetic effects between OSA candidate SNPs and other clinical codes and laboratory values available in the EHR. RESULTS: Most previously reported OSA candidate SNPs showed minimal to no evidence for associations with OSA diagnosis or severity in the EHR-derived datasets. Three SNPs in LEPR, MMP-9, and GABBR1 validated for an association with OSA diagnosis in European Americans; the SNP in GABBR1 was associated following meta-analysis of results from both clinical populations. The GABBR1 and LEPR SNPs, and one additional SNP, were associated with OSA severity measures in European Americans from Geisinger. Three additional candidate OSA SNPs were not associated with OSA-related traits but instead with hyperlipidemia and autoimmune diseases of the thyroid. CONCLUSIONS: To our knowledge, this is one of the largest candidate gene studies and one of the first phenome-wide association studies of OSA genomic variation. Results validate genetic associates with OSA in the LEPR, MMP-9 and GABBR1 genes, but suggest that the majority of previously identified genetic associations with OSA may be false positives. Phenome-wide analyses provide evidence of mediated pleiotropy. Future well-powered genome-wide association analyses of OSA risk and severity across populations with diverse ancestral backgrounds are needed. The comprehensive nature of the analyses represents a platform for informing future work focused on understanding how genetic data can be useful to informing treatment of OSA and related comorbidities.

Author Correction: Rare variants in drug target genes contributing to complex diseases, phenome-wide.

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

Healthcare Utilization and Patients' Perspectives After Receiving a Positive Genetic Test for Familial Hypercholesterolemia.

BACKGROUND: The MyCode Community Health Initiative (MyCode) is returning actionable results from whole exome sequencing. Familial hypercholesterolemia (FH) is an inherited condition characterized by premature cardiovascular disease. METHODS: We used multiple methods to assess care in 28 MyCode participants who received FH results. Chart reviews were conducted on 23 individuals in the sample and 7 individuals participated semistructured interviews. RESULTS: Chart reviews for 23 individuals with a Geisinger primary care provider found that 4 individuals (17% of 23) were at LDL-C (low-density lipoprotein cholesterol) goal (of either LDL-C <100 mg/dL for primary prevention and LDL-C <70 mg/dL for secondary prevention) and 17 individuals (74% of 23) were prescribed lipid-lowering therapy before genetic result disclosure. After disclosure of the genetic test result, 5 individuals (22% of 23) met their LDL-C goal and 18 individuals (78% of 23) were prescribed lipid-lowering therapy. Follow-up care about this result was not documented for 4 individuals (17% of 23). Changes to intensity of medication management were made for 8 individuals (47% of 17 individuals previously prescribed lipid-lowering therapy). Interviewed individuals (n=7) were not surprised by their result as all knew they had high cholesterol; however, individuals did not seem to discern FH as a separate condition from their high cholesterol. CONCLUSIONS: Among individuals receiving genetic diagnosis of FH, >25% had no changes to lipid-lowering therapy, despite not being at LDL-C goal and learning their high cholesterol is related to a genetic condition requiring more aggressive treatment. Individuals and clinicians may have an inadequate understanding of FH as a distinct condition requiring enhanced medical management.

Rare variants in drug target genes contributing to complex diseases, phenome-wide.

The DrugBank database consists of ~800 genes that are well characterized drug targets. This list of genes is a useful resource for association testing. For example, loss of function (LOF) genetic variation has the potential to mimic the effect of drugs, and high impact variation in these genes can impact downstream traits. Identifying novel associations between genetic variation in these genes and a range of diseases can also uncover new uses for the drugs that target these genes. Phenome Wide Association Studies (PheWAS) have been successful in identifying genetic associations across hundreds of thousands of diseases. We have conducted a novel gene based PheWAS to test the effect of rare variants in DrugBank genes, evaluating associations between these genes and more than 500 quantitative and dichotomous phenotypes. We used whole exome sequencing data from 38,568 samples in Geisinger MyCode Community Health Initiative. We evaluated the results of this study when binning rare variants using various filters based on potential functional impact. We identified multiple novel associations, and the majority of the significant associations were driven by functionally annotated variation. Overall, this study provides a sweeping exploration of rare variant associations within functionally relevant genes across a wide range of diagnoses.