Hassan Foroughi Asl

Prediction of Causal Candidate Genes in Coronary Artery Disease Loci

Objective—Genome-wide association studies have to date identified 159 significant and suggestive loci for coronary artery disease (CAD). We now report comprehensive bioinformatics analyses of sequence variation in these loci to predict candidate causal genes. Approach and Results—All annotated genes in the loci were evaluated with respect to protein-coding single-nucleotide polymorphism and gene expression parameters. The latter included expression quantitative trait loci, tissue specificity, and miRNA binding. High priority candidate genes were further identified based on literature searches and our experimental data. We conclude that the great majority of causal variations affecting CAD risk occur in noncoding regions, with 41% affecting gene expression robustly versus 6% leading to amino acid changes. Many of these genes differed from the traditionally annotated genes, which was usually based on proximity to the lead single-nucleotide polymorphism. Indeed, we obtained evidence that genetic variants at CAD loci affect 98 genes which had not been linked to CAD previously. Conclusions—Our results substantially revise the list of likely candidates for CAD and suggest that genome-wide association studies efforts in other diseases may benefit from similar bioinformatics analyses.

Cross-Tissue Regulatory Gene Networks in Coronary Artery Disease

SUMMARY Inferring molecular networks can reveal how genetic perturbations interact with environmental factors to cause common complex diseases. We analyzed genetic and gene expression data from seven tissues relevant to coronary artery disease (CAD) and identified regulatory gene networks (RGNs) and their key drivers. By integrating data from genome-wide association studies, we identified 30 CAD-causal RGNs interconnected in vascular and metabolic tissues, and we validated them with corresponding data from the Hybrid Mouse Diversity Panel. As proof of concept, by targeting the key drivers AIP, DRAP1, POLR2I, and PQBP1 in a cross-species-validated, arterial-wall RGN involving RNA-processing genes, we re-identified this RGN in THP-1 foam cells and independent data from CAD macrophages and carotid lesions. This characterization of the molecular landscape in CAD will help better define the regulation of CAD candidate genes identified by genome-wide association studies and is a first step toward achieving the goals of precision medicine.

Plasma cholesterol-induced lesion networks activated before regression of early, mature, and advanced atherosclerosis.

Plasma cholesterol lowering (PCL) slows and sometimes prevents progression of atherosclerosis and may even lead to regression. Little is known about how molecular processes in the atherosclerotic arterial wall respond to PCL and modify responses to atherosclerosis regression. We studied atherosclerosis regression and global gene expression responses to PCL (≥80%) and to atherosclerosis regression itself in early, mature, and advanced lesions. In atherosclerotic aortic wall from Ldlr−/−Apob 100/100 Mttp flox/floxMx1-Cre mice, atherosclerosis regressed after PCL regardless of lesion stage. However, near-complete regression was observed only in mice with early lesions; mice with mature and advanced lesions were left with regression-resistant, relatively unstable plaque remnants. Atherosclerosis genes responding to PCL before regression, unlike those responding to the regression itself, were enriched in inherited risk for coronary artery disease and myocardial infarction, indicating causality. Inference of transcription factor (TF) regulatory networks of these PCL-responsive gene sets revealed largely different networks in early, mature, and advanced lesions. In early lesions, PPARG was identified as a specific master regulator of the PCL-responsive atherosclerosis TF-regulatory network, whereas in mature and advanced lesions, the specific master regulators were MLL5 and SRSF10/XRN2, respectively. In a THP-1 foam cell model of atherosclerosis regression, siRNA targeting of these master regulators activated the time-point-specific TF-regulatory networks and altered the accumulation of cholesterol esters. We conclude that PCL leads to complete atherosclerosis regression only in mice with early lesions. Identified master regulators and related PCL-responsive TF-regulatory networks will be interesting targets to enhance PCL-mediated regression of mature and advanced atherosclerotic lesions.

kruX: matrix-based non-parametric eQTL discovery

BackgroundThe Kruskal-Wallis test is a popular non-parametric statistical test for identifying expression quantitative trait loci (eQTLs) from genome-wide data due to its robustness against variations in the underlying genetic model and expression trait distribution, but testing billions of marker-trait combinations one-by-one can become computationally prohibitive.ResultsWe developed kruX, an algorithm implemented in Matlab, Python and R that uses matrix multiplications to simultaneously calculate the Kruskal-Wallis test statistic for several millions of marker-trait combinations at once. KruX is more than ten thousand times faster than computing associations one-by-one on a typical human dataset. We used kruX and a dataset of more than 500k SNPs and 20k expression traits measured in 102 human blood samples to compare eQTLs detected by the Kruskal-Wallis test to eQTLs detected by the parametric ANOVA and linear model methods. We found that the Kruskal-Wallis test is more robust against data outliers and heterogeneous genotype group sizes and detects a higher proportion of non-linear associations, but is more conservative for calling additive linear associations.ConclusionkruX enables the use of robust non-parametric methods for massive eQTL mapping without the need for a high-performance computing infrastructure and is freely available from http://krux.googlecode.com.

Expression quantitative trait Loci acting across multiple tissues are enriched in inherited risk for coronary artery disease.

Background—Despite recent discoveries of new genetic risk factors, the majority of risk for coronary artery disease (CAD) remains elusive. As the most proximal sensor of DNA variation, RNA abundance can help identify subpopulations of genetic variants active in and across tissues mediating CAD risk through gene expression. Methods and Results—By generating new genomic data on DNA and RNA samples from the Stockholm Atherosclerosis Gene Expression (STAGE) study, 8156 cis-acting expression quantitative trait loci (eQTLs) for 6450 genes across 7 CAD-relevant tissues were detected. The inherited risk enrichments of tissue-defined sets of these eQTLs were assessed using 2 independent genome-wide association data sets. eQTLs acting across increasing numbers of tissues were found increasingly enriched for CAD risk and resided at regulatory hot spots. The risk enrichment of 42 eQTLs acting across 5 to 6 tissues was particularly high (⩽7.3-fold) and confirmed in the combined genome-wide association data from Coronary Artery Disease Genome Wide Replication And Meta-Analysis Consortium. Sixteen of the 42 eQTLs associated with 19 master regulatory genes and 29 downstream gene sets (n>30) were further risk enriched comparable to that of the 153 genome-wide association risk single-nucleotide polymorphisms established for CAD (8.4-fold versus 10-fold). Three gene sets, governed by the master regulators FLYWCH1, PSORSIC3, and G3BP1, segregated the STAGE patients according to extent of CAD, and small interfering RNA targeting of these master regulators affected cholesterol-ester accumulation in foam cells of the THP1 monocytic cell line. Conclusions—eQTLs acting across multiple tissues are significant carriers of inherited risk for CAD. FLYWCH1, PSORSIC3, and G3BP1 are novel master regulatory genes in CAD that may be suitable targets.

Human Validation of Genes Associated With a Murine Atherosclerotic Phenotype

Objective— The genetically modified mouse is the most commonly used animal model for studying the pathogenesis of atherosclerotic disease. We aimed to assess if mice atherosclerosis-related genes could be validated in human disease through examination of results from genome-wide association studies. Approach and Results— We performed a systematic review to identify atherosclerosis-causing genes in mice and carried out gene-based association tests of their human orthologs for an association with human coronary artery disease and human large artery ischemic stroke. Moreover, we investigated the association of these genes with human atherosclerotic plaque characteristics. In addition, we assessed the presence of tissue-specific cis-acting expression quantitative trait loci for these genes in humans. Finally, using pathway analyses we show that the putative atherosclerosis-causing genes revealed few associations with human coronary artery disease, large artery ischemic stroke, or atherosclerotic plaque characteristics, despite the fact that the majority of these genes have cis-acting expression quantitative trait loci. Conclusions— A role for genes that has been observed in mice for atherosclerotic lesion development could scarcely be confirmed by studying associations of disease development with common human genetic variants. The value of murine atherosclerotic models for selection of therapeutic targets in human disease remains unclear.

Systematic analysis of chromatin interactions at disease associated loci links novel candidate genes to inflammatory bowel disease

BackgroundGenome-wide association studies (GWAS) have revealed many susceptibility loci for complex genetic diseases. For most loci, the causal genes have not been identified. Currently, the identification of candidate genes is predominantly based on genes that localize close to or within identified loci. We have recently shown that 92 of the 163 inflammatory bowel disease (IBD)-loci co-localize with non-coding DNA regulatory elements (DREs). Mutations in DREs can contribute to IBD pathogenesis through dysregulation of gene expression. Consequently, genes that are regulated by these 92 DREs are to be considered as candidate genes. This study uses circular chromosome conformation capture-sequencing (4C-seq) to systematically analyze chromatin-interactions at IBD susceptibility loci that localize to regulatory DNA.ResultsUsing 4C-seq, we identify genomic regions that physically interact with the 92 DRE that were found at IBD susceptibility loci. Since the activity of regulatory elements is cell-type specific, 4C-seq was performed in monocytes, lymphocytes, and intestinal epithelial cells. Altogether, we identified 902 novel IBD candidate genes. These include genes specific for IBD-subtypes and many noteworthy genes including ATG9A and IL10RA. We show that expression of many novel candidate genes is genotype-dependent and that these genes are upregulated during intestinal inflammation in IBD. Furthermore, we identify HNF4α as a potential key upstream regulator of IBD candidate genes.ConclusionsWe reveal many novel and relevant IBD candidate genes, pathways, and regulators. Our approach complements classical candidate gene identification, links novel genes to IBD and can be applied to any existing GWAS data.

Lim Domain Binding 2 A Key Driver of Transendothelial Migration of Leukocytes and Atherosclerosis

Objective— Using a multi-tissue, genome-wide gene expression approach, we recently identified a gene module linked to the extent of human atherosclerosis. This atherosclerosis module was enriched with inherited risk for coronary and carotid artery disease (CAD) and overlapped with genes in the transendothelial migration of leukocyte (TEML) pathway. Among the atherosclerosis module genes, the transcription cofactor Lim domain binding 2 (LDB2) was the most connected in a CAD vascular wall regulatory gene network. Here, we used human genomics and atherosclerosis-prone mice to evaluate the possible role of LDB2 in TEML and atherosclerosis. Approach and Results— mRNA profiles generated from blood macrophages in patients with CAD were used to infer transcription factor regulatory gene networks; Ldlr –/– Apob 100/100 mice were used to study the effects of Ldb2 deficiency on TEML activity and atherogenesis. LDB2 was the most connected gene in a transcription factor regulatory network inferred from TEML and atherosclerosis module genes in CAD macrophages. In Ldlr –/– Apob 100/100 mice, loss of Ldb2 increased atherosclerotic lesion size ≈2-fold and decreased plaque stability. The exacerbated atherosclerosis was caused by increased TEML activity, as demonstrated in air-pouch and retinal vasculature models in vivo, by ex vivo perfusion of primary leukocytes, and by leukocyte migration in vitro. In THP1 cells, migration was increased by overexpression and decreased by small interfering RNA inhibition of LDB2. A functional LDB2 variant (rs10939673) was associated with the risk and extent of CAD across several cohorts. Conclusions— As a key driver of the TEML pathway in CAD macrophages, LDB2 is a novel candidate to target CAD by inhibiting the overall activity of TEML.Objective— Using a multi-tissue, genome-wide gene expression approach, we recently identified a gene module linked to the extent of human atherosclerosis. This atherosclerosis module was enriched with inherited risk for coronary and carotid artery disease (CAD) and overlapped with genes in the transendothelial migration of leukocyte (TEML) pathway. Among the atherosclerosis module genes, the transcription cofactor Lim domain binding 2 (LDB2) was the most connected in a CAD vascular wall regulatory gene network. Here, we used human genomics and atherosclerosis-prone mice to evaluate the possible role of LDB2 in TEML and atherosclerosis. Approach and Results— mRNA profiles generated from blood macrophages in patients with CAD were used to infer transcription factor regulatory gene networks; Ldlr –/– Apob 100/100 mice were used to study the effects of Ldb2 deficiency on TEML activity and atherogenesis. LDB2 was the most connected gene in a transcription factor regulatory network inferred from TEML and atherosclerosis module genes in CAD macrophages. In Ldlr –/– Apob 100/100 mice, loss of Ldb2 increased atherosclerotic lesion size ≈2-fold and decreased plaque stability. The exacerbated atherosclerosis was caused by increased TEML activity, as demonstrated in air-pouch and retinal vasculature models in vivo, by ex vivo perfusion of primary leukocytes, and by leukocyte migration in vitro. In THP1 cells, migration was increased by overexpression and decreased by small interfering RNA inhibition of LDB2 . A functional LDB2 variant (rs10939673) was associated with the risk and extent of CAD across several cohorts. Conclusions— As a key driver of the TEML pathway in CAD macrophages, LDB2 is a novel candidate to target CAD by inhibiting the overall activity of TEML. # Significance {#article-title-34}

Network analysis of coronary artery disease risk genes elucidates disease mechanisms and druggable targets

Genome-wide association studies (GWAS) have identified over two hundred chromosomal loci that modulate risk of coronary artery disease (CAD). The genes affected by variants at these loci are largely unknown and an untapped resource to improve our understanding of CAD pathophysiology and identify potential therapeutic targets. Here, we prioritized 68 genes as the most likely causal genes at genome-wide significant loci identified by GWAS of CAD and examined their regulatory roles in 286 metabolic and vascular tissue gene-protein sub-networks (“modules”). The modules and genes within were scored for CAD druggability potential. The scoring enriched for targets of cardiometabolic drugs currently in clinical use and in-depth analysis of the top-scoring modules validated established and revealed novel target tissues, biological processes, and druggable targets. This study provides an unprecedented resource of tissue-defined gene–protein interactions directly affected by genetic variance in CAD risk loci.

Lim Domain Binding 2

Objective— Using a multi-tissue, genome-wide gene expression approach, we recently identified a gene module linked to the extent of human atherosclerosis. This atherosclerosis module was enriched with inherited risk for coronary and carotid artery disease (CAD) and overlapped with genes in the transendothelial migration of leukocyte (TEML) pathway. Among the atherosclerosis module genes, the transcription cofactor Lim domain binding 2 (LDB2) was the most connected in a CAD vascular wall regulatory gene network. Here, we used human genomics and atherosclerosis-prone mice to evaluate the possible role of LDB2 in TEML and atherosclerosis. Approach and Results— mRNA profiles generated from blood macrophages in patients with CAD were used to infer transcription factor regulatory gene networks; Ldlr –/– Apob 100/100 mice were used to study the effects of Ldb2 deficiency on TEML activity and atherogenesis. LDB2 was the most connected gene in a transcription factor regulatory network inferred from TEML and atherosclerosis module genes in CAD macrophages. In Ldlr –/– Apob 100/100 mice, loss of Ldb2 increased atherosclerotic lesion size ≈2-fold and decreased plaque stability. The exacerbated atherosclerosis was caused by increased TEML activity, as demonstrated in air-pouch and retinal vasculature models in vivo, by ex vivo perfusion of primary leukocytes, and by leukocyte migration in vitro. In THP1 cells, migration was increased by overexpression and decreased by small interfering RNA inhibition of LDB2. A functional LDB2 variant (rs10939673) was associated with the risk and extent of CAD across several cohorts. Conclusions— As a key driver of the TEML pathway in CAD macrophages, LDB2 is a novel candidate to target CAD by inhibiting the overall activity of TEML.Objective— Using a multi-tissue, genome-wide gene expression approach, we recently identified a gene module linked to the extent of human atherosclerosis. This atherosclerosis module was enriched with inherited risk for coronary and carotid artery disease (CAD) and overlapped with genes in the transendothelial migration of leukocyte (TEML) pathway. Among the atherosclerosis module genes, the transcription cofactor Lim domain binding 2 (LDB2) was the most connected in a CAD vascular wall regulatory gene network. Here, we used human genomics and atherosclerosis-prone mice to evaluate the possible role of LDB2 in TEML and atherosclerosis. Approach and Results— mRNA profiles generated from blood macrophages in patients with CAD were used to infer transcription factor regulatory gene networks; Ldlr –/– Apob 100/100 mice were used to study the effects of Ldb2 deficiency on TEML activity and atherogenesis. LDB2 was the most connected gene in a transcription factor regulatory network inferred from TEML and atherosclerosis module genes in CAD macrophages. In Ldlr –/– Apob 100/100 mice, loss of Ldb2 increased atherosclerotic lesion size ≈2-fold and decreased plaque stability. The exacerbated atherosclerosis was caused by increased TEML activity, as demonstrated in air-pouch and retinal vasculature models in vivo, by ex vivo perfusion of primary leukocytes, and by leukocyte migration in vitro. In THP1 cells, migration was increased by overexpression and decreased by small interfering RNA inhibition of LDB2 . A functional LDB2 variant (rs10939673) was associated with the risk and extent of CAD across several cohorts. Conclusions— As a key driver of the TEML pathway in CAD macrophages, LDB2 is a novel candidate to target CAD by inhibiting the overall activity of TEML. # Significance {#article-title-34}