Search

Supplemental Digital Content is available in the text., Background: Atrial fibrillation (AF) is associated with substantial morbidity, especially when it goes undetected. If new-onset AF could be predicted, targeted screening could be used to find it early. We hypothesized that a deep neural network could predict new-onset AF from the resting 12-lead ECG and that this prediction may help identify those at risk of AF-related stroke. Methods: We used 1.6 M resting 12-lead digital ECG traces from 430 000 patients collected from 1984 to 2019. Deep neural networks were trained to predict new-onset AF (within 1 year) in patients without a history of AF. Performance was evaluated using areas under the receiver operating characteristic curve and precision-recall curve. We performed an incidence-free survival analysis for a period of 30 years following the ECG stratified by model predictions. To simulate real-world deployment, we trained a separate model using all ECGs before 2010 and evaluated model performance on a test set of ECGs from 2010 through 2014 that were linked to our stroke registry. We identified the patients at risk for AF-related stroke among those predicted to be high risk for AF by the model at different prediction thresholds. Results: The area under the receiver operating characteristic curve and area under the precision-recall curve were 0.85 and 0.22, respectively, for predicting new-onset AF within 1 year of an ECG. The hazard ratio for the predicted high- versus low-risk groups over a 30-year span was 7.2 (95% CI, 6.9–7.6). In a simulated deployment scenario, the model predicted new-onset AF at 1 year with a sensitivity of 69% and specificity of 81%. The number needed to screen to find 1 new case of AF was 9. This model predicted patients at high risk for new-onset AF in 62% of all patients who experienced an AF-related stroke within 3 years of the index ECG. Conclusions: Deep learning can predict new-onset AF from the 12-lead ECG in patients with no previous history of AF. This prediction may help identify patients at risk for AF-related strokes.

Background Heart failure is a prevalent, costly disease for which new value-based payment models demand optimized population management strategies. Objectives This study sought to generate a strategy for managing populations of patients with heart failure by leveraging large clinical datasets and machine learning. Methods Geisinger electronic health record data were used to train machine learning models to predict 1-year all-cause mortality in 26,971 patients with heart failure who underwent 276,819 clinical episodes. There were 26 clinical variables (demographics, laboratory test results, medications), 90 diagnostic codes, 41 electrocardiogram measurements and patterns, 44 echocardiographic measurements, and 8 evidence-based “care gaps”: flu vaccine, blood pressure of Results Machine learning models achieved areas under the receiver-operating characteristic curve (AUCs) of 0.74 to 0.77 in a split-by-year training/test scheme, with the nonlinear XGBoost model (AUC: 0.77) outperforming linear logistic regression (AUC: 0.74). Out of 13,238 currently living patients, 2,844 were predicted to die within a year, and closing all care gaps was predicted to save 231 of these lives. Prioritizing patients for intervention by using the predicted reduction in 1-year mortality risk outperformed all other priority rankings (e.g., random selection or Seattle Heart Failure risk score). Conclusions Machine learning can be used to priority-rank patients most likely to benefit from interventions to optimize evidence-based therapies. This approach may prove useful for optimizing heart failure population health management teams within value-based payment models.

Background: Atrial fibrillation (AF) risk estimation using clinical factors with or without genetic information may identify AF screening candidates more accurately than the guideline-based age threshold of ≥65 years. Methods: We analyzed 4 samples across the United States and Europe (derivation: UK Biobank; validation: FINRISK, Geisinger MyCode Initiative, and Framingham Heart Study). We estimated AF risk using the CHARGE-AF (Cohorts for Heart and Aging Research in Genomic Epidemiology AF) score and a combination of CHARGE-AF and a 1168-variant polygenic score (Predict-AF). We compared the utility of age, CHARGE-AF, and Predict-AF for predicting 5-year AF by quantifying discrimination and calibration. Results: Among 543 093 individuals, 8940 developed AF within 5 years. In the validation sets, CHARGE-AF (C index range, 0.720–0.824) and Predict-AF (0.749–0.831) had largely comparable discrimination, both favorable to continuous age (0.675–0.801). Calibration was similar using CHARGE-AF (slope range, 0.67–0.87) and Predict-AF (0.65–0.83). Net reclassification improvement using Predict-AF versus CHARGE-AF was modest (net reclassification improvement range, 0.024–0.057) but more favorable among individuals aged Conclusions: AF risk estimation using clinical factors may prioritize individuals for AF screening more precisely than the age threshold endorsed in current guidelines. The additional value of genetic predisposition is modest but greatest among younger individuals.

Background Mutations in PKD1 or PKD2 cause typical autosomal dominant polycystic kidney disease (ADPKD), the most common monogenic kidney disease. Dominantly inherited polycystic kidney and liver diseases on the ADPKD spectrum are also caused by mutations in at least six other genes required for protein biogenesis in the endoplasmic reticulum, the loss of which results in defective production of the PKD1 gene product, the membrane protein polycystin-1 (PC1). Methods We used whole-exome sequencing in a cohort of 122 patients with genetically unresolved clinical diagnosis of ADPKD or polycystic liver disease to identify a candidate gene, ALG9, and in vitro cell-based assays of PC1 protein maturation to functionally validate it. For further validation, we identified carriers of ALG9 loss-of-function mutations and noncarrier matched controls in a large exome-sequenced population-based cohort and evaluated the occurrence of polycystic phenotypes in both groups. Results Two patients in the clinically defined cohort had rare loss-of-function variants in ALG9, which encodes a protein required for addition of specific mannose molecules to the assembling N-glycan precursors in the endoplasmic reticulum lumen. In vitro assays showed that inactivation of Alg9 results in impaired maturation and defective glycosylation of PC1. Seven of the eight (88%) cases selected from the population-based cohort based on ALG9 mutation carrier state who had abdominal imaging after age 50; seven (88%) had at least four kidney cysts, compared with none in matched controls without ALG9 mutations. Conclusions ALG9 is a novel disease gene in the genetically heterogeneous ADPKD spectrum. This study supports the utility of phenotype characterization in genetically-defined cohorts to validate novel disease genes, and provide much-needed genotype-phenotype correlations.

Objectives The goal of this study was to use machine learning to more accurately predict survival after echocardiography. Background Predicting patient outcomes (e.g., survival) following echocardiography is primarily based on ejection fraction (EF) and comorbidities. However, there may be significant predictive information within additional echocardiography-derived measurements combined with clinical electronic health record data. Methods Mortality was studied in 171,510 unselected patients who underwent 331,317 echocardiograms in a large regional health system. We investigated the predictive performance of nonlinear machine learning models compared with that of linear logistic regression models using 3 different inputs: 1) clinical variables, including 90 cardiovascular-relevant International Classification of Diseases, Tenth Revision, codes, and age, sex, height, weight, heart rate, blood pressures, low-density lipoprotein, high-density lipoprotein, and smoking; 2) clinical variables plus physician-reported EF; and 3) clinical variables and EF, plus 57 additional echocardiographic measurements. Missing data were imputed with a multivariate imputation by using a chained equations algorithm (MICE). We compared models versus each other and baseline clinical scoring systems by using a mean area under the curve (AUC) over 10 cross-validation folds and across 10 survival durations (6 to 60 months). Results Machine learning models achieved significantly higher prediction accuracy (all AUC >0.82) over common clinical risk scores (AUC = 0.61 to 0.79), with the nonlinear random forest models outperforming logistic regression (p Conclusions Machine learning can fully utilize large combinations of disparate input variables to predict survival after echocardiography with superior accuracy.

Machine learning promises to assist physicians with predictions of mortality and of other future clinical events by learning complex patterns from historical data, such as longitudinal electronic health records. Here we show that a convolutional neural network trained on raw pixel data in 812,278 echocardiographic videos from 34,362 individuals provides superior predictions of one-year all-cause mortality. The model’s predictions outperformed the widely used pooled cohort equations, the Seattle Heart Failure score (measured in an independent dataset of 2,404 patients with heart failure who underwent 3,384 echocardiograms), and a machine learning model involving 58 human-derived variables from echocardiograms and 100 clinical variables derived from electronic health records. We also show that cardiologists assisted by the model substantially improved the sensitivity of their predictions of one-year all-cause mortality by 13% while maintaining prediction specificity. Large unstructured datasets may enable deep learning to improve a wide range of clinical prediction models. A deep learning model trained on raw pixel data in hundreds of thousands of echocardiographic videos for the prediction of one-year all-cause mortality outperforms clinical scores and improves predictions by cardiologists.

Aims We investigated the relationship between clinically assessed left ventricular ejection fraction (LVEF) and survival in a large, heterogeneous clinical cohort. Methods and results Physician-reported LVEF on 403 977 echocardiograms from 203 135 patients were linked to all-cause mortality using electronic health records (1998–2018) from US regional healthcare system. Cox proportional hazards regression was used for analyses while adjusting for many patient characteristics including age, sex, and relevant comorbidities. A dataset including 45 531 echocardiograms and 35 976 patients from New Zealand was used to provide independent validation of analyses. During follow-up of the US cohort, 46 258 (23%) patients who had undergone 108 578 (27%) echocardiograms died. Overall, adjusted hazard ratios (HR) for mortality showed a u-shaped relationship for LVEF with a nadir of risk at an LVEF of 60–65%, a HR of 1.71 [95% confidence interval (CI) 1.64–1.77] when ≥70% and a HR of 1.73 (95% CI 1.66–1.80) at LVEF of 35–40%. Similar relationships with a nadir at 60–65% were observed in the validation dataset as well as for each age group and both sexes. The results were similar after further adjustments for conditions associated with an elevated LVEF, including mitral regurgitation, increased wall thickness, and anaemia and when restricted to patients reported to have heart failure at the time of the echocardiogram. Conclusion Deviation of LVEF from 60% to 65% is associated with poorer survival regardless of age, sex, or other relevant comorbidities such as heart failure. These results may herald the recognition of a new phenotype characterized by supra-normal LVEF.

Type 2 diabetes mellitus (T2DM) is a major health and economic burden because of the seriousness of the disease and its complications. Improvements in short- and long-term glycemic control is the goal of diabetes treatment. To investigate the longitudinal management of T2DM at Geisinger, we interrogated the electronic health record (EHR) information and identified a T2DM cohort including 125,477 patients using the Electronic Medical Records and Genomics Network (eMERGE) T2DM phenotyping algorithm. We investigated the annual anti-diabetic medication usage and the overall glycemic control using hemoglobin A1c (HbA1c). Metformin remains the most frequently medication despite the availability of the new classes of anti-diabetic medications. Median value of HbA1c decreased to 7% in 2002 and since remained stable, indicating a good glycemic management in Geisinger population. Using metformin as a pilot study, we identified three groups of patients with distinct HbA1c trajectories after metformin treatment. The variabilities in metformin response is mainly explained by the baseline HbA1c. The pharmacogenomic analysis of metformin identified a missense variant rs75740279 (Leu/Val) for STAU2 associated with the metformin response. This strategy can be applied to study other anti-diabeticmedications. Such research will facilitate the translational healthcare for better T2DM management.

Population genomic screening has been demonstrated to detect at-risk individuals who would not be clinically identified otherwise. However, there are concerns about the increased utilization of unnecessary services and the associated increase in costs. The objectives of this study are twofold: (1) determine whether there is a difference in healthcare utilization and costs following disclosure of a pathogenic/likely pathogenic (P/LP) BRCA1/2 variant via a genomic screening program, and (2) measure the post-disclosure uptake of National Comprehensive Cancer Network (NCCN) guideline-recommended risk management. We retrospectively reviewed electronic health record (EHR) and billing data from a female population of BRCA1/2 P/LP variant carriers without a personal history of breast or ovarian cancer enrolled in Geisinger&rsquo, s MyCode genomic screening program with at least a one-year post-disclosure observation period. We identified 59 women for the study cohort out of 50,726 MyCode participants. We found no statistically significant differences in inpatient and outpatient utilization and average total costs between one-year pre- and one-year post-disclosure periods ($18,821 vs. $19,359, p = 0.76). During the first year post-disclosure, 49.2% of women had a genetic counseling visit, 45.8% had a mammography and 32.2% had an MRI. The uptake of mastectomy and oophorectomy was 3.5% and 11.8%, respectively, and 5% of patients received chemoprevention.

BackgroundPolycystic ovary syndrome (PCOS) is the most common endocrine disorder affecting women of reproductive age. Previous studies have identified genetic variants associated with PCOS identified by different diagnostic criteria. The Rotterdam Criteria is the broadest and able to identify the most PCOS cases.ObjectivesTo identify novel associated genetic variants, we extracted PCOS cases and controls from the electronic health records (EHR) based on the Rotterdam Criteria and performed a genome-wide association study (GWAS).Study DesignWe developed a PCOS phenotyping algorithm based on the Rotterdam criteria and applied it to three EHR-linked biobanks to identify cases and controls for genetic study. In discovery phase, we performed individual GWAS using the Geisinger’s MyCode and the eMERGE cohorts, which were then meta-analyzed. We attempted validation of the significantly association loci (P−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-analyses. Additionally, we examined the top variants to evaluate their associations with each criterion in the phenotyping algorithm. We used STRING to identify protein-protein interaction network.ResultsWe identified 2,995 PCOS cases and 53,599 controls in total (2,742cases and 51,438 controls from the discovery phase; 253 cases and 2,161 controls in the validation phase). GWAS identified one novel genome-wide significant variant rs17186366 (OR=1.37 [1.23,1.54], P=2.8×10−8) located nearSOD2. Additionally, two loci with suggestive association were also identified: rs113168128 (OR=1.72 [1.42,2.10], P=5.2 x10−8), an intronic variant ofERBB4that 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 inWWTR1. In the further association tests of the top 3 SNPs with each criterion in the PCOS algorithm, we found that rs17186366 was associated with polycystic and hyperandrogenism, while rs11316812 and rs144248326 were mainly associated with oligomenorrhea or infertility. Besides ERBB4, we also validated the association withDENND1A1.ConclusionThrough a discovery-validation GWAS on PCOS cases and controls identified from EHR using an algorithm based on Rotterdam criteria, we identified and validated a novel association with variants withinERBB4. We also identified novel associations nearbySOD2andWWTR1. These results suggest the eGFR and Hippo pathways in the disease etiology. With previously identified PCOS-associated lociYAP1, theERBB4-YAP1-WWTR1network implicates the epidermal growth factor receptor and the Hippo pathway in the multifactorial etiology of PCOS.

Background: Arrhythmogenic right ventricular cardiomyopathy (ARVC) is associated with variants in desmosome genes. Secondary findings of pathogenic/likely pathogenic variants, primarily loss-of-function (LOF) variants, are recommended for clinical reporting; however, their prevalence and associated phenotype in a general clinical population are not fully characterized. Methods: From whole-exome sequencing of 61 019 individuals in the DiscovEHR cohort, we screened for putative loss-of-function variants in PKP2 , DSC2 , DSG2 , and DSP . We evaluated measures from prior clinical ECG and echocardiograms, manually over-read to evaluate ARVC diagnostic criteria, and performed a PheWAS (phenome-wide association study). Finally, we estimated expected penetrance using Bayesian inference. Results: One hundred forty individuals (0.23%; 59±18 years old at last encounter; 33% male) had an ARVC variant (G + ). None had an existing diagnosis of ARVC in the electronic health record, nor significant differences in prior ECG or echocardiogram findings compared with matched controls without variants. Several G + individuals satisfied major repolarization (n=4) and ventricular function (n=5) criteria, but this prevalence matched controls. PheWAS showed no significant associations of other heart disease diagnoses. Combining our best genetic and disease prevalence estimates yields an estimated penetrance of 6.0%. Conclusions: The prevalence of ARVC loss-of-function variants is ≈1:435 in a general clinical population of predominantly European descent, but with limited electronic health record-based evidence of phenotypic association in our population, consistent with a low penetrance estimate. Prospective deep phenotyping and longitudinal follow-up of a large sequenced cohort is needed to determine the true clinical relevance of an incidentally identified ARVC loss-of-function variant.

Deletion of glutathione S-transferase µ1 (GSTM1) is common in populations and has been asserted to associate with chronic kidney disease progression in some research studies. The association needs to be validated. We estimated GSTM1 copy number using whole exome sequencing data in the DiscovEHR cohort. Kidney failure was defined as requiring dialysis or receiving kidney transplant using data from the electronic health record and linkage to the United States Renal Data System, or the most recent eGFR < 15 ml/min/1.73m2. In a cohort of 46,983 unrelated participants, 28.8% of blacks and 52.1% of whites had 0 copies of GSTM1. Over a mean of 9.2 years follow-up, 645 kidney failure events were observed in 46,187 white participants, and 28 in 796 black participants. No significant association was observed between GSTM1 copy number and kidney failure in Cox regression adjusting for age, sex, BMI, smoking status, genetic principal components, or co-morbid conditions (hypertension, diabetes, heart failure, coronary artery disease, and stroke), whether using a genotypic, dominant, or recessive model. In sensitivity analyses, GSTM1 copy number was not associated with kidney failure in participants that were 45 years or older at baseline, had baseline eGFR < 60 ml/min per 1.73 m2, or with baseline year between 1996-2002. In conclusion, we found no association between GSTM1 copy number and kidney failure in a large cohort study.Translational StatementDeletion of GSTM1 has been shown to be associated with higher risk of kidney failure. However, inconsistent results have been reported. We used electronic health record and whole exome sequencing data of a large cohort from a single healthcare system to evaluate the association between GSTM1 copy number and risk of kidney failure. We found no significant association between GSTM1 copy number and risk of kidney failure overall, or in multiple sensitivity and subgroup analyses.

Background: Truncating variants in the Titin gene (TTNtvs) are common in individuals with idiopathic dilated cardiomyopathy (DCM). However, a comprehensive genomics-first evaluation of the impact of TTNtvs in different clinical contexts, and the evaluation of modifiers such as genetic ancestry, has not been performed. Methods: We reviewed whole exome sequence data for >71 000 individuals (61 040 from the Geisinger MyCode Community Health Initiative (2007 to present) and 10 273 from the PennMedicine BioBank (2013 to present) to identify anyone with TTNtvs. We further selected individuals with TTNtvs in exons highly expressed in the heart (proportion spliced in [PSI] >0.9). Using linked electronic health records, we evaluated associations of TTNtvs with diagnoses and quantitative echocardiographic measures, including subanalyses for individuals with and without DCM diagnoses. We also reviewed data from the Jackson Heart Study to validate specific analyses for individuals of African ancestry. Results: Identified with a TTNtv in a highly expressed exon (hiPSI) were 1.2% individuals in PennMedicine BioBank and 0.6% at Geisinger. The presence of a hiPSI TTNtv was associated with increased odds of DCM in individuals of European ancestry (odds ratio [95% CI]: 18.7 [9.1–39.4] {PennMedicine BioBank} and 10.8 [7.0–16.0] {Geisinger}). hiPSI TTNtvs were not associated with DCM in individuals of African ancestry, despite a high DCM prevalence (odds ratio, 1.8 [0.2–13.7]; P =0.57). Among 244 individuals of European ancestry with DCM in PennMedicine BioBank, hiPSI TTNtv carriers had lower left ventricular ejection fraction (β=–12%, P =3×10 –7 ), and increased left ventricular diameter (β=0.65 cm, P =9×10 –3 ). In the Geisinger cohort, hiPSI TTNtv carriers without a cardiomyopathy diagnosis had more atrial fibrillation (odds ratio, 2.4 [1.6–3.6]) and heart failure (odds ratio, 3.8 [2.4–6.0]), and lower left ventricular ejection fraction (β=–3.4%, P =1×10 –7 ). Conclusions: Individuals of European ancestry with hiPSI TTNtv have an abnormal cardiac phenotype characterized by lower left ventricular ejection fraction, irrespective of the clinical manifestation of cardiomyopathy. Associations with arrhythmias, including atrial fibrillation, were observed even when controlling for cardiomyopathy diagnosis. In contrast, no association between hiPSI TTNtvs and DCM was discerned among individuals of African ancestry. Given these findings, clinical identification of hiPSI TTNtv carriers may alter clinical management strategies.

The electrocardiogram (ECG) is a widely used medical test, consisting of voltage versus time traces collected from surface recordings over the heart1. Here we hypothesized that a deep neural network (DNN) can predict an important future clinical event, 1-year all-cause mortality, from ECG voltage–time traces. By using ECGs collected over a 34-year period in a large regional health system, we trained a DNN with 1,169,662 12-lead resting ECGs obtained from 253,397 patients, in which 99,371 events occurred. The model achieved an area under the curve (AUC) of 0.88 on a held-out test set of 168,914 patients, in which 14,207 events occurred. Even within the large subset of patients (n = 45,285) with ECGs interpreted as ‘normal’ by a physician, the performance of the model in predicting 1-year mortality remained high (AUC = 0.85). A blinded survey of cardiologists demonstrated that many of the discriminating features of these normal ECGs were not apparent to expert reviewers. Finally, a Cox proportional-hazard model revealed a hazard ratio of 9.5 (P

Background Empirical data on conditions that increase risk of coronavirus disease 2019 (COVID-19) progression are needed to identify high risk individuals. We performed a comprehensive quantitative assessment of pre-existing clinical phenotypes associated with COVID-19-related hospitalization. Methods Phenome-wide association study (PheWAS) of SARS-CoV-2-positive patients from an integrated health system (Geisinger) with system-level outpatient/inpatient COVID-19 testing capacity and retrospective electronic health record (EHR) data to assess pre-COVID-19 pandemic clinical phenotypes associated with hospital admission (hospitalization). Results Of 12,971 individuals tested for SARS-CoV-2 with sufficient pre-COVID-19 pandemic EHR data at Geisinger, 1604 were SARS-CoV-2 positive and 354 required hospitalization. We identified 21 clinical phenotypes in 5 disease categories meeting phenome-wide significance (P-4), including: six kidney phenotypes, e.g. end stage renal disease or stage 5 CKD (OR = 11.07, p = 1.96x10-8), six cardiovascular phenotypes, e.g. congestive heart failure (OR = 3.8, p = 3.24x10-5), five respiratory phenotypes, e.g. chronic airway obstruction (OR = 2.54, p = 3.71x10-5), and three metabolic phenotypes, e.g. type 2 diabetes (OR = 1.80, p = 7.51x10-5). Additional analyses defining CKD based on estimated glomerular filtration rate, confirmed high risk of hospitalization associated with pre-existing stage 4 CKD (OR 2.90, 95% CI: 1.47, 5.74), stage 5 CKD/dialysis (OR 8.83, 95% CI: 2.76, 28.27), and kidney transplant (OR 14.98, 95% CI: 2.77, 80.8) but not stage 3 CKD (OR 1.03, 95% CI: 0.71, 1.48). Conclusions This study provides quantitative estimates of the contribution of pre-existing clinical phenotypes to COVID-19 hospitalization and highlights kidney disorders as the strongest factors associated with hospitalization in an integrated US healthcare system.

Background: Heterologous functional validation studies of putative long-QT syndrome subtype 2–associated variants clarify their pathological potential and identify disease mechanism(s) for most variants studied. The purpose of this study is to clarify the pathological potential for rare nonsynonymous KCNH2 variants seemingly associated with sudden infant death syndrome. Methods: Genetic testing of 292 sudden infant death syndrome cases identified 9 KCNH2 variants: E90K, R181Q, A190T, G294V, R791W, P967L, R1005W, R1047L, and Q1068R. Previous studies show R181Q-, P967L-, and R1047L-Kv11.1 channels function similar to wild-type Kv11.1 channels, whereas Q1068R-Kv11.1 channels accelerate inactivation gating. We studied the biochemical and biophysical properties for E90K-, G294V-, R791W-, and R1005W-Kv11.1 channels expressed in human embryonic kidney 293 cells; examined the electronic health records of patients who were genotype positive for the sudden infant death syndrome–linked KCNH2 variants; and simulated their functional impact using computational models of the human ventricular action potential. RESULTS: Western blot and voltage-clamping analyses of cells expressing E90K-, G294V-, R791W-, and R1005W-Kv11.1 channels demonstrated these variants express and generate peak Kv11.1 current levels similar to cells expressing wild-type-Kv11.1 channels, but R791W- and R1005W-Kv11.1 channels accelerated deactivation and activation gating, respectively. Electronic health records of patients with the sudden infant death syndrome–linked KCNH2 variants showed that the patients had median heart rate–corrected QT intervals Conclusions: We conclude that these rare Kv11.1 missense variants are not long-QT syndrome subtype 2–causative variants and therefore do not represent the pathogenic substrate for sudden infant death syndrome in the variant-positive infants.

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.

Importance The activity of lipoprotein lipase (LPL) is the rate-determining step in clearing triglyceride-rich lipoproteins from the circulation. Mutations that damage the LPL gene ( LPL ) lead to lifelong deficiency in enzymatic activity and can provide insight into the relationship of LPL to human disease. Objective To determine whether rare and/or common variants in LPL are associated with early-onset coronary artery disease (CAD). Design, Setting, and Participants In a cross-sectional study, LPL was sequenced in 10 CAD case-control cohorts of the multinational Myocardial Infarction Genetics Consortium and a nested CAD case-control cohort of the Geisinger Health System DiscovEHR cohort between 2010 and 2015. Common variants were genotyped in up to 305 699 individuals of the Global Lipids Genetics Consortium and up to 120 600 individuals of the CARDIoGRAM Exome Consortium between 2012 and 2014. Study-specific estimates were pooled via meta-analysis. Exposures Rare damaging mutations in LPL included loss-of-function variants and missense variants annotated as pathogenic in a human genetics database or predicted to be damaging by computer prediction algorithms trained to identify mutations that impair protein function. Common variants in the LPL gene region included those independently associated with circulating triglyceride levels. Main Outcomes and Measures Circulating lipid levels and CAD. Results Among 46 891 individuals with LPL gene sequencing data available, the mean (SD) age was 50 (12.6) years and 51% were female. A total of 188 participants (0.40%; 95% CI, 0.35%-0.46%) carried a damaging mutation in LPL , including 105 of 32 646 control participants (0.32%) and 83 of 14 245 participants with early-onset CAD (0.58%). Compared with 46 703 noncarriers, the 188 heterozygous carriers of an LPL damaging mutation displayed higher plasma triglyceride levels (19.6 mg/dL; 95% CI, 4.6-34.6 mg/dL) and higher odds of CAD (odds ratio = 1.84; 95% CI, 1.35-2.51; P LPL variants resulted in an odds ratio for CAD of 1.51 (95% CI, 1.39-1.64; P = 1.1 × 10 −22 ) per 1-SD increase in triglycerides. Conclusions and Relevance The presence of rare damaging mutations in LPL was significantly associated with higher triglyceride levels and presence of coronary artery disease. However, further research is needed to assess whether there are causal mechanisms by which heterozygous lipoprotein lipase deficiency could lead to coronary artery disease.

Unleashing the power of precision medicine Precision medicine promises the ability to identify risks and treat patients on the basis of pathogenic genetic variation. Two studies combined exome sequencing results for over 50,000 people with their electronic health records. Dewey et al. found that ∼3.5% of individuals in their cohort had clinically actionable genetic variants. Many of these variants affected blood lipid levels that could influence cardiovascular health. Abul-Husn et al. extended these findings to investigate the genetics and treatment of familial hypercholesterolemia, a risk factor for cardiovascular disease, within their patient pool. Genetic screening helped identify at-risk patients who could benefit from increased treatment. Science , this issue p. 10.1126/science.aaf6814 , p. 10.1126/science.aaf7000

PurposeArrhythmogenic right ventricular cardiomyopathy (ARVC) is an inherited heart disease. Clinical follow-up of incidental findings in ARVC-associated genes is recommended. We aimed to determine the prevalence of disease thus ascertained.MethodsIndividuals (n = 30,716) underwent exome sequencing. Variants in PKP2, DSG2, DSC2, DSP, JUP, TMEM43, or TGFβ3 that were database-listed as pathogenic or likely pathogenic were identified and evidence-reviewed. For subjects with putative loss-of-function (pLOF) variants or variants of uncertain significance (VUS), electronic health records (EHR) were reviewed for ARVC diagnosis, diagnostic criteria, and International Classification of Diseases (ICD-9) codes.ResultsEighteen subjects had pLOF variants; none of these had an EHR diagnosis of ARVC. Of 14 patients with an electrocardiogram, one had a minor diagnostic criterion; the rest were normal. A total of 184 subjects had VUS, none of whom had an ARVC diagnosis. The proportion of subjects with VUS with major (4%) or minor (13%) electrocardiogram diagnostic criteria did not differ from that of variant-negative controls. ICD-9 codes showed no difference in defibrillator use, electrophysiologic abnormalities or nonischemic cardiomyopathies in patients with pLOF or VUSs compared with controls.ConclusionpLOF variants in an unselected cohort were not associated with ARVC phenotypes based on EHR review. The negative predictive value of EHR review remains uncertain.

Unleashing the power of precision medicine Precision medicine promises the ability to identify risks and treat patients on the basis of pathogenic genetic variation. Two studies combined exome sequencing results for over 50,000 people with their electronic health records. Dewey et al. found that ∼3.5% of individuals in their cohort had clinically actionable genetic variants. Many of these variants affected blood lipid levels that could influence cardiovascular health. Abul-Husn et al. extended these findings to investigate the genetics and treatment of familial hypercholesterolemia, a risk factor for cardiovascular disease, within their patient pool. Genetic screening helped identify at-risk patients who could benefit from increased treatment. Science , this issue p. 10.1126/science.aaf6814 , p. 10.1126/science.aaf7000

Key Points Question Can population-level genomic screening identify those at risk for disease? Findings In this cross-sectional study of an unselected population cohort of 50 726 adults who underwent exome sequencing, pathogenic and likely pathogenic BRCA1 and BRCA2 variants were found in a higher proportion of patients than was previously reported. Meaning Current methods to identify BRCA1/2 variant carriers may not be sufficient as a screening tool; population genomic screening for hereditary breast and ovarian cancer may better identify patients at high risk and provide an intervention opportunity to reduce mortality and morbidity., Importance Detection of disease-associated variants in the BRCA1 and BRCA2 (BRCA1/2) genes allows for cancer prevention and early diagnosis in high-risk individuals. Objectives To identify pathogenic and likely pathogenic (P/LP) BRCA1/2 variants in an unselected research cohort, and to characterize the features associated with P/LP variants. Design, Setting, and Participants This is a cross-sectional study of adult volunteers (n = 50 726) who underwent exome sequencing at a single health care system (Geisinger Health System, Danville, Pennsylvania) from January 1, 2014, to March 1, 2016. Participants are part of the DiscovEHR cohort and were identified through the Geisinger MyCode Community Health Initiative. They consented to a research protocol that included sequencing and return of actionable test results. Clinical data from electronic health records and clinical visits were correlated with variants. Comparisons were made between those with (cases) and those without (controls) P/LP variants in BRCA1/2. Main Outcomes Prevalence of P/LP BRCA1/2 variants in cohort, proportion of variant carriers not previously ascertained through clinical testing, and personal and family history of relevant cancers among BRCA1/2 variant carriers and noncarriers. Results Of the 50 726 health system patients who underwent exome sequencing, 50 459 (99.5%) had no expected pathogenic BRCA1/2 variants and 267 (0.5%) were BRCA1/2 carriers. Of the 267 cases (148 [55.4%] were women and 119 [44.6%] were men with a mean [range] age of 58.9 [23-90] years), 183 (68.5%) received clinically confirmed results in their electronic health record. Among the 267 participants with P/LP BRCA1/2 variants, 219 (82.0%) had no prior clinical testing, 95 (35.6%) had BRCA1 variants, and 172 (64.4%) had BRCA2 variants. Syndromic cancer diagnoses were present in 11 (47.8%) of the 23 deceased BRCA1/2 carriers and in 56 (20.9%) of all 267 BRCA1/2 carriers. Among women, 31 (20.9%) of 148 variant carriers had a personal history of breast cancer, compared with 1554 (5.2%) of 29 880 noncarriers (odds ratio [OR], 5.95; 95% CI, 3.88-9.13; P, This cross-sectional study investigates the ascertainment of pathogenic and likely pathogenic BRCA1/2 variants among US adults enrolled in the MyCode Community Health Initiative.