Ritchie, Marylyn D.

Apart from ancestry, personal or environmental covariates may contribute to differences in polygenic score (PGS) performance. We analyzed the effects of covariate stratification and interaction on body mass index (BMI) PGS (PGS BMI ) across four cohorts of European (N = 491,111) and African (N = 21,612) ancestry. Stratifying on binary covariates and quintiles for continuous covariates, 18/62 covariates had significant and replicable R 2 differences among strata. Covariates with the largest differences included age, sex, blood lipids, physical activity, and alcohol consumption, with R 2 being nearly double between best- and worst-performing quintiles for certain covariates. Twenty-eight covariates had significant PGS BMI –covariate interaction effects, modifying PGS BMI effects by nearly 20% per standard deviation change. We observed overlap between covariates that had significant R 2 differences among strata and interaction effects – across all covariates, their main effects on BMI were correlated with their maximum R 2 differences and interaction effects (0.56 and 0.58, respectively), suggesting high-PGS BMI individuals have highest R 2 and increase in PGS effect. Using quantile regression, we show the effect of PGS BMI increases as BMI itself increases, and that these differences in effects are directly related to differences in R 2 when stratifying by different covariates. Given significant and replicable evidence for context-specific PGS BMI performance and effects, we investigated ways to increase model performance taking into account nonlinear effects. Machine learning models (neural networks) increased relative model R 2 (mean 23%) across datasets. Finally, creating PGS BMI directly from GxAge genome-wide association studies effects increased relative R 2 by 7.8%. These results demonstrate that certain covariates, especially those most associated with BMI, significantly affect both PGS BMI performance and effects across diverse cohorts and ancestries, and we provide avenues to improve model performance that consider these effects.

Abstract Pharmacogenomics (PGx) is focused on the relationship between an individual's genetic makeup and their response to medications, with the overarching aim of guiding prescribing decisions to improve drug efficacy and reduce adverse events. The PGx and genomic medicine communities have worked independently for over 2 decades, developing separate standards and terminology, making implementation of PGx across all areas of genomic medicine difficult. To address this issue, the Clinical Genome Resource (ClinGen) Pharmacogenomics Working Group (PGxWG) was established by the National Institutes of Health (NIH)-funded ClinGen to initially create frameworks for evaluating gene–drug response clinical validity and actionability aligned with the ClinGen frameworks for evaluating monogenic gene–disease relationships, and a framework for classifying germline PGx variants similar to the American College of Medical Genetics (ACMG) and Association of Molecular Pathology (AMP) system for interpretation of disease-causing variants. These frameworks will leverage decades of work from well-established PGx resources facilitating buy-in among PGx stakeholders. In this report, we describe the background and major activities of the ClinGen PGxWG, and how this initiative will facilitate the critical inclusion of PGx into the larger context of genomic medicine.

Apart from ancestry, personal or environmental covariates may contribute to differences in polygenic score (PGS) performance. We analyzed the effects of covariate stratification and interaction on body mass index (BMI) PGS (PGS BMI ) across four cohorts of European (N = 491,111) and African (N = 21,612) ancestry. Stratifying on binary covariates and quintiles for continuous covariates, 18/62 covariates had significant and replicable R 2 differences among strata. Covariates with the largest differences included age, sex, blood lipids, physical activity, and alcohol consumption, with R 2 being nearly double between best- and worst-performing quintiles for certain covariates. Twenty-eight covariates had significant PGS BMI –covariate interaction effects, modifying PGS BMI effects by nearly 20% per standard deviation change. We observed overlap between covariates that had significant R 2 differences among strata and interaction effects – across all covariates, their main effects on BMI were correlated with their maximum R 2 differences and interaction effects (0.56 and 0.58, respectively), suggesting high-PGS BMI individuals have highest R 2 and increase in PGS effect. Using quantile regression, we show the effect of PGS BMI increases as BMI itself increases, and that these differences in effects are directly related to differences in R 2 when stratifying by different covariates. Given significant and replicable evidence for context-specific PGS BMI performance and effects, we investigated ways to increase model performance taking into account nonlinear effects. Machine learning models (neural networks) increased relative model R 2 (mean 23%) across datasets. Finally, creating PGS BMI directly from GxAge genome-wide association studies effects increased relative R 2 by 7.8%. These results demonstrate that certain covariates, especially those most associated with BMI, significantly affect both PGS BMI performance and effects across diverse cohorts and ancestries, and we provide avenues to improve model performance that consider these effects.

Abstract Objective The Genomic Medicine Working Group of the National Advisory Council for Human Genome Research virtually hosted its 13th genomic medicine meeting titled “Developing a Clinical Genomic Informatics Research Agenda”. The meeting’s goal was to articulate a research strategy to develop Genomics-based Clinical Informatics Tools and Resources (GCIT) to improve the detection, treatment, and reporting of genetic disorders in clinical settings. Materials and Methods Experts from government agencies, the private sector, and academia in genomic medicine and clinical informatics were invited to address the meeting's goals. Invitees were also asked to complete a survey to assess important considerations needed to develop a genomic-based clinical informatics research strategy. Results Outcomes from the meeting included identifying short-term research needs, such as designing and implementing standards-based interfaces between laboratory information systems and electronic health records, as well as long-term projects, such as identifying and addressing barriers related to the establishment and implementation of genomic data exchange systems that, in turn, the research community could help address. Discussion Discussions centered on identifying gaps and barriers that impede the use of GCIT in genomic medicine. Emergent themes from the meeting included developing an implementation science framework, defining a value proposition for all stakeholders, fostering engagement with patients and partners to develop applications under patient control, promoting the use of relevant clinical workflows in research, and lowering related barriers to regulatory processes. Another key theme was recognizing pervasive biases in data and information systems, algorithms, access, value, and knowledge repositories and identifying ways to resolve them.

Abstract Coronary artery disease (CAD) is a leading cause of death, yet its genetic determinants are not fully elucidated. We report a multi-ethnic genome-wide association study of CAD involving nearly a quarter of a million cases, incorporating the largest cohorts to date of Whites, Blacks, and Hispanics from the Million Veteran Program with existing studies including CARDIoGRAMplusC4D, UK Biobank, and Biobank Japan. We verify substantial and nearly equivalent heritability of CAD across multiple ancestral groups, discover 107 novel loci including the first nine on the X-chromosome, identify the first eight genome-wide significant loci among Blacks and Hispanics, and demonstrate that two common haplotypes are largely responsible for the risk stratification at the well-known 9p21 locus in most populations except those of African origin where both haplotypes are virtually absent. We identify 15 loci for angiographically derived burden of coronary atherosclerosis, which robustly overlap with the strongest and earliest loci reported to date for clinical CAD. Phenome-wide association analyses of novel loci and externally validated polygenic risk scores (PRS) augment signals from the insulin resistance cluster of risk factors and consequences, extend previously established pleiotropic associations of loci with traditional risk factors to include smoking and family history, and confirm a substantially reduced transferability of existing PRS to Blacks. Downstream integrative genomic analyses reinforce the critical role of endothelial, fibroblast, and smooth muscle cells within the coronary vessel wall in CAD susceptibility. Our study highlights the value of a multi-ethnic design in efficiently characterizing the genetic architecture of CAD across all human populations.

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.

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.

Background: Proteomic approaches allow measurement of thousands of proteins in a single specimen, which can accelerate biomarker discovery. However, applying these technologies to massive biobanks is not currently feasible because of the practical barriers and costs of implementing such assays at scale. To overcome these challenges, we used a “virtual proteomic” approach, linking genetically predicted protein levels to clinical diagnoses in >40 000 individuals. Methods: We used genome-wide association data from the Framingham Heart Study (n=759) to construct genetic predictors for 1129 plasma protein levels. We validated the genetic predictors for 268 proteins and used them to compute predicted protein levels in 41 288 genotyped individuals in the Electronic Medical Records and Genomics (eMERGE) cohort. We tested associations for each predicted protein with 1128 clinical phenotypes. Lead associations were validated with directly measured protein levels and either low-density lipoprotein cholesterol or subclinical atherosclerosis in the MDCS (Malmö Diet and Cancer Study; n=651). Results: In the virtual proteomic analysis in eMERGE, 55 proteins were associated with 89 distinct diagnoses at a false discovery rate q<0.1. Among these, 13 associations involved lipid (n=7) or atherosclerosis (n=6) phenotypes. We tested each association for validation in MDCS using directly measured protein levels. At Bonferroni-adjusted significance thresholds, levels of apolipoprotein E isoforms were associated with hyperlipidemia, and circulating C-type lectin domain family 1 member B and platelet-derived growth factor receptor-β predicted subclinical atherosclerosis. Odds ratios for carotid atherosclerosis were 1.31 (95% CI, 1.08–1.58; P =0.006) per 1-SD increment in C-type lectin domain family 1 member B and 0.79 (0.66–0.94; P =0.008) per 1-SD increment in platelet-derived growth factor receptor-β. Conclusions: We demonstrate a biomarker discovery paradigm to identify candidate biomarkers of cardiovascular and other diseases.