Development of electronic health record based algorithms to identify individuals with diabetic retinopathy.

OBJECTIVES: To develop, validate, and implement algorithms to identify diabetic retinopathy (DR) cases and controls from electronic health care records (EHRs). MATERIALS AND METHODS: We developed and validated electronic health record (EHR)-based algorithms to identify DR cases and individuals with type I or II diabetes without DR (controls) in 3 independent EHR systems: Vanderbilt University Medical Center Synthetic Derivative (VUMC), the VA Northeast Ohio Healthcare System (VANEOHS), and Massachusetts General Brigham (MGB). Cases were required to meet 1 of the following 3 criteria: (1) 2 or more dates with any DR ICD-9/10 code documented in the EHR, (2) at least one affirmative health-factor or EPIC code for DR along with an ICD9/10 code for DR on a different day, or (3) at least one ICD-9/10 code for any DR occurring within 24 hours of an ophthalmology examination. Criteria for controls included affirmative evidence for diabetes as well as an ophthalmology examination. RESULTS: The algorithms, developed and evaluated in VUMC through manual chart review, resulted in a positive predictive value (PPV) of 0.93 for cases and negative predictive value (NPV) of 0.91 for controls. Implementation of algorithms yielded similar metrics in VANEOHS (PPV = 0.94; NPV = 0.86) and lower in MGB (PPV = 0.84; NPV = 0.76). In comparison, the algorithm for DR implemented in Phenome-wide association study (PheWAS) in VUMC yielded similar PPV (0.92) but substantially reduced NPV (0.48). Implementation of the algorithms to the Million Veteran Program identified over 62 000 DR cases with genetic data including 14 549 African Americans and 6209 Hispanics with DR. CONCLUSIONS/DISCUSSION: We demonstrate the robustness of the algorithms at 3 separate healthcare centers, with a minimum PPV of 0.84 and substantially improved NPV than existing automated methods. We strongly encourage independent validation and incorporation of features unique to each EHR to enhance algorithm performance for DR cases and controls.

Genetic predictors of blood pressure traits are associated with preeclampsia.

Preeclampsia, a pregnancy complication characterized by hypertension after 20 gestational weeks, is a major cause of maternal and neonatal morbidity and mortality. Mechanisms leading to preeclampsia are unclear; however, there is evidence of high heritability. We evaluated the association of polygenic scores (PGS) for blood pressure traits and preeclampsia to assess whether there is shared genetic architecture. Non-Hispanic Black and White reproductive age females with pregnancy indications and genotypes were obtained from Vanderbilt University's BioVU, Electronic Medical Records and Genomics network, and Penn Medicine Biobank. Preeclampsia was defined by ICD codes. Summary statistics for diastolic blood pressure (DBP), systolic blood pressure (SBP), and pulse pressure (PP) PGS were acquired from Giri et al. Associations between preeclampsia and each PGS were evaluated separately by race and data source before subsequent meta-analysis. Ten-fold cross validation was used for prediction modeling. In 3504 Black and 5009 White included individuals, the rate of preeclampsia was 15.49%. In cross-ancestry meta-analysis, all PGSs were associated with preeclampsia (ORDBP = 1.10, 95% CI 1.02-1.17, p = 7.68 × 10-3; ORSBP = 1.16, 95% CI 1.09-1.23, p = 2.23 × 10-6; ORPP = 1.14, 95% CI 1.07-1.27, p = 9.86 × 10-5). Addition of PGSs to clinical prediction models did not improve predictive performance. Genetic factors contributing to blood pressure regulation in the general population also predispose to preeclampsia.

Genomic Insights into Gestational Weight Gain: Uncovering Tissue-Specific Mechanisms and Pathways.

Increasing gestational weight gain (GWG) is linked to adverse outcomes in pregnant persons and their children. The Early Growth Genetics (EGG) Consortium identified previously genetic variants that could contribute to early, late, and total GWG from fetal and maternal genomes. However, the biologic mechanisms and tissue-Specificity of these variants in GWG is unknown. We evaluated the association between genetically predicted gene expression in five relevant maternal (subcutaneous and visceral adipose, breast, uterus, and whole blood) from GTEx (v7) and fetal (placenta) tissues and early, late, and total GWG using S-PrediXcan. We tested enrichment of pre-defined biological pathways for nominally (P < 0.05) significant associations using the GENE2FUNC module from Functional Mapping and Annotation of Genome-Wide Association Studies. After multiple testing correction, we did not find significant associations between maternal and fetal gene expression and early, late, or total GWG. There was significant enrichment of several biological pathways, including metabolic processes, secretion, and intracellular transport, among nominally significant genes from the maternal analyses (false discovery rate p-values: 0.016 to 9.37×10). Enriched biological pathways varied across pregnancy. Though additional research is necessary, these results indicate that diverse biological pathways are likely to impact GWG, with their influence varying by tissue and weeks of gestation.

Soluble Glycoprotein VI Predicts Abdominal Aortic Aneurysm Growth Rate and is a Novel Therapeutic Target.

A common feature in patients with abdominal aortic aneurysms (AAA) is the formation of a nonocclusive intraluminal thrombus (ILT) in regions of aortic dilation. Platelets are known to maintain hemostasis and propagate thrombosis through several redundant activation mechanisms, yet the role of platelet activation in the pathogenesis of AAA associated ILT is still poorly understood. Thus, we sought to investigate how platelet activation impacts the pathogenesis of AAA. Using RNA-sequencing, we identify that the platelet-associated transcripts are significantly enriched in the ILT compared to the adjacent aneurysm wall and healthy control aortas. We found that the platelet specific receptor glycoprotein VI (GPVI) is among the top enriched genes in AAA ILT and is increased on the platelet surface of AAA patients. Examination of a specific indicator of platelet activity, soluble GPVI (sGPVI), in two independent AAA patient cohorts is highly predictive of a AAA diagnosis and associates more strongly with aneurysm growth rate when compared to D-dimer in humans. Finally, intervention with the anti-GPVI antibody (JAQ1) in mice with established aneurysms blunted the progression of AAA in two independent mouse models. In conclusion, we show that levels of sGPVI in humans can predict a diagnosis of AAA and AAA growth rate, which may be critical in the identification of high-risk patients. We also identify GPVI as a novel platelet-specific AAA therapeutic target, with minimal risk of adverse bleeding complications, where none currently exist.

Adaptive selection at G6PD and disparities in diabetes complications.

Diabetes complications occur at higher rates in individuals of African ancestry. Glucose-6-phosphate dehydrogenase deficiency (G6PDdef), common in some African populations, confers malaria resistance, and reduces hemoglobin A1c (HbA1c) levels by shortening erythrocyte lifespan. In a combined-ancestry genome-wide association study of diabetic retinopathy, we identified nine loci including a G6PDdef causal variant, rs1050828 -T (Val98Met), which was also associated with increased risk of other diabetes complications. The effect of rs1050828 -T on retinopathy was fully mediated by glucose levels. In the years preceding diabetes diagnosis and insulin prescription, glucose levels were significantly higher and HbA1c significantly lower in those with versus without G6PDdef. In the Action to Control Cardiovascular Risk in Diabetes (ACCORD) trial, participants with G6PDdef had significantly higher hazards of incident retinopathy and neuropathy. At the same HbA1c levels, G6PDdef participants in both ACCORD and the Million Veteran Program had significantly increased risk of retinopathy. We estimate that 12% and 9% of diabetic retinopathy and neuropathy cases, respectively, in participants of African ancestry are due to this exposure. Across continentally defined ancestral populations, the differences in frequency of rs1050828 -T and other G6PDdef alleles contribute to disparities in diabetes complications. Diabetes management guided by glucose or potentially genotype-adjusted HbA1c levels could lead to more timely diagnoses and appropriate intensification of therapy, decreasing the risk of diabetes complications in patients with G6PDdef alleles.

Overcome the Limitation of Phenome-Wide Association Studies (PheWAS): Extension of PheWAS to Efficient and Robust Large-Scale ICD Codes Analysis.

The Phenome-wide association studies (PheWAS) have become widely used for efficient, high-throughput evaluation of relationship between a genetic factor and a large number of disease phenotypes, typically extracted from a DNA biobank linked with electronic medical records (EMR). Phecodes, billing code-derived disease case-control status, are usually used as outcome variables in PheWAS and logistic regression has been the standard choice of analysis method. Since the clinical diagnoses in EMR are often inaccurate with errors which can lead to biases in the odds ratio estimates, much effort has been put to accurately define the cases and controls to ensure an accurate analysis. Specifically in order to correctly classify controls in the population, an exclusion criteria list for each Phecode was manually compiled to obtain unbiased odds ratios. However, the accuracy of the list cannot be guaranteed without extensive data curation process. The costly curation process limits the efficiency of large-scale analyses that take full advantage of all structured phenotypic information available in EMR. Here, we proposed to estimate relative risks (RR) instead. We first demonstrated the desired nature of RR that overcomes the inaccuracy in the controls via theoretical formula. With simulation and real data application, we further confirmed that RR is unbiased without compiling exclusion criteria lists. With RR as estimates, we are able to efficiently extend PheWAS to a larger-scale, phenome construction agnostic analysis of phenotypes, using ICD 9/10 codes, which preserve much more disease-related clinical information than Phecodes.

A new test for trait mean and variance detects unreported loci for blood-pressure variation.

Variability in quantitative traits has clinical, ecological, and evolutionary significance. Most genetic variants identified for complex quantitative traits have only a detectable effect on the mean of trait. We have developed the mean-variance test (MVtest) to simultaneously model the mean and log-variance of a quantitative trait as functions of genotypes and covariates by using estimating equations. The advantages of MVtest include the facts that it can detect effect modification, that multiple testing can follow conventional thresholds, that it is robust to non-normal outcomes, and that association statistics can be meta-analyzed. In simulations, we show control of type I error of MVtest over several alternatives. We identified 51 and 37 previously unreported associations for effects on blood-pressure variance and mean, respectively, in the UK Biobank. Transcriptome-wide association studies revealed 633 significant unique gene associations with blood-pressure mean variance. MVtest is broadly applicable to studies of complex quantitative traits and provides an important opportunity to detect novel loci.

X-chromosome and kidney function: evidence from a multi-trait genetic analysis of 908,697 individuals reveals sex-specific and sex-differential findings in genes regulated by androgen response elements.

X-chromosomal genetic variants are understudied but can yield valuable insights into sexually dimorphic human traits and diseases. We performed a sex-stratified cross-ancestry X-chromosome-wide association meta-analysis of seven kidney-related traits (n = 908,697), identifying 23 loci genome-wide significantly associated with two of the traits: 7 for uric acid and 16 for estimated glomerular filtration rate (eGFR), including four novel eGFR loci containing the functionally plausible prioritized genes ACSL4, CLDN2, TSPAN6 and the female-specific DRP2. Further, we identified five novel sex-interactions, comprising male-specific effects at FAM9B and AR/EDA2R, and three sex-differential findings with larger genetic effect sizes in males at DCAF12L1 and MST4 and larger effect sizes in females at HPRT1. All prioritized genes in loci showing significant sex-interactions were located next to androgen response elements (ARE). Five ARE genes showed sex-differential expressions. This study contributes new insights into sex-dimorphisms of kidney traits along with new prioritized gene targets for further molecular research.

Session Introduction: Drug-repurposing and discovery in the era of "big" real-world data: how the incorporation of observational data, genetics, and other -omic technologies can move us forward.

This PSB 2024 session discusses the many broad biological, computational, and statistical approaches currently being used for therapeutic drug target identification and repurposing of existing treatments. Drug repurposing efforts have the potential to dramatically improve the treatment landscape by more rapidly identifying drug targets and alternative strategies for untreated or poorly managed diseases. The overarching theme for this session is the use and integration of real-world data to identify drug-disease pairs with potential therapeutic use. These drug-disease pairs may be identified through genomic, proteomic, biomarkers, protein interaction analyses, electronic health records, and chemical profiling. Taken together, this session combines novel applications of methods and innovative modeling strategies with diverse real-world data to suggest new pharmaceutical treatments for human diseases.

EVALUATING THE RELATIONSHIPS BETWEEN GENETIC ANCESTRY AND THE CLINICAL PHENOME.

There is a desire in research to move away from the concept of race as a clinical factor because it is a societal construct used as an imprecise proxy for geographic ancestry. In this study, we leverage the biobank from Vanderbilt University Medical Center, BioVU, to investigate relationships between genetic ancestry proportion and the clinical phenome. For all samples in BioVU, we calculated six ancestry proportions based on 1000 Genomes references: eastern African (EAFR), western African (WAFR), northern European (NEUR), southern European (SEUR), eastern Asian (EAS), and southern Asian (SAS). From PheWAS, we found phecode categories significantly enriched neoplasms for EAFR, WAFR, and SEUR, and pregnancy complication in SEUR, NEUR, SAS, and EAS (p < 0.003). We then selected phenotypes hypertension (HTN) and atrial fibrillation (AFib) to further investigate the relationships between these phenotypes and EAFR, WAFR, SEUR, and NEUR using logistic regression modeling and non-linear restricted cubic spline modeling (RCS). For EAS and SAS, we chose renal failure (RF) for further modeling. The relationships between HTN and AFib and the ancestries EAFR, WAFR, and SEUR were best fit by the linear model (beta p < 1x10-4 for all) while the relationships with NEUR were best fit with RCS (HTN ANOVA p = 0.001, AFib ANOVA p < 1x10-4). For RF, the relationship with SAS was best fit with a linear model (beta p < 1x10-4) while RCS model was a better fit for EAS (ANOVA p < 1x10-4). In this study, we identify relationships between genetic ancestry and phenotypes that are best fit with non-linear modeling techniques. The assumption of linearity for regression modeling is integral for proper fitting of a model and there is no knowing a priori to modeling if the relationship is truly linear.

Evidence of recent and ongoing admixture in the U.S. and influences on health and disparities.

Many researchers in genetics and social science incorporate information about race in their work. However, migrations (historical and forced) and social mobility have brought formerly separated populations of humans together, creating younger generations of individuals who have more complex and diverse ancestry and race profiles than older age groups. Here, we sought to better understand how temporal changes in genetic admixture influence levels of heterozygosity and impact health outcomes. We evaluated variation in genetic ancestry over 100 birth years in a cohort of 35,842 individuals with electronic health record (EHR) information in the Southeastern United States. Using the software STRUCTURE, we analyzed 2,678 ancestrally informative markers relative to three ancestral clusters (African, East Asian, and European) and observed rising levels of admixture for all clinically-defined race groups since 1990. Most race groups also exhibited increases in heterozygosity and long-range linkage disequilibrium over time, further supporting the finding of increasing admixture in young individuals in our cohort. These data are consistent with United States Census information from broader geographic areas and highlight the changing demography of the population. This increased diversity challenges classic approaches to studies of genotype-phenotype relationships which motivated us to explore the relationship between heterozygosity and disease diagnosis. Using a phenome-wide association study approach, we explored the relationship between admixture and disease risk and found that increased admixture resulted in protective associations with female reproductive disorders and increased risk for diseases with links to autoimmune dysfunction. These data suggest that tendencies in the United States population are increasing ancestral complexity over time. Further, these observations imply that, because both prevalence and severity of many diseases vary by race groups, complexity of ancestral origins influences health and disparities.

Whole genome sequence analysis of apparent treatment resistant hypertension status in participants from the Trans-Omics for Precision Medicine program.

Introduction: Apparent treatment-resistant hypertension (aTRH) is characterized by the use of four or more antihypertensive (AHT) classes to achieve blood pressure (BP) control. In the current study, we conducted single-variant and gene-based analyses of aTRH among individuals from 12 Trans-Omics for Precision Medicine cohorts with whole-genome sequencing data. Methods: Cases were defined as individuals treated for hypertension (HTN) taking three different AHT classes, with average systolic BP ≥ 140 or diastolic BP ≥ 90 mmHg, or four or more medications regardless of BP (n = 1,705). A normotensive control group was defined as individuals with BP < 140/90 mmHg (n = 22,079), not on AHT medication. A second control group comprised individuals who were treatment responsive on one AHT medication with BP < 140/ 90 mmHg (n = 5,424). Logistic regression with kinship adjustment using the Scalable and Accurate Implementation of Generalized mixed models (SAIGE) was performed, adjusting for age, sex, and genetic ancestry. We assessed variants using SKAT-O in rare-variant analyses. Single-variant and gene-based tests were conducted in a pooled multi-ethnicity stratum, as well as self-reported ethnic/racial strata (European and African American). Results: One variant in the known HTN locus, KCNK3, was a top finding in the multi-ethnic analysis (p = 8.23E-07) for the normotensive control group [rs12476527, odds ratio (95% confidence interval) = 0.80 (0.74-0.88)]. This variant was replicated in the Vanderbilt University Medical Center's DNA repository data. Aggregate gene-based signals included the genes AGTPBP, MYL4, PDCD4, BBS9, ERG, and IER3. Discussion: Additional work validating these loci in larger, more diverse populations, is warranted to determine whether these regions influence the pathobiology of aTRH.