Comparison of phenomic profiles in the All of Us Research Program against the US general population and the UK Biobank.

IMPORTANCE: Knowledge gained from cohort studies has dramatically advanced both public and precision health. The All of Us Research Program seeks to enroll 1 million diverse participants who share multiple sources of data, providing unique opportunities for research. It is important to understand the phenomic profiles of its participants to conduct research in this cohort. OBJECTIVES: More than 280 000 participants have shared their electronic health records (EHRs) in the All of Us Research Program. We aim to understand the phenomic profiles of this cohort through comparisons with those in the US general population and a well-established nation-wide cohort, UK Biobank, and to test whether association results of selected commonly studied diseases in the All of Us cohort were comparable to those in UK Biobank. MATERIALS AND METHODS: We included participants with EHRs in All of Us and participants with health records from UK Biobank. The estimates of prevalence of diseases in the US general population were obtained from the Global Burden of Diseases (GBD) study. We conducted phenome-wide association studies (PheWAS) of 9 commonly studied diseases in both cohorts. RESULTS: This study included 287 012 participants from the All of Us EHR cohort and 502 477 participants from the UK Biobank. A total of 314 diseases curated by the GBD were evaluated in All of Us, 80.9% (N = 254) of which were more common in All of Us than in the US general population [prevalence ratio (PR) >1.1, P < 2 × 10-5]. Among 2515 diseases and phenotypes evaluated in both All of Us and UK Biobank, 85.6% (N = 2152) were more common in All of Us (PR >1.1, P < 2 × 10-5). The Pearson correlation coefficients of effect sizes from PheWAS between All of Us and UK Biobank were 0.61, 0.50, 0.60, 0.57, 0.40, 0.53, 0.46, 0.47, and 0.24 for ischemic heart diseases, lung cancer, chronic obstructive pulmonary disease, dementia, colorectal cancer, lower back pain, multiple sclerosis, lupus, and cystic fibrosis, respectively. DISCUSSION: Despite the differences in prevalence of diseases in All of Us compared to the US general population or the UK Biobank, our study supports that All of Us can facilitate rapid investigation of a broad range of diseases. CONCLUSION: Most diseases were more common in All of Us than in the general US population or the UK Biobank. Results of disease-disease association tests from All of Us are comparable to those estimated in another well-studied national cohort.

Bronchopulmonary dysplasia is associated with polyhydramnios in a scan for novel perinatal risk factors.

BACKGROUND: The pathogenesis of bronchopulmonary dysplasia (BPD) is multifactorial, and there are limited data about prenatal exposures and risk of BPD. STUDY DESIGN: Our study performed parallel analyses using a logistic regression model in a cohort of 4527 infants with data from a curated registry and using a phenome wide association study (PheWAS) based on ICD9/10-based phecodes. We examined 20 prenatal exposures from a neonatal intensive care unit (NICU) curated registry database related to pregnancy and maternal health as well as 94 maternal diagnosis phecodes with a PheWAS analysis. RESULT: In both the curated registry and PheWAS analyses, polyhydramnios was associated with an increased risk of BPD (OR 5.70, 95% CI 2.78-11.44, p = 1.37 × 10-6). CONCLUSION: Our data suggest that polyhydramnios may be a clinical indicator of premature infants at increased risk for bronchopulmonary dysplasia. Combining curated registry data with PheWAS analysis creates a valuable tool to generate hypotheses. IMPACT: Polyhydramnios was significantly associated with bronchopulmonary dysplasia in both a curated registry and by ICD coding analysis with a phenome wide association study (PheWAS). Preterm polyhydramnios may be a clinical indicator of infants at increased risk for developing bronchopulmonary dysplasia after preterm birth. Combining curated registry with PheWAS analysis creates a valuable tool to generate hypotheses about perinatal risk factors and morbidities associated with preterm birth.

Association of Pathogenic Variants in Hereditary Cancer Genes With Multiple Diseases.

Importance: Knowledge about the spectrum of diseases associated with hereditary cancer syndromes may improve disease diagnosis and management for patients and help to identify high-risk individuals. Objective: To identify phenotypes associated with hereditary cancer genes through a phenome-wide association study. Design, Setting, and Participants: This phenome-wide association study used health data from participants in 3 cohorts. The Electronic Medical Records and Genomics Sequencing (eMERGEseq) data set recruited predominantly healthy individuals from 10 US medical centers from July 16, 2016, through February 18, 2018, with a mean follow-up through electronic health records (EHRs) of 12.7 (7.4) years. The UK Biobank (UKB) cohort recruited participants from March 15, 2006, through August 1, 2010, with a mean (SD) follow-up of 12.4 (1.0) years. The Hereditary Cancer Registry (HCR) recruited patients undergoing clinical genetic testing at Vanderbilt University Medical Center from May 1, 2012, through December 31, 2019, with a mean (SD) follow-up through EHRs of 8.8 (6.5) years. Exposures: Germline variants in 23 hereditary cancer genes. Pathogenic and likely pathogenic variants for each gene were aggregated for association analyses. Main Outcomes and Measures: Phenotypes in the eMERGEseq and HCR cohorts were derived from the linked EHRs. Phenotypes in UKB were from multiple sources of health-related data. Results: A total of 214 020 participants were identified, including 23 544 in eMERGEseq cohort (mean [SD] age, 47.8 [23.7] years; 12 611 women [53.6%]), 187 234 in the UKB cohort (mean [SD] age, 56.7 [8.1] years; 104 055 [55.6%] women), and 3242 in the HCR cohort (mean [SD] age, 52.5 [15.5] years; 2851 [87.9%] women). All 38 established gene-cancer associations were replicated, and 19 new associations were identified. These included the following 7 associations with neoplasms: CHEK2 with leukemia (odds ratio [OR], 3.81 [95% CI, 2.64-5.48]) and plasma cell neoplasms (OR, 3.12 [95% CI, 1.84-5.28]), ATM with gastric cancer (OR, 4.27 [95% CI, 2.35-7.44]) and pancreatic cancer (OR, 4.44 [95% CI, 2.66-7.40]), MUTYH (biallelic) with kidney cancer (OR, 32.28 [95% CI, 6.40-162.73]), MSH6 with bladder cancer (OR, 5.63 [95% CI, 2.75-11.49]), and APC with benign liver/intrahepatic bile duct tumors (OR, 52.01 [95% CI, 14.29-189.29]). The remaining 12 associations with nonneoplastic diseases included BRCA1/2 with ovarian cysts (OR, 3.15 [95% CI, 2.22-4.46] and 3.12 [95% CI, 2.36-4.12], respectively), MEN1 with acute pancreatitis (OR, 33.45 [95% CI, 9.25-121.02]), APC with gastritis and duodenitis (OR, 4.66 [95% CI, 2.61-8.33]), and PTEN with chronic gastritis (OR, 15.68 [95% CI, 6.01-40.92]). Conclusions and Relevance: The findings of this genetic association study analyzing the EHRs of 3 large cohorts suggest that these new phenotypes associated with hereditary cancer genes may facilitate early detection and better management of cancers. This study highlights the potential benefits of using EHR data in genomic medicine.

LabWAS: Novel findings and study design recommendations from a meta-analysis of clinical labs in two independent biobanks.

Phenotypes extracted from Electronic Health Records (EHRs) are increasingly prevalent in genetic studies. EHRs contain hundreds of distinct clinical laboratory test results, providing a trove of health data beyond diagnoses. Such lab data is complex and lacks a ubiquitous coding scheme, making it more challenging than diagnosis data. Here we describe the first large-scale cross-health system genome-wide association study (GWAS) of EHR-based quantitative laboratory-derived phenotypes. We meta-analyzed 70 lab traits matched between the BioVU cohort from the Vanderbilt University Health System and the Michigan Genomics Initiative (MGI) cohort from Michigan Medicine. We show high replication of known association for these traits, validating EHR-based measurements as high-quality phenotypes for genetic analysis. Notably, our analysis provides the first replication for 699 previous GWAS associations across 46 different traits. We discovered 31 novel associations at genome-wide significance for 22 distinct traits, including the first reported associations for two lab-based traits. We replicated 22 of these novel associations in an independent tranche of BioVU samples. The summary statistics for all association tests are freely available to benefit other researchers. Finally, we performed mirrored analyses in BioVU and MGI to assess competing analytic practices for EHR lab traits. We find that using the mean of all available lab measurements provides a robust summary value, but alternate summarizations can improve power in certain circumstances. This study provides a proof-of-principle for cross health system GWAS and is a framework for future studies of quantitative EHR lab traits.

Efficiently controlling for case-control imbalance and sample relatedness in large-scale genetic association studies.

In genome-wide association studies (GWAS) for thousands of phenotypes in large biobanks, most binary traits have substantially fewer cases than controls. Both of the widely used approaches, the linear mixed model and the recently proposed logistic mixed model, perform poorly; they produce large type I error rates when used to analyze unbalanced case-control phenotypes. Here we propose a scalable and accurate generalized mixed model association test that uses the saddlepoint approximation to calibrate the distribution of score test statistics. This method, SAIGE (Scalable and Accurate Implementation of GEneralized mixed model), provides accurate P values even when case-control ratios are extremely unbalanced. SAIGE uses state-of-art optimization strategies to reduce computational costs; hence, it is applicable to GWAS for thousands of phenotypes by large biobanks. Through the analysis of UK Biobank data of 408,961 samples from white British participants with European ancestry for > 1,400 binary phenotypes, we show that SAIGE can efficiently analyze large sample data, controlling for unbalanced case-control ratios and sample relatedness.

PregOMICS-Leveraging systems biology and bioinformatics for drug repurposing in maternal-child health.

Obstetric diseases remain underserved and understudied. Drug repurposing-utilization of a drug whose use is accepted in one condition for a different condition-could represent a rapid and low-cost way to identify new therapies that are known to be safe. In diseases of pregnancy, the known safety profile is a strong additional incentive. We describe the techniques and steps used in the use of 'omics data for drug repurposing. We illustrate these techniques using case studies of published drug repurposing projects. We provide a set of available databases with low barriers to entry which investigators can use to perform their own projects. The promise of 'omics techniques is unbiased screening, either of all drug targets or of all patients using particular drugs to find which are likely to alter disease risk or progression. However, we caution that reproducibility across the underlying studies, and thus the drugs suggested for repurposing, can be poor. We suggest that improved nosology, for example correlating patient clinical conditions with placental pathology, could yield more robust associations. We conclude that 'omics-driven drug repurposing represents a potential fruitful path to discover new, safe treatments of obstetric diseases.

Calcium channel blockers as drug repurposing candidates for gestational diabetes: Mining large scale genomic and electronic health records data to repurpose medications.

New therapeutic approaches are needed for gestational diabetes mellitus (GDM), but must show safety and efficacy in a historically understudied population. We studied associations between electronic medical record (EMR) phenotypes and genetic variants to uncover drugs currently considered safe in pregnancy that could treat or prevent GDM. We identified 129 systemically active drugs considered safe in pregnancy targeting the proteins produced from 196 genes. We tested for associations between GDM and/or type 2 diabetes (DM2) and 306 SNPs in 130 genes represented on the Illumina Infinium Human Exome Bead Chip (DM2 was included due to shared pathophysiological features with GDM). In parallel, we tested the association between drugs and glucose tolerance during pregnancy as measured by the glucose recorded during a routine 50-g glucose tolerance test (GTT). We found an association between GDM/DM2 and the genes targeted by 11 drug classes. In the EMR analysis, 6 drug classes were associated with changes in GTT. Two classes were identified in both analyses. L-type calcium channel blocking antihypertensives (CCBs), were associated with a 3.18 mg/dL (95% CI -6.18 to -0.18) decrease in glucose during GTT, and serotonin receptor type 3 (5HT-3) antagonist antinausea medications were associated with a 3.54 mg/dL (95% CI 1.86-5.23) increase in glucose during GTT. CCBs were identified as a class of drugs considered safe in pregnancy could have efficacy in treating or preventing GDM. 5HT-3 antagonists may be associated with worse glucose tolerance.

Using Human 'Experiments of Nature' to Predict Drug Safety Issues: An Example with PCSK9 Inhibitors.

INTRODUCTION: When a new drug enters the market, its full array of side effects remains to be defined. Current surveillance approaches targeting these effects remain largely reactive. There is a need for development of methods to predict specific safety events that should be sought for a given new drug during development and postmarketing activities. OBJECTIVE: We present here a safety signal identification approach applied to a new set of drug entities, inhibitors of the serine protease proprotein convertase subtilisin/kexin type 9 (PCSK9). METHODS: Using phenome-wide association study (PheWAS) methods, we analyzed available genotype and clinical data from 29,722 patients, leveraging the known effects of changes in PCSK9 to identify novel phenotypes in which this protein and its inhibitors may have impact. RESULTS: PheWAS revealed a significantly reduced risk of hypercholesterolemia (odds ratio [OR] 0.68, p = 7.6 × 10-4) in association with a known loss-of-function variant in PCSK9, R46L. Similarly, laboratory data indicated significantly reduced beta mean low-density lipoprotein cholesterol (- 14.47 mg/dL, p = 2.58 × 10-23) in individuals carrying the R46L variant. The R46L variant was also associated with an increased risk of spina bifida (OR 5.90, p = 2.7 × 10-4), suggesting that further investigation of potential connections between inhibition of PCSK9 and neural tube defects may be warranted. CONCLUSION: This novel methodology provides an opportunity to put in place new mechanisms to assess the safety and long-term tolerability of PCSK9 inhibitors specifically, and other new agents in general, as they move into human testing and expanded clinical use.

Evaluating phecodes, clinical classification software, and ICD-9-CM codes for phenome-wide association studies in the electronic health record.

OBJECTIVE: To compare three groupings of Electronic Health Record (EHR) billing codes for their ability to represent clinically meaningful phenotypes and to replicate known genetic associations. The three tested coding systems were the International Classification of Diseases, Ninth Revision, Clinical Modification (ICD-9-CM) codes, the Agency for Healthcare Research and Quality Clinical Classification Software for ICD-9-CM (CCS), and manually curated "phecodes" designed to facilitate phenome-wide association studies (PheWAS) in EHRs. METHODS AND MATERIALS: We selected 100 disease phenotypes and compared the ability of each coding system to accurately represent them without performing additional groupings. The 100 phenotypes included 25 randomly-chosen clinical phenotypes pursued in prior genome-wide association studies (GWAS) and another 75 common disease phenotypes mentioned across free-text problem lists from 189,289 individuals. We then evaluated the performance of each coding system to replicate known associations for 440 SNP-phenotype pairs. RESULTS: Out of the 100 tested clinical phenotypes, phecodes exactly matched 83, compared to 53 for ICD-9-CM and 32 for CCS. ICD-9-CM codes were typically too detailed (requiring custom groupings) while CCS codes were often not granular enough. Among 440 tested known SNP-phenotype associations, use of phecodes replicated 153 SNP-phenotype pairs compared to 143 for ICD-9-CM and 139 for CCS. Phecodes also generally produced stronger odds ratios and lower p-values for known associations than ICD-9-CM and CCS. Finally, evaluation of several SNPs via PheWAS identified novel potential signals, some seen in only using the phecode approach. Among them, rs7318369 in PEPD was associated with gastrointestinal hemorrhage. CONCLUSION: Our results suggest that the phecode groupings better align with clinical diseases mentioned in clinical practice or for genomic studies. ICD-9-CM, CCS, and phecode groupings all worked for PheWAS-type studies, though the phecode groupings produced superior results.

A common deletion in the haptoglobin gene associated with blood cholesterol levels among Chinese women.

Haptoglobin (HP) protein plays a critical role in binding and removing free hemoglobin from blood. A deletion in the HP gene affects the protein structure and function. A recent study developed a novel method to impute this variant and discovered significant association of this variant with low-density lipoprotein (LDL) and total cholesterol levels among European descendants. In the present study, we investigated this variant among 3608 Chinese women. Consistent with findings from Europeans, we found significant associations between the deletion with lower cholesterol levels; women homozygous for the deletion allele (HP1-HP1), had a lower level of total cholesterol (-4.24 mg dl-1, P=0.02) and LDL cholesterol (-3.43 mg dl-1, P=0.03) than those not carrying the deletion allele (HP2-HP2). Especially, women carrying the HP1S-HP1S, had an even lower level of total cholesterol (-5.59 mg dl-1, P=7.0 × 10-3) and LDL cholesterol (-4.68 mg dl-1, P=8.0 × 10-3) compared to those carrying HP2-HP2. These associations remained significant after an adjustment for an established cholesterol level-related variant, rs2000999. Our study extends the previous findings regarding the association of HP structure variant with blood cholesterol levels to East Asians and affirms the validity of the new methodology for assessing HP structure variation.

Evaluating electronic health record data sources and algorithmic approaches to identify hypertensive individuals.

OBJECTIVE: Phenotyping algorithms applied to electronic health record (EHR) data enable investigators to identify large cohorts for clinical and genomic research. Algorithm development is often iterative, depends on fallible investigator intuition, and is time- and labor-intensive. We developed and evaluated 4 types of phenotyping algorithms and categories of EHR information to identify hypertensive individuals and controls and provide a portable module for implementation at other sites. MATERIALS AND METHODS: We reviewed the EHRs of 631 individuals followed at Vanderbilt for hypertension status. We developed features and phenotyping algorithms of increasing complexity. Input categories included International Classification of Diseases, Ninth Revision (ICD9) codes, medications, vital signs, narrative-text search results, and Unified Medical Language System (UMLS) concepts extracted using natural language processing (NLP). We developed a module and tested portability by replicating 10 of the best-performing algorithms at the Marshfield Clinic. RESULTS: Random forests using billing codes, medications, vitals, and concepts had the best performance with a median area under the receiver operator characteristic curve (AUC) of 0.976. Normalized sums of all 4 categories also performed well (0.959 AUC). The best non-NLP algorithm combined normalized ICD9 codes, medications, and blood pressure readings with a median AUC of 0.948. Blood pressure cutoffs or ICD9 code counts alone had AUCs of 0.854 and 0.908, respectively. Marshfield Clinic results were similar. CONCLUSION: This work shows that billing codes or blood pressure readings alone yield good hypertension classification performance. However, even simple combinations of input categories improve performance. The most complex algorithms classified hypertension with excellent recall and precision.