Genome-wide Association Studies Risk Variants Research Articles

Cross-Tissue Single-Nucleus RNA Sequencing Discovers Tissue-Resident Adipocytes Involved in Propanoate Metabolism in the Human Heart.

Adipocytes are crucial regulators of cardiovascular health. However, not much is known about gene expression profiles of adipocytes residing in nonfat cardiovascular tissues, their genetic regulation, and contribution to coronary artery disease. Here, we investigated whether and how the gene expression profiles of adipocytes in the subcutaneous adipose tissue differ from adipocytes residing in the heart. We used single-nucleus RNA-sequencing data sets of subcutaneous adipose tissue and heart and performed in-depth analysis of tissue-resident adipocytes and their cell-cell interactions. We first discovered tissue-specific features of tissue-resident adipocytes, identified functional pathways involved in their tissue specificity, and found genes with cell type-specific expression enrichment in tissue-resident adipocytes. By following up these results, we discovered the propanoate metabolism pathway as a novel distinct characteristic of the heart-resident adipocytes and found a significant enrichment of coronary artery disease genome-wide association study risk variants among the right atrium-specific adipocyte marker genes. Our cell-cell communication analysis identified 22 specific heart adipocyte-associated ligand-receptor pairs and signaling pathways, including THBS (thrombospondin) and EPHA (ephrin type-A), further supporting the distinct tissue-resident role of heart adipocytes. Our results also suggest chamber-level coordination of heart adipocyte expression profiles as we observed a consistently larger number of adipocyte-associated ligand-receptor interactions and functional pathways in the atriums than ventricles. Overall, we introduce a new function and genetic link to coronary artery disease for the previously unexplored heart-resident adipocytes.

Splicing annotation of endometrial cancer GWAS risk loci reveals potentially causal variants and supports a role for NF1 and SKAP1 as susceptibility genes.

Alternative splicing contributes to cancer development. Indeed, splicing analysis of cancer genome-wide association study (GWAS) risk variants has revealed likely causal variants. To systematically assess GWAS variants for splicing effects, we developed a prioritization workflow using a combination of splicing prediction tools, alternative transcript isoforms, and splicing quantitative trait locus (sQTL) annotations. Application of this workflow to candidate causal variants from 16 endometrial cancer GWAS risk loci highlighted single-nucleotide polymorphisms (SNPs) that were predicted to upregulate alternative transcripts. For two variants, sQTL data supported the predicted impact on splicing. At the 17q11.2 locus, the protective allele for rs7502834 was associated with increased splicing of an exon in a NF1 alternative transcript encoding a truncated protein in adipose tissue and is consistent with an endometrial cancer transcriptome-wide association study (TWAS) finding in adipose tissue. Notably, NF1 haploinsufficiency is protective for obesity, a well-established risk factor for endometrial cancer. At the 17q21.32 locus, the rs2278868 risk allele was predicted to upregulate a SKAP1 transcript that is subject to nonsense-mediated decay, concordant with a corresponding sQTL in lymphocytes. This is consistent with a TWAS finding that indicates decreased SKAP1 expression in blood increases endometrial cancer risk. As SKAP1 is involved in Tcell immune responses, decreased SKAP1 expression may impact endometrial tumor immunosurveillance. In summary, our analysis has identified potentially causal endometrial cancer GWAS risk variants with plausible biological mechanisms and provides a splicing annotation workflow to aid interpretation of other GWAS datasets.

Transcriptomics insights into interpreting AMD-GWAS discoveries for biological and clinical applications

Genome-wide association studies (GWAS) have been successful in identifying genetic risk factors for a large number of complex diseases, including age-related macular degeneration (AMD), which is a highly heritable complex disease affecting millions of elderly individuals. However, the progress of elucidating the functional relevance of genetic findings in AMD has been slow, as most risk factors are non-coding, and we have little insight into the causal genes and disease mechanisms. In the last few years, gene expression regulation is emerging as a dominant mechanism through which GWAS risk variants lead to the disease. The purpose of this review is to provide an overview of how transcriptome studies can help in identifying the genes, pathways and therapeutic targets underlying GWAS discoveries in AMD. These approaches help pave the road for mechanistic understanding of GWAS findings and drive translational advances that will lead to improved AMD management and treatment.

Allele-specific epigenetic activity in prostate cancer and normal prostate tissue implicates prostate cancer risk mechanisms

Genome-wide association studies (GWASs) of prostate cancer have identified >250 significant risk loci, but the causal variants and mechanisms for these loci remain largely unknown. Here, we sought to identify and characterize risk-harboring regulatory elements by integrating epigenomes from primary prostate tumor and normal tissues of 27 individuals across the H3K27ac, H3K4me3, and H3K4me2 histone marks and FOXA1 and HOXB13 transcription factors. We identified 7,371 peaks with significant allele specificity (allele-specific quantitative trait locus [asQTL] peaks). Showcasing their relevance to prostate cancer risk, H3K27ac T-asQTL peaks were the single annotation most enriched for prostate cancer GWAS heritability (40×), significantly higher than corresponding non-asQTL H3K27ac peaks (14×) or coding regions (14×). Surprisingly, fine-mapped GWAS risk variants were most significantly enriched for asQTL peaks observed in tumors, including asQTL peaks that were differentially imbalanced with respect to tumor-normal states. These data pinpointed putative causal regulatory elements at 20 GWAS loci, of which 11 were detected only in the tumor samples. More broadly, tumor-specific asQTLs were enriched for expression QTLs in benign tissues as well as accessible regions found in stem cells, supporting a hypothesis where some germline variants become reactivated during or after transformation and can be captured by epigenomic profiling of the tumor. Our study demonstrates the power of allele specificity in chromatin signals to uncover GWAS mechanisms, highlights the relevance of tumor-specific regulation in the context of cancer risk, and prioritizes multiple loci for experimental follow-up.

A comprehensive integrated post-GWAS analysis of Type 1 diabetes reveals enhancer-based immune dysregulation.

Type 1 diabetes (T1D) is an organ-specific autoimmune disease, whereby immune cell-mediated killing leads to loss of the insulin-producing β cells in the pancreas. Genome-wide association studies (GWAS) have identified over 200 genetic variants associated with risk for T1D. The majority of the GWAS risk variants reside in the non-coding regions of the genome, suggesting that gene regulatory changes substantially contribute to T1D. However, identification of causal regulatory variants associated with T1D risk and their affected genes is challenging due to incomplete knowledge of non-coding regulatory elements and the cellular states and processes in which they function. Here, we performed a comprehensive integrated post-GWAS analysis of T1D to identify functional regulatory variants in enhancers and their cognate target genes. Starting with 1,817 candidate T1D SNPs defined from the GWAS catalog and LDlink databases, we conducted functional annotation analysis using genomic data from various public databases. These include 1) Roadmap Epigenomics, ENCODE, and RegulomeDB for epigenome data; 2) GTEx for tissue-specific gene expression and expression quantitative trait loci data; and 3) lncRNASNP2 for long non-coding RNA data. Our results indicated a prevalent enhancer-based immune dysregulation in T1D pathogenesis. We identified 26 high-probability causal enhancer SNPs associated with T1D, and 64 predicted target genes. The majority of the target genes play major roles in antigen presentation and immune response and are regulated through complex transcriptional regulatory circuits, including those in HLA (6p21) and non-HLA (16p11.2) loci. These candidate causal enhancer SNPs are supported by strong evidence and warrant functional follow-up studies.

Hearing difficulty is linked to Alzheimer\u2019s disease by common genetic vulnerability, not shared genetic architecture

Age-related hearing loss was recently established as the largest modifiable risk factor for Alzheimer’s disease (AD), however, the reasons for this link remain unclear. We investigate shared underlying genetic associations using results from recent large genome-wide association studies (GWAS) on adult hearing difficulty and AD. Genetic correlation and Mendelian randomization (MR) analysis do not support a genetic correlation between the disorders, but suggest a direct causal link from AD genetic risk to hearing difficulty, driven by APOE. Systematic MR analyses on the effect of other traits revealed shared effects of glutamine, gamma-glutamylglutamine, and citrate levels on reduced risk of both hearing difficulty and AD. In addition, pathway analysis on GWAS risk variants suggests shared function in neuronal signalling pathways as well as etiology of diabetes and cardiovascular disease. However, after multiple testing corrections, neither analysis led to statistically significant associations. Altogether, our genetic-driven analysis suggests hearing difficulty and AD are linked by a shared vulnerability in molecular pathways rather than by a shared genetic architecture.

Sex-specific nicotine sensitization and imprinting of self-administration in rats inform GWAS findings on human addiction phenotypes

Repeated nicotine exposure leads to sensitization (SST) and enhances self-administration (SA) in rodents. However, the molecular basis of nicotine SST and SA and their biological relevance to the mounting genome-wide association study (GWAS) loci of human addictive behaviors are poorly understood. Considering a gateway drug role of nicotine, we modeled nicotine SST and SA in F1 progeny of inbred rats (F344/BN) and conducted integrative genomics analyses. We unexpectedly observed male-specific nicotine SST and a parental effect of SA only present in paternal F344 crosses. Transcriptional profiling in the ventral tegmental area (VTA) and nucleus accumbens (NAc) core and shell further revealed sex- and brain region-specific transcriptomic signatures of SST and SA. We found that genes associated with SST and SA were enriched for those related to synaptic processes, myelin sheath, and tobacco use disorder or chemdependency. Interestingly, SST-associated genes were often downregulated in male VTA but upregulated in female VTA, and strongly enriched for smoking GWAS risk variants, possibly explaining the male-specific SST. For SA, we found widespread region-specific allelic imbalance of expression (AIE), of which genes showing AIE bias toward paternal F344 alleles in NAc core were strongly enriched for SA-associated genes and for GWAS risk variants of smoking initiation, likely contributing to the parental effect of SA. Our study suggests a mechanistic link between transcriptional changes underlying the NIC SST and SA and human nicotine addiction, providing a resource for understanding the neurobiology basis of the GWAS findings on human smoking and other addictive phenotypes.

MCF-7 as a Model for Functional Analysis of Breast Cancer Risk Variants.

Breast cancer genetic predisposition is governed by more than 142 loci as revealed by genome-wide association studies (GWAS). The functional contribution of these risk loci to breast cancer remains unclear, and additional post-GWAS analyses are required. We identified active regulatory elements (enhancers, promoters, and chromatin organizing elements) by histone H3K27 acetylation and CTCF occupancy and determined the enrichment of risk variants at these sites. We compared these results with previously published data and for other cell lines, including human mammary epithelial cells, and related these data to gene expression. In terms of mapping accuracy and resolution, our data augment previous annotations of the MCF-7 epigenome. After intersection with GWAS risk variants, we found 39 enhancers and 15 CTCF occupancy sites that, between them, overlapped 96 breast cancer credible risk variants at 42 loci. These risk enhancers likely regulate the expression of dozens of genes, which are enriched for GO categories, including estrogen and prolactin signaling. Ten (of 142) breast cancer risk loci likely function via enhancers that are active in MCF-7 and are well suited to targeted manipulation in this system. In contrast, risk loci cannot be mapped to specific CTCF-binding sites, and the genes linked to risk CTCF sites did not show functional enrichment. The identity of risk enhancers and their associated genes suggests that some risk may function during later processes in cancer progression. Here, we report how the ER+ cell line MCF-7 can be used to dissect risk mechanisms for breast cancer.

Transcriptome-wide association study of multiple myeloma identifies candidate susceptibility genes

BackgroundWhile genome-wide association studies (GWAS) of multiple myeloma (MM) have identified variants at 23 regions influencing risk, the genes underlying these associations are largely unknown. To identify candidate causal genes at these regions and search for novel risk regions, we performed a multi-tissue transcriptome-wide association study (TWAS).ResultsGWAS data on 7319 MM cases and 234,385 controls was integrated with Genotype-Tissue Expression Project (GTEx) data assayed in 48 tissues (sample sizes, N = 80–491), including lymphocyte cell lines and whole blood, to predict gene expression. We identified 108 genes at 13 independent regions associated with MM risk, all of which were in 1 Mb of known MM GWAS risk variants. Of these, 94 genes, located in eight regions, had not previously been considered as a candidate gene for that locus.ConclusionsOur findings highlight the value of leveraging expression data from multiple tissues to identify candidate genes responsible for GWAS associations which provide insight into MM tumorigenesis. Among the genes identified, a number have plausible roles in MM biology, notably APOBEC3C, APOBEC3H, APOBEC3D, APOBEC3F, APOBEC3G, or have been previously implicated in other malignancies. The genes identified in this TWAS can be explored for follow-up and validation to further understand their role in MM biology.

SP140 regulates the expression of immune-related genes associated with multiple sclerosis and other autoimmune diseases by NF-κB inhibition

SP140 locus has been associated with multiple sclerosis (MS) as well as other autoimmune diseases by genome-wide association studies (GWAS). The causal variant of these associations (rs28445040-T) alters the splicing of the SP140 gene transcripts reducing the protein expression. We aimed to understand why the reduction of SP140 expression produced by the risk variant can increase the susceptibility to MS. To this end, we determined by RNA sequencing (RNA-seq) analysis the differentially expressed genes after SP140 silencing in lymphoblastoid cell lines (LCLs). We analyzed these genes by gene ontology (GO), comparative transcriptome profiles, enrichment of transcription factors (TFs) in the promoters of these genes and colocalization with GWAS risk variants. We also monitored the activity of the nuclear factor kappa-light-chain-enhancer of activated B cells (NF-κB) in SP140-silenced cells by luciferase reporter system. We identified 100 genes that were up-regulated and 22 genes down-regulated in SP140-silenced LCLs. GO analysis revealed that genes affected by SP140 were involved in regulation of cytokine production, inflammatory response and cell-cell adhesion. We observed enrichment of NF-κB TF in the promoter of up-regulated genes and NF-κB-increased activity in SP140-silenced cell lines. We showed enrichment of genes regulated by SP140 in GWAS-detected risk loci for MS (14.63 folds), Crohn's disease (4.82 folds) and inflammatory bowel disease (4.47 folds), not observed in other unrelated immune diseases. Our findings showed that SP140 is an important repressor of genes implicated in inflammation, suggesting that decreased expression of SP140, promoted by the rs28445040-T risk variant, may lead to up-regulation of these genes by means of NF-κB inhibition in B cells.

Abstract 1221: Chronic lymphocytic leukemia risk variants are associated with alteration of epigenetic states in regulatory regions

Abstract Chronic lymphocytic leukemia risk variants are associated with alteration of epigenetic state in regulatory regions The vast majority of disease-associated single-nucleotide polymorphisms (SNPs) identified by genome-wide association studies (GWASs) reside in noncoding genomic regions. Functional impact of noncoding casual variants on target gene expression often involves regulatory mechanisms, typically through altering local chromatin accessibility, histone modifications, transcription factor (TF) binding affinity and chromatin interactions. In chronic lymphocytic leukemia (CLL), we and others have identified 41 GWAS risk variants, whose mechanism of action in CLL etiology remains largely unknown. In this study, we aim to decipher the regulatory potential of these variants and to better understand the mechanisms contributing to dysregulation of genes implicated in CLL pathogenesis. We performed ChIP-seq for H3K4me3, H3K4me1, H3K27ac and H3K27me3 in CLL. In addition, we also analyzed over 250 public ChIP-seq, ATAC-seq, DNase-seq, FAIRE-seq, and chromatin interaction data in CLL, B cells and/or lymphoblastoid cell lines. Comparing the 41 index SNPs and the nearby SNPs they associated with in linkage disequilibrium (R2 at least 0.5 in the EUR ethnic group in 1000 Genomes Project) to the above epigenetic data sets revealed overlaps with regulatory regions in 38 of the risk loci, with 28 of the loci overlapping super-enhancers. We used two approaches to elucidate the mechanisms of potential SNP causality. First, we used GATK UnifiedGenotyper to infer genotype for the SNPs that overlapped open chromatin regions and ChIP-seq peaks, by which evidence of imbalance between reference and alternative allele was identified at heterozygous sites. Second, for SNPs located within open chromatin regions or ChIP-seq peaks, we used FIMO program to scan the 100bp sequences centered on the SNPs for significant matches (p=1e-05) to the position weight matrices in five TF binding motif databases; SNPs that had TF motifs altered by the alternative alleles were retained. Together, these two approaches revealed evidence for allele imbalance or TF binding motif alterations at 36 of the 38 CLL risk loci, suggesting the generality of epigenetic regulation in CLL pathogenesis. Further analysis of chromatin interaction, gene expression and expression quantitative trait loci data will likely link the casual SNPs to their target genes. Citation Format: Huihuang Yan, Shulan Tian, Geffen Kleinstern, Susan Slager. Chronic lymphocytic leukemia risk variants are associated with alteration of epigenetic states in regulatory regions [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2018; 2018 Apr 14-18; Chicago, IL. Philadelphia (PA): AACR; Cancer Res 2018;78(13 Suppl):Abstract nr 1221.

Integration of Enhancer-Promoter Interactions with GWAS Summary Results Identifies Novel Schizophrenia-Associated Genes and Pathways.

It remains challenging to boost statistical power of genome-wide association studies (GWASs) to identify more risk variants or loci that can account for "missing heritability." Furthermore, since most identified variants are not in gene-coding regions, a biological interpretation of their function is largely lacking. On the other hand, recent biotechnological advances have made it feasible to experimentally measure the three-dimensional organization of the genome, including enhancer-promoter interactions in high resolutions. Due to the well-known critical roles of enhancer-promoter interactions in regulating gene expression programs, such data have been applied to link GWAS risk variants to their putative target genes, gaining insights into underlying biological mechanisms. However, their direct use in GWAS association testing is yet to be exploited. Here we propose integrating enhancer-promoter interactions into GWAS association analysis to both boost statistical power and enhance interpretability. We demonstrate that through an application to two large-scale schizophrenia (SCZ) GWAS summary data sets, the proposed method could identify some novel SCZ-associated genes and pathways (containing no significant SNPs). For example, after the Bonferroni correction, for the larger SCZ data set with 36,989 cases and 113,075 controls, our method applied to the gene body and enhancer regions identified 27 novel genes and 11 novel KEGG pathways to be significant, all missed by the transcriptome-wide association study (TWAS) approach. We conclude that our proposed method is potentially useful and is complementary to TWAS and other standard gene- and pathway-based methods.

P.1.a.016 Investigating schizophrenia genome-wide association studies risk variants in cognitive and neuroimaging defects in the GENUS Consortium collection

P.1.a.016 Investigating schizophrenia genome-wide association studies risk variants in cognitive and neuroimaging defects in the GENUS Consortium collection

Top associated SNPs in prostate cancer are significantly enriched in cis-expression quantitative trait loci and at transcription factor binding sites

While genome-wide association studies (GWAS) have revealed thousands of disease risk single nucleotide polymorphisms (SNPs), their functions remain largely unknown. Recent studies have suggested the regulatory roles of GWAS risk variants in several common diseases; however, the complex regulatory structure in prostate cancer is unclear. We investigated the potential regulatory roles of risk variants in two prostate cancer GWAS datasets by their interactions with expression quantitative trait loci (eQTL) and/or transcription factor binding sites (TFBSs) in three populations. Our results indicated that the moderately associated GWAS SNPs were significantly enriched with cis-eQTLs and TFBSs in Caucasians (CEU), but not in African Americans (AA) or Japanese (JPT); this was also observed in an independent pan-cancer related SNPs from the GWAS Catalog. We found that the eQTL enrichment in the CEU population was tissue-specific to eQTLs from CEU lymphoblastoid cell lines. Importantly, we pinpointed two SNPs, rs2861405 and rs4766642, by overlapping results from cis-eQTL and TFBS as applied to the CEU data. These results suggested that prostate cancer associated SNPs and pan-cancer associated SNPs are likely to play regulatory roles in CEU. However, the negative enrichment results in AA or JPT and the potential mechanisms remain to be elucidated in additional samples.

Genome wide association identified colorectal cancer susceptibility loci and colorectal cancer risk in Lynch syndrome

Background Recent genome-wide association studies (GWAS) have identified several common low-risk variants for colorectal cancer. Although some of the GWAS were enriched for young onset and family history positive colorectal cancer cases, it is not clear if these variants modify colorectal cancer risk for people with Lynch syndrome. In a case-unaffected sibling analysis of population and clinic based sibships, a study by Poynter et al. (2007) found a stronger association for risk variants in the 8q24 region (rs10505477 and rs6983627) with microsatellite instability (MSI)-high tumors. MSI-high tumors are characteristic of Lynch syndrome. Therefore we hypothesize that these variants may also influence colorectal cancer risk in Lynch syndrome. Recently Wijnen et al. (2009) have reported on some of the GWAS risk variants in Dutch Lynch syndrome families and although they did not find an association between rs6983627 and colorectal cancer risk, other risk variants in 8q24 and 9p24 have yet to be characterized in Lynch syndrome.