Detection of pleiotropy through a Phenome-wide association study (PheWAS) of epidemiologic data as part of the Environmental Architecture for Genes Linked to Environment (EAGLE) study

The interplay of sex and genotype in disease associations: a comprehensive network analysis in the UK Biobank

BACKGROUND

Disease comorbidities and longer-term complications, arising from biologically related associations across phenotypes, can lead to increased risk of severe health outcomes. Given that many diseases exhibit sex-specific differences in their genetics, our objective was to determine whether genotype-by-sex (GxS) interactions similarly influence cross-phenotype associations. Through comparison of sex-stratified disease-disease networks (DDNs)-where nodes represent diseases and edges represent their relationships-we investigate sex differences in patterns of polygenicity and pleiotropy between diseases.

RESULTS

Using UK Biobank summary statistics, we built male- and female-specific DDNs for 103 diseases. This revealed that male and female diseasomes have similar topology and central diseases (e.g., hypertensive, chronic respiratory, and thyroid-based disorders), yet some phenotypes exhibit sex-specific influence in cross-phenotype associations. Multiple sclerosis and osteoarthritis are central only in the female DDN, while cardiometabolic diseases and skin cancer are more prominent in the male DDN. Edge comparison indicated similar shared genetics between the two graphs relative to a random model of disease association, though notable discrepancies in embedding distances and clustering patterns imply a more expansive genetic influence on multimorbidity risk for females than males. Analysis of pleiotropic contributions of two sexually-dimorphic single-nucleotide polymorphisms related to thyroid disorders further validated a distinct genetic architecture across sexes that influences associations, confirmed through examination of corresponding gene expression profiles from the GTEx Portal.

CONCLUSIONS

Our analysis affirms the presence of GxS interactions in cross-phenotype associations, emphasizing the need to investigate the role of sex in disease onset and its importance in biomedical discovery and precision medicine research.

Genetic variants for Alzheimer's disease and comorbid conditions

BACKGROUND

Alzheimer's disease and related dementias (ADRD) frequently co-occur with comorbidities such as diabetes and cardiovascular diseases in elderly populations.

OBJECTIVE

Utilize a life-course approach to identify genetic variants that are associated with the co-occurrence of ADRD and another comorbid condition.

METHODS

Research data from African American participants of the Indianapolis-Ibadan Dementia Project (IIDP) linked with electronic medical record (EMR) data and genome-wide association study (GWAS) data were utilized. The age of onset for ADRD was obtained from longitudinal follow-up of the IIDP study. Age of onset for comorbid conditions was obtained from EMR. The analysis included 1177 African Americans, among whom 174 were diagnosed with ADRD. A semi-parametric marginal bivariate survival model was used to examine the influence of single nucleotide polymorphisms (SNPs) on dual time-to-event outcomes while adjusting for sex, years of education, and the first principal component of GWAS data.

RESULTS

Targeted analysis of 20 SNPs that were reported to be associated with ADRD revealed that six were significantly associated with dual-disease outcomes, specifically congestive heart failure and cancer. In addition, eight novel SNPs were identified for associations with both ADRD and a comorbid condition.

CONCLUSIONS

Using a bivariate survival model approach, we identified genetic variants associated not only with ADRD, but also with comorbid conditions. Our utilization of dual-disease models represents a novel analytic strategy for uncovering shared genetic variants for multiple disease phenotypes.

Differential effects of environmental exposures on clinically relevant endophenotypes between sexes

Sex and gender differences play a crucial role in health and disease outcomes. This study used data from the National Health and Nutrition Examination Survey to explore how environmental exposures affect health-related traits differently in males and females. We utilized a sex-stratified phenomic environment-wide association study (PheEWAS), which allowed the identification of associations across a wide range of phenotypes and environmental exposures. We examined associations between 272 environmental exposures, including smoking-related exposures such as cotinine levels and smoking habits, and 58 clinically relevant blood phenotypes, such as serum albumin and homocysteine levels. Our analysis identified 119 sex-specific associations. For example, smoking-related exposures had a stronger impact on increasing homocysteine, hemoglobin, and hematocrit levels in females while reducing serum albumin and bilirubin levels and increasing c-reactive protein levels more significantly in males. These findings suggest mechanisms by which smoking exposure may pose higher cardiovascular risks and greater induced hypoxia for women, and greater inflammatory and immune responses in men. The results highlight the importance of considering sex differences in biomedical research. Understanding these differences can help develop more personalized and effective health interventions and improve clinical outcomes for both men and women.

Uncovering genetic associations in the human diseasome using an endophenotype-augmented disease network

MOTIVATION

Many diseases, particularly cardiometabolic disorders, exhibit complex multimorbidities with one another. An intuitive way to model the connections between phenotypes is with a disease-disease network (DDN), where nodes represent diseases and edges represent associations, such as shared single-nucleotide polymorphisms (SNPs), between pairs of diseases. To gain further genetic understanding of molecular contributors to disease associations, we propose a novel version of the shared-SNP DDN (ssDDN), denoted as ssDDN+, which includes connections between diseases derived from genetic correlations with intermediate endophenotypes. We hypothesize that a ssDDN+ can provide complementary information to the disease connections in a ssDDN, yielding insight into the role of clinical laboratory measurements in disease interactions.

RESULTS

Using PheWAS summary statistics from the UK Biobank, we constructed a ssDDN+ revealing hundreds of genetic correlations between diseases and quantitative traits. Our augmented network uncovers genetic associations across different disease categories, connects relevant cardiometabolic diseases, and highlights specific biomarkers that are associated with cross-phenotype associations. Out of the 31 clinical measurements under consideration, HDL-C connects the greatest number of diseases and is strongly associated with both type 2 diabetes and heart failure. Triglycerides, another blood lipid with known genetic causes in non-mendelian diseases, also adds a substantial number of edges to the ssDDN. This work demonstrates how association with clinical biomarkers can better explain the shared genetics between cardiometabolic disorders. Our study can facilitate future network-based investigations of cross-phenotype associations involving pleiotropy and genetic heterogeneity, potentially uncovering sources of missing heritability in multimorbidities.

AVAILABILITY AND IMPLEMENTATION

The generated ssDDN+ can be explored at https://hdpm.biomedinfolab.com/ddn/biomarkerDDN.

Causal genetic regulation of DNA replication on immune microenvironment in colorectal tumorigenesis: Evidenced by an integrated approach of trans-omics and GWAS

The interplay between DNA replication stress and immune microenvironment alterations is known to play a crucial role in colorectal tumorigenesis, but a comprehensive understanding of their association with and relevant biomarkers involved in colorectal tumorigenesis is lacking. To address this gap, we conducted a study aiming to investigate this association and identify relevant biomarkers. We analyzed transcriptomic and proteomic profiles of 904 colorectal tumor tissues and 342 normal tissues to examine pathway enrichment, biological activity, and the immune microenvironment. Additionally, we evaluated genetic effects of single variants and genes on colorectal cancer susceptibility using data from genome-wide association studies (GWASs) involving both East Asian (7062 cases and 195745 controls) and European (24476 cases and 23073 controls) populations. We employed mediation analysis to infer the causal pathway, and applied multiplex immunofluorescence to visualize colocalized biomarkers in colorectal tumors and immune cells. Our findings revealed that both DNA replication activity and the flap structure-specific endonuclease 1 ( FEN1) gene were significantly enriched in colorectal tumor tissues, compared with normal tissues. Moreover, a genetic variant rs4246215 G>T in FEN1 was associated with a decreased risk of colorectal cancer (odds ratio = 0.94, 95% confidence interval: 0.90-0.97, P meta = 4.70 × 10 -9). Importantly, we identified basophils and eosinophils that both exhibited a significantly decreased infiltration in colorectal tumors, and were regulated by rs4246215 through causal pathways involving both FEN1 and DNA replication. In conclusion, this trans-omics incorporating GWAS data provides insights into a plausible pathway connecting DNA replication and immunity, expanding biological knowledge of colorectal tumorigenesis and therapeutic targets.

Uncovering genetic associations in the human diseasome using an endophenotype-augmented disease network

Many diseases exhibit complex multimorbidities with one another. An intuitive way to model the connections between phenotypes is with a disease-disease network (DDN), where nodes represent diseases and edges represent associations, such as shared single-nucleotide polymorphisms (SNPs), between pairs of diseases. To gain further genetic understanding of molecular contributors to disease associations, we propose a novel version of the shared-SNP DDN (ssDDN), denoted as ssDDN+, which includes connections between diseases derived from genetic correlations with endophenotypes. We hypothesize that a ssDDN+ can provide complementary information to the disease connections in a ssDDN, yielding insight into the role of clinical laboratory measurements in disease interactions. Using PheWAS summary statistics from the UK Biobank, we constructed a ssDDN+ revealing hundreds of genetic correlations between disease phenotypes and quantitative traits. Our augmented network uncovers genetic associations across different disease categories, connects relevant cardiometabolic diseases, and highlights specific biomarkers that are associated with cross-phenotype associations. Out of the 31 clinical measurements under consideration, HDL-C connects the greatest number of diseases and is strongly associated with both type 2 diabetes and diabetic retinopathy. Triglycerides, another blood lipid with known genetics causes in non-mendelian diseases, also adds a substantial number of edges to the ssDDN. Our study can facilitate future network-based investigations of cross-phenotype associations involving pleiotropy and genetic heterogeneity, potentially uncovering sources of missing heritability in multimorbidities.

Genetic Factors Associated with Response to Vitamin E Treatment in NAFLD

Non-alcoholic fatty liver disease (NAFLD) is becoming the predominant liver disease worldwide, and vitamin E has been clinically shown to improve histological parameters in a subset of patients. In this narrative review, we investigate whether genetic factors may help to explain why some patients show histological improvements upon high-dose alpha-tocopherol (αT) treatment while others do not. In summary, we identified two factors that are associated with treatment response, including genetic variations in haptoglobin as well as fatty acid desaturase 1/2 (FADS1/FADS2). Other genetic variants such as in alpha-tocopherol transfer protein (αTTP), tocopherol associated protein (TAP), transmembrane 6 superfamily 2 (TM6SF2), cluster of differentiation 36 (CD36), and proteins involved in lipoprotein metabolism may also play a role, but have not yet been investigated in a clinical context. We propose to further validate these associations in larger populations, to then use them as a clinical tool to identify the subset of patients that will benefit the most from vitamin E supplementation.

Large-scale genomic analyses reveal insights into pleiotropy across circulatory system diseases and nervous system disorders

Clinical and epidemiological studies have shown that circulatory system diseases and nervous system disorders often co-occur in patients. However, genetic susceptibility factors shared between these disease categories remain largely unknown. Here, we characterized pleiotropy across 107 circulatory system and 40 nervous system traits using an ensemble of methods in the eMERGE Network and UK Biobank. Using a formal test of pleiotropy, five genomic loci demonstrated statistically significant evidence of pleiotropy. We observed region-specific patterns of direction of genetic effects for the two disease categories, suggesting potential antagonistic and synergistic pleiotropy. Our findings provide insights into the relationship between circulatory system diseases and nervous system disorders which can provide context for future prevention and treatment strategies.

Genetic variation in ALDH4A1 is associated with muscle health over the lifespan and across species

The influence of genetic variation on the aging process, including the incidence and severity of age-related diseases, is complex. Here, we define the evolutionarily conserved mitochondrial enzyme ALH-6/ALDH4A1 as a predictive biomarker for age-related changes in muscle health by combining Caenorhabditis elegans genetics and a gene-wide association scanning (GeneWAS) from older human participants of the US Health and Retirement Study (HRS). In a screen for mutations that activate oxidative stress responses, specifically in the muscle of C. elegans, we identified 96 independent genetic mutants harboring loss-of-function alleles of alh-6, exclusively. Each of these genetic mutations mapped to the ALH-6 polypeptide and led to the age-dependent loss of muscle health. Intriguingly, genetic variants in ALDH4A1 show associations with age-related muscle-related function in humans. Taken together, our work uncovers mitochondrial alh-6/ALDH4A1 as a critical component to impact normal muscle aging across species and a predictive biomarker for muscle health over the lifespan.

Ageing is inevitable, but what makes one person ‘age well’ and another decline more quickly remains largely unknown. While many aspects of ageing are clearly linked to genetics, the specific genes involved often remain unidentified.

Sarcopenia is an age-related condition affecting the muscles. It involves a gradual loss of muscle mass that becomes faster with age, and is associated with loss of mobility, decreased quality of life, and increased risk of death. Around half of all people aged 80 and over suffer from sarcopenia. Several lifestyle factors, especially poor diet and lack of exercise, are associated with the condition, but genetics is also involved: the condition accelerates more quickly in some people than others, and even fit, physically active individuals can be affected.

To study the genetics of conditions like sarcopenia, researchers often use animals like flies or worms, which have short generation times but share genetic similarities with humans. For example, the worm Caenorhabditis elegans has equivalents of several human muscle genes, including the gene alh-6. In worms, alh-6 is important for maintaining energy supply to the muscles, and mutating it not only leads to muscle damage but also to premature ageing. Given this insight, Villa, Stuhr, Yen et al. wanted to determine if variation in the human version of alh-6, ALDH4A1, also contributes to individual differences in muscle ageing and decline in humans.

Evaluating variation in this gene required a large amount of genetic data from older adults. These were taken from a continuous study that follows >35,000 older adults. Importantly, the study collects not only information on gene sequences but also measures of muscle health and performance over time for each individual. Analysis of these genetic data revealed specific small variations in the DNA of ALDH4A1, all of which associated with reduced muscle health.

Follow-up experiments in worms used genetic engineering techniques to test how variation in the worm alh-6 gene could influence age-related health. The resulting mutant worms developed muscle problems much earlier than their normal counterparts, supporting the role of alh-6/ALDH4A1 in determining muscle health across the lifespan of both worms and humans.

These results have identified a key influencer of muscle health during ageing in worms, and emphasize the importance of validating effects of genetic variation among humans during this process. Villa, Stuhr, Yen et al. hope that this study will help researchers find more genetic ‘markers’ of muscle health, and ultimately allow us to predict an individual’s risk of sarcopenia based on their genetic make-up.

Manifestations of Alzheimer's disease genetic risk in the blood are evident in a multiomic analysis in healthy adults aged 18 to 90

Genetics play an important role in late-onset Alzheimer's Disease (AD) etiology and dozens of genetic variants have been implicated in AD risk through large-scale GWAS meta-analyses. However, the precise mechanistic effects of most of these variants have yet to be determined. Deeply phenotyped cohort data can reveal physiological changes associated with genetic risk for AD across an age spectrum that may provide clues to the biology of the disease. We utilized over 2000 high-quality quantitative measurements obtained from blood of 2831 cognitively normal adult clients of a consumer-based scientific wellness company, each with CLIA-certified whole-genome sequencing data. Measurements included: clinical laboratory blood tests, targeted chip-based proteomics, and metabolomics. We performed a phenome-wide association study utilizing this diverse blood marker data and 25 known AD genetic variants and an AD-specific polygenic risk score (PGRS), adjusting for sex, age, vendor (for clinical labs), and the first four genetic principal components; sex-SNP interactions were also assessed. We observed statistically significant SNP-analyte associations for five genetic variants after correction for multiple testing (for SNPs in or near NYAP1, ABCA7, INPP5D, and APOE), with effects detectable from early adulthood. The ABCA7 SNP and the APOE2 and APOE4 encoding alleles were associated with lipid variability, as seen in previous studies; in addition, six novel proteins were associated with the e2 allele. The most statistically significant finding was between the NYAP1 variant and PILRA and PILRB protein levels, supporting previous functional genomic studies in the identification of a putative causal variant within the PILRA gene. We did not observe associations between the PGRS and any analyte. Sex modified the effects of four genetic variants, with multiple interrelated immune-modulating effects associated with the PICALM variant. In post-hoc analysis, sex-stratified GWAS results from an independent AD case-control meta-analysis supported sex-specific disease effects of the PICALM variant, highlighting the importance of sex as a biological variable. Known AD genetic variation influenced lipid metabolism and immune response systems in a population of non-AD individuals, with associations observed from early adulthood onward. Further research is needed to determine whether and how these effects are implicated in early-stage biological pathways to AD. These analyses aim to complement ongoing work on the functional interpretation of AD-associated genetic variants.

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.

Methodology in phenome-wide association studies: a systematic review

Phenome-wide association study (PheWAS) has been increasingly used to identify novel genetic associations across a wide spectrum of phenotypes. This systematic review aims to summarise the PheWAS methodology, discuss the advantages and challenges of PheWAS, and provide potential implications for future PheWAS studies. Medical Literature Analysis and Retrieval System Online (MEDLINE) and Excerpta Medica Database (EMBASE) databases were searched to identify all published PheWAS studies up until 24 April 2021. The PheWAS methodology incorporating how to perform PheWAS analysis and which software/tool could be used, were summarised based on the extracted information. A total of 1035 studies were identified and 195 eligible articles were finally included. Among them, 137 (77.0%) contained 10 000 or more study participants, 164 (92.1%) defined the phenome based on electronic medical records data, 140 (78.7%) used genetic variants as predictors, and 73 (41.0%) conducted replication analysis to validate PheWAS findings and almost all of them (94.5%) received consistent results. The methodology applied in these PheWAS studies was dissected into several critical steps, including quality control of the phenome, selecting predictors, phenotyping, statistical analysis, interpretation and visualisation of PheWAS results, and the workflow for performing a PheWAS was established with detailed instructions on each step. This study provides a comprehensive overview of PheWAS methodology to help practitioners achieve a better understanding of the PheWAS design, to detect understudied or overstudied outcomes, and to direct their research by applying the most appropriate software and online tools for their study data structure.

Integrative omics provide biological and clinical insights into acute respiratory distress syndrome

PURPOSE

Acute respiratory distress syndrome (ARDS) is accompanied by a dysfunctional immune-inflammatory response following lung injury, including during coronavirus disease 2019 (COVID-19). Limited causal biomarkers exist for ARDS development. We sought to identify novel genetic susceptibility targets for ARDS to focus further investigation on their biological mechanism and therapeutic potential.

METHODS

Meta-analyses of ARDS genome-wide association studies were performed with 1250 cases and 1583 controls in Europeans, and 387 cases and 387 controls in African Americans. The functionality of novel loci was determined in silico using multiple omics approaches. The causality of 114 factors potentially involved in ARDS development was assessed using Mendelian Randomization analysis.

RESULTS

There was distinct genetic heterogeneity in ARDS between Europeans and African Americans. rs7967111 at 12p13.2 was functionally associated with ARDS susceptibility in Europeans (odds ratio = 1.38; P = 2.15 × 10^-8). Expression of two genes annotated at this locus, BORCS5 and DUSP16, was dynamic but ultimately decreased during ARDS development, as well as downregulated in immune cells alongside COVID-19 severity. Causal inference implied that comorbidity of inflammatory bowel disease and elevated levels of C-reactive protein and interleukin-10 causally increased ARDS risk, while vitamin D supplementation and vasodilator use ameliorated risk.

CONCLUSION

Our findings suggest a novel susceptibility locus in ARDS pathophysiology that implicates BORCS5 and DUSP16 as potentially acting in immune-inflammatory processes. This locus warrants further investigation to inform the development of therapeutic targets and clinical care strategies for ARDS, including those induced by COVID-19.

Statistical Impact of Sample Size and Imbalance on Multivariate Analysis in silico and A Case Study in the UK Biobank

Large-scale biobank cohorts coupled with electronic health records offer unprecedented opportunities to study genotype-phenotype relationships. Genome-wide association studies uncovered disease-associated loci through univariate methods, with the focus on one trait at a time. With genetic variants being identifiedfor thousands of traits, researchers found that 90% of human genetic loci are associated with more than one trait, highlighting the ubiquity of pleiotropy. Recently, multivariate methods have been proposed to effectively identify pleiotropy. However, the statistical performance in natural biomedical data, which often have unbalanced case-control sample sizes, is largely known. In this work, we designed 21 scenarios of real-data informed simulations to thoroughly evaluate the statistical characteristics of univariate and multivariate methods. Our results can serve as a reference guide for the application of multivariate methods. We also investigated potential pleiotropy across type II diabetes, Alzheimer's disease, atherosclerosis of arteries, depression, and atherosclerotic heart disease in the UK Biobank.

Functional Haplotype of LIPC Induces Triglyceride-Mediated Suppression of HDL-C Levels According to Genome-Wide Association Studies

Hepatic lipase (encoded by LIPC) is a glycoprotein in the triacylglycerol lipase family and mainly synthesized in and secreted from the liver. Previous studies demonstrated that hepatic lipase is crucial for reverse cholesterol transport and modulating metabolism and the plasma levels of several lipoproteins. This study was conducted to investigate the suppression effect of high-density lipoprotein cholesterol (HDL-C) levels in a genome-wide association study and explore the possible mechanisms linking triglyceride (TG) to LIPC variants and HDL-C. Genome-wide association data for TG and HDL-C were available for 4657 Taiwan-biobank participants. The prevalence of haplotypes in the LIPC promoter region and their effects were calculated. The cloned constructs of the haplotypes were expressed transiently in HepG2 cells and evaluated in a luciferase reporter assay. Genome-wide association analysis revealed that HDL-C was significantly associated with variations in LIPC after adjusting for TG. Three haplotypes (H1: TCG, H2: CTA and H3: CCA) in LIPC were identified. H2: CTA was significantly associated with HDL-C levels and H1: TCG suppressed HDL-C levels when a third factor, TG, was included in mediation analysis. The luciferase reporter assay further showed that the H2: CTA haplotype significantly inhibited luciferase activity compared with the H1: TCG haplotype. In conclusion, we identified a suppressive role for TG in the genome-wide association between LIPC and HDL-C. A functional haplotype of hepatic lipase may reduce HDL-C levels and is suppressed by TG.

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.

Mendelian randomization study of serum uric acid levels and diabetes risk: evidence from the Dongfeng-Tongji cohort

OBJECTIVE

Limited Mendelian randomization (MR) studies have assessed the causal relationship between serum uric acid levels and diabetes risk. Here we investigated causality between the serum uric acid concentration and diabetes risk in Chinese population.

RESEARCH DESIGN AND METHODS

The observational analysis, based on the Dongfeng-Tongji prospective cohort (n=15 195) we tested the association of serum uric acid levels with incident diabetes risk. In the instrumental variable analysis, we examined the association of the genetic risk score (GRS) of serum uric acid with diabetes risk in case-control design (2539 cases and 4595 controls) via MR analysis.

RESULTS

During a mean (SD) follow-up of 4.5 (0.5) years, 1156 incident diabetes cases were identified. Compared with those in the lowest quintile of serum uric acid levels, the HRs of incident diabetes were 1.19 (95% CI 0.96 to 1.48), 1.12 (95% CI 0.90 to 1.40), 1.38 (95% CI 1.12 to 1.70), and 1.51 (95% CI 1.23 to 1.87) for Q2, Q3, Q4 and Q5, respectively (P-trend <0.001). The GRS was strongly associated with serum uric acid levels (β=0.17, 95% CI 0.15 to 0.19; P=2.81×10^-67). However, no significant association was observed between the GRS and diabetes risk (OR=1.01, 95 CI 0.95 to 1.06; P=0.75).

CONCLUSIONS

Even though serum uric acid levels were significantly associated with increased incident diabetes risk, the results did not provide evidence for a causal relationship between them.

Large-Scale Simultaneous Testing of Cross-Covariance Matrices with Applications to PheWAS

Motivated by applications in phenome-wide association studies (PheWAS), we consider in this paper simultaneous testing of columns of high-dimensional cross-covariance matrices and develop a multiple testing procedure with theoretical guarantees. It is shown that the proposed testing procedure maintains a desired false discovery rate (FDR) and false discovery proportion (FDP) under mild regularity conditions. We also provide results on the magnitudes of the signals that can be detected with high power. Simulation studies demonstrate that the proposed procedure can be substantially more powerful than existing FDR controlling procedures in the presence of correlation of unknown structure. The proposed multiple testing procedure is applied to a PheWAS of two auto-immune genetic markers using a rheumatoid arthritis patient cohort constructed from the electronic medical records of Partners Healthcare System.

Genetic contributions to NAFLD: leveraging shared genetics to uncover systems biology

Nonalcoholic fatty liver disease (NAFLD) affects around a quarter of the global population, paralleling worldwide increases in obesity and metabolic syndrome. NAFLD arises in the context of systemic metabolic dysfunction that concomitantly amplifies the risk of cardiovascular disease and diabetes. These interrelated conditions have long been recognized to have a heritable component, and advances using unbiased association studies followed by functional characterization have created a paradigm for unravelling the genetic architecture of these conditions. A novel perspective is to characterize the shared genetic basis of NAFLD and other related disorders. This information on shared genetic risks and their biological overlap should in future enable the development of precision medicine approaches through better patient stratification, and enable the identification of preventive and therapeutic strategies. In this Review, we discuss current knowledge of the genetic basis of NAFLD and of possible pleiotropy between NAFLD and other liver diseases as well as other related metabolic disorders. We also discuss evidence of causality in NAFLD and other related diseases and the translational significance of such evidence, and future challenges from the study of genetic pleiotropy.

Fine mapping and identification of serum urate loci in American Indians: The Strong Heart Family Study

While studies have reported genetic loci affecting serum urate (SU) concentrations, few studies have been conducted in minority populations. Our objective for this study was to identify genetic loci regulating SU in a multigenerational family-based cohort of American Indians, the Strong Heart Family Study (SHFS). We genotyped 162,718 single nucleotide polymorphisms (SNPs) in 2000 SHFS participants using an Illumina MetaboChip array. A genome-wide association analysis of SU was conducted using measured genotype analysis approach accounting for kinships in SOLAR, and meta-analysis in METAL. Our results showed strong association of SU with rs4481233, rs9998811, rs7696092 and rs13145758 (minor allele frequency (MAF) = 25-44%; P < 3 × 10-14) of solute carrier family 2, member 9 (SLC2A9) and rs41481455, rs2231142 and rs1481012 (MAF = 29%; p < 3 × 10-9) of ATP-binding cassette protein, subfamily G, member 2 (ABCG2). Carriers of G alleles of rs9998811, rs4148155 and rs1481012 and A alleles of rs4481233, rs7696092 and rs13145758 and rs2231142 had lower SU concentrations as compared to non-carriers. Genetic analysis of SU conditional on significant SLC2A9 and ABCG2 SNPs revealed new loci, nucleobindin 1 (NUCB1) and neuronal PAS domain protein 4 (NPAS4) (p <6× 10-6). To identify American Indian-specific SNPs, we conducted targeted sequencing of key regions of SLC2A9. A total of 233 SNPs were identified of which 89 were strongly associated with SU (p < 7.1 × 10-10) and 117 were American Indian specific. Analysis of key SNPs in cohorts of Mexican-mestizos, European, Indian and East Asian ancestries showed replication of common SNPs, including our lead SNPs. Our results demonstrate the association of SU with uric acid transporters in a minority population of American Indians and potential novel associations of SU with neuronal-related genes which warrant further investigation.

Pleiotropic association of LIPC variants with lipid and urinary 8-hydroxy deoxyguanosine levels in a Taiwanese population

BACKGROUND

Hepatic lipase (HL, encoded by LIPC) is a glycoprotein primarily synthesized and secreted by hepatocytes. Previous studies had demonstrated that HL is crucial for reverse cholesterol transport and affects the metabolism, composition, and level of several lipoproteins. In current study, we investigated the association of LIPC (Lipase C, Hepatic Type) variants with circulating and urinary biomarker levels by using subgroup and mediation analyses.

METHODS

A total of 572 participants from Taiwan were genotyped for three LIPC single nucleotide polymorphisms (SNPs) by using TaqMan assay. Fasting levels of glucose, lipid profile, inflammation markers, urine creatinine and 8-hydroxy deoxyguanosine (8-OHdG) were measured. The chi-square test, 2-sample t test and Analysis of variance (ANOVA) were used to examine differences among variables and genotype frequencies.

RESULTS

SNPs rs2043085 and rs1532085 were significantly associated with urinary 8-OHdG levels, whereas all three SNPs were more significantly associated with Triglycerides (TG) or HDL-cholesterol (HDL-C) levels after additional adjustment for HDL-C or TG levels, respectively. Subgroup analyses revealed that the association of the LIPC SNPs with the levels of serum TG, HDL-C, and urinary 8-OHdG were predominantly observed in the men but not in the women. Differential associations of the LIPC SNPs with various lipid levels were observed in participants with different adiposity statuses. Mediation analyses indicated that TG levels acted as a suppressor masking the association of the LIPC genotypes with HDL-C levels, particularly in the men (Sobel test, all P < 0.01).

CONCLUSION

Our data revealed that interaction and suppression effects mediated the pleiotropic association of the LIPC variants. The effects of the LIPC SNPs depended on sex, adiposity status, and TG levels. Thus, our findings can provide a method for identifying high-risk populations of cardiovascular diseases for clinical diagnosis.

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.

An exploratory phenome wide association study linking asthma and liver disease genetic variants to electronic health records from the Estonian Biobank

The Estonian Biobank, governed by the Institute of Genomics at the University of Tartu (Biobank), has stored genetic material/DNA and continuously collected data since 2002 on a total of 52,274 individuals representing ~5% of the Estonian adult population and is increasing. To explore the utility of data available in the Biobank, we conducted a phenome-wide association study (PheWAS) in two areas of interest to healthcare researchers; asthma and liver disease. We used 11 asthma and 13 liver disease-associated single nucleotide polymorphisms (SNPs), identified from published genome-wide association studies, to test our ability to detect established associations. We confirmed 2 asthma and 5 liver disease associated variants at nominal significance and directionally consistent with published results. We found 2 associations that were opposite to what was published before (rs4374383:AA increases risk of NASH/NAFLD, rs11597086 increases ALT level). Three SNP-diagnosis pairs passed the phenome-wide significance threshold: rs9273349 and E06 (thyroiditis, p = 5.50x10-8); rs9273349 and E10 (type-1 diabetes, p = 2.60x10-7); and rs2281135 and K76 (non-alcoholic liver diseases, including NAFLD, p = 4.10x10-7). We have validated our approach and confirmed the quality of the data for these conditions. Importantly, we demonstrate that the extensive amount of genetic and medical information from the Estonian Biobank can be successfully utilized for scientific research.

GRIK5 Genetically Regulated Expression Associated with Eye and Vascular Phenomes: Discovery through Iteration among Biobanks, Electronic Health Records, and Zebrafish

Although the use of model systems for studying the mechanism of mutations that have a large effect is common, we highlight here the ways that zebrafish-model-system studies of a gene, GRIK5, that contributes to the polygenic liability to develop eye diseases have helped to illuminate a mechanism that implicates vascular biology in eye disease. A gene-expression prediction derived from a reference transcriptome panel applied to BioVU, a large electronic health record (EHR)-linked biobank at Vanderbilt University Medical Center, implicated reduced GRIK5 expression in diverse eye diseases. We tested the function of GRIK5 by depletion of its ortholog in zebrafish, and we observed reduced blood vessel numbers and integrity in the eye and increased vascular permeability. Analyses of EHRs in >2.6 million Vanderbilt subjects revealed significant comorbidity of eye and vascular diseases (relative risks 2-15); this comorbidity was confirmed in 150 million individuals from a large insurance claims dataset. Subsequent studies in >60,000 genotyped BioVU participants confirmed the association of reduced genetically predicted expression of GRIK5 with comorbid vascular and eye diseases. Our studies pioneer an approach that allows a rapid iteration of the discovery of gene-phenotype relationships to the primary genetic mechanism contributing to the pathophysiology of human disease. Our findings also add dimension to the understanding of the biology driven by glutamate receptors such as GRIK5 (also referred to as GLUK5 in protein form) and to mechanisms contributing to human eye diseases.

A phenome-wide association study to discover pleiotropic effects of PCSK9, APOB, and LDLR

We conducted an electronic health record (EHR)-based phenome-wide association study (PheWAS) to discover pleiotropic effects of variants in three lipoprotein metabolism genes PCSK9, APOB, and LDLR. Using high-density genotype data, we tested the associations of variants in the three genes with 1232 EHR-derived binary phecodes in 51,700 European-ancestry (EA) individuals and 585 phecodes in 10,276 African-ancestry (AA) individuals; 457 PCSK9, 730 APOB, and 720 LDLR variants were filtered by imputation quality (r 2 > 0.4), minor allele frequency (>1%), linkage disequilibrium (r 2 < 0.3), and association with LDL-C levels, yielding a set of two PCSK9, three APOB, and five LDLR variants in EA but no variants in AA. Cases and controls were defined for each phecode using the PheWAS package in R. Logistic regression assuming an additive genetic model was used with adjustment for age, sex, and the first two principal components. Significant associations were tested in additional cohorts from Vanderbilt University (n = 29,713), the Marshfield Clinic Personalized Medicine Research Project (n = 9562), and UK Biobank (n = 408,455). We identified one PCSK9, two APOB, and two LDLR variants significantly associated with an examined phecode. Only one of the variants was associated with a non-lipid disease phecode, ("myopia") but this association was not significant in the replication cohorts. In this large-scale PheWAS we did not find LDL-C-related variants in PCSK9, APOB, and LDLR to be associated with non-lipid-related phenotypes including diabetes, neurocognitive disorders, or cataracts.

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.

Precision medicine review: rare driver mutations and their biophysical classification

How can biophysical principles help precision medicine identify rare driver mutations? A major tenet of pragmatic approaches to precision oncology and pharmacology is that driver mutations are very frequent. However, frequency is a statistical attribute, not a mechanistic one. Rare mutations can also act through the same mechanism, and as we discuss below, “latent driver” mutations may also follow the same route, with “helper” mutations. Here, we review how biophysics provides mechanistic guidelines that extend precision medicine. We outline principles and strategies, especially focusing on mutations that drive cancer. Biophysics has contributed profoundly to deciphering biological processes. However, driven by data science, precision medicine has skirted some of its major tenets. Data science embodies genomics, tissue- and cell-specific expression levels, making it capable of defining genome- and systems-wide molecular disease signatures. It classifies cancer driver genes/mutations and affected pathways, and its associated protein structural data guide drug discovery. Biophysics complements data science. It considers structures and their heterogeneous ensembles, explains how mutational variants can signal through distinct pathways, and how allo-network drugs can be harnessed. Biophysics clarifies how one mutation—frequent or rare—can affect multiple phenotypic traits by populating conformations that favor interactions with other network modules. It also suggests how to identify such mutations and their signaling consequences. Biophysics offers principles and strategies that can help precision medicine push the boundaries to transform our insight into biological processes and the practice of personalized medicine. By contrast, “phenotypic drug discovery,” which capitalizes on physiological cellular conditions and first-in-class drug discovery, may not capture the proper molecular variant. This is because variants of the same protein can express more than one phenotype, and a phenotype can be encoded by several variants.

Human-Disease Phenotype Map Derived from PheWAS across 38,682 Individuals

Phenome-wide association studies (PheWASs) have been a useful tool for testing associations between genetic variations and multiple complex traits or diagnoses. Linking PheWAS-based associations between phenotypes and a variant or a genomic region into a network provides a new way to investigate cross-phenotype associations, and it might broaden the understanding of genetic architecture that exists between diagnoses, genes, and pleiotropy. We created a network of associations from one of the largest PheWASs on electronic health record (EHR)-derived phenotypes across 38,682 unrelated samples from the Geisinger's biobank; the samples were genotyped through the DiscovEHR project. We computed associations between 632,574 common variants and 541 diagnosis codes. Using these associations, we constructed a "disease-disease" network (DDN) wherein pairs of diseases were connected on the basis of shared associations with a given genetic variant. The DDN provides a landscape of intra-connections within the same disease classes, as well as inter-connections across disease classes. We identified clusters of diseases with known biological connections, such as autoimmune disorders (type 1 diabetes, rheumatoid arthritis, and multiple sclerosis) and cardiovascular disorders. Previously unreported relationships between multiple diseases were identified on the basis of genetic associations as well. The network approach applied in this study can be used to uncover interactions between diseases as a result of their shared, potentially pleiotropic SNPs. Additionally, this approach might advance clinical research and even clinical practice by accelerating our understanding of disease mechanisms on the basis of similar underlying genetic associations.

Detecting potential pleiotropy across cardiovascular and neurological diseases using univariate, bivariate, and multivariate methods on 43,870 individuals from the eMERGE network

The link between cardiovascular diseases and neurological disorders has been widely observed in the aging population. Disease prevention and treatment rely on understanding the potential genetic nexus of multiple diseases in these categories. In this study, we were interested in detecting pleiotropy, or the phenomenon in which a genetic variant influences more than one phenotype. Marker-phenotype association approaches can be grouped into univariate, bivariate, and multivariate categories based on the number of phenotypes considered at one time. Here we applied one statistical method per category followed by an eQTL colocalization analysis to identify potential pleiotropic variants that contribute to the link between cardiovascular and neurological diseases. We performed our analyses on ~530,000 common SNPs coupled with 65 electronic health record (EHR)-based phenotypes in 43,870 unrelated European adults from the Electronic Medical Records and Genomics (eMERGE) network. There were 31 variants identified by all three methods that showed significant associations across late onset cardiac- and neurologic- diseases. We further investigated functional implications of gene expression on the detected "lead SNPs" via colocalization analysis, providing a deeper understanding of the discovered associations. In summary, we present the framework and landscape for detecting potential pleiotropy using univariate, bivariate, multivariate, and colocalization methods. Further exploration of these potentially pleiotropic genetic variants will work toward understanding disease causing mechanisms across cardiovascular and neurological diseases and may assist in considering disease prevention as well as drug repositioning in future research.

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.

MR-PheWAS: exploring the causal effect of SUA level on multiple disease outcomes by using genetic instruments in UK Biobank

OBJECTIVES

We aimed to investigate the role of serum uric acid (SUA) level in a broad spectrum of disease outcomes using data for 120 091 individuals from UK Biobank.

METHODS

We performed a phenome-wide association study (PheWAS) to identify disease outcomes associated with SUA genetic risk loci. We then implemented conventional Mendelianrandomisation (MR) analysis to investigate the causal relevance between SUA level and disease outcomes identified from PheWAS. We next applied MR Egger analysis to detect and account for potential pleiotropy, which conventional MR analysis might mistake for causality, and used the HEIDI (heterogeneity in dependent instruments) test to remove cross-phenotype associations that were likely due to genetic linkage.

RESULTS

Our PheWAS identified 25 disease groups/outcomes associated with SUA genetic risk loci after multiple testing correction (P<8.57e-05). Our conventional MR analysis implicated a causal role of SUA level in three disease groups: inflammatory polyarthropathies (OR=1.22, 95% CI 1.11 to 1.34), hypertensive disease (OR=1.08, 95% CI 1.03 to 1.14) and disorders of metabolism (OR=1.07, 95% CI 1.01 to 1.14); and four disease outcomes: gout (OR=4.88, 95% CI 3.91 to 6.09), essential hypertension (OR=1.08, 95% CI 1.03 to 1.14), myocardial infarction (OR=1.16, 95% CI 1.03 to 1.30) and coeliac disease (OR=1.41, 95% CI 1.05 to 1.89). After balancing pleiotropic effects in MR Egger analysis, only gout and its encompassing disease group of inflammatory polyarthropathies were considered to be causally associated with SUA level. Our analysis highlighted a locus (ATXN2/S2HB3) that may influence SUA level and multiple cardiovascular and autoimmune diseases via pleiotropy.

CONCLUSIONS

Elevated SUA level is convincing to cause gout and inflammatory polyarthropathies, and might act as a marker for the wider range of diseases with which it associates. Our findings support further investigation on the clinical relevance of SUA level with cardiovascular, metabolic, autoimmune and respiratory diseases.

Recent Genetics and Epigenetics Approaches to PTSD

PURPOSE OF REVIEW

Following a life-threatening traumatic exposure, about 10% of those exposed are at considerable risk for developing posttraumatic stress disorder (PTSD), a severe and disabling syndrome characterized by uncontrollable intrusive memories, nightmares, avoidance behaviors, and hyperarousal in addition to impaired cognition and negative emotion symptoms. This review will explore recent genetic and epigenetic approaches to PTSD that explain some of the differential risk following trauma exposure.

RECENT FINDINGS

A substantial portion of the variance explaining differential risk responses to trauma exposure may be explained by differential inherited and acquired genetic and epigenetic risk. This biological risk is complemented by alterations in the functional regulation of genes via environmentally induced epigenetic changes, including prior childhood and adult trauma exposure. This review will cover recent findings from large-scale genome-wide association studies as well as newer epigenome-wide studies. We will also discuss future "phenome-wide" studies utilizing electronic medical records as well as targeted genetic studies focusing on mechanistic ways in which specific genetic or epigenetic alterations regulate the biological risk for PTSD.

PheWAS and Beyond: The Landscape of Associations with Medical Diagnoses and Clinical Measures across 38,662 Individuals from Geisinger

Most phenome-wide association studies (PheWASs) to date have used a small to moderate number of SNPs for association with phenotypic data. We performed a large-scale single-cohort PheWAS, using electronic health record (EHR)-derived case-control status for 541 diagnoses using International Classification of Disease version 9 (ICD-9) codes and 25 median clinical laboratory measures. We calculated associations between these diagnoses and traits with ∼630,000 common frequency SNPs with minor allele frequency > 0.01 for 38,662 individuals. In this landscape PheWAS, we explored results within diseases and traits, comparing results to those previously reported in genome-wide association studies (GWASs), as well as previously published PheWASs. We further leveraged the context of functional impact from protein-coding to regulatory regions, providing a deeper interpretation of these associations. The comprehensive nature of this PheWAS allows for novel hypothesis generation, the identification of phenotypes for further study for future phenotypic algorithm development, and identification of cross-phenotype associations.

A simulation study investigating power estimates in phenome-wide association studies

Background

Phenome-wide association studies (PheWAS) are a high-throughput approach to evaluate comprehensive associations between genetic variants and a wide range of phenotypic measures. PheWAS has varying sample sizes for quantitative traits, and variable numbers of cases and controls for binary traits across the many phenotypes of interest, which can affect the statistical power to detect associations. The motivation of this study is to investigate the various parameters which affect the estimation of statistical power in PheWAS, including sample size, case-control ratio, minor allele frequency, and disease penetrance.

Results

We performed a PheWAS simulation study, where we investigated variations in statistical power based on different parameters, such as overall sample size, number of cases, case-control ratio, minor allele frequency, and disease penetrance. The simulation was performed on both binary and quantitative phenotypic measures. Our simulation on binary traits suggests that the number of cases has more impact on statistical power than the case to control ratio; also, we found that a sample size of 200 cases or more maintains the statistical power to identify associations for common variants. For quantitative traits, a sample size of 1000 or more individuals performed best in the power calculations. We focused on common genetic variants (MAF > 0.01) in this study; however, in future studies, we will be extending this effort to perform similar simulations on rare variants.

Conclusions

This study provides a series of PheWAS simulation analyses that can be used to estimate statistical power for some potential scenarios. These results can be used to provide guidelines for appropriate study design for future PheWAS analyses.

Electronic supplementary material

The online version of this article (10.1186/s12859-018-2135-0) contains supplementary material, which is available to authorized users.

LIPC variants as genetic determinants of adiposity status, visceral adiposity indicators, and triglyceride-glucose (TyG) index-related parameters mediated by serum triglyceride levels

BACKGROUND

Visceral adiposity indicators and the product of triglyceride and fasting plasma glucose (TyG) index-related parameters are effective surrogate markers for insulin resistance (IR) and are predictors of metabolic syndrome and diabetes mellitus. However, their genetic determinants have not been previously reported. Pleiotropic associations of LIPC variants have been observed in lipid profiles and atherosclerotic cardiovascular diseases. We aimed to investigate LIPC polymorphisms as the genetic determinants of adiposity status, visceral adiposity indicators and TyG index-related parameters.

METHODS

A total of 592 participants from Taiwan were genotyped for three LIPC single nucleotide polymorphisms (SNPs).

RESULTS

The LIPC SNPs rs2043085 and rs1532085 were significantly associated with body mass index (BMI), waist circumference (WC), lipid accumulation product, visceral adiposity index, and TyG index-related parameters [including the TyG index, TyG with adiposity status (TyG-BMI), and TyG-WC index], whereas the rs1800588 SNP was only significantly associated with the TyG index. The associations became nonsignificant after further adjustment for serum TG levels. No significant association was observed between any the studied LIPC SNPs and IR status.

CONCLUSION

Our data revealed a pleiotropic association between the LIPC variants and visceral adiposity indicators and TyG index-related parameters, which are mediated by serum TG levels.

Current Scope and Challenges in Phenome-Wide Association Studies

PURPOSE OF REVIEW

Over many decades, researchers have been designing studies to investigate the relationship between genotypes and phenotypes to gain an understanding about the effect of genetics on disease. Recently, a high-throughput approach called phenome-wide associations studies (PheWAS) have been extensively used to identify associations between genetic variants and many diseases and traits simultaneously. In this review, we describe the value of PheWAS along with methodological issues and challenges in interpretation for current applications of PheWAS.

RECENT FINDINGS

PheWAS have uncovered a paradigm to identify new associations for genetic loci across many diseases. The application of PheWAS have been effective with phenotype data from electronic health records, epidemiological studies, and clinical trials data.

SUMMARY

The key strength of a PheWAS is to identify the association of one or more genetic variants with multiple phenotypes, which can showcase interconnections among the phenotypes due to shared genetic associations. While the PheWAS approach appears promising, there are a number of challenges that need to be addressed to provide additional robustness to PheWAS findings.

Embracing Complex Associations in Common Traits: Critical Considerations for Precision Medicine

Genome-wide association studies (GWAS) have identified numerous loci associated with human phenotypes. This approach, however, does not consider the richly diverse and complex environment with which humans interact throughout the life course, nor does it allow for interrelationships between genetic loci and across traits. As we move toward making precision medicine a reality, whereby we make predictions about disease risk based on genomic profiles, we need to identify improved predictive models of the relationship between genome and phenome. Methods that embrace pleiotropy (the effect of one locus on more than one trait), and gene-environment (G×E) and gene-gene (G×G) interactions, will further unveil the impact of alterations in biological pathways and identify genes that are only involved with disease in the context of the environment. This valuable information can be used to assess personal risk and choose the most appropriate medical interventions based on the genotype and environment of an individual, the whole premise of precision medicine.

Comparative analyses of population-scale phenomic data in electronic medical records reveal race-specific disease networks

Motivation: Underrepresentation of racial groups represents an important challenge and major gap in phenomics research. Most of the current human phenomics research is based primarily on European populations; hence it is an important challenge to expand it to consider other population groups. One approach is to utilize data from EMR databases that contain patient data from diverse demographics and ancestries. The implications of this racial underrepresentation of data can be profound regarding effects on the healthcare delivery and actionability. To the best of our knowledge, our work is the first attempt to perform comparative, population-scale analyses of disease networks across three different populations, namely Caucasian (EA), African American (AA) and Hispanic/Latino (HL).

Results: We compared susceptibility profiles and temporal connectivity patterns for 1988 diseases and 37 282 disease pairs represented in a clinical population of 1 025 573 patients. Accordingly, we revealed appreciable differences in disease susceptibility, temporal patterns, network structure and underlying disease connections between EA, AA and HL populations. We found 2158 significantly comorbid diseases for the EA cohort, 3265 for AA and 672 for HL. We further outlined key disease pair associations unique to each population as well as categorical enrichments of these pairs. Finally, we identified 51 key ‘hub’ diseases that are the focal points in the race-centric networks and of particular clinical importance. Incorporating race-specific disease comorbidity patterns will produce a more accurate and complete picture of the disease landscape overall and could support more precise understanding of disease relationships and patient management towards improved clinical outcomes.

Contacts: rong.chen@mssm.edu or joel.dudley@mssm.edu

Supplementary information:Supplementary data are available at Bioinformatics online.

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

OBJECTIVE

To compare three groupings of Electronic Health Record (EHR) billing codes for their ability to represent clinically meaningful phenotypes and to replicate known genetic associations. The three tested coding systems were the International Classification of Diseases, Ninth Revision, Clinical Modification (ICD-9-CM) codes, the Agency for Healthcare Research and Quality Clinical Classification Software for ICD-9-CM (CCS), and manually curated "phecodes" designed to facilitate phenome-wide association studies (PheWAS) in EHRs.

METHODS AND MATERIALS

We selected 100 disease phenotypes and compared the ability of each coding system to accurately represent them without performing additional groupings. The 100 phenotypes included 25 randomly-chosen clinical phenotypes pursued in prior genome-wide association studies (GWAS) and another 75 common disease phenotypes mentioned across free-text problem lists from 189,289 individuals. We then evaluated the performance of each coding system to replicate known associations for 440 SNP-phenotype pairs.

RESULTS

Out of the 100 tested clinical phenotypes, phecodes exactly matched 83, compared to 53 for ICD-9-CM and 32 for CCS. ICD-9-CM codes were typically too detailed (requiring custom groupings) while CCS codes were often not granular enough. Among 440 tested known SNP-phenotype associations, use of phecodes replicated 153 SNP-phenotype pairs compared to 143 for ICD-9-CM and 139 for CCS. Phecodes also generally produced stronger odds ratios and lower p-values for known associations than ICD-9-CM and CCS. Finally, evaluation of several SNPs via PheWAS identified novel potential signals, some seen in only using the phecode approach. Among them, rs7318369 in PEPD was associated with gastrointestinal hemorrhage.

CONCLUSION

Our results suggest that the phecode groupings better align with clinical diseases mentioned in clinical practice or for genomic studies. ICD-9-CM, CCS, and phecode groupings all worked for PheWAS-type studies, though the phecode groupings produced superior results.

Phenome-wide scanning identifies multiple diseases and disease severity phenotypes associated with HLA variants

Although many phenotypes have been associated with variants in human leukocyte antigen (HLA) genes, the full phenotypic impact of HLA variants across all diseases is unknown. We imputed HLA genomic variation from two populations of 28,839 and 8431 European ancestry individuals and tested association of HLA variation with 1368 phenotypes. A total of 104 four-digit and 92 two-digit HLA allele phenotype associations were significant in both discovery and replication cohorts, the strongest being HLA-DQB1*03:02 and type 1 diabetes. Four previously unidentified associations were identified across the spectrum of disease with two- and four-digit HLA alleles and 10 with nonsynonymous variants. Some conditions associated with multiple HLA variants and stronger associations with more severe disease manifestations were identified. A comprehensive, publicly available catalog of clinical phenotypes associated with HLA variation is provided. Examining HLA variant disease associations in this large data set allows comprehensive definition of disease associations to drive further mechanistic insights.

The joint effect of air pollution exposure and copy number variation on risk for autism

Autism spectrum disorder is a complex trait with a high degree of heritability as well as documented susceptibility from environmental factors. In this study the contributions of copy number variation, exposure to air pollutants, and the interaction between the two on autism risk, were evaluated in the population-based case-control Childhood Autism Risks from Genetics and Environment (CHARGE) Study. For the current investigation, we included only those CHARGE children (a) who met criteria for autism or typical development and (b) for whom our team had conducted both genetic evaluation of copy number burden and determination of environmental air pollution exposures based on mapping addresses from the pregnancy and early childhood. This sample consisted of 158 cases of children with autism and 147 controls with typical development. Multiple logistic regression models were fit with and without environmental variable-copy number burden interactions. We found no correlation between average air pollution exposure from conception to age 2 years and the child's CNV burden. We found a significant interaction in which a 1SD increase in duplication burden combined with a 1SD increase in ozone exposure was associated with an elevated autism risk (OR 3.4, P < 0.005) much greater than the increased risks associated with either genomic duplication (OR 1.85, 95% CI 1.25-2.73) or ozone (OR 1.20, 95% CI 0.93-1.54) alone. Similar results were obtained when CNV and ozone were dichotomized to compare those in the top quartile relative to those having a smaller CNV burden and lower exposure to ozone, and when exposures were assessed separately for pregnancy, the first year of life, and the second year of life. No interactions were observed for other air pollutants, even those that demonstrated main effects; ozone tends to be negatively correlated with the other pollutants examined. While earlier work has demonstrated interactions between the presence of a pathogenic CNV and an environmental exposure [Webb et al., 2016], these findings appear to be the first indication that global copy number variation may increase susceptibility to certain environmental factors, and underscore the need to consider both genomics and environmental exposures as well as the mechanisms by which each may amplify the risks for autism associated with the other. Autism Res 2017, 10: 1470-1480. © 2017 International Society for Autism Research, Wiley Periodicals, Inc.

Informatics and Data Analytics to Support Exposome-Based Discovery for Public Health

The complexity of the human exposome-the totality of environmental exposures encountered from birth to death-motivates systematic, high-throughput approaches to discover new environmental determinants of disease. In this review, we describe the state of science in analyzing the human exposome and provide recommendations for the public health community to consider in dealing with analytic challenges of exposome-based biomedical research. We describe extant and novel analytic methods needed to associate the exposome with critical health outcomes and contextualize the data-centered challenges by drawing parallels to other research endeavors such as human genomics research. We discuss efforts for training scientists who can bridge public health, genomics, and biomedicine in informatics and statistics. If an exposome data ecosystem is brought to fruition, it will likely play a role as central as genomic science has had in molding the current and new generations of biomedical researchers, computational scientists, and public health research programs.

Biomechanisms of Comorbidity: Reviewing Integrative Analyses of Multi-omics Datasets and Electronic Health Records

OBJECTIVES

Disease comorbidity is a pervasive phenomenon impacting patients' health outcomes, disease management, and clinical decisions. This review presents past, current and future research directions leveraging both phenotypic and molecular information to uncover disease similarity underpinning the biology and etiology of disease comorbidity.

METHODS

We retrieved ~130 publications and retained 59, ranging from 2006 to 2015, that comprise a minimum number of five diseases and at least one type of biomolecule. We surveyed their methods, disease similarity metrics, and calculation of comorbidities in the electronic health records, if present.

RESULTS

Among the surveyed studies, 44% generated or validated disease similarity metrics in context of comorbidity, with 60% being published in the last two years. As inputs, 87% of studies utilized intragenic loci and proteins while 13% employed RNA (mRNA, LncRNA or miRNA). Network modeling was predominantly used (35%) followed by statistics (28%) to impute similarity between these biomolecules and diseases. Studies with large numbers of biomolecules and diseases used network models or naïve overlap of disease-molecule associations, while machine learning, statistics, and information retrieval were utilized in smaller and moderate sized studies. Multiscale computations comprising shared function, network topology, and phenotypes were performed exclusively on proteins.

CONCLUSION

This review highlighted the growing methods for identifying the molecular mechanisms underpinning comorbidities that leverage multiscale molecular information and patterns from electronic health records. The survey unveiled that intergenic polymorphisms have been overlooked for similarity imputation compared to their intragenic counterparts, offering new opportunities to bridge the mechanistic and similarity gaps of comorbidity.

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

Beginning in the early 2000s, the accumulation of biospecimens linked to electronic health records (EHRs) made possible genome-phenome studies (i.e., comparative analyses of genetic variants and phenotypes) using only data collected as a by-product of typical health care. In addition to disease and trait genetics, EHRs proved a valuable resource for analyzing pharmacogenetic traits and developing reverse genetics approaches such as phenome-wide association studies (PheWASs). PheWASs are designed to survey which of many phenotypes may be associated with a given genetic variant. PheWAS methods have been validated through replication of hundreds of known genotype-phenotype associations, and their use has differentiated between true pleiotropy and clinical comorbidity, added context to genetic discoveries, and helped define disease subtypes, and may also help repurpose medications. PheWAS methods have also proven to be useful with research-collected data. Future efforts that integrate broad, robust collection of phenotype data (e.g., EHR data) with purpose-collected research data in combination with a greater understanding of EHR data will create a rich resource for increasingly more efficient and detailed genome-phenome analysis to usher in new discoveries in precision medicine.

Leveraging electronic health records to study pleiotropic effects on bipolar disorder and medical comorbidities

Patients with bipolar disorder (BD) have a high prevalence of comorbid medical illness. However, the mechanisms underlying these comorbidities with BD are not well known. Certain genetic variants may have pleiotropic effects, increasing the risk of BD and other medical illnesses simultaneously. In this study, we evaluated the association of BD-susceptibility genetic variants with various medical conditions that tend to co-exist with BD, using electronic health records (EHR) data linked to genome-wide single-nucleotide polymorphism (SNP) data. Data from 7316 Caucasian subjects were used to test the association of 19 EHR-derived phenotypes with 34 SNPs that were previously reported to be associated with BD. After Bonferroni multiple testing correction, P<7.7 × 10(-5) was considered statistically significant. The top association findings suggested that the BD risk alleles at SNP rs4765913 in CACNA1C gene and rs7042161 in SVEP1 may be associated with increased risk of 'cardiac dysrhythmias' (odds ratio (OR)=1.1, P=3.4 × 10(-3)) and 'essential hypertension' (OR=1.1, P=3.5 × 10(-3)), respectively. Although these associations are not statistically significant after multiple testing correction, both genes have been previously implicated with cardiovascular phenotypes. Moreover, we present additional evidence supporting these associations, particularly the association of the SVEP1 SNP with hypertension. This study shows the potential for EHR-based analyses of large cohorts to discover pleiotropic effects contributing to complex psychiatric traits and commonly co-occurring medical conditions.

eMERGE Phenome-Wide Association Study (PheWAS) identifies clinical associations and pleiotropy for stop-gain variants

BACKGROUND

We explored premature stop-gain variants to test the hypothesis that variants, which are likely to have a consequence on protein structure and function, will reveal important insights with respect to the phenotypes associated with them. We performed a phenome-wide association study (PheWAS) exploring the association between a selected list of functional stop-gain genetic variants (variation resulting in truncated proteins or in nonsense-mediated decay) and an extensive group of diagnoses to identify novel associations and uncover potential pleiotropy.

RESULTS

In this study, we selected 25 stop-gain variants: 5 stop-gain variants with previously reported phenotypic associations, and a set of 20 putative stop-gain variants identified using dbSNP. For the PheWAS, we used data from the electronic MEdical Records and GEnomics (eMERGE) Network across 9 sites with a total of 41,057 unrelated patients. We divided all these samples into two datasets by equal proportion of eMERGE site, sex, race, and genotyping platform. We calculated single effect associations between these 25 stop-gain variants and ICD-9 defined case-control diagnoses. We also performed stratified analyses for samples of European and African ancestry. Associations were adjusted for sex, site, genotyping platform and the first three principal components to account for global ancestry. We identified previously known associations, such as variants in LPL associated with hyperglyceridemia indicating that our approach was robust. We also found a total of three significant associations with p < 0.01 in both datasets, with the most significant replicating result being LPL SNP rs328 and ICD-9 code 272.1 "Disorder of Lipoid metabolism" (pdiscovery = 2.59x10-6, preplicating = 2.7x10-4). The other two significant replicated associations identified by this study are: variant rs1137617 in KCNH2 gene associated with ICD-9 code category 244 "Acquired Hypothyroidism" (pdiscovery = 5.31x103, preplicating = 1.15x10-3) and variant rs12060879 in DPT gene associated with ICD-9 code category 996 "Complications peculiar to certain specified procedures" (pdiscovery = 8.65x103, preplicating = 4.16x10-3).

CONCLUSION

In conclusion, this PheWAS revealed novel associations of stop-gained variants with interesting phenotypes (ICD-9 codes) along with pleiotropic effects.

Phenome-Wide Association Study to Explore Relationships between Immune System Related Genetic Loci and Complex Traits and Diseases

We performed a Phenome-Wide Association Study (PheWAS) to identify interrelationships between the immune system genetic architecture and a wide array of phenotypes from two de-identified electronic health record (EHR) biorepositories. We selected variants within genes encoding critical factors in the immune system and variants with known associations with autoimmunity. To define case/control status for EHR diagnoses, we used International Classification of Diseases, Ninth Revision (ICD-9) diagnosis codes from 3,024 Geisinger Clinic MyCode^® subjects (470 diagnoses) and 2,899 Vanderbilt University Medical Center BioVU biorepository subjects (380 diagnoses). A pooled-analysis was also carried out for the replicating results of the two data sets. We identified new associations with potential biological relevance including SNPs in tumor necrosis factor (TNF) and ankyrin-related genes associated with acute and chronic sinusitis and acute respiratory tract infection. The two most significant associations identified were for the C6orf10 SNP rs6910071 and “rheumatoid arthritis” (ICD-9 code category 714) (p_METAL = 2.58 x 10⁻⁹) and the ATN1 SNP rs2239167 and “diabetes mellitus, type 2” (ICD-9 code category 250) (p_METAL = 6.39 x 10⁻⁹). This study highlights the utility of using PheWAS in conjunction with EHRs to discover new genotypic-phenotypic associations for immune-system related genetic loci.

Integrating electronic health record genotype and phenotype datasets to transform patient care

The Health Information Technology for Economic and Clinical Health (HITECH) Act of 2009 mandates the development and implementation of electronic health record (EHR) systems across the country. While a primary goal is to improve the care of individual patients, EHRs are also key enabling resources for a vision of individualized (or personalized or precision) medicine: the aggregation of multiple EHRs within or across healthcare systems should allow discovery of patient subsets that have unusual and definable clinical trajectories that deviate importantly from the expected response in a "typical" patient. The spectrum of such personalized care can then extend from prevention to choice of medication to intensity or nature of follow-up.