Denny, Joshua C.

AbstractObjectiveThe aim of this study was to analyze a publicly available sample of rule-based phenotype definitions to characterize and evaluate the variability of logical constructs used.Materials and MethodsA sample of 33 preexisting phenotype definitions used in research that are represented using Fast Healthcare Interoperability Resources and Clinical Quality Language (CQL) was analyzed using automated analysis of the computable representation of the CQL libraries.ResultsMost of the phenotype definitions include narrative descriptions and flowcharts, while few provide pseudocode or executable artifacts. Most use 4 or fewer medical terminologies. The number of codes used ranges from 5 to 6865, and value sets from 1 to 19. We found that the most common expressions used were literal, data, and logical expressions. Aggregate and arithmetic expressions are the least common. Expression depth ranges from 4 to 27.DiscussionDespite the range of conditions, we found that all of the phenotype definitions consisted of logical criteria, representing both clinical and operational logic, and tabular data, consisting of codes from standard terminologies and keywords for natural language processing. The total number and variety of expressions are low, which may be to simplify implementation, or authors may limit complexity due to data availability constraints.ConclusionsThe phenotype definitions analyzed show significant variation in specific logical, arithmetic, and other operators but are all composed of the same high-level components, namely tabular data and logical expressions. A standard representation for phenotype definitions should support these formats and be modular to support localization and shared logic.

AbstractObjectiveHigh BMI is associated with many comorbidities and mortality. This study aimed to elucidate the overall clinical risk of obesity using a genome‐ and phenome‐wide approach.MethodsThis study performed a phenome‐wide association study of BMI using a clinical cohort of 736,726 adults. This was followed by genetic association studies using two separate cohorts: one consisting of 65,174 adults in the Electronic Medical Records and Genomics (eMERGE) Network and another with 405,432 participants in the UK Biobank.ResultsClass 3 obesity was associated with 433 phenotypes, representing 59.3% of all billing codes in individuals with severe obesity. A genome‐wide polygenic risk score for BMI, accounting for 7.5% of variance in BMI, was associated with 296 clinical diseases, including strong associations with type 2 diabetes, sleep apnea, hypertension, and chronic liver disease. In all three cohorts, 199 phenotypes were associated with class 3 obesity and polygenic risk for obesity, including novel associations such as increased risk of renal failure, venous insufficiency, and gastroesophageal reflux.ConclusionsThis combined genomic and phenomic systematic approach demonstrated that obesity has a strong genetic predisposition and is associated with a considerable burden of disease across all disease classes.

Background: Sequencing Mendelian arrhythmia genes in individuals without an indication for arrhythmia genetic testing can identify carriers of pathogenic or likely pathogenic (P/LP) variants. However, the extent to which these variants are associated with clinically meaningful phenotypes before or after return of variant results is unclear. In addition, the majority of discovered variants are currently classified as variants of uncertain significance, limiting clinical actionability. Methods: The eMERGE-III study (Electronic Medical Records and Genomics Phase III) is a multicenter prospective cohort that included 21 846 participants without previous indication for cardiac genetic testing. Participants were sequenced for 109 Mendelian disease genes, including 10 linked to arrhythmia syndromes. Variant carriers were assessed with electronic health record–derived phenotypes and follow-up clinical examination. Selected variants of uncertain significance (n=50) were characterized in vitro with automated electrophysiology experiments in HEK293 cells. Results: As previously reported, 3.0% of participants had P/LP variants in the 109 genes. Herein, we report 120 participants (0.6%) with P/LP arrhythmia variants. Compared with noncarriers, arrhythmia P/LP carriers had a significantly higher burden of arrhythmia phenotypes in their electronic health records. Fifty-four participants had variant results returned. Nineteen of these 54 participants had inherited arrhythmia syndrome diagnoses (primarily long-QT syndrome), and 12 of these 19 diagnoses were made only after variant results were returned (0.05%). After in vitro functional evaluation of 50 variants of uncertain significance, we reclassified 11 variants: 3 to likely benign and 8 to P/LP. Conclusions: Genome sequencing in a large population without indication for arrhythmia genetic testing identified phenotype-positive carriers of variants in congenital arrhythmia syndrome disease genes. As the genomes of large numbers of people are sequenced, the disease risk from rare variants in arrhythmia genes can be assessed by integrating genomic screening, electronic health record phenotypes, and in vitro functional studies. Registration: URL: https://www.clinicaltrials.gov ; Unique identifier; NCT03394859.

Abstract Background With limited severe acute respiratory syndrome coronavirus (SARS-CoV-2) testing capacity in the United States at the start of the epidemic (January–March 2020), testing was focused on symptomatic patients with a travel history throughout February, obscuring the picture of SARS-CoV-2 seeding and community transmission. We sought to identify individuals with SARS-CoV-2 antibodies in the early weeks of the US epidemic. Methods All of Us study participants in all 50 US states provided blood specimens during study visits from 2 January to 18 March 2020. Participants were considered seropositive if they tested positive for SARS-CoV-2 immunoglobulin G (IgG) antibodies with the Abbott Architect SARS-CoV-2 IgG enzyme-linked immunosorbent assay (ELISA) and the EUROIMMUN SARS-CoV-2 ELISA in a sequential testing algorithm. The sensitivity and specificity of these ELISAs and the net sensitivity and specificity of the sequential testing algorithm were estimated, along with 95% confidence intervals (CIs). Results The estimated sensitivities of the Abbott and EUROIMMUN assays were 100% (107 of 107 [95% CI: 96.6%–100%]) and 90.7% (97 of 107 [83.5%–95.4%]), respectively, and the estimated specificities were 99.5% (995 of 1000 [98.8%–99.8%]) and 99.7% (997 of 1000 [99.1%–99.9%]), respectively. The net sensitivity and specificity of our sequential testing algorithm were 90.7% (97 of 107 [95% CI: 83.5%–95.4%]) and 100.0% (1000 of 1000 [99.6%–100%]), respectively. Of the 24 079 study participants with blood specimens from 2 January to 18 March 2020, 9 were seropositive, 7 before the first confirmed case in the states of Illinois, Massachusetts, Wisconsin, Pennsylvania, and Mississippi. Conclusions Our findings identified SARS-CoV-2 infections weeks before the first recognized cases in 5 US states.

Background: Return of incidental genetic findings is recommended for pathogenic/likely pathogenic (P/LP) variants in Mendelian arrhythmia genes. The extent to which these variants are associated with arrhythmia phenotypes in unselected populations is unknown. Objective: Assess the impact of return of results (RoR) of P/LP variants in 4 arrhythmia genes in research participants not selected for cardiovascular disease. Methods: The cohort included 24,410 participants from the Electronic Medical Records and Genomics Network. Rare variants (minor allele frequency <1%) in KCNE1 , KCNH2 , KCNQ1 , and SCN5A were identified and classified according to ACMG guidelines. Arrhythmia phenotypes extracted from electronic health records (EHRs) included atrial fibrillation/flutter (AF/AFL), conduction system disease, long QT syndrome (LQTS), ventricular tachycardia/fibrillation (VT/VF), and premature ventricular contractions (PVCs). P/LP carriers’ phenotypes were compared to non-carriers using logistic regression adjusted for demographic covariates and site. Participants with P/LP variants were informed of their results and the impact of this RoR was assessed after 1 year of follow-up clinical care. Results: The participants included 53.7% females and 72.5% white individuals, with a median age of 57 years. 71 participants had a heterozygous P/LP variant in one of the 4 genes. LQTS diagnoses were associated with P/LP carrier status for all 4 genes ( KCNE1 OR 22.8, p=1.0e-4; KCNH2 OR 35.1, p=1.8e-7; KCNQ1 OR 27.3, p=3.6e-18; SCN5A OR 5.6, p=0.02). Carriers of SCN5A P/LP variants were also more likely to have PVCs (OR 4.3, p=0.003). Sensitivity analyses showed the associations predated the RoR. EHRs included a diagnosis of LQTS or Brugada Syndrome in 11/71 (16%) pre-RoR, and 16 more were diagnosed after RoR (27/71, 38%). A total of 47/71 (66%) participants had a documented arrhythmia phenotype 1 year after RoR, although 20 of these participants were not diagnosed with an inherited arrhythmia syndrome. Conclusion: Inherited arrhythmia diagnoses could be established in nearly half of participants with P/LP variants. These data indicate that return of incidental P/LP variants in these genes may facilitate the early diagnosis of arrhythmias.

Background: Abdominal aortic aneurysm (AAA) is an important cause of cardiovascular mortality; however, its genetic determinants remain incompletely defined. In total, 10 previously identified risk loci explain a small fraction of AAA heritability. Methods: We performed a genome-wide association study in the Million Veteran Program testing ≈18 million DNA sequence variants with AAA (7642 cases and 172 172 controls) in veterans of European ancestry with independent replication in up to 4972 cases and 99 858 controls. We then used mendelian randomization to examine the causal effects of blood pressure on AAA. We examined the association of AAA risk variants with aneurysms in the lower extremity, cerebral, and iliac arterial beds, and derived a genome-wide polygenic risk score (PRS) to identify a subset of the population at greater risk for disease. Results: Through a genome-wide association study, we identified 14 novel loci, bringing the total number of known significant AAA loci to 24. In our mendelian randomization analysis, we demonstrate that a genetic increase of 10 mm Hg in diastolic blood pressure (odds ratio, 1.43 [95% CI, 1.24–1.66]; P =1.6×10 −6 ), as opposed to systolic blood pressure (odds ratio, 1.06 [95% CI, 0.97–1.15]; P =0.2), likely has a causal relationship with AAA development. We observed that 19 of 24 AAA risk variants associate with aneurysms in at least 1 other vascular territory. A 29-variant PRS was strongly associated with AAA (odds ratio PRS , 1.26 [95% CI, 1.18–1.36]; P PRS =2.7×10 −11 per SD increase in PRS), independent of family history and smoking risk factors (odds ratio PRS+family history+smoking , 1.24 [95% CI, 1.14–1.35]; P PRS =1.27×10 −6 ). Using this PRS, we identified a subset of the population with AAA prevalence greater than that observed in screening trials informing current guidelines. Conclusions: We identify novel AAA genetic associations with therapeutic implications and identify a subset of the population at significantly increased genetic risk of AAA independent of family history. Our data suggest that extending current screening guidelines to include testing to identify those with high polygenic AAA risk, once the cost of genotyping becomes comparable with that of screening ultrasound, would significantly increase the yield of current screening at reasonable cost.

Abstract Background: Results from epidemiologic studies examining polyunsaturated fatty acids (PUFA) and colorectal cancer risk are inconsistent. Mendelian randomization may strengthen causal inference from observational studies. Given their shared metabolic pathway, examining the combined effects of aspirin/NSAID use with PUFAs could help elucidate an association between PUFAs and colorectal cancer risk. Methods: Information was leveraged from genome-wide association studies (GWAS) regarding PUFA-associated SNPs to create weighted genetic scores (wGS) representing genetically predicted circulating blood PUFAs for 11,016 non-Hispanic white colorectal cancer cases and 13,732 controls in the Genetics and Epidemiology of Colorectal Cancer Consortium (GECCO). Associations per SD increase in the wGS were estimated using unconditional logistic regression. Interactions between PUFA wGSs and aspirin/NSAID use on colorectal cancer risk were also examined. Results: Modest colorectal cancer risk reductions were observed per SD increase in circulating linoleic acid [ORLA = 0.96; 95% confidence interval (CI) = 0.93–0.98; P = 5.2 × 10−4] and α-linolenic acid (ORALA = 0.95; 95% CI = 0.92–0.97; P = 5.4 × 10−5), whereas modest increased risks were observed for arachidonic (ORAA = 1.06; 95% CI = 1.03–1.08; P = 3.3 × 10−5), eicosapentaenoic (OREPA = 1.04; 95% CI = 1.01–1.07; P = 2.5 × 10−3), and docosapentaenoic acids (ORDPA = 1.03; 95% CI = 1.01–1.06; P = 1.2 × 10−2). Each of these effects was stronger among aspirin/NSAID nonusers in the stratified analyses. Conclusions: Our study suggests that higher circulating shorter-chain PUFAs (i.e., LA and ALA) were associated with reduced colorectal cancer risk, whereas longer-chain PUFAs (i.e., AA, EPA, and DPA) were associated with an increased colorectal cancer risk. Impact: The interaction of PUFAs with aspirin/NSAID use indicates a shared colorectal cancer inflammatory pathway. Future research should continue to improve PUFA genetic instruments to elucidate the independent effects of PUFAs on colorectal cancer.

Background The incidence of VTE in African-Americans (AAs) is similar to or higher than in Americans of European ancestry. However, the carrier frequencies of inherited thrombophilias common in whites (i.e., Factor V Leiden, Prothrombin G20210A) are very low in AAs, suggesting that other inherited thrombophilias may be associated with VTE in AAs. Objective To identify potentially novel non-hypothesis-driven single nucleotide polymorphisms (SNPs) associated with VTE in AAs. Methods We used the resources of the eMERGE Network to perform a GWAS of VTE in AAs. The eMERGE Network (funded by NHGRI) is comprised of nine sites each with DNA biobanks linked to electronic health records (EHRs). Approximately 39,206 unique DNA samples have been genotyped using either Affymetrix or Illumina genome-wide SNP arrays. Led by the Coordinating Center and eMERGE genomics workgroup, an imputation pipeline was developed for merging genomic data across the different SNP arrays used by the eMERGE sites, to maximize sample size and the power to detect associations. The imputation was performed using the 1000 Genomes Cosmopolitan reference panel which includes 1092 individuals and over 36 million SNPs. From our previously-identified cohort of Olmsted County, MN residents with incident or recurrent VTE, 1996-2005, we derived and validated an EHR-driven phenotype extraction algorithm that leveraged structured data based on ICD-9-CM codes and unstructured data from clinical notes via natural language processing to identify VTE cases and controls with 100% and 94% positive and negative predictive values, respectively. We tested for an association between each SNP and VTE among AAs using unconditional logistic regression, adjusting for age, sex and eMERGE site. Results Among 294 AA VTE cases and 3,661 AA controls (total n=3,955; females, n=2,512), the Factor V Leiden (F5 rs6025) was not analyzed due to a very low minor allele frequency (MAF=0.0036). The prothrombin G20210A (F2 rs1799963) and ABO blood type O (ABO rs8176719) SNPs were not genotyped in any of the arrays and could not be imputed. Among SNPs with an imputation score >0.8, the most significant SNPs associated with VTE were ITPR3 (inositol 1,4,5-triphosphate receptor type 3) rs2229637 (OR=1.65; p=3.61E-07; MAF=0.19) and CLEC7A (C-type lectin domain family 7, member A) rs59819090 (OR=2.16; p=1.06E-06; MAF=0.056). ITPR3 SNPs have been associated with coronary artery aneurysm in Kawasaki disease, type 1 diabetes mellitus and other autoimmune disorders. CLEC7A rs59819090 encodes for a serine82 to leucine nonsynonymous amino acid change in dectin-1. Dectin-1 is a transmembrane innate immune pattern recognition receptor on myeloid cells that upon binding it’s agonist, β-glucans, stimulates cytokine production and the respiratory burst, a prerequisite for formation of neutrophil extracelluar traps (NETs). NETs have been associated with VTE (van Montfoort, et al. Ateriosclero Thromb Vasc Biol 2013;33:147-51). Dectin-1 serine82 is not evolutionarily conserved and the serine82 to leucine amino acid change is predicted to be tolerated, with a moderate effect on protein function. Conclusions ITPR3 rs2229637 and CLAC7A rs59819090 may be associated with VTE in African-Americans. These observations require replication and functional studies toward understanding how they may lie on the causal pathway to VTE. Disclosures: No relevant conflicts of interest to declare.

AbstractBackgroundThe extent to which obesity and genetics determine postoperative complications is incompletely understood.MethodsWe performed a retrospective study using two population cohorts with electronic health record (EHR) data. The first included 736,726 adults with body mass index (BMI) recorded between 1990 and 2017 at Vanderbilt University Medical Center. The second cohort consisted of 65,174 individuals from 12 institutions contributing EHR and genome‐wide genotyping data to the Electronic Medical Records and Genomics (eMERGE) Network. Pairwise logistic regression analyses were used to measure the association of BMI categories with postoperative complications derived from International Classification of Disease‐9 codes, including postoperative infection, incisional hernia, and intestinal obstruction. A genetic risk score was constructed from 97 obesity‐risk single‐nucleotide polymorphisms for a Mendelian randomization study to determine the association of genetic risk of obesity on postoperative complications. Logistic regression analyses were adjusted for sex, age, site, and race/principal components.ResultsIndividuals with overweight or obese BMI (≥25 kg/m2) had increased risk of incisional hernia (odds ratio [OR] 1.7–5.5, p < 3.1 × 10−20), and people with obesity (BMI ≥ 30 kg/m2) had increased risk of postoperative infection (OR 1.2–2.3, p < 2.5 × 10−5). In the eMERGE cohort, genetically predicted BMI was associated with incisional hernia (OR 2.1 [95% CI 1.8–2.5], p = 1.4 × 10−6) and postoperative infection (OR 1.6 [95% CI 1.4–1.9], p = 3.1 × 10−6). Association findings were similar after limitation of the cohorts to those who underwent abdominal procedures.ConclusionsClinical and Mendelian randomization studies suggest that obesity, as measured by BMI, is associated with the development of postoperative incisional hernia and infection.

PURPOSEDrug development is becoming increasingly expensive and time consuming. Drug repurposing is one potential solution to accelerate drug discovery. However, limited research exists on the use of electronic health record (EHR) data for drug repurposing, and most published studies have been conducted in a hypothesis-driven manner that requires a predefined hypothesis about drugs and new indications. Whether EHRs can be used to detect drug repurposing signals is not clear. We want to demonstrate the feasibility of mining large, longitudinal EHRs for drug repurposing by detecting candidate noncancer drugs that can potentially be used for the treatment of cancer.PATIENTS AND METHODSBy linking cancer registry data to EHRs, we identified 43,310 patients with cancer treated at Vanderbilt University Medical Center (VUMC) and 98,366 treated at the Mayo Clinic. We assessed the effect of 146 noncancer drugs on cancer survival using VUMC EHR data and sought to replicate significant associations (false discovery rate < .1) using the identical approach with Mayo Clinic EHR data. To evaluate replicated signals further, we reviewed the biomedical literature and clinical trials on cancers for corroborating evidence.RESULTSWe identified 22 drugs from six drug classes (statins, proton pump inhibitors, angiotensin-converting enzyme inhibitors, β-blockers, nonsteroidal anti-inflammatory drugs, and α-1 blockers) associated with improved overall cancer survival (false discovery rate < .1) from VUMC; nine of the 22 drug associations were replicated at the Mayo Clinic. Literature and cancer clinical trial evaluations also showed very strong evidence to support the repurposing signals from EHRs.CONCLUSIONMining of EHRs for drug exposure–mediated survival signals is feasible and identifies potential candidates for antineoplastic repurposing. This study sets up a new model of mining EHRs for drug repurposing signals.