Phenome-Wide Association Studies as a Tool to Advance Precision Medicine

Clinical Phenotypes Associated With Cerebral Small Vessel Disease: A Study of 45,013 UK Biobank Participants

BACKGROUND AND OBJECTIVES

Cerebral small vessel disease (cSVD) is the most common pathology underlying vascular cognitive impairment. Although other clinical features of cSVD are increasingly recognized, it is likely that certain symptoms are being overlooked. A comprehensive description of cSVD associations with clinical phenotypes at scale is lacking. The objective of this study was to conduct a large-scale, hypothesis-free study of associations between cSVD and clinical phenotypes in UK Biobank (UKB).

METHODS

We included participants from the UKB imaging study who had available information on total volume of white matter hyperintensities (WMHs), the most common cSVD neuroimaging feature. We included various UKB variables describing clinical phenotypes, defined as observable signs and symptoms of individuals with concurrent neuroimaging evidence of cSVD. We conducted a phenome scan using the open-source PHESANT software package. Total volume of WMHs was introduced as the independent variable and clinical phenotypes as the dependent variables in the regression model. The association of each phenotype with total volume of WMHs was tested using one of several regression analyses (all age at recruitment and sex-adjusted). All associations were corrected for multiple comparisons using the false discovery rate (FDR) correction method.

RESULTS

We included 45,013 participants in the analysis (mean age = 54.97 years, SD = 7.55). We confirm previously reported associations with depression (odds ratio [OR] = 1.07 [95% CI 1.05-1.10]), apathy (OR = 1.11 [95% CI 1.08-1.14]), falls (OR = 1.11 [95% CI 1.09-1.13]), respiratory problems (OR = 1.14 [95% CI 1.04-1.25]), and sleep disturbance (OR = 1.07 [95% CI 1.04-1.09], all FDR-adjusted p < 0.001). We further identified associations with all-cause dental issues (OR = 0.94 [95% CI 0.96-0.92]), hearing problems (OR = 1.06 [95% CI 1.03-1.08]), and eye problems (OR = 0.93 [95% CI 0.91-0.95], all FDR-adjusted p < 0.001).

DISCUSSION

Our findings suggest that presence of cSVD associates with concurrent clinical phenotypes across several body systems. We have corroborated established associations of cSVD and present novel ones. While our results do not provide causality or direction of association because of the cross-sectional nature of our study, they support the need for a more holistic view of cSVD in research, practice, and policy.

Clinical associations with treatment resistance in depression: An electronic health record study

Treatment resistance is common in major depressive disorder (MDD), yet clinical risk factors are not well understood. Using a discovery-replication design, we conducted phenome-wide association studies (PheWASs) of MDD treatment resistance in two electronic health record (EHR)-linked biobanks. The PheWAS included participants with an MDD diagnosis in the EHR and at least one antidepressant (AD) prescription. Participant lifetime diagnoses were mapped to phecodes. PheWASs were conducted for three treatment resistance outcomes based on AD prescription data: number of unique ADs prescribed, ≥1 and ≥2 CE switches. Of the 180 phecodes significantly associated with these outcomes in the discovery cohort (n = 12,558), 71 replicated (n = 8,206). In addition to identifying known clinical factors for treatment resistance in MDD, the total unique AD prescriptions was associated with additional clinical variables including irritable bowel syndrome, gastroesophageal reflux disease, symptomatic menopause, and spondylosis. We calculated polygenic risk of specific-associated conditions and tested their association with AD outcomes revealing that genetic risk for many of these conditions is also associated with the total unique AD prescriptions. The number of unique ADs prescribed, which is easily assessed in EHRs, provides a more nuanced measure of treatment resistance, and may facilitate future research and clinical application in this area.

Genetic association studies using disease liabilities from deep neural networks

The case-control study is a widely used method for investigating the genetic underpinnings of binary traits. However, long-term, prospective cohort studies often grapple with absent or evolving health-related outcomes. Here, we propose two methods, liability and meta , for conducting genome-wide association study (GWAS) that leverage disease liabilities calculated from deep patient phenotyping. Analyzing 38 common traits in ∼300,000 UK Biobank participants, we identified an increased number of loci compared to the conventional case-control approach, with high replication rates in larger external GWAS. Further analyses confirmed the disease-specificity of the genetic architecture with the meta method demonstrating higher robustness when phenotypes were imputed with low accuracy. Additionally, polygenic risk scores based on disease liabilities more effectively predicted newly diagnosed cases in the 2022 dataset, which were controls in the earlier 2019 dataset. Our findings demonstrate that integrating high-dimensional phenotypic data into deep neural networks enhances genetic association studies while capturing disease-relevant genetic architecture.

Identifying potential causal effects of Parkinson's disease: A polygenic risk score-based phenome-wide association and mendelian randomization study in UK Biobank

There is considerable uncertainty regarding the associations between various risk factors and Parkinson's Disease (PD). This study systematically screened and validated a wide range of potential PD risk factors from 502,364 participants in the UK Biobank. Baseline data for 1851 factors across 11 categories were analyzed through a phenome-wide association study (PheWAS). Polygenic risk scores (PRS) for PD were used to diagnose Parkinson's Disease and identify factors associated with PD diagnosis through PheWAS. Two-sample Mendelian randomization (MR) analysis was employed to assess causal relationships. PheWAS results revealed 267 risk factors significantly associated with PD-PRS among the 1851 factors, and of these, 27 factors showed causal evidence from MR analysis. Compelling evidence suggests that fluid intelligence score, age at first sexual intercourse, cereal intake, dried fruit intake, and average total household income before tax have emerged as newly identified risk factors for PD. Conversely, maternal smoking around birth, playing computer games, salt added to food, and time spent watching television have been identified as novel protective factors against PD. The integration of phenotypic and genomic data may help to identify risk factors and prevention targets for PD.

Clinical and genetic contributions to medical comorbidity in bipolar disorder: a study using electronic health records-linked biobank data

Bipolar disorder is a chronic and complex polygenic disease with high rates of comorbidity. However, the independent contribution of either diagnosis or genetic risk of bipolar disorder to the medical comorbidity profile of individuals with the disease remains unresolved. Here, we conducted a multi-step phenome-wide association study (PheWAS) of bipolar disorder using phenomes derived from the electronic health records of participants enrolled in the Mayo Clinic Biobank and the Mayo Clinic Bipolar Disorder Biobank. First, we explored the conditions associated with a diagnosis of bipolar disorder by conducting a phenotype-based PheWAS followed by LASSO-penalized regression to account for correlations within the phenome. Then, we explored the conditions associated with bipolar disorder polygenic risk score (BD-PRS) using a PRS-based PheWAS with a sequential exclusion approach to account for the possibility that diagnosis, instead of genetic risk, may drive such associations. 53,386 participants (58.7% women) with a mean age at analysis of 67.8 years (SD = 15.6) were included. A bipolar disorder diagnosis (n = 1479) was associated with higher rates of psychiatric conditions, injuries and poisonings, endocrine/metabolic and neurological conditions, viral hepatitis C, and asthma. BD-PRS was associated with psychiatric comorbidities but, in contrast, had no positive associations with general medical conditions. While our findings warrant confirmation with longitudinal-prospective studies, the limited associations between bipolar disorder genetics and medical conditions suggest that shared environmental effects or environmental consequences of diagnosis may have a greater impact on the general medical comorbidity profile of individuals with bipolar disorder than its genetic risk.

Identifying behaviour-related and physiological risk factors for suicide attempts in the UK Biobank

Suicide is a global public health challenge, yet considerable uncertainty remains regarding the associations of both behaviour-related and physiological factors with suicide attempts (SA). Here we first estimated polygenic risk scores (PRS) for SA in 334,706 UK Biobank participants and conducted phenome-wide association analyses considering 2,291 factors. We identified 246 (63.07%) behaviour-related and 200 (10.41%, encompassing neuroimaging, blood and metabolic biomarkers, and proteins) physiological factors significantly associated with SA-PRS, with robust associations observed in lifestyle factors and mental health. Further case-control analyses involving 3,558 SA cases and 149,976 controls mirrored behaviour-related associations observed with SA-PRS. Moreover, Mendelian randomization analyses supported a potential causal effect of liability to 58 factors on SA, such as age at first intercourse, neuroticism, smoking, overall health rating and depression. Notably, machine-learning classification models based on behaviour-related factors exhibited high discriminative accuracy in distinguishing those with and without SA (area under the receiver operating characteristic curve 0.909 ± 0.006). This study provides comprehensive insights into diverse risk factors for SA, shedding light on potential avenues for targeted prevention and intervention strategies.

Autoencoder to Identify Sex-Specific Sub-phenotypes in Alzheimer's Disease Progression Using Longitudinal Electronic Health Records

Alzheimer's Disease (AD) is a complex neurodegenerative disorder significantly influenced by sex differences, with approximately two-thirds of AD patients being women. Characterizing the sex-specific AD progression and identifying its progression trajectory is a crucial step to developing effective risk stratification and prevention strategies. In this study, we developed an autoencoder to uncover sex-specific sub-phenotypes in AD progression leveraging longitudinal electronic health record (EHR) data from OneFlorida+ Clinical Research Consortium. Specifically, we first constructed temporal patient representation using longitudinal EHRs from a sex-stratified AD cohort. We used a long short-term memory (LSTM)-based autoencoder to extract and generate latent representation embeddings from sequential clinical records of patients. We then applied hierarchical agglomerative clustering to the learned representations, grouping patients based on their progression sub-phenotypes. The experimental results show we successfully identified five primary sex-based AD sub-phenotypes with corresponding progression pathways with high confidence. These sex-specific sub-phenotypes not only illustrated distinct AD progression patterns but also revealed differences in clinical characteristics and comorbidities between females and males in AD development. These findings could provide valuable insights for advancing personalized AD intervention and treatment strategies.

Using electronic health records for clinical pharmacology research: Challenges and considerations

Electronic health records (EHRs) contain a vast array of phenotypic data on large numbers of individuals, often collected over decades. Due to the wealth of information, EHR data have emerged as a powerful resource to make first discoveries and identify disparities in our healthcare system. While the number of EHR-based studies has exploded in recent years, most of these studies are directed at associations with disease rather than pharmacotherapeutic outcomes, such as drug response or adverse drug reactions. This is largely due to challenges specific to deriving drug-related phenotypes from the EHR. There is great potential for EHR-based discovery in clinical pharmacology research, and there is a critical need to address specific challenges related to accurate and reproducible derivation of drug-related phenotypes from the EHR. This review provides a detailed evaluation of challenges and considerations for deriving drug-related data from EHRs. We provide an examination of EHR-based computable phenotypes and discuss cutting-edge approaches to map medication information for clinical pharmacology research, including medication-based computable phenotypes and natural language processing. We also discuss additional considerations such as data structure, heterogeneity and missing data, rare phenotypes, and diversity within the EHR. By further understanding the complexities associated with conducting clinical pharmacology research using EHR-based data, investigators will be better equipped to design thoughtful studies with more reproducible results. Progress in utilizing EHRs for clinical pharmacology research should lead to significant advances in our ability to understand differential drug response and predict adverse drug reactions.

Multi-ancestry meta-analysis of tobacco use disorder identifies 461 potential risk genes and reveals associations with multiple health outcomes

Tobacco use disorder (TUD) is the most prevalent substance use disorder in the world. Genetic factors influence smoking behaviours and although strides have been made using genome-wide association studies to identify risk variants, most variants identified have been for nicotine consumption, rather than TUD. Here we leveraged four US biobanks to perform a multi-ancestral meta-analysis of TUD (derived via electronic health records) in 653,790 individuals (495,005 European, 114,420 African American and 44,365 Latin American) and data from UK Biobank (ncombined = 898,680). We identified 88 independent risk loci; integration with functional genomic tools uncovered 461 potential risk genes, primarily expressed in the brain. TUD was genetically correlated with smoking and psychiatric traits from traditionally ascertained cohorts, externalizing behaviours in children and hundreds of medical outcomes, including HIV infection, heart disease and pain. This work furthers our biological understanding of TUD and establishes electronic health records as a source of phenotypic information for studying the genetics of TUD.

Prevalence and associations of multiple hypnotic prescriptions in a clinical sample

STUDY OBJECTIVES

We examined the prevalence of multiple hypnotic prescriptions and its association with clinical and demographic characteristics from the electronic health record (EHR) in the Mayo Clinic Biobank.

METHODS

Adult participants enrolled in the Mayo Clinic Biobank with an EHR number of ≥ 1 year were included (n = 52,940). Clinical and demographic characteristics were compared between participants who were and were not prescribed any hypnotic approved for insomnia by the US Food and Drug Administration and/or trazodone and in those prescribed a single vs multiple (≥ 2) hypnotics. A phenotype-based, phenome-wide association study (PheWAS) examining associations between hypnotic prescriptions and diagnoses across the EHR was performed adjusting for demographic and other confounders.

RESULTS

A total of 17,662 (33%) participants were prescribed at least 1 hypnotic and 5,331 (10%) received ≥ 2 hypnotics. Participants who were prescribed a hypnotic were more likely to be older, female, White, with a longer EHR, and a greater number of diagnostic codes (all P < .001). Those with multiple hypnotic prescriptions were more likely to be younger, female, with a longer EHR, and a greater number of diagnostic codes (all P < .001) compared with those prescribed a single hypnotic. The PheWAS revealed that participants with multiple hypnotic prescriptions had higher rates of mood disorders, anxiety disorders, suicidal ideation, restless legs syndrome, and chronic pain (all P < 1 e-10).

CONCLUSIONS

Receiving multiple hypnotic prescriptions is common and associated with a greater prevalence of psychiatric, chronic pain, and sleep-related movement disorders. Future studies should examine potential genetic associations with multiple hypnotic prescriptions to personalize treatments for chronic insomnia.

CITATION

Kolla BP, Mansukhani MP, Chakravorty S, Frank JA, Coombes BJ. Prevalence and associations of multiple hypnotic prescriptions in a clinical sample. J Clin Sleep Med. 2024;20(5):793-800.

Phenotyping Down syndrome: discovery and predictive modelling with electronic medical records

BACKGROUND

Individuals with Down syndrome (DS) have a heightened risk for various co-occurring health conditions, including congenital heart disease (CHD). In this two-part study, electronic medical records (EMRs) were leveraged to examine co-occurring health conditions among individuals with DS (Study 1) and to investigate health conditions linked to surgical intervention among DS cases with CHD (Study 2).

METHODS

De-identified EMRs were acquired from Vanderbilt University Medical Center and facilitated creating a cohort of N = 2282 DS cases (55% females), along with comparison groups for each study. In Study 1, DS cases were one-by-two sex and age matched with samples of case-controls and of individuals with other intellectual and developmental difficulties (IDDs). The phenome-disease association study (PheDAS) strategy was employed to reveal co-occurring health conditions in DS versus comparison groups, which were then ranked for how often they are discussed in relation to DS using the PubMed database and Novelty Finding Index. In Study 2, a subset of DS individuals with CHD [N = 1098 (48%)] were identified to create longitudinal data for N = 204 cases with surgical intervention (19%) versus 204 case-controls. Data were included in predictive models and assessed which model-based health conditions, when more prevalent, would increase the likelihood of surgical intervention.

RESULTS

In Study 1, relative to case-controls and those with other IDDs, co-occurring health conditions among individuals with DS were confirmed to include heart failure, pulmonary heart disease, atrioventricular block, heart transplant/surgery and primary pulmonary hypertension (circulatory); hypothyroidism (endocrine/metabolic); and speech and language disorder and Alzheimer's disease (neurological/mental). Findings also revealed more versus less prevalent co-occurring health conditions in individuals with DS when comparing with those with other IDDs. Findings with high Novelty Finding Index were abnormal electrocardiogram, non-rheumatic aortic valve disorders and heart failure (circulatory); acid-base balance disorder (endocrine/metabolism); and abnormal blood chemistry (symptoms). In Study 2, the predictive models revealed that among individuals with DS and CHD, presence of health conditions such as congestive heart failure (circulatory), valvular heart disease and cardiac shunt (congenital), and pleural effusion and pulmonary collapse (respiratory) were associated with increased likelihood of surgical intervention.

CONCLUSIONS

Research efforts using EMRs and rigorous statistical methods could shed light on the complexity in health profile among individuals with DS and other IDDs and motivate precision-care development.

A novel method for multiple phenotype association studies based on genotype and phenotype network

Joint analysis of multiple correlated phenotypes for genome-wide association studies (GWAS) can identify and interpret pleiotropic loci which are essential to understand pleiotropy in diseases and complex traits. Meanwhile, constructing a network based on associations between phenotypes and genotypes provides a new insight to analyze multiple phenotypes, which can explore whether phenotypes and genotypes might be related to each other at a higher level of cellular and organismal organization. In this paper, we first develop a bipartite signed network by linking phenotypes and genotypes into a Genotype and Phenotype Network (GPN). The GPN can be constructed by a mixture of quantitative and qualitative phenotypes and is applicable to binary phenotypes with extremely unbalanced case-control ratios in large-scale biobank datasets. We then apply a powerful community detection method to partition phenotypes into disjoint network modules based on GPN. Finally, we jointly test the association between multiple phenotypes in a network module and a single nucleotide polymorphism (SNP). Simulations and analyses of 72 complex traits in the UK Biobank show that multiple phenotype association tests based on network modules detected by GPN are much more powerful than those without considering network modules. The newly proposed GPN provides a new insight to investigate the genetic architecture among different types of phenotypes. Multiple phenotypes association studies based on GPN are improved by incorporating the genetic information into the phenotype clustering. Notably, it might broaden the understanding of genetic architecture that exists between diagnoses, genes, and pleiotropy.

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 R R that overcomes the inaccuracy in the controls via theoretical formula. With simulation and real data application, we further confirmed that R R is unbiased without compiling exclusion criteria lists. With R R 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.

Novel strategies for drug repurposing

Synthetic biology, precision medicine, and nanobiotechnology are the three main emerging areas that drive translational innovation toward commercialization. There are several strategies used in precision medicine and drug repurposing is one of the key approaches as it addresses the challenges in drug discovery (high cost and time). Here, we provide a perspective on various new approaches to drug repurposing for cancer precision medicine. We report here our optimized wound healing methodology that can be used to validate drug sensitivity and drug repurposing. Using HeLa as our benchmark, we demonstrated that the assay can be applied to identify drugs that limit cell proliferation. From a future perspective, this assay can be expanded to ex vivo culturing of solid tumors in 2D culture and leukemia in 3D culture.

Unbiased characterization of atrial fibrillation phenotypic architecture provides insight to genetic liability and clinically relevant outcomes

Background

Atrial Fibrillation (AF) is a common and clinically heterogeneous arrythmia. Machine learning (ML) algorithms can define data-driven disease subtypes in an unbiased fashion, but whether the AF subgroups defined in this way align with underlying mechanisms, such as high polygenic liability to AF or inflammation, and associate with clinical outcomes is unclear.

Methods

We identified individuals with AF in a large biobank linked to electronic health records (EHR) and genome-wide genotyping. The phenotypic architecture in the AF cohort was defined using principal component analysis of 35 expertly curated and uncorrelated clinical features. We applied an unsupervised co-clustering machine learning algorithm to the 35 features to identify distinct phenotypic AF clusters. The clinical inflammatory status of the clusters was defined using measured biomarkers (CRP, ESR, WBC, Neutrophil %, Platelet count, RDW) within 6 months of first AF mention in the EHR. Polygenic risk scores (PRS) for AF and cytokine levels were used to assess genetic liability of clusters to AF and inflammation, respectively. Clinical outcomes were collected from EHR up to the last medical contact.

Results

The analysis included 23,271 subjects with AF, of which 6,023 had available genome-wide genotyping. The machine learning algorithm identified 3 phenotypic clusters that were distinguished by increasing prevalence of comorbidities, particularly renal dysfunction, and coronary artery disease. Polygenic liability to AF across clusters was highest in the low comorbidity cluster. Clinically measured inflammatory biomarkers were highest in the high comorbid cluster, while there was no difference between groups in genetically predicted levels of inflammatory biomarkers. Subgroup assignment was associated with multiple clinical outcomes including mortality, stroke, bleeding, and use of cardiac implantable electronic devices after AF diagnosis.

Conclusion

Patient subgroups identified by unsupervised clustering were distinguished by comorbidity burden and associated with risk of clinically important outcomes. Polygenic liability to AF across clusters was greatest in the low comorbidity subgroup. Clinical inflammation, as reflected by measured biomarkers, was lowest in the subgroup with lowest comorbidities. However, there were no differences in genetically predicted levels of inflammatory biomarkers, suggesting associations between AF and inflammation is driven by acquired comorbidities rather than genetic predisposition.

Associations of a multidimensional polygenic sleep health score and a sleep lifestyle index on health outcomes and their interaction in a clinical biobank

Sleep is a complex behavior regulated by genetic and environmental factors, and is known to influence health outcomes. However, the effect of multidimensional sleep encompassing several sleep dimensions on diseases has yet to be fully elucidated. Using the Mass General Brigham Biobank, we aimed to examine the association of multidimensional sleep with health outcomes and investigate whether sleep behaviors modulate genetic predisposition to unfavorable sleep on mental health outcomes. First, we generated a Polygenic Sleep Health Score using previously identified single nucleotide polymorphisms for sleep health and constructed a Sleep Lifestyle Index using data from self-reported sleep questions and electronic health records; second, we performed phenome-wide association analyses between these indexes and clinical phenotypes; and third, we analyzed the interaction between the indexes on prevalent mental health outcomes. Fifteen thousand eight hundred and eighty-four participants were included in the analysis (mean age 54.4; 58.6% female). The Polygenic Sleep Health Score was associated with the Sleep Lifestyle Index (β=0.050, 95%CI=0.032, 0.068) and with 114 disease outcomes spanning 12 disease groups, including obesity, sleep, and substance use disease outcomes (p<3.3×10-5). The Sleep Lifestyle Index was associated with 458 disease outcomes spanning 17 groups, including sleep, mood, and anxiety disease outcomes (p<5.1×10-5). No interactions were found between the indexes on prevalent mental health outcomes. These findings suggest that favorable sleep behaviors and genetic predisposition to healthy sleep may independently be protective of disease outcomes. This work provides novel insights into the role of multidimensional sleep on population health and highlights the need to develop prevention strategies focused on healthy sleep habits.

Genome-wide association analysis of composite sleep health scores in 413,904 individuals

Recent genome-wide association studies (GWASs) of several individual sleep traits have identified hundreds of genetic loci, suggesting diverse mechanisms. Moreover, sleep traits are moderately correlated, and together may provide a more complete picture of sleep health, while also illuminating distinct domains. Here we construct novel sleep health scores (SHSs) incorporating five core self-report measures: sleep duration, insomnia symptoms, chronotype, snoring, and daytime sleepiness, using additive (SHS-ADD) and five principal components-based (SHS-PCs) approaches. GWASs of these six SHSs identify 28 significant novel loci adjusting for multiple testing on six traits (p<8.3e-9), along with 341 previously reported loci (p<5e-08). The heritability of the first three SHS-PCs equals or exceeds that of SHS-ADD (SNP-h2=0.094), while revealing sleep-domain-specific genetic discoveries. Significant loci enrich in multiple brain tissues and in metabolic and neuronal pathways. Post GWAS analyses uncover novel genetic mechanisms underlying sleep health and reveal connections to behavioral, psychological, and cardiometabolic traits.

Alzheimer's disease phenotype based upon the carrier status of the apolipoprotein E ɛ4 allele

The apolipoprotein E ɛ4 allele (APOE4) is universally acknowledged as the most potent genetic risk factor for Alzheimer's disease (AD). APOE4 promotes the initiation and progression of AD. Although the underlying mechanisms are unclearly understood, differences in lipid-bound affinity among the three APOE isoforms may constitute the basis. The protein APOE4 isoform has a high affinity with triglycerides and cholesterol. A distinction in lipid metabolism extensively impacts neurons, microglia, and astrocytes. APOE4 carriers exhibit phenotypic differences from non-carriers in clinical examinations and respond differently to multiple treatments. Therefore, we hypothesized that phenotypic classification of AD patients according to the status of APOE4 carrier will help specify research and promote its use in diagnosing and treating AD. Recent reviews have mainly evaluated the differences between APOE4 allele carriers and non-carriers from gene to protein structures, clinical features, neuroimaging, pathology, the neural network, and the response to various treatments, and have provided the feasibility of phenotypic group classification based on APOE4 carrier status. This review will facilitate the application of APOE phenomics concept in clinical practice and promote further medical research on AD.

Phenome-wide association study and precision medicine of cardiovascular diseases in the post-COVID-19 era

SARS-CoV-2 infection causes injuries of not only the lungs but also the heart and endothelial cells in vasculature of multiple organs, and induces systemic inflammation and immune over-reactions, which makes COVID-19 a disease phenome that simultaneously affects multiple systems. Cardiovascular diseases (CVD) are intrinsic risk and causative factors for severe COVID-19 comorbidities and death. The wide-spread infection and reinfection of SARS-CoV-2 variants and the long-COVID may become a new common threat to human health and propose unprecedented impact on the risk factors, pathophysiology, and pharmacology of many diseases including CVD for a long time. COVID-19 has highlighted the urgent demand for precision medicine which needs new knowledge network to innovate disease taxonomy for more precise diagnosis, therapy, and prevention of disease. A deeper understanding of CVD in the setting of COVID-19 phenome requires a paradigm shift from the current phenotypic study that focuses on the virus or individual symptoms to phenomics of COVID-19 that addresses the inter-connectedness of clinical phenotypes, i.e., clinical phenome. Here, we summarize the CVD manifestations in the full clinical spectrum of COVID-19, and the phenome-wide association study of CVD interrelated to COVID-19. We discuss the underlying biology for CVD in the COVID-19 phenome and the concept of precision medicine with new phenomic taxonomy that addresses the overall pathophysiological responses of the body to the SARS-CoV-2 infection. We also briefly discuss the unique taxonomy of disease as Zheng-hou patterns in traditional Chinese medicine, and their potential implications in precision medicine of CVD in the post-COVID-19 era.

Multi-ancestry meta-analysis of tobacco use disorder prioritizes novel candidate risk genes and reveals associations with numerous health outcomes

Tobacco use disorder (TUD) is the most prevalent substance use disorder in the world. Genetic factors influence smoking behaviors, and although strides have been made using genome-wide association studies (GWAS) to identify risk variants, the majority of variants identified have been for nicotine consumption, rather than TUD. We leveraged five biobanks to perform a multi-ancestral meta-analysis of TUD (derived via electronic health records, EHR) in 898,680 individuals (739,895 European, 114,420 African American, 44,365 Latin American). We identified 88 independent risk loci; integration with functional genomic tools uncovered 461 potential risk genes, primarily expressed in the brain. TUD was genetically correlated with smoking and psychiatric traits from traditionally ascertained cohorts, externalizing behaviors in children, and hundreds of medical outcomes, including HIV infection, heart disease, and pain. This work furthers our biological understanding of TUD and establishes EHR as a source of phenotypic information for studying the genetics of TUD.

Association Between Fuchs Endothelial Corneal Dystrophy, Diabetes Mellitus, and Multimorbidity

PURPOSE

The aim of this study was to assess risk for demographic variables and other health conditions that are associated with Fuchs endothelial corneal dystrophy (FECD).

METHODS

We developed a FECD case-control algorithm based on structured electronic health record data and confirmed accuracy by individual review of charts at 3 Veterans Affairs (VA) Medical Centers. This algorithm was applied to the Department of VA Million Veteran Program cohort from whom sex, genetic ancestry, comorbidities, diagnostic phecodes, and laboratory values were extracted. Single-variable and multiple variable logistic regression models were used to determine the association of these risk factors with FECD diagnosis.

RESULTS

Being a FECD case was associated with female sex, European genetic ancestry, and a greater number of comorbidities. Of 1417 diagnostic phecodes evaluated, 213 had a significant association with FECD, falling in both ocular and nonocular conditions, including diabetes mellitus (DM). Five of 69 laboratory values were associated with FECD, with the direction of change for 4 being consistent with DM. Insulin dependency and type 1 DM raised risk to a greater degree than type 2 DM, like other microvascular diabetic complications.

CONCLUSIONS

Female sex, European ancestry, and multimorbidity increased FECD risk. Endocrine/metabolic clinic encounter codes and altered patterns of laboratory values support DM increasing FECD risk. Our results evoke a threshold model in which the FECD phenotype is intensified by DM and potentially other health conditions that alter corneal physiology. Further studies to better understand the relationship between FECD and DM are indicated and may help identify opportunities for slowing FECD progression.

Therapeutic outcomes wide association scan of different antipsychotics in patients with schizophrenia: Randomized clinical trials and multi-ancestry validation

AIM

This study identified discrepant therapeutic outcomes of antipsychotics.

METHODS

A total of 5191 patients with schizophrenia were enrolled, 3030 as discovery cohort, 1395 as validation cohort, and 766 as multi-ancestry validation cohort. Therapeutic Outcomes Wide Association Scan was conducted. Types of antipsychotics (one antipsychotic vs other antipsychotics) were dependent variables, therapeutic outcomes including efficacy and safety were independent variables.

RESULTS

In discovery cohort, olanzapine related to higher risk of weight gain (AIWG, OR: 2.21-2.86), liver dysfunction (OR: 1.75-2.33), sedation (OR: 1.76-2.86), increased lipid level (OR: 2.04-2.12), and lower risk of extrapyramidal syndrome (EPS, OR: 0.14-0.46); risperidone related to higher risk of hyperprolactinemia (OR: 12.45-20.53); quetiapine related to higher risk of sedation (OR = 1.73), palpitation (OR = 2.87), increased lipid level (OR = 1.69), lower risk of hyperprolactinemia (OR: 0.09-0.11), and EPS (OR: 0.15-0.44); aripiprazole related to lower risk of hyperprolactinemia (OR: 0.09-0.14), AIWG (OR = 0.44), sedation (OR: 0.33-0.47), and QTc prolongation (β = -2.17); ziprasidone related to higher risk of increased QT interval (β range: 3.11-3.22), nausea (OR: 3.22-3.91), lower risk of AIWG (OR: 0.27-0.46), liver dysfunction (OR: 0.41-0.38), and increased lipid level (OR: 0.41-0.55); haloperidol related to higher risk of EPS (OR: 2.64-6.29), hyperprolactinemia (OR: 5.45-9.44), and increased salivation (OR: 3.50-3.68). Perphenazine related to higher risk of EPS (OR: 1.89-2.54). Higher risk of liver dysfunction in olanzapine and lower risk of hyperprolactinemia in aripiprazole were confirmed in validation cohort, and higher risk of AIWG in olanzapine and hyperprolactinemia in risperidone were confirmed in multi-ancestry validation cohort.

CONCLUSION

Future precision medicine should focus on personalized side-effects.

A phenome-wide association study (PheWAS) to identify the health impacts of 4-cresol sulfate in the Nagahama Study

Gut-microbiota derived metabolites are important regulators of host biology and metabolism. To understand the impacts of the microbial metabolite 4-cresol sulfate (4-CS) on four chronic diseases [type 2 diabetes mellitus, metabolic syndrome (MetS), non-alcoholic fatty liver disease, and chronic kidney disease (CKD)], we conducted association analyses of plasma 4-CS quantified by liquid chromatography coupled to mass spectrometry (LC-MS) in 3641 participants of the Nagahama study. Our results validated the elevation of 4-CS in CKD and identified a reducing trend in MetS. To delineate the holistic effects of 4-CS, we performed a phenome-wide association analysis (PheWAS) with 937 intermediate biological and behavioral traits. We detected associations between 4-CS and 39 phenotypes related to blood pressure regulation, hepatic and renal functions, hematology, sleep quality, intraocular pressure, ion regulation, ketone and fatty acid metabolisms, disease history and dietary habits. Among them, 19 PheWAS significant traits, including fatty acids and 14 blood pressure indices, were correlated with MetS, suggesting that 4-CS is a potential biomarker for MetS. Consistent associations of this gut microbial-derived metabolite on multiple endophenotypes underlying distinct etiopathogenesis support its role in the overall host health, with prospects of probiotic-based therapeutic solutions in chronic diseases.

Brief Report: The Characterization of Medical Comorbidity Prior to Autism Diagnosis in Children Before Age Two

In autism spectrum disorder (ASD), medical conditions in infancy could be predictive markers for later ASD diagnosis. In this study, electronic medical records of 579 autistic individuals and 1897 matched controls prior to age 2 were analyzed for potential predictive conditions. Using a novel tool, the relative association of each condition in the autistic group was compared to the control group using logistic regressions across medical records. Generalized convulsive epilepsy, nystagmus, lack of normal physiological development, delayed milestones, and strabismus were more likely in those later diagnosed with ASD while perinatal jaundice was less likely to be associated. Lesser-known conditions, such as strabismus and nystagmus, may point to novel predictive co-occurring condition profiles which could improve screening practices for ASD.

Schizophrenia and co-morbidity risk: Evidence from a data driven phenomewide association study

Schizophrenia is a chronic debilitating psychiatric disorder with significant morbidity and mortality. In this study, we used information from 337,484 UK Biobank participants and performed PheWAS using schizophrenia genetic risk score on 1135 disease outcomes. Signals that passed the false discovery rate threshold were further analyzed for evidence on the causality of the association. We extended the analysis to 30 serum, four urine, and six neuroimaging biomarkers to identify biomarkers that could be affected by schizophrenia. Schizophrenia GRS was associated with 54 (39 distinct) disease outcomes including schizophrenia in the PheWAS analysis. Of these, a causal association were found with 10 distinct diseases in the MR analysis. Schizophrenia causally linked with higher odds of anxiety (OR = 1.41, 95%CI 1.12 to 1.21), bipolar disorder (OR = 1.52, 95%CI 1.36 to 1.70), major depressive disorder (OR = 1.12, 95%CI 1.08 to 1.16) and suicidal ideation (OR = 1.30, 95%CI 1.19 to 1.42). Lower odds were found for several diseases including type 1 diabetes, coronary atherosclerosis and some musculoskeletal disorders. In analyses using biomarkers, schizophrenia was associated with lower serum 25(OH)D, gamma glutamyltransferase, cystatin C, serum creatinine. However, we did not find association with any of the brain imaging markers. Our analyses confirmed the co-existence of schizophrenia with other mental health disorders but did not otherwise suggest strong effects on disease risk. Biomarker analyses reflected associations which could be explained by unhealthy lifestyles, suggesting patients with schizophrenia may benefit from screening for and managing broader health aspects.

Ontologizing health systems data at scale: making translational discovery a reality

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.

The Human Extracellular Matrix Diseasome Reveals Genotype-Phenotype Associations with Clinical Implications for Age-Related Diseases

The extracellular matrix (ECM) is earning an increasingly relevant role in many disease states and aging. The analysis of these disease states is possible with the GWAS and PheWAS methodologies, and through our analysis, we aimed to explore the relationships between polymorphisms in the compendium of ECM genes (i.e., matrisome genes) in various disease states. A significant contribution on the part of ECM polymorphisms is evident in various types of disease, particularly those in the core-matrisome genes. Our results confirm previous links to connective-tissue disorders but also unearth new and underexplored relationships with neurological, psychiatric, and age-related disease states. Through our analysis of the drug indications for gene-disease relationships, we identify numerous targets that may be repurposed for age-related pathologies. The identification of ECM polymorphisms and their contributions to disease will play an integral role in future therapeutic developments, drug repurposing, precision medicine, and personalized care.

Joint analysis of multiple phenotypes for extremely unbalanced case-control association studies

In genome-wide association studies (GWAS) for thousands of phenotypes in biobanks, most binary phenotypes have substantially fewer cases than controls. Many widely used approaches for joint analysis of multiple phenotypes produce inflated type I error rates for such extremely unbalanced case-control phenotypes. In this research, we develop a method to jointly analyze multiple unbalanced case-control phenotypes to circumvent this issue. We first group multiple phenotypes into different clusters based on a hierarchical clustering method, then we merge phenotypes in each cluster into a single phenotype. In each cluster, we use the saddlepoint approximation to estimate the p value of an association test between the merged phenotype and a single nucleotide polymorphism (SNP) which eliminates the issue of inflated type I error rate of the test for extremely unbalanced case-control phenotypes. Finally, we use the Cauchy combination method to obtain an integrated p value for all clusters to test the association between multiple phenotypes and a SNP. We use extensive simulation studies to evaluate the performance of the proposed approach. The results show that the proposed approach can control type I error rate very well and is more powerful than other available methods. We also apply the proposed approach to phenotypes in category IX (diseases of the circulatory system) in the UK Biobank. We find that the proposed approach can identify more significant SNPs than the other viable methods we compared with.

Using human genetics to improve safety assessment of therapeutics

Human genetics research has discovered thousands of proteins associated with complex and rare diseases. Genome-wide association studies (GWAS) and studies of Mendelian disease have resulted in an increased understanding of the role of gene function and regulation in human conditions. Although the application of human genetics has been explored primarily as a method to identify potential drug targets and support their relevance to disease in humans, there is increasing interest in using genetic data to identify potential safety liabilities of modulating a given target. Human genetic variants can be used as a model to anticipate the effect of lifelong modulation of therapeutic targets and identify the potential risk for on-target adverse events. This approach is particularly useful for non-clinical safety evaluation of novel therapeutics that lack pharmacologically relevant animal models and can contribute to the intrinsic safety profile of a drug target. This Review illustrates applications of human genetics to safety studies during drug discovery and development, including assessing the potential for on- and off-target associated adverse events, carcinogenicity risk assessment, and guiding translational safety study designs and monitoring strategies. A summary of available human genetic resources and recommended best practices is provided. The challenges and future perspectives of translating human genetic information to identify risks for potential drug effects in preclinical and clinical development are discussed.

New insights from the last decade of research in psychiatric genetics: discoveries, challenges and clinical implications

Psychiatric genetics has made substantial progress in the last decade, providing new insights into the genetic etiology of psychiatric disorders, and paving the way for precision psychiatry, in which individual genetic profiles may be used to personalize risk assessment and inform clinical decision-making. Long recognized to be heritable, recent evidence shows that psychiatric disorders are influenced by thousands of genetic variants acting together. Most of these variants are commonly occurring, meaning that every individual has a genetic risk to each psychiatric disorder, from low to high. A series of large-scale genetic studies have discovered an increasing number of common and rare genetic variants robustly associated with major psychiatric disorders. The most convincing biological interpretation of the genetic findings implicates altered synaptic function in autism spectrum disorder and schizophrenia. However, the mechanistic understanding is still incomplete. In line with their extensive clinical and epidemiological overlap, psychiatric disorders appear to exist on genetic continua and share a large degree of genetic risk with one another. This provides further support to the notion that current psychiatric diagnoses do not represent distinct pathogenic entities, which may inform ongoing attempts to reconceptualize psychiatric nosology. Psychiatric disorders also share genetic influences with a range of behavioral and somatic traits and diseases, including brain structures, cognitive function, immunological phenotypes and cardiovascular disease, suggesting shared genetic etiology of potential clinical importance. Current polygenic risk score tools, which predict individual genetic susceptibility to illness, do not yet provide clinically actionable information. However, their precision is likely to improve in the coming years, and they may eventually become part of clinical practice, stressing the need to educate clinicians and patients about their potential use and misuse. This review discusses key recent insights from psychiatric genetics and their possible clinical applications, and suggests future directions.

Scanning the medical phenome to identify new diagnoses after recovery from COVID-19 in a US cohort

OBJECTIVE

COVID-19 survivors are at risk for long-term health effects, but assessing the sequelae of COVID-19 at large scales is challenging. High-throughput methods to efficiently identify new medical problems arising after acute medical events using the electronic health record (EHR) could improve surveillance for long-term consequences of acute medical problems like COVID-19.

MATERIALS AND METHODS

We augmented an existing high-throughput phenotyping method (PheWAS) to identify new diagnoses occurring after an acute temporal event in the EHR. We then used the temporal-informed phenotypes to assess development of new medical problems among COVID-19 survivors enrolled in an EHR cohort of adults tested for COVID-19 at Vanderbilt University Medical Center.

RESULTS

The study cohort included 186 105 adults tested for COVID-19 from March 5, 2020 to November 1, 2021; of which 30 088 (16.2%) tested positive. Median follow-up after testing was 412 days (IQR 274-528). Our temporal-informed phenotyping was able to distinguish phenotype chapters based on chronicity of their constituent diagnoses. PheWAS with temporal-informed phenotypes identified increased risk for 43 diagnoses among COVID-19 survivors during outpatient follow-up, including multiple new respiratory, cardiovascular, neurological, and pregnancy-related conditions. Findings were robust to sensitivity analyses, and several phenotypic associations were supported by changes in outpatient vital signs or laboratory tests from the pretesting to postrecovery period.

CONCLUSION

Temporal-informed PheWAS identified new diagnoses affecting multiple organ systems among COVID-19 survivors. These findings can inform future efforts to enable longitudinal health surveillance for survivors of COVID-19 and other acute medical conditions using the EHR.

POPDx: an automated framework for patient phenotyping across 392 246 individuals in the UK Biobank study

OBJECTIVE

For the UK Biobank, standardized phenotype codes are associated with patients who have been hospitalized but are missing for many patients who have been treated exclusively in an outpatient setting. We describe a method for phenotype recognition that imputes phenotype codes for all UK Biobank participants.

MATERIALS AND METHODS

POPDx (Population-based Objective Phenotyping by Deep Extrapolation) is a bilinear machine learning framework for simultaneously estimating the probabilities of 1538 phenotype codes. We extracted phenotypic and health-related information of 392 246 individuals from the UK Biobank for POPDx development and evaluation. A total of 12 803 ICD-10 diagnosis codes of the patients were converted to 1538 phecodes as gold standard labels. The POPDx framework was evaluated and compared to other available methods on automated multiphenotype recognition.

RESULTS

POPDx can predict phenotypes that are rare or even unobserved in training. We demonstrate substantial improvement of automated multiphenotype recognition across 22 disease categories, and its application in identifying key epidemiological features associated with each phenotype.

CONCLUSIONS

POPDx helps provide well-defined cohorts for downstream studies. It is a general-purpose method that can be applied to other biobanks with diverse but incomplete data.

Genetically proxied therapeutic prolyl-hydroxylase inhibition and cardiovascular risk

Prolyl hydroxylase (PHD) inhibitors are in clinical development for anaemia in chronic kidney disease. Epidemiological studies have reported conflicting results regarding safety of long-term therapeutic haemoglobin (Hgb) rises through PHD inhibition on risk of cardiovascular disease. Genetic variation in genes encoding PHDs can be used as partial proxies to investigate the potential effects of long-term Hgb rises. We used Mendelian randomization to investigate the effect of long-term Hgb level rises through genetically proxied PHD inhibition on coronary artery disease (CAD: 60 801 cases; 123 504 controls), myocardial infarction (MI: 42 561 cases; 123 504 controls) or stroke (40 585 cases; 406 111 controls). To further characterize long-term effects of Hgb level rises, we performed a phenome-wide association study (PheWAS) in up to 451 099 UK Biobank individuals. Genetically proxied therapeutic PHD inhibition, equivalent to a 1.00 g/dl increase in Hgb levels, was not associated (at P < 0.05) with increased odds of CAD; odd ratio (OR) [95% confidence intervals (CI)] = 1.06 (0.84, 1.35), MI [OR (95% CI) = 1.02 (0.79, 1.33)] or stroke [OR (95% CI) = 0.91 (0.66, 1.24)]. PheWAS revealed associations with blood related phenotypes consistent with EGLN's role, relevant kidney- and liver-related biomarkers like estimated glomerular filtration rate and microalbuminuria, and non-alcoholic fatty liver disease (Bonferroni-adjusted P < 5.42E-05) but these were not clinically meaningful. These findings suggest that long-term alterations in Hgb through PHD inhibition are unlikely to substantially increase cardiovascular disease risk; using large disease genome-wide association study data, we could exclude ORs of 1.35 for cardiovascular risk with a 1.00 g/dl increase in Hgb.

HCLC-FC: A novel statistical method for phenome-wide association studies

The emergence of genetic data coupled to longitudinal electronic medical records (EMRs) offers the possibility of phenome-wide association studies (PheWAS). In PheWAS, the whole phenome can be divided into numerous phenotypic categories according to the genetic architecture across phenotypes. Currently, statistical analyses for PheWAS are mainly univariate analyses, which test the association between one genetic variant and one phenotype at a time. In this article, we derived a novel and powerful multivariate method for PheWAS. The proposed method involves three steps. In the first step, we apply the bottom-up hierarchical clustering method to partition a large number of phenotypes into disjoint clusters within each phenotypic category. In the second step, the clustering linear combination method is used to combine test statistics within each category based on the phenotypic clusters and obtain p-values from each phenotypic category. In the third step, we propose a new false discovery rate (FDR) control approach. We perform extensive simulation studies to compare the performance of our method with that of other existing methods. The results show that our proposed method controls FDR very well and outperforms other methods we compared with. We also apply the proposed approach to a set of EMR-based phenotypes across more than 300,000 samples from the UK Biobank. We find that the proposed approach not only can well-control FDR at a nominal level but also successfully identify 1,244 significant SNPs that are reported to be associated with some phenotypes in the GWAS catalog. Our open-access tools and instructions on how to implement HCLC-FC are available at https://github.com/XiaoyuLiang/HCLCFC.

The Role of Systems Biology in Deciphering Asthma Heterogeneity

Asthma is one of the most common and lifelong and chronic inflammatory diseases characterized by inflammation, bronchial hyperresponsiveness, and airway obstruction episodes. It is a heterogeneous disease of varying and overlapping phenotypes with many confounding factors playing a role in disease susceptibility and management. Such multifactorial disorders will benefit from using systems biology as a strategy to elucidate molecular insights from complex, quantitative, massive clinical, and biological data that will help to understand the underlying disease mechanism, early detection, and treatment planning. Systems biology is an approach that uses the comprehensive understanding of living systems through bioinformatics, mathematical, and computational techniques to model diverse high-throughput molecular, cellular, and the physiologic profiling of healthy and diseased populations to define biological processes. The use of systems biology has helped understand and enrich our knowledge of asthma heterogeneity and molecular basis; however, such methods have their limitations. The translational benefits of these studies are few, and it is recommended to reanalyze the different studies and omics in conjugation with one another which may help understand the reasons for this variation and help overcome the limitations of understanding the heterogeneity in asthma pathology. In this review, we aim to show the different factors that play a role in asthma heterogeneity and how systems biology may aid in understanding and deciphering the molecular basis of asthma.

An integrated genome and phenome-wide association study approach to understanding Alzheimer's disease predisposition

Genome-wide association studies (GWAS) have identified common single nucleotide polymorphisms (SNPs) that increase late-onset Alzheimer's disease (LOAD) risk. To identify additional LOAD-associated variants and provide insight into underlying disease biology, we performed a phenome-wide association study on 23 known LOAD-associated SNPs and 4:1 matched control SNPs using UK Biobank data. LOAD-associated SNPs were significantly enriched for associations with 8/778 queried traits, including 3 platelet traits. The strongest enrichment was for platelet distribution width (PDW) (p = 1.2 × 10-5), but increased PDW was not associated with LOAD susceptibility in Mendelian randomization analysis. Of 384 PDW-associated SNPs identified by prior GWAS, 36 were nominally associated with LOAD risk (17,008 cases; 37,154 controls) and 5 survived false-discovery rate correction. Associations confirmed known LOAD risk loci near PICALM, CD2AP, SPI1, and NDUFAF6, and identified a novel risk locus in epidermal growth factor receptor. Integrating GWAS and phenome-wide association study data reveals substantial pleiotropy between genetic determinants of LOAD and of platelet morphology, and for the first time implicates epidermal growth factor receptor - a mediator of β-amyloid toxicity - in Alzheimer's disease susceptibility.

Using the Phecode System to Identify the Preoperative Clinical Phenotypes Associated with Surgical Site Infection in Patients Undergoing Primary Total Knee Arthroplasty: The Sex Differences

Sex-related differences among comorbid conditions associated with surgical site infection (SSI) after total knee arthroplasty (TKA) are unclear. This population-based cohort study used a novel approach with a Phecode system to evaluate preoperative clinical phenotypes (i.e., comorbid conditions) associated with SSI after TKA and delineate sex-related differences in phenotypes. Using the Taiwan National Health Insurance Research Database (2014-2018), 83,870 patients who underwent TKA were identified. Demographic and SSI data during the 90-day postoperative follow-up were obtained. Comorbidities identified by the International Classification of Diseases within 1 year before TKA were recorded and mapped into Phecodes representing phenotypes. The overall rate of 90-day SSI was 1.3%. In total, 1663 phenotypes were identified among 83,870 patients-1585 and 1458 phenotypes for female (n = 62,018) and male (n = 21,852) patients, respectively. According to multivariate logistic regression analysis, the SSI odds ratio significantly increased with the presence of each of the 16 phenotypes. Subgroup analysis revealed that the presence of 10 and 4 phenotypes significantly increased SSI risk in both sexes; only one phenotype was common to both sexes. Therefore, comorbid conditions and sex should be considered in preoperative SSI risk evaluation in patients undergoing primary TKA. These findings provide new perspectives on susceptibility, prevention, and treatment in these patients.

Questionnaire-based exposome-wide association studies (ExWAS) reveal expected and novel risk factors associated with cardiovascular outcomes in the Personalized Environment and Genes Study

While multiple factors are associated with cardiovascular disease (CVD), many environmental exposures that may contribute to CVD have not been examined. To understand environmental effects on cardiovascular health, we performed an exposome-wide association study (ExWAS), a hypothesis-free approach, using survey data on endogenous and exogenous exposures at home and work and data from health and medical histories from the North Carolina-based Personalized Environment and Genes Study (PEGS) (n = 5015). We performed ExWAS analyses separately on six cardiovascular outcomes (cardiac arrhythmia, congestive heart failure, coronary artery disease, heart attack, stroke, and a combined atherogenic-related outcome comprising angina, angioplasty, atherosclerosis, coronary artery disease, heart attack, and stroke) using logistic regression and a false discovery rate of 5%. For each CVD outcome, we tested 502 single exposures and built multi-exposure models using the deletion-substitution-addition (DSA) algorithm. To evaluate complex nonlinear relationships, we employed the knockoff boosted tree (KOBT) algorithm. We adjusted all analyses for age, sex, race, BMI, and annual household income. ExWAS analyses revealed novel associations that include blood type A (Rh-) with heart attack (OR[95%CI] = 8.2[2.2:29.7]); paint exposures with stroke (paint related chemicals: 6.1[2.2:16.0], acrylic paint: 8.1[2.6:22.9], primer: 6.7[2.2:18.6]); biohazardous materials exposure with arrhythmia (1.8[1.5:2.3]); and higher paternal education level with reduced risk of multiple CVD outcomes (stroke, heart attack, coronary artery disease, and combined atherogenic outcome). In multi-exposure models, trouble sleeping and smoking remained important risk factors. KOBT identified significant nonlinear effects of sleep disorder, regular intake of grapefruit, and a family history of blood clotting problems for multiple CVD outcomes (combined atherogenic outcome, congestive heart failure, and coronary artery disease). In conclusion, using statistics and machine learning, these findings identify novel potential risk factors for CVD, enable hypothesis generation, provide insights into the complex relationships between risk factors and CVD, and highlight the importance of considering multiple exposures when examining CVD outcomes.

Challenges and opportunities associated with rare-variant pharmacogenomics

Recent advances in next-generation sequencing (NGS) have resulted in the identification of tens of thousands of rare pharmacogenetic variations with unknown functional effects. However, although such pharmacogenetic variations have been estimated to account for a considerable amount of the heritable variability in drug response and toxicity, accurate interpretation at the level of the individual patient remains challenging. We discuss emerging strategies and concepts to close this translational gap. We illustrate how massively parallel experimental assays, artificial intelligence (AI), and machine learning can synergize with population-scale biobank projects to facilitate the interpretation of NGS data to individualize clinical decision-making and personalized medicine.

Predictive Modelling in Clinical Bioinformatics: Key Concepts for Startups

Clinical bioinformatics is a newly emerging field that applies bioinformatics techniques for facilitating the identification of diseases, discovery of biomarkers, and therapy decision. Mathematical modelling is part of bioinformatics analysis pipelines and a fundamental step to extract clinical insights from genomes, transcriptomes and proteomes of patients. Often, the chosen modelling techniques relies on either statistical, machine learning or deterministic approaches. Research that combines bioinformatics with modelling techniques have been generating innovative biomedical technology, algorithms and models with biotech applications, attracting private investment to develop new business; however, startups that emerge from these technologies have been facing difficulties to implement clinical bioinformatics pipelines, protect their technology and generate profit. In this commentary, we discuss the main concepts that startups should know for enabling a successful application of predictive modelling in clinical bioinformatics. Here we will focus on key modelling concepts, provide some successful examples and briefly discuss the modelling framework choice. We also highlight some aspects to be taken into account for a successful implementation of cost-effective bioinformatics from a business perspective.

The causes and consequences of Alzheimer's disease: phenome-wide evidence from Mendelian randomization

Alzheimer's disease (AD) has no proven causal and modifiable risk factors, or effective interventions. We report a phenome-wide association study (PheWAS) of genetic liability for AD in 334,968 participants of the UK Biobank study, stratified by age. We also examined the effects of AD genetic liability on previously implicated risk factors. We replicated these analyses in the HUNT study. PheWAS hits and previously implicated risk factors were followed up in a Mendelian randomization (MR) framework to identify the causal effect of each risk factor on AD risk. A higher genetic liability for AD was associated with medical history and cognitive, lifestyle, physical and blood-based measures as early as 39 years of age. These effects were largely driven by the APOE gene. The follow-up MR analyses were primarily null, implying that most of these associations are likely to be a consequence of prodromal disease or selection bias, rather than the risk factor causing the disease.

International electronic health record-derived post-acute sequelae profiles of COVID-19 patients

The risk profiles of post-acute sequelae of COVID-19 (PASC) have not been well characterized in multi-national settings with appropriate controls. We leveraged electronic health record (EHR) data from 277 international hospitals representing 414,602 patients with COVID-19, 2.3 million control patients without COVID-19 in the inpatient and outpatient settings, and over 221 million diagnosis codes to systematically identify new-onset conditions enriched among patients with COVID-19 during the post-acute period. Compared to inpatient controls, inpatient COVID-19 cases were at significant risk for angina pectoris (RR 1.30, 95% CI 1.09–1.55), heart failure (RR 1.22, 95% CI 1.10–1.35), cognitive dysfunctions (RR 1.18, 95% CI 1.07–1.31), and fatigue (RR 1.18, 95% CI 1.07–1.30). Relative to outpatient controls, outpatient COVID-19 cases were at risk for pulmonary embolism (RR 2.10, 95% CI 1.58–2.76), venous embolism (RR 1.34, 95% CI 1.17–1.54), atrial fibrillation (RR 1.30, 95% CI 1.13–1.50), type 2 diabetes (RR 1.26, 95% CI 1.16–1.36) and vitamin D deficiency (RR 1.19, 95% CI 1.09–1.30). Outpatient COVID-19 cases were also at risk for loss of smell and taste (RR 2.42, 95% CI 1.90–3.06), inflammatory neuropathy (RR 1.66, 95% CI 1.21–2.27), and cognitive dysfunction (RR 1.18, 95% CI 1.04–1.33). The incidence of post-acute cardiovascular and pulmonary conditions decreased across time among inpatient cases while the incidence of cardiovascular, digestive, and metabolic conditions increased among outpatient cases. Our study, based on a federated international network, systematically identified robust conditions associated with PASC compared to control groups, underscoring the multifaceted cardiovascular and neurological phenotype profiles of PASC.

Are phenome-wide association studies feasible in a developing country?

Phenome-wide association studies (PheWASs), a powerful approach that examines phenotypes associated with a genetic marker, have been used extensively in highly developed countries. Although there may be a clear need for PheWAS in a developing country such as the Philippines, limitations related to resources and practicality would make conducting them a challenge.

A Process Mining Pipeline to Characterize COVID-19 Patients' Trajectories and Identify Relevant Temporal Phenotypes From EHR Data

The impact of the COVID-19 pandemic involved the disruption of the processes of care and the need for immediately effective re-organizational procedures. In the context of digital health, it is of paramount importance to determine how a specific patients' population reflects into the healthcare dynamics of the hospital, to investigate how patients' sub-group/strata respond to the different care processes, in order to generate novel hypotheses regarding the most effective healthcare strategies. We present an analysis pipeline based on the heterogeneous collected data aimed at identifying the most frequent healthcare processes patterns, jointly analyzing them with demographic and physiological disease trajectories, and stratify the observed cohort on the basis of the mined patterns. This is a process-oriented pipeline which integrates process mining algorithms, and trajectory mining by topological data analyses and pseudo time approaches. Data was collected for 1,179 COVID-19 positive patients, hospitalized at the Italian Hospital "Istituti Clinici Salvatore Maugeri" in Lombardy, integrating different sources including text admission letters, EHR and hospital infrastructure data. We identified five temporal phenotypes, from laboratory values trajectories, which are characterized by statistically significant different death risk estimates. The process mining algorithms allowed splitting the data in sub-cohorts as function of the pandemic waves and of the temporal trajectories showing statistically significant differences in terms of events characteristics.

Cox regression is robust to inaccurate EHR-extracted event time: an application to EHR-based GWAS

MOTIVATION

Logistic regression models are used in genomic studies to analyze the genetic data linked to electronic health records (EHRs), and do not take full usage of the time-to-event information available in EHRs. Previous work has shown that Cox regression, which can account for left truncation and right censoring in EHRs, increased the power to detect genotype-phenotype associations compared to logistic regression. We extend this to evaluate the relative performance of Cox regression and various logistic regression models in the presence of positive errors in event time (delayed event time), relating to recorded event time accuracy.

RESULTS

One Cox model and three logistic regression models were considered under different scenarios of delayed event time. Extensive simulations and a genomic study application were used to evaluate the impact of delayed event time. While logistic regression does not model the time-to-event directly, various logistic regression models used in the literature were more sensitive to delayed event time than Cox regression. Results highlighted the importance to identify and exclude the patients diagnosed before entry time. Cox regression had similar or modest improvement in statistical power over various logistic regression models at controlled type I error. This was supported by the empirical data, where the Cox models steadily had the highest sensitivity to detect known genotype-phenotype associations under all scenarios of delayed event time.

AVAILABILITY AND IMPLEMENTATION

Access to individual-level EHR and genotype data is restricted by the IRB. Simulation code and R script for data process are at: https://github.com/QingxiaCindyChen/CoxRobustEHR.git.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Large-scale real-world data analysis identifies comorbidity patterns in schizophrenia

Schizophrenia affects >3.2 million people in the USA. However, its comorbidity patterns have not been systematically characterized in real-world populations. To address this gap, we conducted an observational study using a cohort of 86 million patients in a nationwide health insurance dataset. We identified participants with schizophrenia and those without schizophrenia matched by age, sex, and the first three digits of zip code. For each phenotype encoded in phecodes, we compared their prevalence in schizophrenia patients and the matched non-schizophrenic participants, and we performed subgroup analyses stratified by age and sex. Results show that anxiety, posttraumatic stress disorder, and substance abuse commonly occur in adolescents and young adults prior to schizophrenia diagnoses. Patients aged 60 and above are at higher risks of developing delirium, alcoholism, dementia, pelvic fracture, and osteomyelitis than their matched controls. Type 2 diabetes, sleep apnea, and eating disorders were more prevalent in women prior to schizophrenia diagnosis, whereas acute renal failure, rhabdomyolysis, and developmental delays were found at higher rates in men. Anxiety and obesity are more commonly seen in patients with schizoaffective disorders compared to patients with other types of schizophrenia. Leveraging a large-scale insurance claims dataset, this study identified less-known comorbidity patterns of schizophrenia and confirmed known ones. These comorbidity profiles can guide clinicians and researchers to take heed of early signs of co-occurring diseases.

Pharmacogenomics: A Step forward Precision Medicine in Childhood Asthma

Personalized medicine, an approach to care in which individual characteristics are used for targeting interventions and maximizing health outcomes, is rapidly becoming a reality for many diseases. Childhood asthma is a heterogeneous disease and many children have uncontrolled symptoms. Therefore, an individualized approach is needed for improving asthma outcomes in children. The rapidly evolving fields of genomics and pharmacogenomics may provide a way to achieve asthma control and reduce future risks in children with asthma. In particular, pharmacogenomics can provide tools for identifying novel molecular mechanisms and biomarkers to guide treatment. Emergent high-throughput technologies, along with patient pheno-endotypization, will increase our knowledge of several molecular mechanisms involved in asthma pathophysiology and contribute to selecting and stratifying appropriate treatment for each patient.

Association of Preinjury Medical Diagnoses With Pediatric Persistent Postconcussion Symptoms in Electronic Health Records

OBJECTIVE

To identify risk factors and generate hypotheses for pediatric persistent postconcussion symptoms (PPCS).

SETTING

A regional healthcare system in the Southeastern United States.

PARTICIPANTS

An electronic health record-based algorithm was developed and validated to identify PPCS cases and controls from an institutional database of more than 2.8 million patients. PPCS cases (n = 274) were patients aged 5 to 18 years with PPCS-related diagnostic codes or with PPCS key words identified by natural language processing of clinical notes. Age, sex, and year of index event-matched controls (n = 1096) were patients with mild traumatic brain injury codes only. Patients with moderate or severe traumatic brain injury were excluded. All patients used our healthcare system at least 3 times 180 days before their injury.

DESIGN

Case-control study.

MAIN MEASURES

The outcome was algorithmic classification of PPCS. Exposures were all preinjury medical diagnoses assigned at least 180 days before the injury.

RESULTS

Cases and controls both had a mean of more than 9 years of healthcare system use preinjury. Of 221 preinjury medical diagnoses, headache disorder was associated with PPCS after accounting for multiple testing (odds ratio [OR] = 2.9; 95% confidence interval [CI]: 1.6-5.0; P = 2.1e-4). Six diagnoses were associated with PPCS at a suggestive threshold for statistical significance (false discovery rate P < .10): gastritis/duodenitis (OR = 2.8; 95% CI: 1.6-5.1; P = 5.0e-4), sleep disorders (OR = 2.3; 95% CI: 1.4-3.7; P = 7.4e-4), abdominal pain (OR = 1.6; 95% CI: 1.2-2.2; P = 9.2e-4), chronic sinusitis (OR = 2.8; 95% CI: 1.5-5.2; P = 1.3e-3), congenital anomalies of the skin (OR = 2.9; 95% CI: 1.5-5.5; P = 1.9e-3), and chronic pharyngitis/nasopharyngitis (OR = 2.4; 95% CI: 1.4-4.3; P = 2.5e-3).

CONCLUSIONS

These results support the strong association of preinjury headache disorders with PPCS. An association of PPCS with prior gastritis/duodenitis, sinusitis, and pharyngitis/nasopharyngitis suggests a role for chronic inflammation in PPCS pathophysiology and risk, although results could equally be attributable to a higher likelihood of somatization among PPCS cases. Identified risk factors should be investigated further and potentially considered during the management of pediatric mild traumatic brain injury cases.

Leveraging deep phenotyping from health check-up cohort with 10,000 Korean individuals for phenome-wide association study of 136 traits

The expanding use of the phenome-wide association study (PheWAS) faces challenges in the context of using International Classification of Diseases billing codes for phenotype definition, imbalanced study population ethnicity, and constrained application of the results in research. We performed a PheWAS utilizing 136 deep phenotypes corroborated by comprehensive health check-ups in a Korean population, along with trans-ethnic comparisons through using the UK Biobank and Biobank Japan Project. Meta-analysis with Korean and Japanese population was done. The PheWAS associated 65 phenotypes with 14,101 significant variants (P < 4.92 × 10-10). Network analysis, visualization of cross-phenotype mapping, and causal inference mapping with Mendelian randomization were conducted. Among phenotype pairs from the genotype-driven cross-phenotype associations, we evaluated penetrance in correlation analysis using a clinical database. We focused on the application of PheWAS in order to make it robust and to aid the derivation of biological meaning post-PheWAS. This comprehensive analysis of PheWAS results based on a health check-up database will provide researchers and clinicians with a panoramic overview of the networks among multiple phenotypes and genetic variants, laying groundwork for the practical application of precision medicine.