Allele-specific Enhancer Activity Research Articles

A functional variant rs912304 for late-onset T1D risk contributes to islet dysfunction by regulating proinflammatory cytokine-responsive gene STXBP6 expression

BackgroundOur previous genome‑wide association studies (GWAS) have suggested rs912304 in 14q12 as a suggestive risk variant for type 1 diabetes (T1D). However, the association between this risk region and T1D subgroups and related clinical risk features, the underlying causal functional variant(s), putative candidate gene(s), and related mechanisms are yet to be elucidated.MethodsWe assessed the association between variant rs912304 and T1D, as well as islet autoimmunity and islet function, stratified by the diagnosed age of 12. We used epigenome bioinformatics analyses, dual luciferase reporter assays, and expression quantitative trait loci (eQTL) analyses to prioritize the most likely functional variant and potential causal gene. We also performed functional experiments to evaluate the role of the causal gene on islet function and its related mechanisms.ResultsWe identified rs912304 as a risk variant for T1D subgroups with diagnosed age ≥ 12 but not < 12. This variant is associated with residual islet function but not islet-specific autoantibody positivity in T1D individuals. Bioinformatics analysis indicated that rs912304 is a functional variant exhibiting spatial overlaps with enhancer active histone marks (H3K27ac and H3K4me1) and open chromatin status (ATAC-seq) in the human pancreas and islet tissues. Luciferase reporter gene assays and eQTL analyses demonstrated that the biallelic sites of rs912304 had differential allele-specific enhancer activity in beta cell lines and regulated STXBP6 expression, which was defined as the most putative causal gene based on Open Targets Genetics, GTEx v8 and Tiger database. Moreover, Stxbp6 was upregulated by T1D-related proinflammatory cytokines but not high glucose/fat. Notably, Stxbp6 over-expressed INS-1E cells exhibited decreasing insulin secretion and increasing cell apoptosis through Glut1 and Gadd45β, respectively.ConclusionsThis study expanded the genomic landscape regarding late-onset T1D risk and supported islet function mechanistically connected to T1D pathogenesis.

Uromodulin and progression of IgA nephropathy.

This study investigates the link between genetic variants associated with kidney function and immunoglobulin A (IgA) nephropathy (IgAN) progression. We recruited 961 biopsy-proven IgAN patients and 651 non-IgAN end-stage renal disease (ESRD) patients from Ruijin Hospital. Clinical and renal pathological data were collected. The primary outcome was the time to ESRD. A healthy population was defined as estimated glomerular filtration rate >60mL/min/1.73 m2 without albuminuria or hematuria. Fifteen single-nucleotide polymorphisms (SNPs) were selected from a genome-wide association study of kidney function and genotyped by the SNaPshot. Immunohistochemistry in renal tissue and ELISA in urine samples were performed to explore the potential functions of genetic variations. The rs77924615-G was independently associated with an increased risk for ESRD in IgAN patients after adjustments for clinical and pathologic indices, and treatment (adjusted hazard ratio 2.10; 95% confidence interval 1.14-3.88). No significant differences in ESRD-free survival time were found among different genotypes in non-IgAN ESRD patients (log-rank, P=.480). Moreover, rs77924615 exhibited allele-specific enhancer activity by dual-luciferase reporter assay. Accordingly, the urinary uromodulin-creatinine ratio (uUCR) was significantly higher in healthy individuals with rs77924615 AG or GG than in individuals with AA. Furthermore, uromodulin expression in tubular epithelial cells was higher in patients with rs77924615 AG or GG. Finally, we confirmed that an increased uUCR (P=.009) was associated with faster IgAN progression. The SNP rs77924615, which modulates the enhancer activity of the UMOD gene, is associated with renal function deterioration in IgAN patients by increasing uromodulin levels in both the renal tubular epithelium and urine.

KLF2 enhancer variant rs4808485 increases lupus risk by modulating inflammasome machinery and cellular homoeostasis

ObjectiveA recent genome-wide association study linked KLF2 as a novel Asian-specific locus for systemic lupus erythematosus (SLE) susceptibility. However, the underlying causal functional variant(s), cognate target gene(s) and genetic mechanisms...

A regulatory variant at 19p13.3 is associated with primary biliary cholangitis risk and ARID3A expression

Genome-wide association studies have identified 19p13.3 locus associated with primary biliary cholangitis (PBC). Here we aim to identify causative variant(s) and initiate efforts to define the mechanism by which the 19p13.3 locus variant(s) contributes to the pathogenesis of PBC. A genome-wide meta-analysis of 1931 PBC subjects and 7852 controls in two Han Chinese cohorts confirms the strong association between 19p13.3 locus and PBC. By integrating functional annotations, luciferase reporter assay and allele-specific chromatin immunoprecipitation, we prioritize rs2238574, an AT-Rich Interaction Domain 3A (ARID3A) intronic variant, as a potential causal variant at 19p13.3 locus. The risk allele of rs2238574 shows higher binding affinity of transcription factors, leading to an increased enhancer activity in myeloid cells. Genome-editing demonstrates the regulatory effect of rs2238574 on ARID3A expression through allele-specific enhancer activity. Furthermore, knock-down of ARID3A inhibits myeloid differentiation and activation pathway, and overexpression of the gene has the opposite effect. Finally, we find ARID3A expression and rs2238574 genotypes linked to disease severity in PBC. Our work provides several lines of evidence that a non-coding variant regulates ARID3A expression, presenting a mechanistic basis for association of 19p13.3 locus with the susceptibility to PBC.

Regulatory de novo mutations underlying intellectual disability.

The genetic aetiology of a major fraction of patients with intellectual disability (ID) remains unknown. De novo mutations (DNMs) in protein-coding genes explain up to 40% of cases, but the potential role of regulatory DNMs is still poorly understood. We sequenced 63 whole genomes from 21 ID probands and their unaffected parents. In addition, we analysed 30 previously sequenced genomes from exome-negative ID probands. We found that regulatory DNMs were selectively enriched in fetal brain-specific enhancers as compared with adult brain enhancers. DNM-containing enhancers were associated with genes that show preferential expression in the prefrontal cortex. Furthermore, we identified recurrently mutated enhancer clusters that regulate genes involved in nervous system development (CSMD1, OLFM1, and POU3F3). Most of the DNMs from ID probands showed allele-specific enhancer activity when tested using luciferase assay. Using CRISPR-mediated mutation and editing of epigenomic marks, we show that DNMs at regulatory elements affect the expression of putative target genes. Our results, therefore, provide new evidence to indicate that DNMs in fetal brain-specific enhancers play an essential role in the aetiology of ID.

Functional characterization of 5p15.33 risk locus in uveal melanoma reveals rs452384 as a functional variant and NKX2.4 as an allele-specific interactor

The TERT/CLPTM1L risk locus on chromosome 5p15.33 is a pleiotropic cancer risk locus in which multiple independent risk alleles have been identified, across well over ten cancer types. We previously conducted a genome-wide association study in uveal melanoma (UM), which uncovered a role for the TERT/CLPTM1L risk locus in this intraocular tumor and identified multiple highly correlated risk alleles. Aiming to unravel the biological mechanisms in UM of this locus, which contains a domain enriched in active chromatin marks and enhancer elements, we demonstrated the allele-specific enhancer activity of this risk region using reporter assays. In UM, we identified the functional variant rs452384, of which the C risk allele is associated with higher gene expression, increased CLPTM1L expression in UM tumors, and a longer telomere length in peripheral blood mononuclear cells. Electrophoretic mobility shift assays and quantitative mass spectrometry identified NKX2.4 as an rs452384-T-specific binding protein, whereas GATA4 preferentially interacted with rs452384-C. Knockdown of NKX2.4 but not GATA4 resulted in increased TERT and CLPTM1L expression. In summary, the UM risk conferred by the 5p locus is at least partly due to rs452384, for which NKX2.4 presents strong differential binding activity and regulates CLPTM1L and TERT expression. Altogether, our work unraveled some of the complex regulatory mechanisms at the 5p15.33 susceptibility region in UM, and this might also shed light on shared mechanisms with other tumor types affected by this susceptibility region.

Haplotype-specific chromatin looping reveals genetic interactions of regulatory regions modulating gene expression in 8p23.1.

A major goal of genetics research is to elucidate mechanisms explaining how genetic variation contributes to phenotypic variation. The genetic variants identified in genome-wide association studies (GWASs) generally explain only a small proportion of heritability of phenotypic traits, the so-called missing heritability problem. Recent evidence suggests that additional common variants beyond lead GWAS variants contribute to phenotypic variation; however, their mechanistic underpinnings generally remain unexplored. Herein, we undertake a study of haplotype-specific mechanisms of gene regulation at 8p23.1 in the human genome, a region associated with a number of complex diseases. The FAM167A-BLK locus in this region has been consistently found in the genome-wide association studies (GWASs) of systemic lupus erythematosus (SLE) in all major ancestries. Our haplotype-specific chromatin interaction (Hi-C) experiments, allele-specific enhancer activity measurements, genetic analyses, and epigenome editing experiments revealed that: 1) haplotype-specific long-range chromatin interactions are prevalent in 8p23.1; 2) BLK promoter and cis-regulatory elements cooperatively interact with haplotype-specificity; 3) genetic variants at distal regulatory elements are allele-specific modifiers of the promoter variants at FAM167A-BLK; 4) the BLK promoter interacts with and, as an enhancer-like promoter, regulates FAM167A expression and 5) local allele-specific enhancer activities are influenced by global haplotype structure due to chromatin looping. Although systemic lupus erythematosus causal variants at the FAM167A-BLK locus are thought to reside in the BLK promoter region, our results reveal that genetic variants at distal regulatory elements modulate promoter activity, changing BLK and FAM167A gene expression and disease risk. Our results suggest that global haplotype-specific 3-dimensional chromatin looping architecture has a strong influence on local allelic BLK and FAM167A gene expression, providing mechanistic details for how regional variants controlling the BLK promoter may influence disease risk.

A risk variant for Barrett's esophagus and esophageal adenocarcinoma at chr8p23.1 affects enhancer activity and implicates multiple gene targets.

Nineteen genetic susceptibility loci for esophageal adenocarcinoma (EAC) and its precursor Barrett's esophagus (BE) have been identified through genome-wide association studies (GWAS). Clinical translation of such discoveries, however, has been hindered by the slow pace of discovery of functional/causal variants and gene targets at these loci. We previously developed a systematic informatics pipeline to prioritize candidate functional variants using functional potential scores, applied the pipeline to select high-scoring BE/EAC risk loci and validated a functional variant at chr19p13.11 (rs10423674). Here, we selected two additional prioritized loci for experimental interrogation: chr3p13/rs1522552 and chr8p23.1/rs55896564. Candidate enhancer regions encompassing these variants were evaluated using luciferase reporter assays in two EAC cell lines. One of the two regions tested exhibited allele-specific enhancer activity - 8p23.1/rs55896564. CRISPR-mediated deletion of the putative enhancer in EAC cell lines correlated with reduced expression of three candidate gene targets: B lymphocyte kinase (BLK), nei like DNA glycosylase 2 (NEIL2) and cathepsin B (CTSB). Expression quantitative trait locus (eQTL) mapping in normal esophagus and stomach revealed strong associations between the BE/EAC risk allele at rs55896564 (G) and lower expression of CTSB, a protease gene implicated in epithelial wound repair. These results further support the utility of functional potential scores for GWAS variant prioritization, and provide the first experimental evidence of a functional variant and risk enhancer at the 8p23.1 GWAS locus. Identification of CTSB, BLK and NEIL2 as candidate gene targets suggests that altered expression of these genes may underlie the genetic risk association at 8p23.1 with BE/EAC.

An enhancer variant at 16q22.1 predisposes to hepatocellular carcinoma via regulating PRMT7 expression

Most cancer causal variants are found in gene regulatory elements, e.g., enhancers. However, enhancer variants predisposing to hepatocellular carcinoma (HCC) remain unreported. Here we conduct a genome-wide survey of HCC-susceptible enhancer variants through a three-stage association study in 11,958 individuals and identify rs73613962 (T > G) within the intronic region of PRMT7 at 16q22.1 as a susceptibility locus of HCC (OR = 1.41, P = 6.02 × 10−10). An enhancer dual-luciferase assay indicates that the rs73613962-harboring region has allele-specific enhancer activity. CRISPR-Cas9/dCas9 experiments further support the enhancer activity of this region to regulate PRMT7 expression. Mechanistically, transcription factor HNF4A binds to this enhancer region, with preference to the risk allele G, to promote PRMT7 expression. PRMT7 upregulation contributes to in vitro, in vivo, and clinical HCC-associated phenotypes, possibly by affecting the p53 signaling pathway. This concept of HCC pathogenesis may open a promising window for HCC prevention/treatment.

Identification of potential functional variants and genes at 18q21.1 associated with the carcinogenesis of colorectal cancer.

Genome-wide association studies (GWAS) have identified more than 160 susceptibility loci for colorectal cancer (CRC). The effects of these variants, particularly their mechanisms, however, remain unclear. In this study, a comprehensive functional annotation of CRC-related GWAS signals was firstly conducted to identify the potential causal variants. We found that the SNP rs7229639 in intron 3 of SMAD7 at 18q21.1 might serve as a putative functional variant in CRC. The SNP rs7229639 is located in a region with evidence of regulatory potential. Dual-luciferase reporter assays revealed that three other SNPs (rs77544449, rs60385309 and rs72917785), in strong linkage disequilibrium (LD) with rs7229639, exhibited allele-specific enhancer activity, of which one of the target genes may conceivably be LIPG, as suggested by eQTL association data and Hi-C data. We also verified that LIPG promoted malignancy of CRC cells in vitro, with supporting clinical data indicating that LIPG is upregulated and correlated with a poor prognosis in CRC. Finally, pitavastatin was observed to exhibit an anti-CRC activity and modest inhibition of LIPG mRNA levels. Collectively, our data suggest that these functional variants at 18q21.1 are involved in the pathogenesis of CRC by modulating enhancer activity, and possibly LIPG expression, thus indicating a promising therapeutic target for CRC. The results of functional annotation in our investigation could also serve as an inventory for CRC susceptibility SNPs and offer guides for post-GWAS downstream functional studies.

Identifying the lungs as a susceptible site for allele-specific regulatory changes associated with type 1 diabetes risk

Type 1 diabetes (T1D) etiology is complex. We developed a machine learning approach that ranked the tissue-specific transcription regulatory effects for T1D SNPs and estimated their relative contributions to conversion to T1D by integrating case and control genotypes (Wellcome Trust Case Control Consortium and UK Biobank) with tissue-specific expression quantitative trait loci (eQTL) data. Here we show an eQTL (rs6679677) associated with changes to AP4B1-AS1 transcript levels in lung tissue makes the largest gene regulatory contribution to the risk of T1D development. Luciferase reporter assays confirmed allele-specific enhancer activity for the rs6679677 tagged locus in lung epithelial cells (i.e. A549 cells; C > A reduces expression, p = 0.005). Our results identify tissue-specific eQTLs for SNPs associated with T1D. The strongest tissue-specific eQTL effects were in the lung and may help explain associations between respiratory infections and risk of islet autoantibody seroconversion in young children.

Prioritization and functional analysis of GWAS risk loci for Barrett's esophagus and esophageal adenocarcinoma.

Genome-wide association studies (GWAS) have identified ~20 genetic susceptibility loci for esophageal adenocarcinoma (EAC), and its precursor, Barrett's esophagus (BE). Despite such advances, functional/causal variants and gene targets at these loci remain undefined, hindering clinical translation. A key challenge is that most causal variants map to non-coding regulatory regions such as enhancers, and typically, numerous potential candidate variants at GWAS loci require testing. We developed a systematic informatics pipeline for prioritizing candidate functional variants via integrative functional potential scores (FPS) consolidated from multi-omics annotations, and used this pipeline to identify two high-scoring variants for experimental interrogation: chr9q22.32/rs11789015 and chr19p13.11/rs10423674. Minimal candidate enhancer regions spanning these variants were evaluated using luciferase reporter assays in two EAC cell lines. One of the two variants tested (rs10423674) exhibited allele-specific enhancer activity. CRISPR-mediated deletion of the putative enhancer region in EAC cell lines correlated with reduced expression of two genes-CREB-regulated transcription coactivator 1 (CRTC1) and Cartilage oligomeric matrix protein (COMP); expression of five other genes remained unchanged (CRLF1, KLHL26, TMEM59L, UBA52, RFXANK). Expression quantitative trait locus mapping indicated that rs10423674 genotype correlated with CRTC1 and COMP expression in normal esophagus. This study represents the first experimental effort to bridge GWAS associations to biology in BE/EAC and supports the utility of FPS to guide variant prioritization. Our findings reveal a functional variant and candidate risk enhancer at chr19p13.11 and implicate CRTC1 and COMP as putative gene targets, suggesting that altered expression of these genes may underlie the BE/EAC risk association.

Integrative analysis of liver-specific non-coding regulatory SNPs associated with the risk of coronary artery disease

Genetic factors underlying coronary artery disease (CAD) have been widely studied using genome-wide association studies (GWASs). However, the functional understanding of the CAD loci has been limited by the fact that a majority of GWAS variants are located within non-coding regions with no functional role. High cholesterol and dysregulation of the liver metabolism such as non-alcoholic fatty liver disease confer an increased risk of CAD. Here, we studied the function of non-coding single-nucleotide polymorphisms in CAD GWAS loci located within liver-specific enhancer elements by identifying their potential target genes using liver cis-eQTL analysis and promoter Capture Hi-C in HepG2 cells. Altogether, 734 target genes were identified of which 121 exhibited correlations to liver-related traits. To identify potentially causal regulatory SNPs, the allele-specific enhancer activity was analyzed by (1) sequence-based computational predictions, (2) quantification of allele-specific transcription factor binding, and (3) STARR-seq massively parallel reporter assay. Altogether, our analysis identified 1,277 unique SNPs that display allele-specific regulatory activity. Among these, susceptibility enhancers near important cholesterol homeostasis genes (APOB, APOC1, APOE, and LIPA) were identified, suggesting that altered gene regulatory activity could represent another way by which genetic variation regulates serum lipoprotein levels. Using CRISPR-based perturbation, we demonstrate how the deletion/activation of a single enhancer leads to changes in the expression of many target genes located in a shared chromatin interaction domain. Our integrative genomics approach represents a comprehensive effort in identifying putative causal regulatory regions and target genes that could predispose to clinical manifestation of CAD by affecting liver function.

Mechanistic Characterization of RASGRP1 Variants Identifies an hnRNP-K-Regulated Transcriptional Enhancer Contributing to SLE Susceptibility.

Systemic lupus erythematosus (SLE) is an autoimmune disease with a strong genetic component. We recently identified a novel SLE susceptibility locus near RASGRP1, which governs the ERK/MAPK kinase cascade and B-/T-cell differentiation and development. However, precise causal RASGRP1 functional variant(s) and their mechanisms of action in SLE pathogenesis remain undefined. Our goal was to fine-map this locus, prioritize genetic variants likely to be functional, experimentally validate their biochemical mechanisms, and determine the contribution of these SNPs to SLE risk. We performed a meta-analysis across six Asian and European cohorts (9,529 cases; 22,462 controls), followed by in silico bioinformatic and epigenetic analyses to prioritize potentially functional SNPs. We experimentally validated the functional significance and mechanism of action of three SNPs in cultured T-cells. Meta-analysis identified 18 genome-wide significant (p < 5 × 10−8) SNPs, mostly concentrated in two haplotype blocks, one intronic and the other intergenic. Epigenetic fine-mapping, allelic, eQTL, and imbalance analyses predicted three transcriptional regulatory regions with four SNPs (rs7170151, rs11631591-rs7173565, and rs9920715) prioritized for functional validation. Luciferase reporter assays indicated significant allele-specific enhancer activity for intronic rs7170151 and rs11631591-rs7173565 in T-lymphoid (Jurkat) cells, but not in HEK293 cells. Following up with EMSA, mass spectrometry, and ChIP-qPCR, we detected allele-dependent interactions between heterogeneous nuclear ribonucleoprotein K (hnRNP-K) and rs11631591. Furthermore, inhibition of hnRNP-K in Jurkat and primary T-cells downregulated RASGRP1 and ERK/MAPK signaling. Comprehensive association, bioinformatics, and epigenetic analyses yielded putative functional variants of RASGRP1, which were experimentally validated. Notably, intronic variant (rs11631591) is located in a cell type-specific enhancer sequence, where its risk allele binds to the hnRNP-K protein and modulates RASGRP1 expression in Jurkat and primary T-cells. As risk allele dosage of rs11631591 correlates with increased RASGRP1 expression and ERK activity, we suggest that this SNP may underlie SLE risk at this locus.

Abstract LB-232: Inherited determinants of early recurrent somatic mutations in prostate cancer

Abstract Background: Prostate cancer (PCa) is a clinically heterogeneous hormone regulated disease with a strong genetic component and characterized by distinct molecular subclasses. Since hormone levels and related transcription factors vary during individual’s lifetime, we hypothesized that inherited variants within regulatory elements of the human genome might act as facilitators of specific early somatic events in PCa in age/hormone dependent manner. Methods: We developed a computational approach to nominate heritable facilitators of somatic genomic events in the context of the Androgen Receptor (AR) signaling. Germline variants within functionally active regions of the genome (N=21,364) defined as per ENCODE data were ranked by the so called Trigger Score, which quantifies the fraction of the transcriptome putatively modulated by each variant leveraging individuals’ genotypes and transcript levels from a cohort of human benign prostate tissues. Results: We identified 300 polymorphic regulatory elements that associate with DNA repair and hormone regulated transcript levels. In a cohort of 539 PCa and more than 3,000 controls, a variant on 7p14.3 strongly associated with an early recurrent prostate cancer specific somatic mutation in the Speckle-Type POZ protein (SPOP) gene (OR=5.54, P=1.22e-08). Integrated computational analysis suggested a broad transcriptional regulatory impact of the variant. In vitro luciferase assay in PC-3 cells demonstrated allele-specific enhancer activity of the locus that is concomitantly modulated by AR and by CCAAT/Enhancer Binding Protein beta (CEBPB). ChIP assays experiments showed AR and CEBPB binding to the locus. Upon deletion of the genomic locus via CRISPR-Cas9 followed by transcriptome sequencing, we observed differential transcript levels along chromosome 7 between edited versus control cells, but not in edited versus edited or control versus control, and nominated genes that were also predicted to physically interact with the 7p14.3 locus by previously generated Hi-C data from prostate cells. Conclusions: Our study suggests a relevant role in cancer predisposition of non-coding variants that lead to allele specific transactivation with age dependent and tissue specific effect manifesting in early somatic genomic events. Citation Format: Alessandro Romanel, Sonia Iolanda Garritano, Blerta Stringa, Mirjam Blattner, Davide Dalfovo, Dimple Chakravarty, Kellie Cotter, Priyanka Dhingra, Paola Gasperini, Olivier Elemento, Andrea Sboner, Alberto Inga, Ekta Khurana, MarK Rubin, Francesca Demichelis. Inherited determinants of early recurrent somatic mutations in prostate cancer [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2017; 2017 Apr 1-5; Washington, DC. Philadelphia (PA): AACR; Cancer Res 2017;77(13 Suppl):Abstract nr LB-232. doi:10.1158/1538-7445.AM2017-LB-232

Abstract 97: Genotype to Phenotype: Function of Common Noncoding and Rare Coding Variants In ANGPTL3

Human genetics studies have demonstrated a strong link between ANGPTL3 , which encodes lipoprotein lipase inhibitor Angiopoietin-like 3, and blood lipid phenotypes. Rare nonsense ANGPTL3 mutations were identified in patients with familial combined hypolipidemia, while common variants at the ANGPTL3 locus have been found by genome-wide association studies (GWASs) to associate with lower triglycerides (TGs) and low-density lipoprotein cholesterol. In light of the seemingly favorable clinical consequences of ANGPTL3 deficiency, we established an experimental framework to identify (1) causal common variants that regulate ANGPTL3 expression and (2) rare missense mutations that disrupt ANGPTL3 function. Using massively parallel reporter assays, we profiled the regulatory activity of all the common variants linked ( r 2 ≥ 0.5) to the lead GWAS SNP in the ANGPTL3 locus and found that rs10889356 demonstrated significant allele-specific enhancer activity. To validate this finding, we used CRISPR-Cas9 to alter the SNP in a human pluripotent stem cell line. When differentiated into hepatocytes, altered cells displayed a 67% increase in ANGPTL3 expression ( n = 4 wild-type and 4 mutant clones, P = 0.007). CRISPR interference using each of three guide RNAs targeting the SNP in HepG2 cells also substantially increased ANGPTL3 expression. These findings support rs10889356- ANGPTL3 as a causal SNP-gene set. Next, we examined the coding regions of ANGPTL3 in 20,000 sequenced individuals and sought to experimentally define rare missense variants using a mouse model. We used CRISPR-Cas9 to generate Angptl3 knockout mice, which exhibited decreased TG (61%, P &lt; 0.001) and decreased cholesterol (31%, P &lt; 0.002). We reconstituted the knockout mice to normal expression levels with adenoviruses expressing either wild-type ANGPTL3 or missense variant ANGPTL3 . So far we have assessed 28 rare missense variants computationally predicted to be deleterious, of which only 10—D42N, K58E, S117P, P264S, Q286H, L315S, L360Q, T383I, T383S, and Y417C—were validated as loss-of-function (conferring &lt;25% of wild-type activity as assessed by changes in both TG and cholesterol levels), underscoring the need for functional characterization of variants of uncertain significance.

Identification and validation of regulatory SNPs that modulate transcription factor chromatin binding and gene expression in prostate cancer.

Prostate cancer (PCa) is the second most common solid tumor for cancer related deaths in American men. Genome wide association studies (GWAS) have identified single nucleotide polymorphisms (SNPs) associated with the increased risk of PCa. Because most of the susceptibility SNPs are located in noncoding regions, little is known about their functional mechanisms. We hypothesize that functional SNPs reside in cell type-specific regulatory elements that mediate the binding of critical transcription factors (TFs), which in turn result in changes in target gene expression. Using PCa-specific functional genomics data, here we identify 38 regulatory candidate SNPs and their target genes in PCa. Through risk analysis by incorporating gene expression and clinical data, we identify 6 target genes (ZG16B, ANKRD5, RERE, FAM96B, NAALADL2 and GTPBP10) as significant predictors of PCa biochemical recurrence. In addition, 5 SNPs (rs2659051, rs10936845, rs9925556, rs6057110 and rs2742624) are selected for experimental validation using Chromatin immunoprecipitation (ChIP), dual-luciferase reporter assay in LNCaP cells, showing allele-specific enhancer activity. Furthermore, we delete the rs2742624-containing region using CRISPR/Cas9 genome editing and observe the drastic downregulation of its target gene UPK3A. Taken together, our results illustrate that this new methodology can be applied to identify regulatory SNPs and their target genes that likely impact PCa risk. We suggest that similar studies can be performed to characterize regulatory variants in other diseases.

Regulation of the imprinted Dlk1-Dio3 locus by allele-specific enhancer activity.

Genomic imprinting is a critical developmental process characteristic of parent of origin-specific gene expression. It is well accepted that differentially DNA-methylated regions (DMRs) and enhancers are two major classes of cis-elements determining parent of origin-specific gene expression, with each recruiting different sets of transcription factors. Previously, we identified the AF4/FMR2 (AFF) family protein AFF3 within the transcription elongation complex SEC-L3. Here, we report that AFF3 can specifically bind both gametic DMRs (gDMRs) and enhancers within imprinted loci in an allele-specific manner. We identify the molecular regulators involved in the recruitment of AFF3 to gDMRs and provide mechanistic insight into the requirement of AFF3 at an enhancer for the expression of an ∼200-kb polycistronic transcript within the imprinted Dlk1-Dio3 locus. Our data suggest that the heterochromatic environment at the gDMR reinforces silencing of its related enhancer by controlling the binding and activity of AFF3 in an allele-specific manner. In summary, this study provides molecular details about the regulation of dosage-critical imprinted gene expression through the regulated binding of the transcription elongation factor AFF3 between a DMR and an enhancer.

In-silico identification and functional validation of allele-dependent AR enhancers.

Androgen Receptor (AR) and Estrogen Receptors (ERs) are key nuclear receptors that can cooperate in orchestrating gene expression programs in multiple tissues and diseases, targeting binding elements in promoters and distant enhancers. We report the unbiased identification of enhancer elements bound by AR and ER-α whose activity can be allele-specific depending on the status of nearby Single Nucleotide Polymorphisms (SNP). ENCODE data were computationally mined to nominate genomic loci with: (i) chromatin signature of enhancer activity from activation histone marks, (ii) binding evidence by AR and ER-α, (iii) presence of a SNP. Forty-one loci were identified and two, on 1q21.3 and 13q34, selected for characterization by gene reporter, Chromatin immunoprecipitation (ChIP) and RT-qPCR assays in breast (MCF7) and prostate (PC-3) cancer-derived cell lines. We observed allele-specific enhancer activity, responsiveness to ligand-bound AR, and potentially influence on the transcription of closely located genes (RAB20, ING1, ARHGEF7, ADAM15). The 1q21.3 variant, rs2242193, showed impact on AR binding in MCF7 cells that are heterozygous for the SNP. Our unbiased genome-wide search proved to be an efficient methodology to discover new functional polymorphic regulatory regions (PRR) potentially acting as risk modifiers in hormone-driven cancers and overall nominated SNPs in PRR across 136 transcription factors.

An 8q24 gene desert variant associated with prostate cancer risk confers differential in vivo activity to a MYC enhancer

Genome-wide association studies (GWAS) routinely identify risk variants in noncoding DNA, as exemplified by reports of multiple single nucleotide polymorphisms (SNPs) associated with prostate cancer in five independent regions in a gene desert on 8q24. Two of these regions also have been associated with breast and colorectal cancer. These findings implicate functional variation within long-range cis-regulatory elements in disease etiology. We used an in vivo bacterial artificial chromosome (BAC) enhancer-trapping strategy in mice to scan a half-megabase of the 8q24 gene desert encompassing the prostate cancer-associated regions for long-range cis-regulatory elements. These BAC assays identified both prostate and mammary gland enhancer activities within the region. We demonstrate that the 8q24 cancer-associated variant rs6983267 lies within an in vivo prostate enhancer whose expression mimics that of the nearby MYC proto-oncogene. Additionally, we show that the cancer risk allele increases prostate enhancer activity in vivo relative to the non-risk allele. This allele-specific enhancer activity is detectable during early prostate development and throughout prostate maturation, raising the possibility that this SNP could assert its influence on prostate cancer risk before tumorigenesis occurs. Our study represents an efficient strategy to build experimentally on GWAS findings with an in vivo method for rapidly scanning large regions of noncoding DNA for functional cis-regulatory sequences harboring variation implicated in complex diseases.