The use of phenome-wide association studies (PheWAS) for exploration of novel genotype-phenotype relationships and pleiotropy discovery

Phenomics in Livestock Research: Bottlenecks and Promises of Digital Phenotyping and Other Quantification Techniques on a Global Scale

Bottlenecks in moving genomics to real-life applications also include phenomics. This is true not only for genomics medicine and public health genomics but also in ecology and livestock phenomics. This expert narrative review explores the intricate relationship between genetic makeup and observable phenotypic traits across various biological levels in the context of livestock research. We unpack and emphasize the significance of precise phenotypic data in selective breeding outcomes and examine the multifaceted applications of phenomics, ranging from improvement to assessing welfare, reproductive traits, and environmental adaptation in livestock. As phenotypic traits exhibit strong correlations, their measurement alongside specific biological outcomes provides insights into performance, overall health, and clinical endpoints like morbidity and disease. In addition, automated assessment of livestock holds potential for monitoring the dynamic phenotypic traits across various species, facilitating a deeper comprehension of how they adapt to their environment and attendant stressors. A key challenge in genetic improvement in livestock is predicting individuals with optimal fitness without direct measurement. Temporal predictions from unmanned aerial systems can surpass genomic predictions, offering in-depth data on livestock. In the near future, digital phenotyping and digital biomarkers may further unravel the genetic intricacies of stress tolerance, adaptation and welfare aspects of animals enabling the selection of climate-resilient and productive livestock. This expert review thus delves into challenges associated with phenotyping and discusses technological advancements shaping the future of biological research concerning livestock.

Dissecting the shared genetic landscape of anxiety, depression, and schizophrenia

BACKGROUND

Numerous studies highlight the genetic underpinnings of mental disorders comorbidity, particularly in anxiety, depression, and schizophrenia. However, their shared genetic loci are not well understood. Our study employs Mendelian randomization (MR) and colocalization analyses, alongside multi-omics data, to uncover potential genetic targets for these conditions, thereby informing therapeutic and drug development strategies.

METHODS

We utilized the Consortium for Linkage Disequilibrium Score Regression (LDSC) and Mendelian Randomization (MR) analysis to investigate genetic correlations among anxiety, depression, and schizophrenia. Utilizing GTEx V8 eQTL and deCODE Genetics pQTL data, we performed a three-step summary-data-based Mendelian randomization (SMR) and protein-protein interaction analysis. This helped assess causal and comorbid loci for these disorders and determine if identified loci share coincidental variations with psychiatric diseases. Additionally, phenome-wide association studies, drug prediction, and molecular docking validated potential drug targets.

RESULTS

We found genetic correlations between anxiety, depression, and schizophrenia, and under a meta-analysis of MR from multiple databases, the causal relationships among these disorders are supported. Based on this, three-step SMR and colocalization analyses identified ITIH3 and CCS as being related to the risk of developing depression, while CTSS and DNPH1 are related to the onset of schizophrenia. BTN3A1, PSMB4, and TIMP4 were identified as comorbidity loci for both disorders. Molecules that could not be determined through colocalization analysis were also presented. Drug prediction and molecular docking showed that some drugs and proteins have good binding affinity and available structural data.

CONCLUSIONS

Our study indicates genetic correlations and shared risk loci between anxiety, depression, and schizophrenia. These findings offer insights into the underlying mechanisms of their comorbidities and aid in drug development.

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

IMPORTANCE

Knowledge gained from cohort studies has dramatically advanced both public and precision health. The All of Us Research Program seeks to enroll 1 million diverse participants who share multiple sources of data, providing unique opportunities for research. It is important to understand the phenomic profiles of its participants to conduct research in this cohort.

OBJECTIVES

More than 280 000 participants have shared their electronic health records (EHRs) in the All of Us Research Program. We aim to understand the phenomic profiles of this cohort through comparisons with those in the US general population and a well-established nation-wide cohort, UK Biobank, and to test whether association results of selected commonly studied diseases in the All of Us cohort were comparable to those in UK Biobank.

MATERIALS AND METHODS

We included participants with EHRs in All of Us and participants with health records from UK Biobank. The estimates of prevalence of diseases in the US general population were obtained from the Global Burden of Diseases (GBD) study. We conducted phenome-wide association studies (PheWAS) of 9 commonly studied diseases in both cohorts.

RESULTS

This study included 287 012 participants from the All of Us EHR cohort and 502 477 participants from the UK Biobank. A total of 314 diseases curated by the GBD were evaluated in All of Us, 80.9% (N = 254) of which were more common in All of Us than in the US general population [prevalence ratio (PR) >1.1, P < 2 × 10-5]. Among 2515 diseases and phenotypes evaluated in both All of Us and UK Biobank, 85.6% (N = 2152) were more common in All of Us (PR >1.1, P < 2 × 10-5). The Pearson correlation coefficients of effect sizes from PheWAS between All of Us and UK Biobank were 0.61, 0.50, 0.60, 0.57, 0.40, 0.53, 0.46, 0.47, and 0.24 for ischemic heart diseases, lung cancer, chronic obstructive pulmonary disease, dementia, colorectal cancer, lower back pain, multiple sclerosis, lupus, and cystic fibrosis, respectively.

DISCUSSION

Despite the differences in prevalence of diseases in All of Us compared to the US general population or the UK Biobank, our study supports that All of Us can facilitate rapid investigation of a broad range of diseases.

CONCLUSION

Most diseases were more common in All of Us than in the general US population or the UK Biobank. Results of disease-disease association tests from All of Us are comparable to those estimated in another well-studied national cohort.

Associations between attention-deficit hyperactivity disorder genetic liability and ICD-10 medical conditions in adults: utilizing electronic health records in a Phenome-Wide Association Study

BACKGROUND

Attention-deficit hyperactivity disorder (ADHD) is often comorbid with other medical conditions in adult patients. However, ADHD is extremely underdiagnosed in adults and little is known about the medical comorbidities in undiagnosed adult individuals with high ADHD liability. In this study we investigated associations between ADHD genetic liability and electronic health record (EHR)-based ICD-10 diagnoses across all diagnostic categories, in individuals without ADHD diagnosis history.

METHODS

We used data from the Estonian Biobank cohort (N = 111 261) and generated polygenic risk scores (PRS) for ADHD (PRSADHD) based on the ADHD genome-wide association study. We performed a phenome-wide association study (PheWAS) to test for associations between standardized PRSADHD and 1515 EHR-based ICD-10 diagnoses in the full and sex-stratified sample. We compared the observed significant ICD-10 associations to associations with (1) ADHD diagnosis and (2) questionnaire-based high ADHD risk analyses.

RESULTS

After Bonferroni correction (p = 3.3 × 10-5) we identified 80 medical conditions associated with PRSADHD. The strongest evidence was seen with chronic obstructive pulmonary disease (OR 1.15, CI 1.11-1.18), obesity (OR 1.13, CI 1.11-1.15), and type 2 diabetes (OR 1.11, CI 1.09-1.14). Sex-stratified analysis generally showed similar associations in males and females. Out of all identified associations, 40% and 78% were also observed using ADHD diagnosis or questionnaire-based ADHD, respectively, as the predictor.

CONCLUSIONS

Overall our findings indicate that ADHD genetic liability is associated with an increased risk of a substantial number of medical conditions in undiagnosed individuals. These results highlight the need for timely detection and improved management of ADHD symptoms in adults.

Phenome-wide association study of ovarian cancer identifies common comorbidities and reveals shared genetics with complex diseases and biomarkers

BACKGROUND

Ovarian cancer (OC) is commonly diagnosed among older women who have comorbidities. This hypothesis-free phenome-wide association study (PheWAS) aimed to identify comorbidities associated with OC, as well as traits that share a genetic architecture with OC.

METHODS

We used data from 181,203 white British female UK Biobank participants and analysed OC and OC subtype-specific genetic risk scores (OC-GRS) for an association with 889 diseases and 43 other traits. We conducted PheWAS and colocalization analyses for individual variants to identify evidence for shared genetic architecture.

RESULTS

The OC-GRS was associated with 10 diseases, and the clear cell OC-GRS was associated with five diseases at the FDR threshold (p = 5.6 × 10-4 ). Mendelian randomizaiton analysis (MR) provided robust evidence for the association of OC with higher risk of "secondary malignant neoplasm of digestive systems" (OR 1.64, 95% CI 1.33, 2.02), "ascites" (1.48, 95% CI 1.17, 1.86), "chronic airway obstruction" (1.17, 95% CI 1.07, 1.29), and "abnormal findings on examination of the lung" (1.51, 95% CI 1.22, 1.87). Analyses of lung spirometry measures provided further support for compromised respiratory function. PheWAS on individual OC variants identified five genetic variants associated with other diseases, and seven variants associated with biomarkers (all, p ≤ 4.5 × 10-8 ). Colocalization analysis identified rs4449583 (from TERT locus) as the shared causal variant for OC and seborrheic keratosis.

CONCLUSIONS

OC is associated with digestive and respiratory comorbidities. Several variants affecting OC risk were associated with other diseases and biomarkers, with this study identifying a novel genetic locus shared between OC and skin conditions.

Gender-dependent multiple cross-phenotype association of interferon lambda genetic variants with peripheral blood profiles in healthy individuals

BACKGROUND

Type III interferons (IFN), also called as lambda IFNs (IFN-λs), are antiviral and immunomodulatory cytokines that are evolutionarily important in humans. Given their central roles in innate immunity, they could be influencing other aspects of human biology. This study aimed to examine the association of genetic variants that control the expression and/or activity of IFN-λ3 and IFN-λ4 with multiple phenotypes in blood profiles of healthy individuals.

METHODS

In a cohort of about 550 self-declared healthy individuals, after applying several exclusion criteria to determine their health status, we measured 30 blood parameters, including cellular, biochemical, and metabolic profiles. We genotyped them at rs12979860 and rs28416813 using competitive allele-specific PCR assays and tested their association with the blood profiles under dominant and recessive models for the minor allele. IFN-λ4 variants rs368234815 and rs117648444 were also genotyped or inferred.

RESULTS

We saw no association in the combined cohort under either of the models for any of the phenotypes. When we stratified the cohort based on gender, we saw a significant association only in males with monocyte (p = 1 × 10-3 ) and SGOT (p = 7 × 10-3 ) levels under the dominant model and with uric acid levels (p = 0.01) under the recessive model. When we tested the IFN-λ4 activity modifying variant within groupings based on absence or presence of one or two copies of IFN-λ4 and on different activity levels of IFN-λ4, we found significant (p < 0.05) association with several phenotypes like monocyte, triglyceride, VLDL, ALP, and uric acid levels, only in males. All the above significant associations did not show any confounding when we tested for the same with up to ten different demographic and lifestyle variables.

CONCLUSIONS

These results show that lambda interferons can have pleiotropic effects. However, gender seems to be an effect modifier, with males being more sensitive than females to the effect.

PheWAS and cross-disorder analysis reveal genetic architecture, pleiotropic loci and phenotypic correlations across 11 autoimmune disorders

Introduction

Autoimmune disorders (ADs) are a group of about 80 disorders that occur when self-attacking autoantibodies are produced due to failure in the self-tolerance mechanisms. ADs are polygenic disorders and associations with genes both in the human leukocyte antigen (HLA) region and outside of it have been described. Previous studies have shown that they are highly comorbid with shared genetic risk factors, while epidemiological studies revealed associations between various lifestyle and health-related phenotypes and ADs.

Methods

Here, for the first time, we performed a comparative polygenic risk score (PRS) - Phenome Wide Association Study (PheWAS) for 11 different ADs (Juvenile Idiopathic Arthritis, Primary Sclerosing Cholangitis, Celiac Disease, Multiple Sclerosis, Rheumatoid Arthritis, Psoriasis, Myasthenia Gravis, Type 1 Diabetes, Systemic Lupus Erythematosus, Vitiligo Late Onset, Vitiligo Early Onset) and 3,254 phenotypes available in the UK Biobank that include a wide range of socio-demographic, lifestyle and health-related outcomes. Additionally, we investigated the genetic relationships of the studied ADs, calculating their genetic correlation and conducting cross-disorder GWAS meta-analyses for the observed AD clusters.

Results

In total, we identified 508 phenotypes significantly associated with at least one AD PRS. 272 phenotypes were significantly associated after excluding variants in the HLA region from the PRS estimation. Through genetic correlation and genetic factor analyses, we identified four genetic factors that run across studied ADs. Cross-trait meta-analyses within each factor revealed pleiotropic genome-wide significant loci.

Discussion

Overall, our study confirms the association of different factors with genetic susceptibility for ADs and reveals novel observations that need to be further explored.

Using Drosophila to identify naturally occurring genetic modifiers of amyloid beta 42- and tau-induced toxicity

Alzheimer's disease is characterized by 2 pathological proteins, amyloid beta 42 and tau. The majority of Alzheimer's disease cases in the population are sporadic and late-onset Alzheimer's disease, which exhibits high levels of heritability. While several genetic risk factors for late-onset Alzheimer's disease have been identified and replicated in independent studies, including the ApoE ε4 allele, the great majority of the heritability of late-onset Alzheimer's disease remains unexplained, likely due to the aggregate effects of a very large number of genes with small effect size, as well as to biases in sample collection and statistical approaches. Here, we present an unbiased forward genetic screen in Drosophila looking for naturally occurring modifiers of amyloid beta 42- and tau-induced ommatidial degeneration. Our results identify 14 significant SNPs, which map to 12 potential genes in 8 unique genomic regions. Our hits that are significant after genome-wide correction identify genes involved in neuronal development, signal transduction, and organismal development. Looking more broadly at suggestive hits (P < 10-5), we see significant enrichment in genes associated with neurogenesis, development, and growth as well as significant enrichment in genes whose orthologs have been identified as significantly or suggestively associated with Alzheimer's disease in human GWAS studies. These latter genes include ones whose orthologs are in close proximity to regions in the human genome that are associated with Alzheimer's disease, but where a causal gene has not been identified. Together, our results illustrate the potential for complementary and convergent evidence provided through multitrait GWAS in Drosophila to supplement and inform human studies, helping to identify the remaining heritability and novel modifiers of complex diseases.

Phenomics as an approach to Comparative Developmental Physiology

The dynamic nature of developing organisms and how they function presents both opportunity and challenge to researchers, with significant advances in understanding possible by adopting innovative approaches to their empirical study. The information content of the phenotype during organismal development is arguably greater than at any other life stage, incorporating change at a broad range of temporal, spatial and functional scales and is of broad relevance to a plethora of research questions. Yet, effectively measuring organismal development, and the ontogeny of physiological regulations and functions, and their responses to the environment, remains a significant challenge. "Phenomics", a global approach to the acquisition of phenotypic data at the scale of the whole organism, is uniquely suited as an approach. In this perspective, we explore the synergies between phenomics and Comparative Developmental Physiology (CDP), a discipline of increasing relevance to understanding sensitivity to drivers of global change. We then identify how organismal development itself provides an excellent model for pushing the boundaries of phenomics, given its inherent complexity, comparably smaller size, relative to adult stages, and the applicability of embryonic development to a broad suite of research questions using a diversity of species. Collection, analysis and interpretation of whole organismal phenotypic data are the largest obstacle to capitalising on phenomics for advancing our understanding of biological systems. We suggest that phenomics within the context of developing organismal form and function could provide an effective scaffold for addressing grand challenges in CDP and phenomics.

International Consensus Statement on Obstructive Sleep Apnea

BACKGROUND

Evaluation and interpretation of the literature on obstructive sleep apnea (OSA) allows for consolidation and determination of the key factors important for clinical management of the adult OSA patient. Toward this goal, an international collaborative of multidisciplinary experts in sleep apnea evaluation and treatment have produced the International Consensus statement on Obstructive Sleep Apnea (ICS:OSA).

METHODS

Using previously defined methodology, focal topics in OSA were assigned as literature review (LR), evidence-based review (EBR), or evidence-based review with recommendations (EBR-R) formats. Each topic incorporated the available and relevant evidence which was summarized and graded on study quality. Each topic and section underwent iterative review and the ICS:OSA was created and reviewed by all authors for consensus.

RESULTS

The ICS:OSA addresses OSA syndrome definitions, pathophysiology, epidemiology, risk factors for disease, screening methods, diagnostic testing types, multiple treatment modalities, and effects of OSA treatment on multiple OSA-associated comorbidities. Specific focus on outcomes with positive airway pressure (PAP) and surgical treatments were evaluated.

CONCLUSION

This review of the literature consolidates the available knowledge and identifies the limitations of the current evidence on OSA. This effort aims to create a resource for OSA evidence-based practice and identify future research needs. Knowledge gaps and research opportunities include improving the metrics of OSA disease, determining the optimal OSA screening paradigms, developing strategies for PAP adherence and longitudinal care, enhancing selection of PAP alternatives and surgery, understanding health risk outcomes, and translating evidence into individualized approaches to therapy.

Future prospects for human genetics and genomics in drug discovery

Evidence from human genetics supporting the therapeutic hypothesis increases the likelihood that a drug will succeed in clinical trials. Rare and common disease genetics yield a wide array of alleles with a range of effect sizes that can proxy for the effect of a drug in disease. Recent advances in large scale population collections and whole genome sequencing approaches have provided a rich resource of human genetic evidence to support drug target selection. As the range of phenotypes profiled increases and ever more alleles are discovered across world-wide populations, these approaches will increasingly influence multiple stages across the lifespan of a drug discovery programme.

Assessment of multi-population polygenic risk scores for lipid traits in African Americans

Polygenic risk scores (PRS) based on genome-wide discoveries are promising predictors or classifiers of disease development, severity, and/or progression for common clinical outcomes. A major limitation of most risk scores is the paucity of genome-wide discoveries in diverse populations, prompting an emphasis to generate these needed data for trans-population and population-specific PRS construction. Given diverse genome-wide discoveries are just now being completed, there has been little opportunity for PRS to be evaluated in diverse populations independent from the discovery efforts. To fill this gap, we leverage here summary data from a recent genome-wide discovery study of lipid traits (HDL-C, LDL-C, triglycerides, and total cholesterol) conducted in diverse populations represented by African Americans, Hispanics, Asians, Native Hawaiians, Native Americans, and others by the Population Architecture using Genomics and Epidemiology (PAGE) Study. We constructed lipid trait PRS using PAGE Study published genetic variants and weights in an independent African American adult patient population linked to de-identified electronic health records and genotypes from the Illumina Metabochip (n = 3,254). Using multi-population lipid trait PRS, we assessed levels of association for their respective lipid traits, clinical outcomes (cardiovascular disease and type 2 diabetes), and common clinical labs. While none of the multi-population PRS were strongly associated with the tested trait or outcome, PRS_LDL-Cwas nominally associated with cardiovascular disease. These data demonstrate the complexity in applying PRS to real-world clinical data even when data from multiple populations are available.

Mapping anorexia nervosa genes to clinical phenotypes

BACKGROUND

Anorexia nervosa (AN) is a psychiatric disorder with complex etiology, with a significant portion of disease risk imparted by genetics. Traditional genome-wide association studies (GWAS) produce principal evidence for the association of genetic variants with disease. Transcriptomic imputation (TI) allows for the translation of those variants into regulatory mechanisms, which can then be used to assess the functional outcome of genetically regulated gene expression (GReX) in a broader setting through the use of phenome-wide association studies (pheWASs) in large and diverse clinical biobank populations with electronic health record phenotypes.

METHODS

Here, we applied TI using S-PrediXcan to translate the most recent PGC-ED AN GWAS findings into AN-GReX. For significant genes, we imputed AN-GReX in the Mount Sinai BioMe™ Biobank and performed pheWASs on over 2000 outcomes to test the clinical consequences of aberrant expression of these genes. We performed a secondary analysis to assess the impact of body mass index (BMI) and sex on AN-GReX clinical associations.

RESULTS

Our S-PrediXcan analysis identified 53 genes associated with AN, including what is, to our knowledge, the first-genetic association of AN with the major histocompatibility complex. AN-GReX was associated with autoimmune, metabolic, and gastrointestinal diagnoses in our biobank cohort, as well as measures of cholesterol, medications, substance use, and pain. Additionally, our analyses showed moderation of AN-GReX associations with measures of cholesterol and substance use by BMI, and moderation of AN-GReX associations with celiac disease by sex.

CONCLUSIONS

Our BMI-stratified results provide potential avenues of functional mechanism for AN-genes to investigate further.

Phenome-wide Association Analysis of Substance Use Disorders in a Deeply Phenotyped Sample

BACKGROUND

Substance use disorders (SUDs) are associated with a variety of co-occurring psychiatric disorders and other SUDs, which partly reflects genetic pleiotropy. Polygenic risk scores (PRSs) and phenome-wide association studies are useful in evaluating pleiotropic effects. However, the comparatively low prevalence of SUDs in population samples and the lack of detailed information available in electronic health records limit these data sets' informativeness for such analyses.

METHODS

We used the deeply phenotyped Yale-Penn sample (n = 10,610 with genetic data; 46.3% African ancestry, 53.7% European ancestry) to examine pleiotropy for 4 major substance-related traits: alcohol use disorder, opioid use disorder, smoking initiation, and lifetime cannabis use. The sample includes both affected and control subjects interviewed using the Semi-Structured Assessment for Drug Dependence and Alcoholism, a comprehensive psychiatric interview.

RESULTS

In African ancestry individuals, PRS for alcohol use disorder, and in European individuals, PRS for alcohol use disorder, opioid use disorder, and smoking initiation were associated with their respective primary DSM diagnoses. These PRSs were also associated with additional phenotypes involving the same substance. Phenome-wide association study analyses of PRS in European individuals identified associations across multiple phenotypic domains, including phenotypes not commonly assessed in phenome-wide association study analyses, such as family environment and early childhood experiences.

CONCLUSIONS

Smaller, deeply phenotyped samples can complement large biobank genetic studies with limited phenotyping by providing greater phenotypic granularity. These efforts allow associations to be identified between specific features of disorders and genetic liability for SUDs, which help to inform our understanding of the pleiotropic pathways underlying them.

Genetic correlates of vitamin D-binding protein and 25-hydroxyvitamin D in neonatal dried blood spots

The vitamin D binding protein (DBP), encoded by the group-specific component (GC) gene, is a component of the vitamin D system. In a genome-wide association study of DBP concentration in 65,589 neonates we identify 26 independent loci, 17 of which are in or close to the GC gene, with fine-mapping identifying 2 missense variants on chromosomes 12 and 17 (within SH2B3 and GSDMA, respectively). When adjusted for GC haplotypes, we find 15 independent loci distributed over 10 chromosomes. Mendelian randomization analyses identify a unidirectional effect of higher DBP concentration and (a) higher 25-hydroxyvitamin D concentration, and (b) a reduced risk of multiple sclerosis and rheumatoid arthritis. A phenome-wide association study confirms that higher DBP concentration is associated with a reduced risk of vitamin D deficiency. Our findings provide valuable insights into the influence of DBP on vitamin D status and a range of health outcomes.

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.

Comorbidities in lichen planus by phenome-wide association study in two biobank population cohorts

BACKGROUND

Lichen planus (LP) is a relatively frequent mucocutaneous inflammatory disease affecting the skin, skin appendages and mucosae, including oral mucosae, and less frequently the anogenital area, conjunctivae, oesophagus or larynx.

OBJECTIVES

To estimate the association of LP, with emphasis on dermatological and gastrointestinal conditions, in two large independent population cohorts.

MATERIALS AND METHODS

We performed a phenome-wide association study (PheWAS) and examined conditions associated with LP in two unrelated cohorts, i.e. the multicentre, community-based UK Biobank (UKB: 501 381 controls; 1130 LP subjects) and the healthcare-associated Penn Medicine BioBank (PMBB; 42 702 controls; 764 LP subjects). The data were analysed in 2021. The 'PheWAS' R package was used to perform the PheWAS analyses and Bonferroni correction was used to adjust for multiple testing. Odds ratios (ORs) were adjusted for age, sex and body mass index.

RESULTS

In the UKB, PheWAS revealed 133 phenome codes (PheCodes) significantly associated with LP and most of them were confirmed in PMBB. Dermatological and digestive PheCodes were the most abundant: 29 and 34 of these disorders, respectively, were significantly overrepresented in LP individuals from both cohorts. The 29 dermatological and 12 oral disorders were often highly enriched, whereas hepatic, gastric, oesophageal and intestinal PheCodes displayed ORs in the range of 1·6-4·5. Several autoimmune disorders also exhibited OR > 5 in both cohorts.

CONCLUSIONS

PheWAS in two large unrelated cohorts identified previously unknown comorbidities and may support clinical counselling of patients with LP. What is already known about this topic? Lichen planus (LP) is known to affect the skin, skin appendages and mucosae, including oral mucosae, and less frequently the anogenital area, conjunctivae, oesophagus or larynx. What does this study add? Our data provide the most comprehensive collection of associated dermatological, digestive and autoimmune disorders to date. Our findings are expected to be useful for the evaluation and management of patients with LP.

Deciphering Pleiotropic Signatures of Regulatory SNPs in Zea mays L. Using Multi-Omics Data and Machine Learning Algorithms

Maize is one of the most widely grown cereals in the world. However, to address the challenges in maize breeding arising from climatic anomalies, there is a need for developing novel strategies to harness the power of multi-omics technologies. In this regard, pleiotropy is an important genetic phenomenon that can be utilized to simultaneously enhance multiple agronomic phenotypes in maize. In addition to pleiotropy, another aspect is the consideration of the regulatory SNPs (rSNPs) that are likely to have causal effects in phenotypic development. By incorporating both aspects in our study, we performed a systematic analysis based on multi-omics data to reveal the novel pleiotropic signatures of rSNPs in a global maize population. For this purpose, we first applied Random Forests and then Markov clustering algorithms to decipher the pleiotropic signatures of rSNPs, based on which hierarchical network models are constructed to elucidate the complex interplay among transcription factors, rSNPs, and phenotypes. The results obtained in our study could help to understand the genetic programs orchestrating multiple phenotypes and thus could provide novel breeding targets for the simultaneous improvement of several agronomic traits.

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.

pyPheWAS: A Phenome-Disease Association Tool for Electronic Medical Record Analysis

Along with the increasing availability of electronic medical record (EMR) data, phenome-wide association studies (PheWAS) and phenome-disease association studies (PheDAS) have become a prominent, first-line method of analysis for uncovering the secrets of EMR. Despite this recent growth, there is a lack of approachable software tools for conducting these analyses on large-scale EMR cohorts. In this article, we introduce pyPheWAS, an open-source python package for conducting PheDAS and related analyses. This toolkit includes 1) data preparation, such as cohort censoring and age-matching; 2) traditional PheDAS analysis of ICD-9 and ICD-10 billing codes; 3) PheDAS analysis applied to a novel EMR phenotype mapping: current procedural terminology (CPT) codes; and 4) novelty analysis of significant disease-phenotype associations found through PheDAS. The pyPheWAS toolkit is approachable and comprehensive, encapsulating data prep through result visualization all within a simple command-line interface. The toolkit is designed for the ever-growing scale of available EMR data, with the ability to analyze cohorts of 100,000 + patients in less than 2 h. Through a case study of Down Syndrome and other intellectual developmental disabilities, we demonstrate the ability of pyPheWAS to discover both known and potentially novel disease-phenotype associations across different experiment designs and disease groups. The software and user documentation are available in open source at https://github.com/MASILab/pyPheWAS .

Maturation and application of phenome-wide association studies

In the past 10 years since its introduction, phenome-wide association studies (PheWAS) have uncovered novel genotype-phenotype relationships. Along the way, PheWAS have evolved in many aspects as a study design with the expanded availability of large data repositories with genome-wide data linked to detailed phenotypic data. Advancement in methods, including algorithms, software, and publicly available integrated resources, makes it feasible to more fully realize the potential of PheWAS, overcoming the previous computational and analytical limitations. We review here the most recent improvements and notable applications of PheWAS since the second half of the decade from its inception. We also note the challenges that remain embedded along the entire PheWAS analytical pipeline that necessitate further development of tools and resources to further advance the understanding of the complex genetic architecture underlying human diseases and traits.

Using Phecodes for Research with the Electronic Health Record: From PheWAS to PheRS

Electronic health records (EHRs) are a rich source of data for researchers, but extracting meaningful information out of this highly complex data source is challenging. Phecodes represent one strategy for defining phenotypes for research using EHR data. They are a high-throughput phenotyping tool based on ICD (International Classification of Diseases) codes that can be used to rapidly define the case/control status of thousands of clinically meaningful diseases and conditions. Phecodes were originally developed to conduct phenome-wide association studies to scan for phenotypic associations with common genetic variants. Since then, phecodes have been used to support a wide range of EHR-based phenotyping methods, including the phenotype risk score. This review aims to comprehensively describe the development, validation, and applications of phecodes and suggest some future directions for phecodes and high-throughput phenotyping.

Salivary metabolites associated with a 5-year tooth loss identified in a population-based setting

BACKGROUND

Periodontitis is among the most common chronic diseases worldwide, and it is one of the main reasons for tooth loss. Comprehensive profiling of the metabolite content of the saliva can enable the identification of novel pathways associated with periodontitis and highlight non-invasive markers to facilitate time and cost-effective screening efforts for the presence of periodontitis and the prediction of tooth loss.

METHODS

We first investigated cross-sectional associations of 13 oral health variables with saliva levels of 562 metabolites, measured by untargeted mass spectrometry among a sub-sample (n = 938) of the Study of Health in Pomerania (SHIP-2) using linear regression models adjusting for common confounders. We took forward any candidate metabolite associated with at least two oral variables, to test for an association with a 5-year tooth loss over and above baseline oral health status using negative binomial regression models.

RESULTS

We identified 84 saliva metabolites that were associated with at least one oral variable cross-sectionally, for a subset of which we observed robust replication in an independent study. Out of 34 metabolites associated with more than two oral variables, baseline saliva levels of nine metabolites were positively associated with a 5-year tooth loss. Across all analyses, the metabolites 2-pyrrolidineacetic acid and butyrylputrescine were the most consistent candidate metabolites, likely reflecting oral dysbiosis. Other candidate metabolites likely reflected tissue destruction and cell proliferation.

CONCLUSIONS

Untargeted metabolic profiling of saliva replicated metabolic signatures of periodontal status and revealed novel metabolites associated with periodontitis and future tooth loss.

Sclerostin Downregulation Globally by Naturally Occurring Genetic Variants, or Locally in Atherosclerotic Plaques, Does Not Associate With Cardiovascular Events in Humans

Inhibition of sclerostin increases bone formation and decreases bone resorption, leading to increased bone mass, bone mineral density, and bone strength and reduced fracture risk. In a clinical study of the sclerostin antibody romosozumab versus alendronate in postmenopausal women (ARCH), an imbalance in adjudicated serious cardiovascular (CV) adverse events driven by an increase in myocardial infarction (MI) and stroke was observed. To explore whether there was a potential mechanistic plausibility that sclerostin expression, or its inhibition, in atherosclerotic (AS) plaques may have contributed to this imbalance, sclerostin was immunostained in human plaques to determine whether it was detected in regions relevant to plaque stability in 94 carotid and 50 femoral AS plaques surgically collected from older female patients (mean age 69.6 ± 10.4 years). Sclerostin staining was absent in most plaques (67%), and when detected, it was of reduced intensity compared with normal aorta and was located in deeper regions of the plaque/wall but was not observed in areas considered relevant to plaque stability (fibrous cap and endothelium). Additionally, genetic variants associated with lifelong reduced sclerostin expression were explored for associations with phenotypes including those related to bone physiology and CV risk factors/events in a population-based phenomewide association study (PheWAS). Natural genetic modulation of sclerostin by variants with a significant positive effect on bone physiology showed no association with lifetime risk of MI or stroke. These data do not support a causal association between the presence of sclerostin, or its inhibition, in the vasculature and increased risk of serious cardiovascular events. © 2021 The Authors. Journal of Bone and Mineral Research published by Wiley Periodicals LLC on behalf of American Society for Bone and Mineral Research (ASBMR).

Neuroimaging PheWAS (Phenome-Wide Association Study): A Free Cloud-Computing Platform for Big-Data, Brain-Wide Imaging Association Studies

Large-scale, case-control genome-wide association studies (GWASs) have revealed genetic variations associated with diverse neurological and psychiatric disorders. Recent advances in neuroimaging and genomic databases of large healthy and diseased cohorts have empowered studies to characterize effects of the discovered genetic factors on brain structure and function, implicating neural pathways and genetic mechanisms in the underlying biology. However, the unprecedented scale and complexity of the imaging and genomic data requires new advanced biomedical data science tools to manage, process and analyze the data. In this work, we introduce Neuroimaging PheWAS (phenome-wide association study): a web-based system for searching over a wide variety of brain-wide imaging phenotypes to discover true system-level gene-brain relationships using a unified genotype-to-phenotype strategy. This design features a user-friendly graphical user interface (GUI) for anonymous data uploading, study definition and management, and interactive result visualizations as well as a cloud-based computational infrastructure and multiple state-of-art methods for statistical association analysis and multiple comparison correction. We demonstrated the potential of Neuroimaging PheWAS with a case study analyzing the influences of the apolipoprotein E (APOE) gene on various brain morphological properties across the brain in the Alzheimer's Disease Neuroimaging Initiative (ADNI) cohort. Benchmark tests were performed to evaluate the system's performance using data from UK Biobank. The Neuroimaging PheWAS system is freely available. It simplifies the execution of PheWAS on neuroimaging data and provides an opportunity for imaging genetics studies to elucidate routes at play for specific genetic variants on diseases in the context of detailed imaging phenotypic data.

Genome-wide association study of circulating interleukin 6 levels identifies novel loci

Interleukin 6 (IL-6) is a multifunctional cytokine with both pro- and anti-inflammatory properties with a heritability estimate of up to 61%. The circulating levels of IL-6 in blood have been associated with an increased risk of complex disease pathogenesis. We conducted a two-staged, discovery and replication meta genome-wide association study (GWAS) of circulating serum IL-6 levels comprising up to 67 428 (ndiscovery = 52 654 and nreplication = 14 774) individuals of European ancestry. The inverse variance fixed effects based discovery meta-analysis, followed by replication led to the identification of two independent loci, IL1F10/IL1RN rs6734238 on chromosome (Chr) 2q14, (Pcombined = 1.8 × 10-11), HLA-DRB1/DRB5 rs660895 on Chr6p21 (Pcombined = 1.5 × 10-10) in the combined meta-analyses of all samples. We also replicated the IL6R rs4537545 locus on Chr1q21 (Pcombined = 1.2 × 10-122). Our study identifies novel loci for circulating IL-6 levels uncovering new immunological and inflammatory pathways that may influence IL-6 pathobiology.

Pleiotropic associations of heterozygosity for the SERPINA1 Z allele in the UK Biobank

Homozygosity for the SERPINA1 Z allele causes α1-antitrypsin deficiency, a rare condition that can cause lung and liver disease. However, the effects of Z allele heterozygosity on nonrespiratory phenotypes, and on lung function in the general population, remain unclear. We conducted a large, population-based study to determine Z allele effects on >2400 phenotypes in the UK Biobank (N=303 353). Z allele heterozygosity was strongly associated with increased height (β=1.02 cm, p=3.91×10-68), and with other nonrespiratory phenotypes including increased risk of gall bladder disease, reduced risk of heart disease and lower blood pressure, reduced risk of osteoarthritis and reduced bone mineral density, increased risk of headache and enlarged prostate, as well as with blood biomarkers of liver function. Heterozygosity was associated with higher height-adjusted forced expiratory volume in 1 s (FEV1) (β=19.36 mL, p=9.21×10-4) and FEV1/forced vital capacity (β=0.0031, p=1.22×10-5) in nonsmokers, whereas in smokers, this protective effect was abolished. Furthermore, we show for the first time that sex modifies the association of the Z allele on lung function. We conclude that Z allele heterozygosity and homozygosity exhibit opposing effects on lung function in the UK population, and that these associations are modified by smoking and sex. In exploratory analyses, heterozygosity for the Z allele also showed pleiotropic associations with nonrespiratory health-related traits and disease risk.

Brain Magnetic Resonance Imaging Phenome-Wide Association Study With Metal Transporter Gene SLC39A8

The SLC39A8 gene encodes a divalent metal transporter, ZIP8. SLC39A8 is associated with pleiotropic effects across multiple tissues, including the brain. We determine the different brain magnetic resonance imaging (MRI) phenotypes associated with SLC39A8. We used a phenome-wide association study approach followed by joint and conditional association analysis. Using the summary statistics datasets from a brain MRI genome-wide association study on adult United Kingdom (UK) Biobank participants, we systematically selected all brain MRI phenotypes associated with single-nucleotide polymorphisms (SNPs) within 500 kb of the SLC39A8 genetic locus. For all significant brain MRI phenotypes, we used GCTA-COJO to determine the number of independent association signals and identify index SNPs for each brain MRI phenotype. Linkage equilibrium for brain phenotypes with multiple independent signals was confirmed by LDpair. We identified 24 brain MRI phenotypes that vary due to MRI type and brain region and contain a SNP associated with the SLC39A8 locus. Missense ZIP8 polymorphism rs13107325 was associated with 22 brain MRI phenotypes. Rare ZIP8 variants present in a published UK Biobank dataset are associated with 6 brain MRI phenotypes also linked to rs13107325. Among the 24 datasets, an additional 4 association signals were identified by GCTA-COJO and confirmed to be in linkage equilibrium with rs13107325 using LDpair. These additional association signals represent new probable causative SNPs in addition to rs13107325. This study provides leads into how genetic variation in SLC39A8, a trace mineral transport gene, is linked to brain structure differences and may affect brain development and nervous system function.

Mendelian randomization for studying the effects of perturbing drug targets

Drugs whose targets have genetic evidence to support efficacy and safety are more likely to be approved after clinical development. In this paper, we provide an overview of how natural sequence variation in the genes that encode drug targets can be used in Mendelian randomization analyses to offer insight into mechanism-based efficacy and adverse effects. Large databases of summary level genetic association data are increasingly available and can be leveraged to identify and validate variants that serve as proxies for drug target perturbation. As with all empirical research, Mendelian randomization has limitations including genetic confounding, its consideration of lifelong effects, and issues related to heterogeneity across different tissues and populations. When appropriately applied, Mendelian randomization provides a useful empirical framework for using population level data to improve the success rates of the drug development pipeline.

Mendelian randomization for studying the effects of perturbing drug targets

Drugs whose targets have genetic evidence to support efficacy and safety are more likely to be approved after clinical development. In this paper, we provide an overview of how natural sequence variation in the genes that encode drug targets can be used in Mendelian randomization analyses to offer insight into mechanism-based efficacy and adverse effects. Large databases of summary level genetic association data are increasingly available and can be leveraged to identify and validate variants that serve as proxies for drug target perturbation. As with all empirical research, Mendelian randomization has limitations including genetic confounding, its consideration of lifelong effects, and issues related to heterogeneity across different tissues and populations. When appropriately applied, Mendelian randomization provides a useful empirical framework for using population level data to improve the success rates of the drug development pipeline.

Genome-Wide Association Studies of CKD and Related Traits

The past few years have seen major advances in genome-wide association studies (GWAS) of CKD and kidney function-related traits in several areas: increases in sample size from >100,000 to >1 million, enabling the discovery of >250 associated genetic loci that are highly reproducible; the inclusion of participants not only of European but also of non-European ancestries; and the use of advanced computational methods to integrate additional genomic and other unbiased, high-dimensional data to characterize the underlying genetic architecture and prioritize potentially causal genes and variants. Together with other large-scale biobank and genetic association studies of complex traits, these GWAS of kidney function-related traits have also provided novel insight into the relationship of kidney function to other diseases with respect to their genetic associations, genetic correlation, and directional relationships. A number of studies also included functional experiments using model organisms or cell lines to validate prioritized potentially causal genes and/or variants. In this review article, we will summarize these recent GWAS of CKD and kidney function-related traits, explain approaches for downstream characterization of associated genetic loci and the value of such computational follow-up analyses, and discuss related challenges along with potential solutions to ultimately enable improved treatment and prevention of kidney diseases through genetics.

Implicit bias of encoded variables: frameworks for addressing structured bias in EHR-GWAS data

The 'discovery' stage of genome-wide association studies required amassing large, homogeneous cohorts. In order to attain clinically useful insights, we must now consider the presentation of disease within our clinics and, by extension, within our medical records. Large-scale use of electronic health record (EHR) data can help to understand phenotypes in a scalable manner, incorporating lifelong and whole-phenome context. However, extending analyses to incorporate EHR and biobank-based analyses will require careful consideration of phenotype definition. Judgements and clinical decisions that occur 'outside' the system inevitably contain some degree of bias and become encoded in EHR data. Any algorithmic approach to phenotypic characterization that assumes non-biased variables will generate compounded biased conclusions. Here, we discuss and illustrate potential biases inherent within EHR analyses, how these may be compounded across time and suggest frameworks for large-scale phenotypic analysis to minimize and uncover encoded bias.

Genome-wide pathogenesis interpretation using a heat diffusion-based systems genetics method and implications for gene function annotation

Background

Genetics is best dedicated to interpreting pathogenesis and revealing gene functions. The past decade has witnessed unprecedented progress in genetics, particularly in genome‐wide identification of disorder variants through Genome‐Wide Association Studies (GWAS) and Phenome‐Wide Association Studies (PheWAS). However, it is still a great challenge to use GWAS/PheWAS‐derived data to elucidate pathogenesis.

Methods

In this study, we used HotNet2, a heat diffusion‐based systems genetics algorithm, to calculate the networks for disease genes obtained from GWAS and PheWAS, with an attempt to get deeper insights into disease pathogenesis at a molecular level.

Results

Through HotNet2 calculation, significant networks for 202 (for GWAS) and 167 (for PheWAS) types of diseases were identified and evaluated, respectively. The GWAS‐derived disease networks exhibit a stronger biomedical relevance than PheWAS counterparts. Therefore, the GWAS‐derived networks were used for pathogenesis interpretation by integrating the accumulated biomedical information. As a result, the pathogenesis for 64 diseases was elucidated in terms of mutation‐caused abnormal transcriptional regulation, and 47 diseases were preliminarily interpreted in terms of mutation‐caused varied protein‐protein interactions. In addition, 3,802 genes (including 46 function‐unknown genes) were assigned with new functions by disease network information, some of which were validated through mice gene knockout experiments.

Conclusions

Systems genetics algorithm HotNet2 can efficiently establish genotype‐phenotype links at the level of biological networks. Compared with original GWAS/PheWAS results, HotNet2‐calculated disease‐gene associations have stronger biomedical significance, hence provide better interpretations for the pathogenesis of genome‐wide variants, and offer new insights into gene functions as well. These results are also helpful in drug development.

Characterization of genetic and phenotypic heterogeneity of obstructive sleep apnea using electronic health records

Background

Obstructive sleep apnea (OSA) is defined by frequent episodes of reduced or complete cessation of airflow during sleep and is linked to negative health outcomes. Understanding the genetic factors influencing expression of OSA may lead to new treatment strategies. Electronic health records (EHRs) can be leveraged to both validate previously reported OSA-associated genomic variation and detect novel relationships between these variants and comorbidities.

Methods

We identified candidate single nucleotide polymorphisms (SNPs) via systematic literature review of existing research. Using datasets available at Geisinger (n = 39,407) and Vanderbilt University Medical Center (n = 24,084), we evaluated associations between 40 previously implicated SNPs and OSA diagnosis, defined using clinical codes. We also evaluated associations between these SNPs and OSA severity measures obtained from sleep reports at Geisinger (n = 6571). Finally, we used a phenome-wide association study approach to help reveal pleiotropic genetic effects between OSA candidate SNPs and other clinical codes and laboratory values available in the EHR.

Results

Most previously reported OSA candidate SNPs showed minimal to no evidence for associations with OSA diagnosis or severity in the EHR-derived datasets. Three SNPs in LEPR, MMP-9, and GABBR1 validated for an association with OSA diagnosis in European Americans; the SNP in GABBR1 was associated following meta-analysis of results from both clinical populations. The GABBR1 and LEPR SNPs, and one additional SNP, were associated with OSA severity measures in European Americans from Geisinger. Three additional candidate OSA SNPs were not associated with OSA-related traits but instead with hyperlipidemia and autoimmune diseases of the thyroid.

Conclusions

To our knowledge, this is one of the largest candidate gene studies and one of the first phenome-wide association studies of OSA genomic variation. Results validate genetic associates with OSA in the LEPR, MMP-9 and GABBR1 genes, but suggest that the majority of previously identified genetic associations with OSA may be false positives. Phenome-wide analyses provide evidence of mediated pleiotropy. Future well-powered genome-wide association analyses of OSA risk and severity across populations with diverse ancestral backgrounds are needed. The comprehensive nature of the analyses represents a platform for informing future work focused on understanding how genetic data can be useful to informing treatment of OSA and related comorbidities.

Exploring Diseases/Traits and Blood Proteins Causally Related to Expression of ACE2, the Putative Receptor of SARS-CoV-2: A Mendelian Randomization Analysis Highlights Tentative Relevance of Diabetes-Related Traits

OBJECTIVE

COVID-19 has become a major public health problem. There is good evidence that ACE2 is a receptor for SARS-CoV-2, and high expression of ACE2 may increase susceptibility to infection. We aimed to explore risk factors affecting susceptibility to infection and prioritize drug repositioning candidates, based on Mendelian randomization (MR) studies on ACE2 lung expression.

RESEARCH DESIGN AND METHODS

We conducted a phenome-wide MR study to prioritize diseases/traits and blood proteins causally linked to ACE2 lung expression in GTEx. We also explored drug candidates whose targets overlapped with the top-ranked proteins in MR, as these drugs may alter ACE2 expression and may be clinically relevant.

RESULTS

The most consistent finding was tentative evidence of an association between diabetes-related traits and increased ACE2 expression. Based on one of the largest genome-wide association studies on type 2 diabetes mellitus (T2DM) to date (N = 898,130), T2DM was causally linked to raised ACE2 expression (P = 2.91E-03; MR-IVW). Significant associations (at nominal level; P < 0.05) with ACE2 expression were observed across multiple diabetes data sets and analytic methods for T1DM, T2DM, and related traits including early start of insulin. Other diseases/traits having nominal significant associations with increased expression included inflammatory bowel disease, (estrogen receptor-positive) breast cancer, lung cancer, asthma, smoking, and elevated alanine aminotransferase. We also identified drugs that may target the top-ranked proteins in MR, such as fostamatinib and zinc.

CONCLUSIONS

Our analysis suggested that diabetes and related traits may increase ACE2 expression, which may influence susceptibility to infection (or more severe infection). However, none of these findings withstood rigorous multiple testing corrections (at false discovery rate <0.05). Proteome-wide MR analyses might help uncover mechanisms underlying ACE2 expression and guide drug repositioning. Further studies are required to verify our findings.

Developing a FHIR-based Framework for Phenome Wide Association Studies: A Case Study with A Pan-Cancer Cohort

Phenome Wide Association Studies (PheWAS) enables phenome-wide scans to discover novel associations between genotype and clinical phenotypes via linking available genomic reports and large-scale Electronic Health Record (EHR). Data heterogeneity from different EHR systems and genetic reports has been a critical challenge that hinders meaningful validation. To address this, we propose an FHIR-based framework to model the PheWAS study in a standard manner. We developed an FHIR-based data model profile to enable the standard representation of data elements from genetic reports and EHR data that are used in the PheWAS study. As a proof-of-concept, we implemented the proposed method using a cohort of 1,595 pan-cancer patients with genetic reports from Foundation Medicine as well as the corresponding lab tests and diagnosis from Mayo EHRs. A PheWAS study is conducted and 81 significant genotype-phenotype associations are identified, in which 36 significant associations for cancers are validated based on a literature review.

Construction, comparison and evolution of networks in life sciences and other disciplines

Network approaches have become pervasive in many research fields. They allow for a more comprehensive understanding of complex relationships between entities as well as their group-level properties and dynamics. Many networks change over time, be it within seconds or millions of years, depending on the nature of the network. Our focus will be on comparative network analyses in life sciences, where deciphering temporal network changes is a core interest of molecular, ecological, neuropsychological and evolutionary biologists. Further, we will take a journey through different disciplines, such as social sciences, finance and computational gastronomy, to present commonalities and differences in how networks change and can be analysed. Finally, we envision how borrowing ideas from these disciplines could enrich the future of life science research.

An atlas of evidence-based phenotypic associations across the mouse phenome

To date, reliable relationships between mammalian phenotypes, based on diagnostic test measurements, have not been reported on a large scale. The purpose of this study was to present a large mouse phenotype-phenotype relationships dataset as a reference resource, alongside detailed evaluation of the resource. We used bias-minimized comprehensive mouse phenotype data and applied association rule mining to a dataset consisting of only binary (normal and abnormal phenotypes) data to determine relationships among phenotypes. We present 3,686 evidence-based significant associations, comprising 345 phenotypes covering 60 biological systems (functions), and evaluate their characteristics in detail. To evaluate the relationships, we defined a set of phenotype-phenotype association pairs (PPAPs) as a module of phenotypic expression for each of the 345 phenotypes. By analyzing each PPAP, we identified phenotype sub-networks consisting of the largest numbers of phenotypes and distinct biological systems. Furthermore, using hierarchical clustering based on phenotype similarities among the 345 PPAPs, we identified seven community types within a putative phenome-wide association network. Moreover, to promote leverage of these data, we developed and published web-application tools. These mouse phenome-wide phenotype-phenotype association data reveal general principles of relationships among mammalian phenotypes and provide a reference resource for biomedical analyses.

A phenome-wide association study (PheWAS) in the Population Architecture using Genomics and Epidemiology (PAGE) study reveals potential pleiotropy in African Americans

We performed a hypothesis-generating phenome-wide association study (PheWAS) to identify and characterize cross-phenotype associations, where one SNP is associated with two or more phenotypes, between thousands of genetic variants assayed on the Metabochip and hundreds of phenotypes in 5,897 African Americans as part of the Population Architecture using Genomics and Epidemiology (PAGE) I study. The PAGE I study was a National Human Genome Research Institute-funded collaboration of four study sites accessing diverse epidemiologic studies genotyped on the Metabochip, a custom genotyping chip that has dense coverage of regions in the genome previously associated with cardio-metabolic traits and outcomes in mostly European-descent populations. Here we focus on identifying novel phenome-genome relationships, where SNPs are associated with more than one phenotype. To do this, we performed a PheWAS, testing each SNP on the Metabochip for an association with up to 273 phenotypes in the participating PAGE I study sites. We identified 133 putative pleiotropic variants, defined as SNPs associated at an empirically derived p-value threshold of p<0.01 in two or more PAGE study sites for two or more phenotype classes. We further annotated these PheWAS-identified variants using publicly available functional data and local genetic ancestry. Amongst our novel findings is SPARC rs4958487, associated with increased glucose levels and hypertension. SPARC has been implicated in the pathogenesis of diabetes and is also known to have a potential role in fibrosis, a common consequence of multiple conditions including hypertension. The SPARC example and others highlight the potential that PheWAS approaches have in improving our understanding of complex disease architecture by identifying novel relationships between genetic variants and an array of common human phenotypes.

Phenotype Algorithm based Big Data Analytics for Cancer Diagnose

Nowadays, Cancer diagnosis is one of the major challenging characteristics for treating cancer. The reality of cancer patients rely on the diagnosis of cancer at the early stages (either in stage 1 or stage 2). If the cancer is diagnosed in stage 3 or later stages means the changes of survival of the patient will become more critical. Normally, single patient records will generate a huge amount of data if the data could be manage and analyze means to solve many problems for identifying the patterns it will leads to diagnose the cancer. Recent work several machine learning algorithms are introduced for the classification of cancer. However still the classification accuracy of machine learning algorithms are reduced because of huge number of samples. So the proposed work introduces a new Hadoop Distributed File System (HDFS) is focused in this work. In this paper, the proposed phenotype techniques are used which handle and classifies the raw EHR (Electronic Health Record) and EMR (Electronic Medical Record). It is based on the HDFS and Two-Phase Map Reduce. Phenotype algorithm uses NLP (National Language Processing) tool which will analyze and classify the cancer patient data like gene mapping, age related data, image and ultrasonic frequency processing, identification and analysis of irregularities, disease and personal histories. In this paper, the three factorized model is used which calculates the mean score values. The values are calculated by disease stage, pain status, etc. This paper focuses big data analytics for cancer diagnosis and the simulation results shows the proposed system produces the highest performance.

Multi-Phenotype Association Decomposition: Unraveling Complex Gene-Phenotype Relationships

Various patterns of multi-phenotype associations (MPAs) exist in the results of Genome Wide Association Studies (GWAS) involving different topologies of single nucleotide polymorphism (SNP)-phenotype associations. These can provide interesting information about the different impacts of a gene on closely related phenotypes or disparate phenotypes (pleiotropy). In this work we present MPA Decomposition, a new network-based approach which decomposes the results of a multi-phenotype GWAS study into three bipartite networks, which, when used together, unravel the multi-phenotype signatures of genes on a genome-wide scale. The decomposition involves the construction of a phenotype powerset space, and subsequent mapping of genes into this new space. Clustering of genes in this powerset space groups genes based on their detailed MPA signatures. We show that this method allows us to find multiple different MPA and pleiotropic signatures within individual genes and to classify and cluster genes based on these SNP-phenotype association topologies. We demonstrate the use of this approach on a GWAS analysis of a large population of 882 Populus trichocarpa genotypes using untargeted metabolomics phenotypes. This method should prove invaluable in the interpretation of large GWAS datasets and aid in future synthetic biology efforts designed to optimize phenotypes of interest.

High-throughput phenotyping for crop improvement in the genomics era

Tremendous progress has been made with continually expanding genomics technologies to unravel and understand crop genomes. However, the impact of genomics data on crop improvement is still far from satisfactory, in large part due to a lack of effective phenotypic data; our capacity to collect useful high quality phenotypic data lags behind the current capacity to generate high-throughput genomics data. Thus, the research bottleneck in plant sciences is shifting from genotyping to phenotyping. This article review the current status of efforts made in the last decade to systematically collect phenotypic data to alleviate this 'phenomics bottlenecks' by recording trait data through sophisticated non-invasive imaging, spectroscopy, image analysis, robotics, high-performance computing facilities and phenomics databases. These modern phenomics platforms and tools aim to record data on traits like plant development, architecture, plant photosynthesis, growth or biomass productivity, on hundreds to thousands of plants in a single day, as a phenomics revolution. It is believed that this revolution will provide plant scientists with the knowledge and tools necessary for unlocking information coded in plant genomes. Efforts have been also made to present the advances made in the last 10 years in phenomics platforms and their use in generating phenotypic data on different traits in several major crops including rice, wheat, barley, and maize. The article also highlights the need for phenomics databases and phenotypic data sharing for crop improvement. The phenomics data generated has been used to identify genes/QTL through QTL mapping, association mapping and genome-wide association studies (GWAS) for genomics-assisted breeding (GAB) for crop improvement.

Bayesian meta-analysis across genome-wide association studies of diverse phenotypes

Genome-wide association studies (GWAS) are a powerful tool for understanding the genetic basis of diseases and traits, but most studies have been conducted in isolation, with a focus on either a single or a set of closely related phenotypes. We describe MetABF, a simple Bayesian framework for performing integrative meta-analysis across multiple GWAS using summary statistics. The approach is applicable across a wide range of study designs and can increase the power by 50% compared with standard frequentist tests when only a subset of studies have a true effect. We demonstrate its utility in a meta-analysis of 20 diverse GWAS which were part of the Wellcome Trust Case Control Consortium 2. The novelty of the approach is its ability to explore, and assess the evidence for a range of possible true patterns of association across studies in a computationally efficient framework.

A PheWAS study of a large observational epidemiological cohort of African Americans from the REGARDS study

BACKGROUND

Cardiovascular disease, diabetes, and kidney disease are among the leading causes of death and disability worldwide. However, knowledge of genetic determinants of those diseases in African Americans remains limited.

RESULTS

In our study, associations between 4956 GWAS catalog reported SNPs and 67 traits were examined among 7726 African Americans from the REasons for Geographic and Racial Differences in Stroke (REGARDS) study, which is focused on identifying factors that increase stroke risk. The prevalent and incident phenotypes studied included inflammation, kidney traits, cardiovascular traits and cognition. Our results validated 29 known associations, of which eight associations were reported for the first time in African Americans.

CONCLUSION

Our cross-racial validation of GWAS findings provide additional evidence for the important roles of these loci in the disease process and may help identify genes especially important for future functional validation.

From high definition precision healthcare to precision public oral health: opportunities and challenges

In anticipation of a major transformation in healthcare, this review provides highlights that anticipate the near future for oral public health (and beyond). Personalized or precision healthcare reflects the expectation that advances in genomics, imaging, and other domains will extend our risk assessment, diagnostic, and prognostic capabilities, and enables more effective prevention and therapeutic options for all Americans. Meanwhile, the current healthcare system does not meet cost, access, or quality criteria for all Americans. It is now an imperative that the success of "smart," quality, and cost-effective high definition precision healthcare requires a public health perspective for several reasons: a) to enhance generalizability, b) to assess methods of implementation, and c) to focus on both risk and prevention in large and small populations, thereby providing a balance between the generation of long-term knowledge and short-term health gains. Sensitivity and resolution, reasonable cost, access to all Americans, coordinated comprehensive care, and advances in whole genome sequencing (WGS) and big data analyses, coupled to other advances in biotechnology and digital/artificial intelligence/machine learning devices, and the behavioral, social, and environmental sciences, offer remarkable opportunities to improve the health and wellness of the American people [genotype + phenotype + environment + behavior = high definition healthcare]. The opportunity is to significantly improve the well-being and life expectancy of all people across the lifespan including the least-advantaged people in our society and potentially increase access, reduce the national costs, and improve health outcomes.

Genomic and Phenomic Research in the 21st Century

The field of human genomics has changed dramatically over time. Initial genomic studies were predominantly restricted to rare disorders in small families. Over the past decade, researchers changed course from family-based studies and instead focused on common diseases and traits in populations of unrelated individuals. With further advancements in biobanking, computer science, electronic health record (EHR) data, and more affordable high-throughput genomics, we are experiencing a new paradigm in human genomic research. Rapidly changing technologies and resources now make it possible to study thousands of diseases simultaneously at the genomic level. This review will focus on these advancements as scientists begin to incorporate phenome-wide strategies in human genomic research to understand the etiology of human diseases and develop new drugs to treat them.

Genome sequencing in the clinic: the past, present, and future of genomic medicine

Genomic sequencing has undergone massive expansion in the past 10 yr, from a rarely used research tool into an approach that has broad applications in a clinical setting. From rare disease to cancer, genomics is transforming our knowledge of biology. The transition from targeted gene sequencing, to whole exome sequencing, to whole genome sequencing has only been made possible due to rapid advancements in technologies and informatics that have plummeted the cost per base of DNA sequencing and analysis. The tools of genomics have resolved the etiology of disease for previously undiagnosable conditions, identified cancer driver gene variants, and have impacted the understanding of pathophysiology for many diseases. However, this expansion of use has also highlighted research's current voids in knowledge. The lack of precise animal models for gene-to-function association, lack of tools for analysis of genomic structural changes, skew in populations used for genetic studies, publication biases, and the "Unknown Proteome" all contribute to voids needing filled for genomics to work in a fast-paced clinical setting. The future will hold the tools to fill in these voids, with new data sets and the continual development of new technologies allowing for expansion of genomic medicine, ushering in the days to come for precision medicine. In this review we highlight these and other points in hopes of advancing and guiding precision medicine into the future for optimal success.