Genome-scale activation screen identifies a lncRNA locus regulating a gene neighbourhood.

Mammalian genomes contain thousands of loci that transcribe long noncoding RNAs (lncRNAs), some of which are known to carry out critical roles in diverse cellular processes through a variety of mechanisms. Although some lncRNA loci encode RNAs that act non-locally (in trans), there is emerging evidence that many lncRNA loci act locally (in cis) to regulate the expression of nearby genes-for example, through functions of the lncRNA promoter, transcription, or transcript itself. Despite their potentially important roles, it remains challenging to identify functional lncRNA loci and distinguish among these and other mechanisms. Here, to address these challenges, we developed a genome-scale CRISPR-Cas9 activation screen that targets more than 10,000 lncRNA transcriptional start sites to identify noncoding loci that influence a phenotype of interest. We found 11 lncRNA loci that, upon recruitment of an activator, mediate resistance to BRAF inhibitors in human melanoma cells. Most candidate loci appear to regulate nearby genes. Detailed analysis of one candidate, termed EMICERI, revealed that its transcriptional activation resulted in dosage-dependent activation of four neighbouring protein-coding genes, one of which confers the resistance phenotype. Our screening and characterization approach provides a CRISPR toolkit with which to systematically discover the functions of noncoding loci and elucidate their diverse roles in gene regulation and cellular function.

Erratum: Genome-scale activation screen identifies a lncRNA locus regulating a gene neighbourhood.

This corrects the article DOI: 10.1038/nature23451.

CRISPR Screens to Discover Functional Noncoding Elements.

A major challenge in genomics is to identify functional elements in the noncoding genome. Recently, pooled clustered regularly interspersed palindromic repeat (CRISPR) mutagenesis screens of noncoding regions have emerged as a novel method for finding elements that impact gene expression and phenotype/disease-relevant biological processes. Here we review and compare different approaches for high-throughput dissection of noncoding elements.

Myc protein is stabilized by suppression of a novel E3 ligase complex in cancer cells.

Rapid Myc protein turnover is critical for maintaining basal levels of Myc activity in normal cells and a prompt response to changing growth signals. We characterize a new Myc-interacting factor, TRPC4AP (transient receptor potential cation channel, subfamily C, member 4-associated protein)/TRUSS (tumor necrosis factor receptor-associated ubiquitous scaffolding and signaling protein), which is the receptor for a DDB1 (damage-specific DNA-binding protein 1)-CUL4 (Cullin 4) E3 ligase complex for selective Myc degradation through the proteasome. TRPC4AP/TRUSS binds specifically to the Myc C terminus and promotes its ubiquitination and destruction through the recognition of evolutionarily conserved domains in the Myc N terminus. TRPC4AP/TRUSS suppresses Myc-mediated transactivation and transformation in a dose-dependent manner. Finally, we found that TRPC4AP/TRUSS expression is strongly down-regulated in most cancer cell lines, leading to Myc protein stabilization. These studies identify a novel pathway targeting Myc degradation that is suppressed in cancer cells.

Distal interleukin-1ß (IL-1ß) response element of human matrix metalloproteinase-13 (MMP-13) binds activator protein 1 (AP-1) transcription factors and regulates gene expression.

The collagenase matrix metalloproteinase-13 (MMP-13) plays an important role in the destruction of cartilage in arthritic joints. MMP-13 expression is strongly up-regulated in arthritis, largely because of stimulation by inflammatory cytokines such as IL-1ß. Treatment of chondrocytes with IL-1ß induces transcription of MMP-13 in vitro. IL-1ß signaling converges upon the activator protein-1 transcription factors, which have been shown to be required for IL-1ß-induced MMP-13 gene expression. Using chromatin immunoprecipitation (ChIP), we detected activator protein-1 binding within an evolutionarily conserved DNA sequence ∼20 kb 5' relative to the MMP-13 transcription start site (TSS). Also using ChIP, we detected histone modifications and binding of RNA polymerase II within this conserved region, all of which are consistent with transcriptional activation. Chromosome conformation capture indicates that chromosome looping brings this region in close proximity with the MMP-13 TSS. Finally, a luciferase reporter construct driven by a component of the conserved region demonstrated an expression pattern similar to that of endogenous MMP-13. These data suggest that a conserved region at 20 kb upstream from the MMP-13 TSS includes a distal transcriptional response element of MMP-13, which contributes to MMP-13 gene expression.

Molecular characterization of precise in vivo targeted gene integration in human cells using AAVHSC15.

Targeted gene integration via precise homologous recombination (HR)-based gene editing has the potential to correct genetic diseases. AAV (adeno-associated virus) can mediate nuclease-free gene integration at a disease-causing locus. Therapeutic application of AAV gene integration requires quantitative molecular characterization of the edited sequence that overcome technical obstacles such as excess episomal vector genomes and lengthy homology arms. Here we describe a novel molecular methodology that utilizes quantitative next-generation sequencing to characterize AAV-mediated targeted insertion and detects the presence of unintended mutations. The methods described here quantify targeted insertion and query the entirety of the target locus for the presence of insertions, deletions, single nucleotide variants (SNVs) and integration of viral components such as inverted terminal repeats (ITR). Using a humanized liver murine model, we demonstrate that hematopoietic stem-cell derived AAVHSC15 mediates in vivo targeted gene integration into human chromosome 12 at the PAH (phenylalanine hydroxylase) locus at 6% frequency, with no sign of co-incident random mutations at or above a lower limit of detection of 0.5% and no ITR sequences at the integration sites. Furthermore, analysis of heterozygous variants across the targeted locus using the methods described shows a pattern of strand cross-over, supportive of an HR mechanism of gene integration with similar efficiencies across two different haplotypes. Rapid advances in the application of AAV-mediated nuclease-free target integration, or gene editing, as a new therapeutic modality requires precise understanding of the efficiency and the nature of the changes being introduced to the target genome at the molecular level. This work provides a framework to be applied to homologous recombination gene editing platforms for assessment of introduced and natural sequence variation across a target site.

Retroviruses hijack chromatin loops to drive oncogene expression and highlight the chromatin architecture around proto-oncogenic loci.

The majority of the genome consists of intergenic and non-coding DNA sequences shown to play a major role in different gene regulatory networks. However, the specific potency of these distal elements as well as how these regions exert function across large genomic distances remains unclear. To address these unresolved issues, we closely examined the chromatin architecture around proto-oncogenic loci in the mouse and human genomes to demonstrate a functional role for chromatin looping in distal gene regulation. Using cell culture models, we show that tumorigenic retroviral integration sites within the mouse genome occur near existing large chromatin loops and that this chromatin architecture is maintained within the human genome as well. Significantly, as mutagenesis screens are not feasible in humans, we demonstrate a way to leverage existing screens in mice to identify disease relevant human enhancers and expose novel disease mechanisms. For instance, we characterize the epigenetic landscape upstream of the human Cyclin D1 locus to find multiple distal interactions that contribute to the complex cis-regulation of this cell cycle gene. Furthermore, we characterize a novel distal interaction upstream of the Cyclin D1 gene which provides mechanistic evidence for the abundant overexpression of Cyclin D1 occurring in multiple myeloma cells harboring a pathogenic translocation event. Through use of mapped retroviral integrations and translocation breakpoints, our studies highlight the importance of chromatin looping in oncogene expression, elucidate the epigenetic mechanisms crucial for distal cis-regulation, and in one particular instance, explain how a translocation event drives tumorigenesis through upregulation of a proto-oncogene.

Breast cancer risk-associated SNPs modulate the affinity of chromatin for FOXA1 and alter gene expression.

Genome-wide association studies (GWAS) have identified thousands of SNPs that are associated with human traits and diseases. But, because the vast majority of these SNPs are located in non-coding regions of the genome, the mechanisms by which they promote disease risk have remained elusive. Employing a new methodology that combines cistromics, epigenomics and genotype imputation, we annotate the non-coding regions of the genome in breast cancer cells and systematically identify the functional nature of SNPs associated with breast cancer risk. Our results show that breast cancer risk-associated SNPs are enriched in the cistromes of FOXA1 and ESR1 and the epigenome of histone H3 lysine 4 monomethylation (H3K4me1) in a cancer- and cell type-specific manner. Furthermore, the majority of the risk-associated SNPs modulate the affinity of chromatin for FOXA1 at distal regulatory elements, thereby resulting in allele-specific gene expression, which is exemplified by the effect of the rs4784227 SNP on the TOX3 gene within the 16q12.1 risk locus.

Upregulation of c-MYC in cis through a large chromatin loop linked to a cancer risk-associated single-nucleotide polymorphism in colorectal cancer cells.

Genome-wide association studies have mapped many single-nucleotide polymorphisms (SNPs) that are linked to cancer risk, but the mechanism by which most SNPs promote cancer remains undefined. The rs6983267 SNP at 8q24 has been associated with many cancers, yet the SNP falls 335 kb from the nearest gene, c-MYC. We show that the beta-catenin-TCF4 transcription factor complex binds preferentially to the cancer risk-associated rs6983267(G) allele in colon cancer cells. We also show that the rs6983267 SNP has enhancer-related histone marks and can form a 335-kb chromatin loop to interact with the c-MYC promoter. Finally, we show that the SNP has no effect on the efficiency of chromatin looping to the c-MYC promoter but that the cancer risk-associated SNP enhances the expression of the linked c-MYC allele. Thus, cancer risk is a direct consequence of elevated c-MYC expression from increased distal enhancer activity and not from reorganization/creation of the large chromatin loop. The findings of these studies support a mechanism for intergenic SNPs that can promote cancer through the regulation of distal genes by utilizing preexisting large chromatin loops.