Search

Background: Limited information is available regarding clinician and participant behaviors after disclosure of genomic risk variants for familial hypercholesterolemia (FH) from a population genomic screening program. Methods: We conducted a retrospective cohort study of MyCode participants with an FH risk variant beginning 2 years before disclosure until January 16, 2019. We analyzed lipid-lowering prescriptions (clinician behavior), medication adherence (participant behavior), and LDL (low-density lipoprotein) cholesterol levels (health outcome impact) pre- and post-disclosure. Data were collected from electronic health records and claims. Results: The cohort included 96 participants of mean age 57 (22–90) years with median follow-up of 14 (range, 3–39) months. Most (90%) had a hypercholesterolemia diagnosis but no specific FH diagnosis before disclosure; 29% had an FH diagnosis post-disclosure. After disclosure, clinicians made 36 prescription changes in 38% of participants, mostly in participants who did not achieve LDL cholesterol goals pre-disclosure (81%). However, clinicians wrote prescriptions for fewer participants post-disclosure (71/96, 74.0%) compared with pre-disclosure (81/96, 84.4%); side effects were documented for most discontinued prescriptions (23/25, 92%). Among the 16 participants with claims data, medication adherence improved (proportion of days covered pre-disclosure of 70% [SD, 24.7%] to post-disclosure of 79.1% [SD, 27.3%]; P =0.05). Among the 52 (54%) participants with LDL cholesterol values both before and after disclosure, average LDL cholesterol decreased from 147 to 132 mg/dL ( P =0.003). Conclusions: Despite disclosure of an FH risk variant, nonprescribing and nonadherence to lipid-lowering therapy remained high. However, when clinicians intensified medication regimens and participants adhered to medications, lipid levels decreased.

Several large trials have employed age or clinical features to select patients for atrial fibrillation (AF) screening to reduce strokes. We hypothesized that a machine learning (ML) model trained to predict AF risk from 12‑lead electrocardiogram (ECG) would be more efficient than criteria based on clinical variables in indicating a population for AF screening to potentially prevent AF-related stroke.We retrospectively included all patients with clinical encounters in Geisinger without a prior history of AF. Incidence of AF within 1 year and AF-related strokes within 3 years of the encounter were identified. AF-related stroke was defined as a stroke where AF was diagnosed at the time of stroke or within a year after the stroke. The efficiency of five methods was evaluated for selecting a cohort for AF screening. The methods were selected from four clinical trials (mSToPS, GUARD-AF, SCREEN-AF and STROKESTOP) and the ECG-based ML model. We simulated patient selection for the five methods between the years 2011 and 2014 and evaluated outcomes for 1 year intervals between 2012 and 2015, resulting in a total of twenty 1-year periods. Patients were considered eligible if they met the criteria before the start of the given 1-year period or within that period. The primary outcomes were numbers needed to screen (NNS) for AF and AF-associated stroke.The clinical trial models indicated large proportions of the population with a prior ECG for AF screening (up to 31%), coinciding with NNS ranging from 14 to 18 for AF and 249-359 for AF-associated stroke. At comparable sensitivity, the ECG ML model indicated a modest number of patients for screening (14%) and had the highest efficiency in NNS for AF (7.3; up to 60% reduction) and AF-associated stroke (223; up to 38% reduction).An ECG-based ML risk prediction model is more efficient than contemporary AF-screening criteria based on age alone or age and clinical features at indicating a population for AF screening to potentially prevent 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.

Introduction: Healthcare cost has increased drastically in the last decade, and over 50% of the cost can be attributed to a small portion (5-10%) of the population. Certain clinical programs, such as home-based care, aim to reduce this utilization but need methods to identify the most appropriate patients to enroll. We hypothesized that machine learning can predict patients with high future utilization with good accuracy. Methods: 683,160 cardiology patients (defined broadly as those with an ECG, echocardiogram or cardiology visit) with ~17 million clinical episodes since 2004 were identified from Geisinger’s electronic health records. Utilization was estimated as total cost of care for outpatient, inpatient and emergency department visits. Patients with the highest 10% utilization in a given year were defined as high utilizers. Machine learning models were used to predict high utilization over the next 3, 6 and 12 months. Input variables (n=191) included age, sex, smoking, 5 vital signs, 21 labs, 18 medications (current and past), 40 ECG and 44 echocardiographic measurements, 43 comorbidities, 7 time / cyclical features, 6 past utilization metrics and 4 social metrics (e.g. distance to healthcare facilities). Results: XGBoost achieved the best performance with areas under the ROC curve (AUC) of 0.82, 0.81 and 0.78 for 3, 6, 12-month models, and average precision scores (AP) of 0.31, 0.36 and 0.37, respectively, while the commonly used Charlson Comorbidity Index had poor performance with AUCs of 0.63 - 0.64 and APs of 0.1 - 0.17. Past utilization was the best predictor of future utilization. Targeting patients with the top 5 and 10% highest risk for utilization achieved sensitivities of 26 and 40% and positive predictive values of 50 and 38% (12-month model, Figure). Conclusions: Machine learning can be used to predict which patients will have high future healthcare utilization. This may help target populations for intervention programs aimed at reducing utilization.

Introduction: Managing heart failure patients relies on mortality prediction models such as the Seattle Heart Failure (SHF) model. We hypothesized that a deep learning model that leverages echocardiography video data could improve prediction performance. Methods: We trained deep neural networks (DNN) to predict 1-year all-cause mortality using echocardiography videos and 158 additional clinical variables (labs, vitals, diagnoses, etc.) acquired from 34,362 patients (812,278 videos). We tested both the DNN and the SHF model on a separate cohort of 2,404 patients with heart failure who underwent 3,384 echocardiograms. We computed the area under the receiver operating characteristic curve (AUC) and its bootstrapped 95% confidence interval (CI) in the test set. Results: The DNN model (AUC of 0.76, 95% CI [0.74, 0.77]) outperformed the SHF model (AUC of 0.70, 95% CI [0.68, 0.71]). This superior performance was observed for patients with both reduced (n=2,026, AUC of 0.76 [0.74, 0.78] vs 0.70 [0.67, 0.72]) and preserved left ventricular ejection fraction (n=1,356, AUC of 0.75 [0.72, 0.78] vs 0.69 [0.66, 0.72]). A Cox Proportional Hazard survival analysis showed that, despite the prediction being for 1-year all-cause mortality, the result held long-term predictive power over the next 9 years with a superior hazard ratio of 2.9 [2.6, 3.2] for the DNN model compared to 2.2 [2, 2.4] for the SHF model. At mid-range operating points, the DNN model also maintained a higher negative predictive value, predicting survival, compared to the SHF model (89% vs 83%, respectively), while maintaining the same positive predictive value of 40%. Conclusion: A deep neural network that automatically analyzes echocardiography videos can outperform traditional risk modeling approaches such as the Seattle Heart Failure Model. This provides additional evidence that deep learning can help improve our ability to make clinically relevant predictions by leveraging complex datasets.

Background: Atrial fibrillation (AF) is associated with stroke, especially when AF goes undetected. Deep neural networks (DNN) can predict incident AF from a 12-lead resting ECG. We hypothesize that use of a DNN to predict new onset AF from an ECG may identify patients at risk of sustaining a potentially preventable AF-related stroke. Methods: We trained a DNN model to predict new-onset AF using 382,604 ECGs prior to 2010. We then evaluated the model performance on a test set of ECGs from 2010 through 2014 linked to patients in an institutional stroke registry. There were 181,969 patients in the test set with at least one ECG and no prior history of AF. Of those patients 3,497 (1.9%) had a stroke following an ECG that did not show AF. Within the set of patients with stroke, 375 had the stroke within 3 years of the ECG and were diagnosed with new AF between -3 and 365 days of the stroke. We considered these potentially preventable AF-related strokes. We report the sensitivity and positive predictive value (PPV) of the model for appropriately risk stratifying these 375 patients who sustained a potentially preventable AF-related stroke. Results: We used F β scores to identify different risk prediction thresholds (operating points) for the model. Operating points chosen by F 0.5 , F 1 , and F 2 scores identified 4, 12, and 21% of the population as high risk for the development of AF within 1 year (Figure 1). Screening 1, 4, 12, and 21% of the overall population resulted in PPV of 28, 21, 15, and 12%, respectively, for identification of new onset AF in one year. Using those same thresholds yielded sensitivities of 4, 17, 45, and 62% for identifying potentially preventable AF-related strokes. The different risk prediction thresholds resulted in a low (120-162) number needed to screen to detect one potentially preventable AF-related stroke at 3 years. Conclusions: Use of a deep learning model to predict new onset AF may identify patients at high risk of sustaining a potentially preventable AF-related stroke.

BackgroundAtrial fibrillation (AF) is associated with substantial morbidity, especially when it goes undetected. If new onset AF could be predicted, targeted population 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 electrocardiogram (ECG) and that this prediction may help prevent AF-related stroke.MethodsWe used 1.6M resting 12-lead ECG voltage-time traces from 430k patients collected from 1984-2019 in this study. 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 (AUROC) and precision-recall curve (AUPRC). 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 prior to 2010 and evaluated model performance on a test set of ECGs from 2010 through 2014 that were linked to our stroke registry. We used standard metrics to explore different prediction thresholds for the model and also calculated how many AF-related strokes might be potentially prevented.ResultsThe AUROC and AUPRC were 0.83 and 0.21, respectively, for predicting new onset AF within 1 year of an ECG. Adding age and sex improved the AUROC to 0.85 and the AUPRC to 0.23. The hazard ratio for the predicted high- vs. low-risk groups over a 30-year span was 7.2 [95% confidence interval: 6.9 – 7.6]. In a simulated deployment scenario, using the F2 score to select the risk prediction threshold, the model predicted new onset AF at 1 year with a sensitivity of 69%, specificity of 81%, and positive predictive value (PPV) of 12%. This model correctly predicted new onset AF in 62% of all patients that experienced an AF-related stroke within 3 years of the ECG.ConclusionsDeep learning can predict new onset AF from the 12-lead ECG in patients with no prior history of AF. This prediction may prove useful in preventing AF-related strokes.

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.

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.

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’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.

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.

PURPOSE: “Genome-first” approaches, in which genetic sequencing is agnostically linked to associated phenotypes, can enhance our understanding of rare variants’ contributions to disease. Loss-of-function variants in LMNA cause a range of rare diseases, including cardiomyopathy. METHODS: We leveraged exome sequencing from 11,451 unselected individuals in the Penn Medicine Biobank to associate rare variants in LMNA with diverse electronic health record (EHR)–derived phenotypes. We used Rare Exome Variant Ensemble Learner (REVEL) to annotate rare missense variants, clustered predicted deleterious and loss-of-function variants into a “gene burden” (N = 72 individuals), and performed a phenome-wide association study (PheWAS). Major findings were replicated in DiscovEHR. RESULTS: The LMNA gene burden was significantly associated with primary cardiomyopathy (p = 1.78E-11) and cardiac conduction disorders (p = 5.27E-07). Most patients had not been clinically diagnosed with LMNA cardiomyopathy. We also noted an association with chronic kidney disease (p = 1.13E-06). Regression analyses on echocardiography and serum labs revealed that LMNA variant carriers had dilated cardiomyopathy and primary renal disease. CONCLUSION: Pathogenic LMNA variants are an underdiagnosed cause of cardiomyopathy. We also find that LMNA loss of function may be a primary cause of renal disease. Finally, we show the value of aggregating rare, annotated variants into a gene burden and using PheWAS to identify novel ontologies for pleiotropic human genes.

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.

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 nonsynonymousGenetic testing of 292 sudden infant death syndrome cases identified 9Western 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-linkedWe 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.

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

Phenome-Wide Association Studies (PheWAS) comprehensively investigate the association between genetic variation and a wide array of outcome traits. Electronic health record (EHR) based PheWAS uses various abstractions of International Classification of Diseases, Ninth Revision (ICD-9) codes to identify case/control status for diagnoses that are used as the phenotypic variables. However, there have not been comparisons within a PheWAS between results from high quality derived phenotypes and high-throughput but potentially inaccurate use of ICD-9 codes for case/control definition. For this study we first developed a group of high quality algorithms for five phenotypes. Next we evaluated the association of these “gold standard” phenotypes and 4,636,178 genetic variants with minor allele frequency > 0.01 and compared the results from high-throughput associations at the 3 digit, 5 digit, and PheWAS codes for defining case/control status. We found that certain diseases contained similar patient populations across phenotyping methods but had differences in PheWAS.

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.

Objective Arrhythmogenic right ventricular cardiomyopathy (ARVC) is a rare myocardial disorder. Genetic etiology is ascribed to variants in 8 genes ( PKP2 , DSP , DSC2 , DSG2 , JUP , TMEM43 , TGFβ2 , RYR2 ). We aimed to establish the prevalence and clinical phenotype of individuals with pathogenic variants (PV) in these genes from an unselected cohort. Methods Whole-exome sequences for 31,036 patients in the MyCode™ Biorepository were evaluated for pathogenic or likely pathogenic (P/LP) ARVC variants reported in ClinVar or the ARVC database. Composite ICD-9 data from the electronic health records of PV individuals vs a non-PV control group were analyzed. Results 57 P/LP variants were identified in 240 individuals (0.77%), aged 61 ± 19 years. Characteristics of this PV group versus a non-PV control group are summarized in the table. Zero patients in the PV group had a documented diagnosis of ARVC. Use of implantable cardioverter defibrillators—the primary treatment for ARVC—was rare. Of individuals ≥55 years, 24% with PV had no history of cardiac disease (vs. 22% for non-PV). However, there was a modest increase in the prevalence of cardiac electrical abnormalities (e.g., tachycardia, fibrillation) in PV vs. non-PV. There was also a trend towards increased diagnoses of 'other primary cardiomyopathy'. Conclusions The prevalence of reported pathogenic ARVC variants (0.77%) greatly exceeds the disease prevalence (0.02 - 0.10%). In this cohort of 240 genotype-positive subjects, there was little evidence of overt ARVC disease, suggesting that the positive-predictive value for classic ARVC in individuals with an incidental positive genetic finding is low. The increased risk of certain clinical findings, not specific to ARVC, supports the notion that the spectrum of phenotypes associated with pathogenic variants in these genes is broader than classically described. Targeted “deep phenotyping” is needed to assess subclinical phenotypes in genomically screened populations.