Joint GWAS Analysis: Comparing similar GWAS at different genomic resolutions identifies novel pathway associations with six complex diseases

HLA-DQ and HLA-DRB1 alleles associated with Henoch-Schönlein purpura nephritis in Finnish pediatric population: a genome-wide association study

BACKGROUND

The pathophysiology of Henoch-Schönlein purpura (HSP) is still unclear, but several findings suggest that genetic factors may influence disease susceptibility. We aimed to perform a genome-wide association study (GWAS) in pediatric HSP patients with an emphasis on severe HSP nephritis.

METHODS

The study included 46 HSP patients, 42 of whom had undergone kidney biopsy. Forty-nine pediatric patients with an inflammatory bowel disease (IBD) served as an autoimmune disease control group while Finnish bone marrow and blood donors represented the general reference population (n = 18,757). GWAS was performed for HSP and IBD samples in a case-control manner against the reference population. The analysis also included imputation of human leukocyte antigen (HLA) alleles.

RESULTS

GWAS analysis in HSP revealed several polymorphisms from the HLA region that surpassed the genome-wide significance level. Three HLA class II alleles were also significantly more frequent in HSP than in the reference population: DQA1*01:01, DQB1*05:01, and DRB1*01:01. Haplotype DQA1*01:01/DQB1*05:01/DRB1*01:01 occurred in 43.5% of HSP patients, whereas its frequency was 8.2% in IBD patients and 15.0% in the reference population. HSP patients with this haplotype showed similar baseline clinical findings and outcome as HSP patients negative for the haplotype. In IBD patients, no polymorphism or HLA allele appeared significant at the genome-wide level.

CONCLUSIONS

Our results suggest that haplotype DQA1*01:01/DQB1*05:01/DRB1*01:01 is associated with susceptibility to HSP, but not with the severity of the kidney involvement. These HLA associations did not occur in IBD patients, suggesting that they are specific to HSP and not related to susceptibility to autoimmune diseases in general.

Coupled mixed model for joint genetic analysis of complex disorders with two independently collected data sets

BACKGROUND

In the last decade, Genome-wide Association studies (GWASs) have contributed to decoding the human genome by uncovering many genetic variations associated with various diseases. Many follow-up investigations involve joint analysis of multiple independently generated GWAS data sets. While most of the computational approaches developed for joint analysis are based on summary statistics, the joint analysis based on individual-level data with consideration of confounding factors remains to be a challenge.

RESULTS

In this study, we propose a method, called Coupled Mixed Model (CMM), that enables a joint GWAS analysis on two independently collected sets of GWAS data with different phenotypes. The CMM method does not require the data sets to have the same phenotypes as it aims to infer the unknown phenotypes using a set of multivariate sparse mixed models. Moreover, CMM addresses the confounding variables due to population stratification, family structures, and cryptic relatedness, as well as those arising during data collection such as batch effects that frequently appear in joint genetic studies. We evaluate the performance of CMM using simulation experiments. In real data analysis, we illustrate the utility of CMM by an application to evaluating common genetic associations for Alzheimer's disease and substance use disorder using datasets independently collected for the two complex human disorders. Comparison of the results with those from previous experiments and analyses supports the utility of our method and provides new insights into the diseases. The software is available at https://github.com/HaohanWang/CMM .

Exploring genetic targets of psoriasis using genome wide association studies (GWAS) for drug repurposing

Psoriasis is a chronic inflammatory disease causing itching in the body and pain in the joints. Currently, no permanent cure is available at a commercial level for this disease. Genome wide association studies (GWAS) provide a deeper insight that helps in better understanding this disease and further possible cure of this disease. The major goal of the present study is to identify potent genetic targets of psoriasis disease using GWAS approach and identify drugs for repurposing. The methods used include GWAS catalogue, GeneAnalytics, canSAR protein annotation tool, VarElect, Drug bank, Proteomics database, ProTox software. By exploring GWAS catalogue, 126 psoriasis associated genes (PAG) were identified. 68 genes found to be druggable were obtained from canSAR protein annotation tool. Localization results depict that maximum genes are present in cytoplasmic cellular components. The superpathways obtained from GeneAnalytics resulted in involvement of these genes in the immune system, Jak/Stat pathway, Th17 and Wnt pathways. Two genes Interleukin 13 (IL13) and POLI are Food and Drug Administration (FDA) approved targets. Small compounds for these targets were analysed for drug-likeliness, toxicity and mutagenecity properties. The FDA approved drug pandel was found to possess desirable properties. The medications used for psoriasis causes mild to severe side effects and does not work well always. Hence we propose drug repurposing strategy to use existing drugs for new therapies. Therefore, the drug pandel could be explored further and repurposed to treat psoriasis.

Whole genome prediction and heritability of childhood asthma phenotypes

Introduction

While whole genome prediction (WGP) methods have recently demonstrated successes in the prediction of complex genetic diseases, they have not yet been applied to asthma and related phenotypes. Longitudinal patterns of lung function differ between asthmatics, but these phenotypes have not been assessed for heritability or predictive ability. Herein, we assess the heritability and genetic predictability of asthma‐related phenotypes.

Methods

We applied several WGP methods to a well‐phenotyped cohort of 832 children with mild‐to‐moderate asthma from CAMP. We assessed narrow‐sense heritability and predictability for airway hyperresponsiveness, serum immunoglobulin E, blood eosinophil count, pre‐ and post‐bronchodilator forced expiratory volume in 1 sec (FEV₁), bronchodilator response, steroid responsiveness, and longitudinal patterns of lung function (normal growth, reduced growth, early decline, and their combinations). Prediction accuracy was evaluated using a training/testing set split of the cohort.

Results

We found that longitudinal lung function phenotypes demonstrated significant narrow‐sense heritability (reduced growth, 95%; normal growth with early decline, 55%). These same phenotypes also showed significant polygenic prediction (areas under the curve [AUCs] 56% to 62%). Including additional demographic covariates in the models increased prediction 4–8%, with reduced growth increasing from 62% to 66% AUC. We found that prediction with a genomic relatedness matrix was improved by filtering available SNPs based on chromatin evidence, and this result extended across cohorts.

Conclusions

Longitudinal reduced lung function growth displayed extremely high heritability. All phenotypes with significant heritability showed significant polygenic prediction. Using SNP‐prioritization increased prediction across cohorts. WGP methods show promise in predicting asthma‐related heritable traits.

Network biology concepts in complex disease comorbidities

The co-occurrence of diseases can inform the underlying network biology of shared and multifunctional genes and pathways. In addition, comorbidities help to elucidate the effects of external exposures, such as diet, lifestyle and patient care. With worldwide health transaction data now often being collected electronically, disease co-occurrences are starting to be quantitatively characterized. Linking network dynamics to the real-life, non-ideal patient in whom diseases co-occur and interact provides a valuable basis for generating hypotheses on molecular disease mechanisms, and provides knowledge that can facilitate drug repurposing and the development of targeted therapeutic strategies.

Harnessing publicly available genetic data to prioritize lipid modifying therapeutic targets for prevention of coronary heart disease based on dysglycemic risk

Therapeutic interventions that lower LDL-cholesterol effectively reduce the risk of coronary artery disease (CAD). However, statins, the most widely prescribed LDL-cholesterol lowering drugs, increase diabetes risk. We used genome-wide association study (GWAS) data in the public domain to investigate the relationship of LDL-C and diabetes and identify loci encoding potential drug targets for LDL-cholesterol modification without causing dysglycemia. We obtained summary-level GWAS data for LDL-C from GLGC, glycemic traits from MAGIC, diabetes from DIAGRAM and CAD from CARDIoGRAMplusC4D consortia. Mendelian randomization analyses identified a one standard deviation (SD) increase in LDL-C caused an increased risk of CAD (odds ratio [OR] 1.63 (95 % confidence interval [CI] 1.55, 1.71), which was not influenced by removing SNPs associated with diabetes. LDL-C/CAD-associated SNPs showed consistent effect directions (binomial P = 6.85 × 10⁻⁵). Conversely, a 1-SD increase in LDL-C was causally protective of diabetes (OR 0.86; 95 % CI 0.81, 0.91), however LDL-cholesterol/diabetes-associated SNPs did not show consistent effect directions (binomial P = 0.15). HMGCR, our positive control, associated with LDL-C, CAD and a glycemic composite (derived from GWAS meta-analysis of four glycemic traits and diabetes). In contrast, PCSK9, APOB, LPA, CETP, PLG, NPC1L1 and ALDH2 were identified as “druggable” loci that alter LDL-C and risk of CAD without displaying associations with dysglycemia. In conclusion, LDL-C increases the risk of CAD and the relationship is independent of any association of LDL-C with diabetes. Loci that encode targets of emerging LDL-C lowering drugs do not associate with dysglycemia, and this provides provisional evidence that new LDL-C lowering drugs (such as PCSK9 inhibitors) may not influence risk of diabetes.

A glutamatergic network mediates lithium response in bipolar disorder as defined by epigenome pathway analysis

AIM

A regulatory network in the human brain mediating lithium response in bipolar patients was revealed by analysis of functional SNPs from genome-wide association studies (GWAS) and published gene association studies, followed by epigenome mapping.

METHODS

An initial set of 23,312 SNPs in linkage disequilibrium with lead SNPs, and sub-threshold GWAS SNPs rescued by pathway analysis, were studied in the same populations. These were assessed using our workflow and annotation by the epigenome roadmap consortium.

RESULTS

Twenty-seven percent of 802 SNPs that were associated with lithium response (13 published studies gene association studies and two GWAS) were shared in common with 1281 SNPs from 18 GWAS examining psychiatric disorders and adverse events associated with lithium treatment. Nineteen SNPs were annotated as active regulatory elements such as enhancers and promoters in a tissue-specific manner. They were located within noncoding regions of ten genes: ANK3, ARNTL, CACNA1C, CACNG2, CDKN1A, CREB1, GRIA2, GSK3B, NR1D1 and SLC1A2. Following gene set enrichment and pathway analysis, these genes were found to be significantly associated (p = 10(-27); Fisher exact test) with an AMPA2 glutamate receptor network in human brain. Our workflow results showed concordance with annotation of regulatory elements from the epigenome roadmap. Analysis of cognate mRNA and enhancer RNA exhibited patterns consistent with an integrated pathway in human brain.

CONCLUSION

This pharmacoepigenomic regulatory pathway is located in the same brain regions that exhibit tissue volume loss in bipolar disorder. Although in silico analysis requires biological validation, the approach provides value for identification of candidate variants that may be used in pharmacogenomic testing to identify bipolar patients likely to respond to lithium.

Epigenomic mapping and effect sizes of noncoding variants associated with psychotropic drug response

AIM

To provide insight into potential regulatory mechanisms of gene expression underlying addiction, analgesia, psychotropic drug response and adverse drug events, genome-wide association studies searching for variants associated with these phenotypes has been undertaken with limited success. We undertook analysis of these results with the aim of applying epigenetic knowledge to aid variant discovery and interpretation.

METHODS

We applied conditional imputation to results from 26 genome-wide association studies and three candidate gene-association studies. The analysis workflow included data from chromatin conformation capture, chromatin state annotation, DNase I hypersensitivity, hypomethylation, anatomical localization and biochronicity. We also made use of chromatin state data from the epigenome roadmap, transcription factor-binding data, spatial maps from published Hi-C datasets and 'guilt by association' methods.

RESULTS

We identified 31 pharmacoepigenomic SNPs from a total of 2024 variants in linkage disequilibrium with lead SNPs, of which only 6% were coding variants. Interrogation of chromatin state using our workflow and the epigenome roadmap showed agreement on 34 of 35 tissue assignments to regulatory elements including enhancers and promoters. Loop boundary domains were inferred by association with CTCF (CCCTC-binding factor) and cohesin, suggesting proximity to topologically associating domain boundaries and enhancer clusters. Spatial interactions between enhancer-promoter pairs detected both known and previously unknown mechanisms. Addiction and analgesia SNPs were common in relevant populations and exhibited large effect sizes, whereas a SNP located in the promoter of the SLC1A2 gene exhibited a moderate effect size for lithium response in bipolar disorder in patients of European ancestry. SNPs associated with drug-induced organ injury were rare but exhibited the largest effect sizes, consistent with the published literature.

CONCLUSION

This work demonstrates that an in silico bioinformatics-based approach using integrative analysis of a diversity of molecular and morphological data types can discover pharmacoepigenomic variants that are suitable candidates for further validation in cell lines, animal models and human clinical trials.

The epigenome, 4D nucleome and next-generation neuropsychiatric pharmacogenomics

The 4D nucleome has the potential to render challenges in neuropsychiatric pharmacogenomics more tractable. The epigenome roadmap consortium has demonstrated the critical role that noncoding regions of the human genome play in determination of human phenotype. Chromosome conformation capture methods have revealed the 4D organization of the nucleus, bringing interactions between distant regulatory elements into close spatial proximity in a periodic manner. These functional interactions have the potential to elucidate mechanisms of CNS drug response and side effects that previously have been unrecognized. This perspective assesses recent advances likely to reveal novel pharmacodynamic regulatory pathways in human brain, charting a future new avenue of pharmacogenomics research, using the spatial and temporal architecture of the human epigenome as its foundation.