ABSTRACT
Variability in drug effectiveness and provider prescribing patterns have been reported in different racial and ethnic populations. We sought to evaluate antihypertensive drug effectiveness and prescribing patterns among self-identified Hispanic/Latino (Hispanic), Non-Hispanic Black (Black), and Non-Hispanic White (White) populations that enrolled in the NIH All of Us Research Program, a US longitudinal cohort. We employed a self-controlled case study method using electronic health record and survey data from 17,718 White, Hispanic, and Black participants who were diagnosed with essential hypertension and prescribed at least one of 19 commonly used antihypertensive medications. Effectiveness was determined by calculating the reduction in systolic blood pressure measurements after 28 or more days of drug exposure. Starting systolic blood pressure and effectiveness for each medication were compared for self-reported Black, Hispanic, and White participants using adjusted linear regressions. Black and Hispanic participants were started on antihypertensive medications at significantly higher SBP than White participants in 13 and 7 out of 19 medications, respectively. More Black participants were prescribed multiple antihypertensive medications (58.46%) than White (52.35%) or Hispanic (49.9%) participants. First-line HTN medications differed by race and ethnicity. Following the 2017 American College of Cardiology and the American Heart Association High Blood Pressure Guideline release, around 64% of Black participants were prescribed a recommended first-line antihypertensive drug compared with 76% of White and 82% of Hispanic participants. Effect sizes suggested that most antihypertensive drugs were less effective in Hispanic and Black, compared with White, participants, and statistical significance was reached in 6 out of 19 drugs. These results indicate that Black and Hispanic populations may benefit from earlier intervention and screening and highlight the potential benefits of personalizing first-line medications.
ABSTRACT
Pharmacogenomics promises improved outcomes through individualized prescribing. However, the lack of diversity in studies impedes clinical translation and equitable application of precision medicine. We evaluated the frequencies of PGx variants, predicted phenotypes, and medication exposures using whole genome sequencing and EHR data from nearly 100k diverse All of Us Research Program participants. We report 100% of participants carried at least one pharmacogenomics variant and nearly all (99.13%) had a predicted phenotype with prescribing recommendations. Clinical impact was high with over 20% having both an actionable phenotype and a prior exposure to an impacted medication with pharmacogenomic prescribing guidance. Importantly, we also report hundreds of alleles and predicted phenotypes that deviate from known frequencies and/or were previously unreported, including within admixed American and African ancestry groups.
ABSTRACT
Clostridioides difficile infection causes pathology that ranges in severity from diarrhea to pseudomembranous colitis. Toxin A and Toxin B are the two primary virulence factors secreted by C. difficile that drive disease severity. The toxins damage intestinal epithelial cells leading to a loss of barrier integrity and induction of a proinflammatory host response. Monoclonal antibodies (mAbs) that neutralize Toxin A and Toxin B, actoxumab and bezlotoxumab, respectively, significantly reduce disease severity in a murine model of C. difficile infection. However, the impact of toxin neutralization on the induction and quality of the innate immune response following infection is unknown. The goal of this study was to define the quality of the host innate immune response in the context of anti-toxin mAbs therapy. At day 2 post-infection, C. difficile-infected, mAbs-treated mice had significantly less disease compared to isotype-treated mice despite remaining colonized with C. difficile. C. difficile-infected mAbs-treated mice still exhibited marked neutrophil infiltration and induction of a subset of proinflammatory cytokines within the intestinal lamina propria following infection that is comparable to isotype-treated mice. Furthermore, both mAbs and isotype-treated mice had an increase in IL-22-producing ILCs in the intestine following infection. MAbs-treated mice exhibited increased infiltration of eosinophils in the intestinal lamina propria, which has been previously reported to promote a protective host response following C. difficile infection. These findings show that activation of host protective mechanisms remain intact in the context of monoclonal antibody-mediated toxin neutralization.
Subject(s)
Antibodies, Monoclonal , Bacterial Toxins , Clostridioides difficile , Clostridium Infections , Immunity, Innate , Animals , Bacterial Toxins/immunology , Antibodies, Monoclonal/immunology , Antibodies, Monoclonal/pharmacology , Clostridioides difficile/immunology , Clostridium Infections/immunology , Clostridium Infections/prevention & control , Clostridium Infections/microbiology , Mice , Enterotoxins/immunology , Disease Models, Animal , Antibodies, Neutralizing/immunology , Bacterial Proteins/immunology , Female , Cytokines/metabolism , Broadly Neutralizing Antibodies/immunology , Mice, Inbred C57BLABSTRACT
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.
ABSTRACT
Summary: With the rapid growth of genetic data linked to electronic health record data in huge cohorts, large-scale phenome-wide association study (PheWAS), have become powerful discovery tools in biomedical research. PheWAS is an analysis method to study phenotype associations utilizing longitudinal electronic health record (EHR) data. Previous PheWAS packages were developed mostly in the days of smaller biobanks and with earlier PheWAS approaches. PheTK was designed to simplify analysis and efficiently handle biobank-scale data. PheTK uses multithreading and supports a full PheWAS workflow including extraction of data from OMOP databases and Hail matrix tables as well as PheWAS analysis for both phecode version 1.2 and phecodeX. Benchmarking results showed PheTK took 64% less time than the R PheWAS package to complete the same workflow. PheTK can be run locally or on cloud platforms such as the All of Us Researcher Workbench ( All of Us ) or the UK Biobank (UKB) Research Analysis Platform (RAP). Availability and implementation: The PheTK package is freely available on the Python Package Index (PyPi) and on GitHub under GNU Public License (GPL-3) at https://github.com/nhgritctran/PheTK . It is implemented in Python and platform independent. The demonstration workspace for All of Us will be made available in the future as a featured workspace. Contact: PheTK@mail.nih.gov.
ABSTRACT
Genome-wide association studies (GWAS) have been instrumental in identifying genetic associations for various diseases and traits. However, uncovering genetic underpinnings among traits beyond univariate phenotype associations remains a challenge. Multi-phenotype associations (MPA), or genetic pleiotropy, offer important insights into shared genes and pathways among traits, enhancing our understanding of genetic architectures of complex diseases. GWAS of biobank-linked electronic health record (EHR) data are increasingly being utilized to identify MPA among various traits and diseases. However, methodologies that can efficiently take advantage of distributed EHR to detect MPA are still lacking. Here, we introduce mixWAS, a novel algorithm that efficiently and losslessly integrates multiple EHRs via summary statistics, allowing the detection of MPA among mixed phenotypes while accounting for heterogeneities across EHRs. Simulations demonstrate that mixWAS outperforms the widely used MPA detection method, Phenome-wide association study (PheWAS), across diverse scenarios. Applying mixWAS to data from seven EHRs in the US, we identified 4,534 MPA among blood lipids, BMI, and circulatory diseases. Validation in an independent EHR data from UK confirmed 97.7% of the associations. mixWAS fundamentally improves the detection of MPA and is available as a free, open-source software.
ABSTRACT
The National Institutes of Health's All of Us Research Program is an accessible platform that hosts genomic and phenotypic data to be collected from 1 million participants in the United States. Its mission is to accelerate medical research and clinical breakthroughs with a special emphasis on diversity.
Subject(s)
Biomedical Research , Population Health , Humans , United States , Data Science , National Institutes of Health (U.S.)ABSTRACT
Clostridioides difficile (C. diff.) infection (CDI) is a leading cause of hospital acquired diarrhea in North America and Europe and a major cause of morbidity and mortality. Known risk factors do not fully explain CDI susceptibility, and genetic susceptibility is suggested by the fact that some patients with colons that are colonized with C. diff. do not develop any infection while others develop severe or recurrent infections. To identify common genetic variants associated with CDI, we performed a genome-wide association analysis in 19,861 participants (1349 cases; 18,512 controls) from the Electronic Medical Records and Genomics (eMERGE) Network. Using logistic regression, we found strong evidence for genetic variation in the DRB locus of the MHC (HLA) II region that predisposes individuals to CDI (P > 1.0 × 10-14; OR 1.56). Altered transcriptional regulation in the HLA region may play a role in conferring susceptibility to this opportunistic enteric pathogen.
Subject(s)
Clostridium Infections , Genome-Wide Association Study , Humans , Clostridium Infections/genetics , Diarrhea , Histocompatibility Antigens , HLA Antigens/genetics , Histocompatibility Antigens Class II , Genetic VariationABSTRACT
Malaria is caused by Plasmodium species and remains a significant cause of morbidity and mortality globally. Gut bacteria can influence the severity of malaria, but the contribution of specific bacteria to the risk of severe malaria is unknown. Here, multiomics approaches demonstrate that specific species of Bacteroides are causally linked to the risk of severe malaria. Plasmodium yoelii hyperparasitemia-resistant mice gavaged with murine-isolated Bacteroides fragilis develop P. yoelii hyperparasitemia. Moreover, Bacteroides are significantly more abundant in Ugandan children with severe malarial anemia than with asymptomatic P. falciparum infection. Human isolates of Bacteroides caccae, Bacteroides uniformis, and Bacteroides ovatus were able to cause susceptibility to severe malaria in mice. While monocolonization of germ-free mice with Bacteroides alone is insufficient to cause susceptibility to hyperparasitemia, meta-analysis across multiple studies support a main role for Bacteroides in susceptibility to severe malaria. Approaches that target gut Bacteroides present an opportunity to prevent severe malaria and associated deaths.
Subject(s)
Anemia , Malaria , Plasmodium yoelii , Child , Humans , Animals , Mice , Microbial Consortia , Bacteroides/genetics , Bacteroides fragilis , Anemia/etiologyABSTRACT
OBJECTIVE: The All of Us Research Program (All of Us) aims to recruit over a million participants to further precision medicine. Essential to the verification of biobanks is a replication of known associations to establish validity. Here, we evaluated how well All of Us data replicated known cigarette smoking associations. MATERIALS AND METHODS: We defined smoking exposure as follows: (1) an EHR Smoking exposure that used International Classification of Disease codes; (2) participant provided information (PPI) Ever Smoking; and, (3) PPI Current Smoking, both from the lifestyle survey. We performed a phenome-wide association study (PheWAS) for each smoking exposure measurement type. For each, we compared the effect sizes derived from the PheWAS to published meta-analyses that studied cigarette smoking from PubMed. We defined two levels of replication of meta-analyses: (1) nominally replicated: which required agreement of direction of effect size, and (2) fully replicated: which required overlap of confidence intervals. RESULTS: PheWASes with EHR Smoking, PPI Ever Smoking, and PPI Current Smoking revealed 736, 492, and 639 phenome-wide significant associations, respectively. We identified 165 meta-analyses representing 99 distinct phenotypes that could be matched to EHR phenotypes. At P < .05, 74 were nominally replicated and 55 were fully replicated. At P < 2.68 × 10-5 (Bonferroni threshold), 58 were nominally replicated and 40 were fully replicated. DISCUSSION: Most phenotypes found in published meta-analyses associated with smoking were nominally replicated in All of Us. Both survey and EHR definitions for smoking produced similar results. CONCLUSION: This study demonstrated the feasibility of studying common exposures using All of Us data.
Subject(s)
Genome-Wide Association Study , Population Health , Humans , Genome-Wide Association Study/methods , Phenotype , Polymorphism, Single Nucleotide , SmokingABSTRACT
Recently, large scale genomic projects such as All of Us and the UK Biobank have introduced a new research paradigm where data are stored centrally in cloud-based Trusted Research Environments (TREs). To characterize the advantages and drawbacks of different TRE attributes in facilitating cross-cohort analysis, we conduct a Genome-Wide Association Study of standard lipid measures using two approaches: meta-analysis and pooled analysis. Comparison of full summary data from both approaches with an external study shows strong correlation of known loci with lipid levels (R2 ~ 83-97%). Importantly, 90 variants meet the significance threshold only in the meta-analysis and 64 variants are significant only in pooled analysis, with approximately 20% of variants in each of those groups being most prevalent in non-European, non-Asian ancestry individuals. These findings have important implications, as technical and policy choices lead to cross-cohort analyses generating similar, but not identical results, particularly for non-European ancestral populations.
Subject(s)
Genome-Wide Association Study , Population Health , Humans , Genomics , Policy , LipidsABSTRACT
The All of Us Research Program's Data and Research Center (DRC) was established to help acquire, curate, and provide access to one of the world's largest and most diverse datasets for precision medicine research. Already, over 500,000 participants are enrolled in All of Us, 80% of whom are underrepresented in biomedical research, and data are being analyzed by a community of over 2,300 researchers. The DRC created this thriving data ecosystem by collaborating with engaged participants, innovative program partners, and empowered researchers. In this review, we first describe how the DRC is organized to meet the needs of this broad group of stakeholders. We then outline guiding principles, common challenges, and innovative approaches used to build the All of Us data ecosystem. Finally, we share lessons learned to help others navigate important decisions and trade-offs in building a modern biomedical data platform.
Subject(s)
Biomedical Research , Population Health , Humans , Ecosystem , Precision MedicineABSTRACT
Population-scale biobanks linked to electronic health record data provide vast opportunities to extend our knowledge of human genetics and discover new phenotype-genotype associations. Given their dense phenotype data, biobanks can also facilitate replication studies on a phenome-wide scale. Here, we introduce the phenotype-genotype reference map (PGRM), a set of 5,879 genetic associations from 523 GWAS publications that can be used for high-throughput replication experiments. PGRM phenotypes are standardized as phecodes, ensuring interoperability between biobanks. We applied the PGRM to five ancestry-specific cohorts from four independent biobanks and found evidence of robust replications across a wide array of phenotypes. We show how the PGRM can be used to detect data corruption and to empirically assess parameters for phenome-wide studies. Finally, we use the PGRM to explore factors associated with replicability of GWAS results.
Subject(s)
Biological Specimen Banks , Data Science , Humans , Phenomics , Phenotype , GenotypeABSTRACT
PURPOSE: Automated use of electronic health records may aid in decreasing the diagnostic delay for rare diseases. The phenotype risk score (PheRS) is a weighted aggregate of syndromically related phenotypes that measures the similarity between an individual's conditions and features of a disease. For some diseases, there are individuals without a diagnosis of that disease who have scores similar to diagnosed patients. These individuals may have that disease but not yet be diagnosed. METHODS: We calculated the PheRS for cystic fibrosis (CF) for 965,626 subjects in the Vanderbilt University Medical Center electronic health record. RESULTS: Of the 400 subjects with the highest PheRS for CF, 248 (62%) had been diagnosed with CF. Twenty-six of the remaining participants, those who were alive and had DNA available in the linked DNA biobank, underwent clinical review and sequencing analysis of CFTR and SERPINA1. This uncovered a potential diagnosis for 2 subjects, 1 with CF and 1 with alpha-1-antitrypsin deficiency. An additional 7 subjects had pathogenic or likely pathogenic variants, 2 in CFTR and 5 in SERPINA1. CONCLUSION: These findings may be clinically actionable for the providers caring for these patients. Importantly, this study highlights feasibility and challenges for future implications of this approach.
Subject(s)
Cystic Fibrosis Transmembrane Conductance Regulator , Cystic Fibrosis , Humans , Cystic Fibrosis Transmembrane Conductance Regulator/genetics , Electronic Health Records , Delayed Diagnosis , Cystic Fibrosis/diagnosis , Cystic Fibrosis/genetics , Cystic Fibrosis/pathology , DNA , MutationABSTRACT
Alpha-1 antitrypsin deficiency (AATD), a relatively common autosomal recessive genetic disorder, is underdiagnosed in symptomatic individuals. We sought to compare the risk of liver transplantation associated with hepatitis C infection with AATD heterozygotes and homozygotes and determine if SERPINA1 sequencing would identify undiagnosed AATD. We performed a retrospective cohort study in a deidentified Electronic Health Record (EHR)-linked DNA biobank with 72,027 individuals genotyped for the M, Z, and S alleles in SERPINA1. We investigated liver transplantation frequency by genotype group and compared with hepatitis C infection. We performed SERPINA1 sequencing in carriers of pathogenic AATD alleles who underwent liver transplantation. Liver transplantation was associated with the Z allele (ZZ: odds ratio [OR] = 1.31, p<2e-16; MZ: OR = 1.02, p = 1.2e-13) and with hepatitis C (OR = 1.20, p<2e-16). For liver transplantation, there was a significant interaction between genotype and hepatitis C (ZZ: interaction OR = 1.23, p = 4.7e-4; MZ: interaction OR = 1.11, p = 6.9e-13). Sequencing uncovered a second, rare, pathogenic SERPINA1 variant in six of 133 individuals with liver transplants and without hepatitis C. Liver transplantation was more common in individuals with AATD risk alleles (including heterozygotes), and AATD and hepatitis C demonstrated evidence of a gene-environment interaction in relation to liver transplantation. The current AATD screening strategy may miss diagnoses whereas SERPINA1 sequencing may increase diagnostic yield for AATD, stratify risk for liver disease, and inform clinical management for individuals with AATD risk alleles and liver disease risk factors.
Subject(s)
Hepatitis C , alpha 1-Antitrypsin Deficiency , Humans , Alleles , Gene-Environment Interaction , Retrospective Studies , alpha 1-Antitrypsin Deficiency/diagnosis , alpha 1-Antitrypsin Deficiency/genetics , Hepatitis C/genetics , Hepacivirus/genetics , Genetics, Population , alpha 1-Antitrypsin/geneticsABSTRACT
Common data models solve many challenges of standardizing electronic health record (EHR) data but are unable to semantically integrate all of the resources needed for deep phenotyping. Open Biological and Biomedical Ontology (OBO) Foundry ontologies provide computable representations of biological knowledge and enable the integration of heterogeneous data. However, mapping EHR data to OBO ontologies requires significant manual curation and domain expertise. We introduce OMOP2OBO, an algorithm for mapping Observational Medical Outcomes Partnership (OMOP) vocabularies to OBO ontologies. Using OMOP2OBO, we produced mappings for 92,367 conditions, 8611 drug ingredients, and 10,673 measurement results, which covered 68-99% of concepts used in clinical practice when examined across 24 hospitals. When used to phenotype rare disease patients, the mappings helped systematically identify undiagnosed patients who might benefit from genetic testing. By aligning OMOP vocabularies to OBO ontologies our algorithm presents new opportunities to advance EHR-based deep phenotyping.
ABSTRACT
Many epidemiologic studies have identified important relationships between leukocyte telomere length (LTL) with genetics and health. Most of these studies have been significantly limited in scope by focusing predominantly on individual diseases or restricted to GWAS analysis. Using two large patient populations derived from Vanderbilt University and Marshfield Clinic biobanks linked to genomic and phenomic data from medical records, we investigated the inter-relationship between LTL, genomics, and human health. Our GWAS confirmed 11 genetic loci previously associated with LTL and two novel loci in SCNN1D and PITPNM1. PheWAS of LTL identified 67 distinct clinical phenotypes associated with both short and long LTL. We demonstrated that several diseases associated with LTL were related to one another but were largely independent from LTL genetics. Age of death was correlated with LTL independent of age. Those with very short LTL (<-1.5 standard deviation [SD]) died 10.4 years (p < 0.0001) younger than those with average LTL (±0.5 SD; mean age of death = 74.2 years). Likewise, those with very long LTL (>1.5 SD) died 1.9 years (p = 0.0175) younger than those with average LTL. This is consistent with the PheWAS results showing diseases associating with both short and long LTL. Finally, we estimated that the genome (12.8%) and age (8.5%) explain the largest proportion of LTL variance, whereas the phenome (1.5%) and sex (0.9%) explained a smaller fraction. In total, 23.7% of LTL variance was explained. These observations provide the rationale for expanded research to understand the multifaceted correlations between TL biology and human health over time, leading to effective LTL usage in medical applications.
Subject(s)
Leukocytes , Telomere , Humans , Aged , Telomere/genetics , Calcium-Binding Proteins/genetics , Eye Proteins/genetics , Membrane Proteins/geneticsABSTRACT
OBJECTIVE: Diverticular disease (DD) is one of the most prevalent conditions encountered by gastroenterologists, affecting ~50% of Americans before the age of 60. Our aim was to identify genetic risk variants and clinical phenotypes associated with DD, leveraging multiple electronic health record (EHR) data sources of 91,166 multi-ancestry participants with a Natural Language Processing (NLP) technique. MATERIALS AND METHODS: We developed a NLP-enriched phenotyping algorithm that incorporated colonoscopy or abdominal imaging reports to identify patients with diverticulosis and diverticulitis from multicenter EHRs. We performed genome-wide association studies (GWAS) of DD in European, African and multi-ancestry participants, followed by phenome-wide association studies (PheWAS) of the risk variants to identify their potential comorbid/pleiotropic effects in clinical phenotypes. RESULTS: Our developed algorithm showed a significant improvement in patient classification performance for DD analysis (algorithm PPVs ≥ 0.94), with up to a 3.5 fold increase in terms of the number of identified patients than the traditional method. Ancestry-stratified analyses of diverticulosis and diverticulitis of the identified subjects replicated the well-established associations between ARHGAP15 loci with DD, showing overall intensified GWAS signals in diverticulitis patients compared to diverticulosis patients. Our PheWAS analyses identified significant associations between the DD GWAS variants and circulatory system, genitourinary, and neoplastic EHR phenotypes. DISCUSSION: As the first multi-ancestry GWAS-PheWAS study, we showcased that heterogenous EHR data can be mapped through an integrative analytical pipeline and reveal significant genotype-phenotype associations with clinical interpretation. CONCLUSION: A systematic framework to process unstructured EHR data with NLP could advance a deep and scalable phenotyping for better patient identification and facilitate etiological investigation of a disease with multilayered data.