CAUSEL: an epigenome- and genome-editing pipeline for establishing function of noncoding GWAS variants

Methods for Scarless, Selection-Free Generation of Human Cells and Allele-Specific Functional Analysis of Disease-Associated SNPs and Variants of Uncertain Significance

With the continued emergence of risk loci from Genome-Wide Association studies and variants of uncertain significance identified from patient sequencing, better methods are required to translate these human genetic findings into improvements in public health. Here we combine CRISPR/Cas9 gene editing with an innovative high-throughput genotyping pipeline utilizing KASP (Kompetitive Allele-Specific PCR) genotyping technology to create scarless isogenic cell models of cancer variants in ~1 month. We successfully modeled two novel variants previously identified by our lab in the PALB2 gene in HEK239 cells, resulting in isogenic cells representing all three genotypes for both variants. We also modeled a known functional risk SNP of colorectal cancer, rs6983267, in HCT-116 cells. Cells with extremely low levels of gene editing could still be identified and isolated using this approach. We also introduce a novel molecular assay, ChIPnQASO (Chromatin Immunoprecipitation and Quantitative Allele-Specific Occupation), which uses the same technology to reveal allele-specific function of these variants at the DNA-protein interaction level. We demonstrated preferential binding of the transcription factor TCF7L2 to the rs6983267 risk allele over the non-risk. Our pipeline provides a platform for functional variant discovery and validation that is accessible and broadly applicable for the progression of efforts towards precision medicine.

A Dementia-Associated Risk Variant near TMEM106B Alters Chromatin Architecture and Gene Expression

Neurodegenerative diseases pose an extraordinary threat to the world's aging population, yet no disease-modifying therapies are available. Although genome-wide association studies (GWASs) have identified hundreds of risk loci for neurodegeneration, the mechanisms by which these loci influence disease risk are largely unknown. Here, we investigated the association between common genetic variants at the 7p21 locus and risk of the neurodegenerative disease frontotemporal lobar degeneration. We showed that variants associated with disease risk correlate with increased expression of the 7p21 gene TMEM106B and no other genes; co-localization analyses implicated a common causal variant underlying both association with disease and association with TMEM106B expression in lymphoblastoid cell lines and human brain. Furthermore, increases in the amount of TMEM106B resulted in increases in abnormal lysosomal phenotypes and cell toxicity in both immortalized cell lines and neurons. We then combined fine-mapping, bioinformatics, and bench-based approaches to functionally characterize all candidate causal variants at this locus. This approach identified a noncoding variant, rs1990620, that differentially recruits CTCF in lymphoblastoid cell lines and human brain to influence CTCF-mediated long-range chromatin-looping interactions between multiple cis-regulatory elements, including the TMEM106B promoter. Our findings thus provide an in-depth analysis of the 7p21 locus linked by GWASs to frontotemporal lobar degeneration, nominating a causal variant and causal mechanism for allele-specific expression and disease association at this locus. Finally, we show that genetic variants associated with risk of neurodegenerative diseases beyond frontotemporal lobar degeneration are enriched in CTCF-binding sites found in brain-relevant tissues, implicating CTCF-mediated gene regulation in risk of neurodegeneration more generally.

Association between genetic polymorphisms of long non-coding RNA PRNCR1 and prostate cancer risk in a sample of the Iranian population

The aim of the present study was to determine whether there is an association between the long non-coding RNA (lncRNA) prostate cancer-associated non-coding RNA 1 (PRNCR1) polymorphisms and prostate cancer (PCa) risk in a sample of the Iranian population. This case-control study was performed on 178 patients with PCa and 180 subjects with benign prostatic hyperplasia (BPH). Genotyping assay was performed by polymerase chain reaction-restriction fragment length polymorphism. The findings indicated that the GG genotype of the rs13252298 A>G variant significantly increased the risk of PCa (odds ratio=3.49, 95% confidence interval: 1.79-6.81, P=0.0001) compared with AA+AG. As regards the rs1456315 G>A polymorphism, the AG genotype and G allele significantly increased the risk of PCa. As regards the rs7841060 T>G variant, the findings demonstrated that this TG genotype and the G allele significantly increased the risk of PCa. The rs7007694 T>C variant was not found to be associated with the risk of PCa. Haplotype analysis indicated that GTGA and GTGG significantly increased the risk of PCa compared with rs1456315A/rs7007694T/rs7841060T/rs13252298G (ATTG). The PRNCR1 variants were not found to be significantly associated with the clinicopathological characteristics of PCa patients. In conclusion, our findings support an association between PRNCR1 variants and the risk of PCa in a sample of the Iranian population.

Functional interrogation of non-coding DNA through CRISPR genome editing

Methodologies to interrogate non-coding regions have lagged behind coding regions despite comprising the vast majority of the genome. However, the rapid evolution of clustered regularly interspaced short palindromic repeats (CRISPR)-based genome editing has provided a multitude of novel techniques for laboratory investigation including significant contributions to the toolbox for studying non-coding DNA. CRISPR-mediated loss-of-function strategies rely on direct disruption of the underlying sequence or repression of transcription without modifying the targeted DNA sequence. CRISPR-mediated gain-of-function approaches similarly benefit from methods to alter the targeted sequence through integration of customized sequence into the genome as well as methods to activate transcription. Here we review CRISPR-based loss- and gain-of-function techniques for the interrogation of non-coding DNA.

Down-regulation of ARID1A is sufficient to initiate neoplastic transformation along with epigenetic reprogramming in non-tumorigenic endometriotic cells

The chromatin remodeler AT-Rich Interactive Domain 1A (ARID1A) is frequently mutated in ovarian clear cell carcinoma (OCCC) and endometriosis precursor lesions. Here, we show that knocking down ARID1A in an immortalized endometriosis cell line is sufficient to induce phenotypic changes indicative of neoplastic transformation as evidenced by higher efficiency of anchorage-independent growth, increased propensity to adhere to collagen, and greater capacity to invade basement membrane extract in vitro. ARID1A knockdown is associated with expression dysregulation of 99 target genes, and many of these expression changes are also observed in primary OCCC tissues. Further, pathway analysis indicates these genes fall within networks highly relevant to tumorigenesis including integrin and paxillin pathways. We demonstrate that the down-regulation of ARID1A does not markedly alter global chromatin accessibility or DNA methylation but unexpectedly, we find strong increases in the active H3K27ac mark in promoter regions and decreases of H3K27ac at potential enhancers. Taken together, these data provide evidence that ARID1A mutation can be an early stage event in the oncogenic transformation of endometriosis cells giving rise to OCCC.

Identification of 12 new susceptibility loci for different histotypes of epithelial ovarian cancer

To identify common alleles associated with different histotypes of epithelial ovarian cancer (EOC), we pooled data from multiple genome-wide genotyping projects totaling 25,509 EOC cases and 40,941 controls. We identified nine new susceptibility loci for different EOC histotypes: six for serous EOC histotypes (3q28, 4q32.3, 8q21.11, 10q24.33, 18q11.2 and 22q12.1), two for mucinous EOC (3q22.3 and 9q31.1) and one for endometrioid EOC (5q12.3). We then performed meta-analysis on the results for high-grade serous ovarian cancer with the results from analysis of 31,448 BRCA1 and BRCA2 mutation carriers, including 3,887 mutation carriers with EOC. This identified three additional susceptibility loci at 2q13, 8q24.1 and 12q24.31. Integrated analyses of genes and regulatory biofeatures at each locus predicted candidate susceptibility genes, including OBFC1, a new candidate susceptibility gene for low-grade and borderline serous EOC.

Modulation of transcription factor binding and epigenetic regulation of the MLH1 CpG island and shore by polymorphism rs1800734 in colorectal cancer

The MLH1 promoter polymorphism rs1800734 is associated with MLH1 CpG island hypermethylation and expression loss in colorectal cancer (CRC). Conversely, variant rs1800734 is associated with MLH1 shore, but not island, hypomethylation in peripheral blood mononuclear cell DNA. To explore these distinct patterns, MLH1 CpG island and shore methylation was assessed in CRC cell lines stratified by rs1800734 genotype. Cell lines containing the variant A allele demonstrated MLH1 shore hypomethylation compared to wild type (GG). There was significant enrichment of transcription factor AP4 at the MLH1 promoter in GG and GA cell lines, but not the AA cell line, by chromatin immunoprecipitation studies. Preferential binding to the G allele was confirmed by sequencing in the GA cell line. The enhancer-associated histone modification H3K4me1 was enriched at the MLH1 shore; however, H3K27ac was not, indicating the shore is an inactive enhancer. These results demonstrate the role of variant rs1800734 in altering transcription factor binding as well as epigenetics at regions beyond the MLH1 CpG island in which it is located.

The OncoArray Consortium: A Network for Understanding the Genetic Architecture of Common Cancers

BACKGROUND

Common cancers develop through a multistep process often including inherited susceptibility. Collaboration among multiple institutions, and funding from multiple sources, has allowed the development of an inexpensive genotyping microarray, the OncoArray. The array includes a genome-wide backbone, comprising 230,000 SNPs tagging most common genetic variants, together with dense mapping of known susceptibility regions, rare variants from sequencing experiments, pharmacogenetic markers, and cancer-related traits.

METHODS

The OncoArray can be genotyped using a novel technology developed by Illumina to facilitate efficient genotyping. The consortium developed standard approaches for selecting SNPs for study, for quality control of markers, and for ancestry analysis. The array was genotyped at selected sites and with prespecified replicate samples to permit evaluation of genotyping accuracy among centers and by ethnic background.

RESULTS

The OncoArray consortium genotyped 447,705 samples. A total of 494,763 SNPs passed quality control steps with a sample success rate of 97% of the samples. Participating sites performed ancestry analysis using a common set of markers and a scoring algorithm based on principal components analysis.

CONCLUSIONS

Results from these analyses will enable researchers to identify new susceptibility loci, perform fine-mapping of new or known loci associated with either single or multiple cancers, assess the degree of overlap in cancer causation and pleiotropic effects of loci that have been identified for disease-specific risk, and jointly model genetic, environmental, and lifestyle-related exposures.

IMPACT

Ongoing analyses will shed light on etiology and risk assessment for many types of cancer. Cancer Epidemiol Biomarkers Prev; 26(1); 126-35. ©2016 AACR.

BNC2 is a putative tumor suppressor gene in high-grade serous ovarian carcinoma and impacts cell survival after oxidative stress

Rs3814113 is the single-nucleotide polymorphism (SNP) showing the strongest association with high-grade serous ovarian carcinoma (HGSOC) incidence and is located in an intergenic region about 44 kb downstream of basonuclin 2 (BNC2) gene. Lifetime number of ovulations is associated with increased risk to develop HGSOC, probably because of cell damage of extrauterine Müllerian epithelium by ovulation-induced oxidative stress. However, the impact of low-penetrance HGSOC risk alleles (e.g. rs3814113) on the damage induced by oxidative stress remains unclear. Therefore, the purpose of this study was to investigate whether rs3814113 genetic interval regulates BNC2 expression and whether BNC2 expression levels impact on cell survival after oxidative stress. To do this, we analyzed gene expression levels of BNC2 first in HGSOC data sets and then in an isogenic cell line that we engineered to carry a 5 kb deletion around rs3814113. Finally, we silenced BNC2 and measured surviving cells after hydrogen peroxide (H2O2) treatment to simulate oxidative stress after ovulation. In this paper, we describe that BNC2 expression levels are reduced in HGSOC samples compared with control samples, and that BNC2 expression levels decrease following oxidative stress and ovulation in vitro and in vivo, respectively. Moreover, deletion of 5 kb surrounding rs3814113 decreases BNC2 expression levels in an isogenic cell line, and silencing of BNC2 expression levels increases cell survival after H2O2 treatment. Altogether, our findings suggest that the intergenic region located around rs3814113 regulates BNC2 expression, which in turn affects cell survival after oxidative stress response. Indeed, HGSOC samples present lower BNC2 expression levels that probably, in the initial phases of oncogenic transformation, conferred resistance to oxidative stress and ultimately reduced the clearance of cells with oxidative-induced damages.

Functional mechanisms underlying pleiotropic risk alleles at the 19p13.1 breast-ovarian cancer susceptibility locus

A locus at 19p13 is associated with breast cancer (BC) and ovarian cancer (OC) risk. Here we analyse 438 SNPs in this region in 46,451 BC and 15,438 OC cases, 15,252 BRCA1 mutation carriers and 73,444 controls and identify 13 candidate causal SNPs associated with serous OC (P=9.2 × 10(-20)), ER-negative BC (P=1.1 × 10(-13)), BRCA1-associated BC (P=7.7 × 10(-16)) and triple negative BC (P-diff=2 × 10(-5)). Genotype-gene expression associations are identified for candidate target genes ANKLE1 (P=2 × 10(-3)) and ABHD8 (P<2 × 10(-3)). Chromosome conformation capture identifies interactions between four candidate SNPs and ABHD8, and luciferase assays indicate six risk alleles increased transactivation of the ADHD8 promoter. Targeted deletion of a region containing risk SNP rs56069439 in a putative enhancer induces ANKLE1 downregulation; and mRNA stability assays indicate functional effects for an ANKLE1 3'-UTR SNP. Altogether, these data suggest that multiple SNPs at 19p13 regulate ABHD8 and perhaps ANKLE1 expression, and indicate common mechanisms underlying breast and ovarian cancer risk.

Identification and function of enhancers in the human genome

The study of gene regulation has rapidly advanced by leveraging next-generation sequencing to identify and characterize the cis and trans elements that are critical for defining cell identity. These advances have paralleled a movement towards whole genome sequencing in clinics. These two tracks have increasingly synergized to underscore the importance of cis-regulatory elements in development as well produce countless studies implicating these elements in human disease. Other studies have emphasized the clinical phenotypes associated with variation or mutations in trans factors, including non-coding RNAs and chromatin regulators. These studies highlight the importance of obtaining a comprehensive understanding of mammalian gene regulation for predicting the impact of genetic variation on patient phenotypes. Currently lagging behind the generation of vast datasets and annotations is our ability to examine these putative elements in the dynamic context of a developing organism.

Every transcription factor deserves its map: Scaling up epitope tagging of proteins to bypass antibody problems

Genome-wide identification of transcription factor binding sites with the ChIP-seq method is an extremely important scientific endeavor - one that should ideally be performed for every transcription factor in as many cell types as possible. A major hurdle on the way to this goal is the necessity for a specific, ChIP-grade antibody for each transcription factor of interest, which is often not available. Here, we describe CETCh-seq, a recently published method utilizing genome engineering with the CRISPR/Cas9 system to circumvent the need for a specific antibody. Using the CETCh-seq method, targeted genomic editing results in an epitope-tagged transcription factor, which is recognized by a well-characterized, standard antibody, efficacious for ChIP-seq. We have used CETCh-seq in human cancer cell lines as well as mouse embryonic stem cells. We find that roughly 60% of transcription factors tagged using CETCh-seq produce a high quality ChIP-seq map, a significant improvement over traditional antibody-based methods.

Decoding transcriptional enhancers: Evolving from annotation to functional interpretation

Deciphering the intricate molecular processes that orchestrate the spatial and temporal regulation of genes has become an increasingly major focus of biological research. The differential expression of genes by diverse cell types with a common genome is a hallmark of complex cellular functions, as well as the basis for multicellular life. Importantly, a more coherent understanding of gene regulation is critical for defining developmental processes, evolutionary principles and disease etiologies. Here we present our current understanding of gene regulation by focusing on the role of enhancer elements in these complex processes. Although functional genomic methods have provided considerable advances to our understanding of gene regulation, these assays, which are usually performed on a genome-wide scale, typically provide correlative observations that lack functional interpretation. Recent innovations in genome editing technologies have placed gene regulatory studies at an exciting crossroads, as systematic, functional evaluation of enhancers and other transcriptional regulatory elements can now be performed in a coordinated, high-throughput manner across the entire genome. This review provides insights on transcriptional enhancer function, their role in development and disease, and catalogues experimental tools commonly used to study these elements. Additionally, we discuss the crucial role of novel techniques in deciphering the complex gene regulatory landscape and how these studies will shape future research.

The Coronary Artery Disease-associated Coding Variant in Zinc Finger C3HC-type Containing 1 (ZC3HC1) Affects Cell Cycle Regulation

Genome-wide association studies have to date identified multiple coronary artery disease (CAD)-associated loci; however, for most of these loci the mechanism by which they affect CAD risk is unclear. The CAD-associated locus 7q32.2 is unusual in that the lead variant, rs11556924, is not in strong linkage disequilibrium with any other variant and introduces a coding change in ZC3HC1, which encodes NIPA. In this study, we show that rs11556924 polymorphism is associated with lower regulatory phosphorylation of NIPA in the risk variant, resulting in NIPA with higher activity. Using a genome-editing approach we show that this causes an effective decrease in cyclin-B1 stability in the nucleus, thereby slowing its nuclear accumulation. By perturbing the rate of nuclear cyclin-B1 accumulation, rs11556924 alters the regulation of mitotic progression resulting in an extended mitosis. This study shows that the CAD-associated coding polymorphism in ZC3HC1 alters the dynamics of cell-cycle regulation by NIPA.

Association between single nucleotide polymorphism in miR-499, miR-196a2, miR-146a and miR-149 and prostate cancer risk in a sample of Iranian population

MicroRNAs (miRNAs) play an important role in regulating gene expression at the post-transcriptional level and are involved in numerous physiological processes. Accumulating evidence suggests that single-nucleotide polymorphisms (SNPs) in human miRNA genes may affect miRNA biogenesis pathway and influence the susceptibility to several diseases such as cancer. The present study aimed to evaluate the impact of miR-499 rs3746444, miR-196a2 rs11614913, miR-149 rs2292832, and miR-146a rs2910164 polymorphisms on prostate cancer (PCa) risk in a sample of Iranian population. This case-control study was done on 169 patients with pathologically confirmed PCa and 182 benign prostatic hyperplasia (BPH). The genotyping assays were done using T-ARMS-PCR or PCR-RFLP methods. The findings indicated that CC genotype of miR-499 rs3746444 polymorphism increased the risk of PCa (OR = 1.76, 95% CI = 1.12-2.79, P = 0.019) compared to TT genotype. No statistically significant association was found between miR-196a2 rs11614913, miR-149 rs2292832, and miR-146a rs2910164 polymorphisms and PCa risk. In summary, the findings indicated that miR-499 rs3746444 polymorphism increased the risk of PCa in an Iranian population. Further studies with larger sample sizes and different ethnicities are necessary to verify the findings of the present study.

NIH working group report-using genomic information to guide weight management: From universal to precision treatment

OBJECTIVE

Precision medicine utilizes genomic and other data to optimize and personalize treatment. Although more than 2,500 genetic tests are currently available, largely for extreme and/or rare phenotypes, the question remains whether this approach can be used for the treatment of common, complex conditions like obesity, inflammation, and insulin resistance, which underlie a host of metabolic diseases.

METHODS

This review, developed from a Trans-NIH Conference titled "Genes, Behaviors, and Response to Weight Loss Interventions," provides an overview of the state of genetic and genomic research in the area of weight change and identifies key areas for future research.

RESULTS

Although many loci have been identified that are associated with cross-sectional measures of obesity/body size, relatively little is known regarding the genes/loci that influence dynamic measures of weight change over time. Although successful short-term weight loss has been achieved using many different strategies, sustainable weight loss has proven elusive for many, and there are important gaps in our understanding of energy balance regulation.

CONCLUSIONS

Elucidating the molecular basis of variability in weight change has the potential to improve treatment outcomes and inform innovative approaches that can simultaneously take into account information from genomic and other sources in devising individualized treatment plans.