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

Effects on the transcriptome upon deletion of a distal element cannot be predicted by the size of the H3K27Ac peak in human cells.

Genome-wide association studies (GWAS) have identified single nucleotide polymorphisms (SNPs) associated with increased risk for colorectal cancer (CRC). A molecular understanding of the functional consequences of this genetic variation is complicated because most GWAS SNPs are located in non-coding regions. We used epigenomic information to identify H3K27Ac peaks in HCT116 colon cancer cells that harbor SNPs associated with an increased risk for CRC. Employing CRISPR/Cas9 nucleases, we deleted a CRC risk-associated H3K27Ac peak from HCT116 cells and observed large-scale changes in gene expression, resulting in decreased expression of many nearby genes. As a comparison, we showed that deletion of a robust H3K27Ac peak not associated with CRC had minimal effects on the transcriptome. Interestingly, although there is no H3K27Ac peak in HEK293 cells in the E7 region, deletion of this region in HEK293 cells decreased expression of several of the same genes that were downregulated in HCT116 cells, including the MYC oncogene. Accordingly, deletion of E7 causes changes in cell culture assays in HCT116 and HEK293 cells. In summary, we show that effects on the transcriptome upon deletion of a distal regulatory element cannot be predicted by the size or presence of an H3K27Ac peak

dCas9-based epigenome editing suggests acquisition of histone methylation is not sufficient for target gene repression.

Distinct epigenomic profiles of histone marks have been associated with gene expression, but questions regarding the causal relationship remain. Here we investigated the activity of a broad collection of genomically targeted epigenetic regulators that could write epigenetic marks associated with a repressed chromatin state (G9A, SUV39H1, Krüppel-associated box (KRAB), DNMT3A as well as the first targetable versions of Ezh2 and Friend of GATA-1 (FOG1)). dCas9 fusions produced target gene repression over a range of 0- to 10-fold that varied by locus and cell type. dCpf1 fusions were unable to repress gene expression. The most persistent gene repression required the action of several effector domains; however, KRAB-dCas9 did not contribute to persistence in contrast to previous reports. A 'direct tethering' strategy attaching the Ezh2 methyltransferase enzyme to dCas9, as well as a 'recruitment' strategy attaching the N-terminal 45 residues of FOG1 to dCas9 to recruit the endogenous nucleosome remodeling and deacetylase complex, were both successful in targeted deposition of H3K27me3. Surprisingly, however, repression was not correlated with deposition of either H3K9me3 or H3K27me3. Our results suggest that so-called repressive histone modifications are not sufficient for gene repression. The easily programmable dCas9 toolkit allowed precise control of epigenetic information and dissection of the relationship between the epigenome and gene regulation

Ezh2-dCas9 and KRAB-dCas9 enable engineering of epigenetic memory in a context-dependent manner.

BackgroundRewriting of the epigenome has risen as a promising alternative to gene editing for precision medicine. In nature, epigenetic silencing can result in complete attenuation of target gene expression over multiple mitotic divisions. However, persistent repression has been difficult to achieve in a predictable manner using targeted systems.ResultsHere, we report that persistent epigenetic memory required both a DNA methyltransferase (DNMT3A-dCas9) and a histone methyltransferase (Ezh2-dCas9 or KRAB-dCas9). We demonstrate that the histone methyltransferase requirement can be locus specific. Co-targeting Ezh2-dCas9, but not KRAB-dCas9, with DNMT3A-dCas9 and DNMT3L induced long-term HER2 repression over at least 50 days (approximately 57 cell divisions) and triggered an epigenetic switch to a heterochromatic environment. An increase in H3K27 trimethylation and DNA methylation was stably maintained and accompanied by a sustained loss of H3K27 acetylation. Interestingly, substitution of Ezh2-dCas9 with KRAB-dCas9 enabled long-term repression at some target genes (e.g., SNURF) but not at HER2, at which H3K9me3 and DNA methylation were transiently acquired and subsequently lost. Off-target DNA hypermethylation occurred at many individual CpG sites but rarely at multiple CpGs in a single promoter, consistent with no detectable effect on transcription at the off-target loci tested. Conversely, robust hypermethylation was observed at HER2. We further demonstrated that Ezh2-dCas9 required full-length DNMT3L for maximal activity and that co-targeting DNMT3L was sufficient for persistent repression by Ezh2-dCas9 or KRAB-dCas9.ConclusionsThese data demonstrate that targeting different combinations of histone and DNA methyltransferases is required to achieve maximal repression at different loci. Fine-tuning of targeting tools is a necessity to engineer epigenetic memory at any given locus in any given cell type

Protein Delivery of an Artificial Transcription Factor Restores Widespread Ube3a Expression in an Angelman Syndrome Mouse Brain.

Angelman syndrome (AS) is a neurological genetic disorder caused by loss of expression of the maternal copy of UBE3A in the brain. Due to brain-specific genetic imprinting at this locus, the paternal UBE3A is silenced by a long antisense transcript. Inhibition of the antisense transcript could lead to unsilencing of paternal UBE3A, thus providing a therapeutic approach for AS. However, widespread delivery of gene regulators to the brain remains challenging. Here, we report an engineered zinc finger-based artificial transcription factor (ATF) that, when injected i.p. or s.c., crossed the blood-brain barrier and increased Ube3a expression in the brain of an adult mouse model of AS. The factor displayed widespread distribution throughout the brain. Immunohistochemistry of both the hippocampus and cerebellum revealed an increase in Ube3a upon treatment. An ATF containing an alternative DNA-binding domain did not activate Ube3a. We believe this to be the first report of an injectable engineered zinc finger protein that can cause widespread activation of an endogenous gene in the brain. These observations have important implications for the study and treatment of AS and other neurological disorders

Genome-Wide Analysis of KAP1 Binding Suggests Autoregulation of KRAB-ZNFs

We performed a genome-scale chromatin immunoprecipitation (ChIP)-chip comparison of two modifications (trimethylation of lysine 9 [H3me3K9] and trimethylation of lysine 27 [H3me3K27]) of histone H3 in Ntera2 testicular carcinoma cells and in three different anatomical sources of primary human fibroblasts. We found that in each of the cell types the two modifications were differentially enriched at the promoters of the two largest classes of transcription factors. Specifically, zinc finger (ZNF) genes were bound by H3me3K9 and homeobox genes were bound by H3me3K27. We have previously shown that the Polycomb repressive complex 2 is responsible for mediating trimethylation of lysine 27 of histone H3 in human cancer cells. In contrast, there is little overlap between H3me3K9 targets and components of the Polycomb repressive complex 2, suggesting that a different histone methyltransferase is responsible for the H3me3K9 modification. Previous studies have shown that SETDB1 can trimethylate H3 on lysine 9, using in vitro or artificial tethering assays. SETDB1 is thought to be recruited to chromatin by complexes containing the KAP1 corepressor. To determine if a KAP1-containing complex mediates trimethylation of the identified H3me3K9 targets, we performed ChIP-chip assays and identified KAP1 target genes using human 5-kb promoter arrays. We found that a large number of genes of ZNF transcription factors were bound by both KAP1 and H3me3K9 in normal and cancer cells. To expand our studies of KAP1, we next performed a complete genomic analysis of KAP1 binding using a 38-array tiling set, identifying ~7,000 KAP1 binding sites. The identified KAP1 targets were highly enriched for C2H2 ZNFs, especially those containing Krüppel-associated box (KRAB) domains. Interestingly, although most KAP1 binding sites were within core promoter regions, the binding sites near ZNF genes were greatly enriched within transcribed regions of the target genes. Because KAP1 is recruited to the DNA via interaction with KRAB-ZNF proteins, we suggest that expression of KRAB-ZNF genes may be controlled via an auto-regulatory mechanism involving KAP1

Artificial escape from XCI by DNA methylation editing of the CDKL5 gene.

A significant number of X-linked genes escape from X chromosome inactivation and are associated with a distinct epigenetic signature. One epigenetic modification that strongly correlates with X-escape is reduced DNA methylation in promoter regions. Here, we created an artificial escape by editing DNA methylation on the promoter of CDKL5, a gene causative for an infantile epilepsy, from the silenced X-chromosomal allele in human neuronal-like cells. We identify that a fusion of the catalytic domain of TET1 to dCas9 targeted to the CDKL5 promoter using three guide RNAs causes significant reactivation of the inactive allele in combination with removal of methyl groups from CpG dinucleotides. Strikingly, we demonstrate that co-expression of TET1 and a VP64 transactivator have a synergistic effect on the reactivation of the inactive allele to levels >60% of the active allele. We further used a multi-omics assessment to determine potential off-targets on the transcriptome and methylome. We find that synergistic delivery of dCas9 effectors is highly selective for the target site. Our findings further elucidate a causal role for reduced DNA methylation associated with escape from X chromosome inactivation. Understanding the epigenetics associated with escape from X chromosome inactivation has potential for those suffering from X-linked disorders

Sole-Search: an integrated analysis program for peak detection and functional annotation using ChIP-seq data

Next-generation sequencing is revolutionizing the identification of transcription factor binding sites throughout the human genome. However, the bioinformatics analysis of large datasets collected using chromatin immunoprecipitation and high-throughput sequencing is often a roadblock that impedes researchers in their attempts to gain biological insights from their experiments. We have developed integrated peak-calling and analysis software (Sole-Search) which is available through a user-friendly interface and (i) converts raw data into a format for visualization on a genome browser, (ii) outputs ranked peak locations using a statistically based method that overcomes the significant problem of false positives, (iii) identifies the gene nearest to each peak, (iv) classifies the location of each peak relative to gene structure, (v) provides information such as the number of binding sites per chromosome and per gene and (vi) allows the user to determine overlap between two different experiments. In addition, the program performs an analysis of amplified and deleted regions of the input genome. This software is web-based and automated, allowing easy and immediate access to all investigators. We demonstrate the utility of our software by collecting, analyzing and comparing ChIP-seq data for six different human transcription factors/cell line combinations