Thousands Of Genomic Loci Research Articles

ImputeHiFI: An Imputation Method for Multiplexed DNA FISH Data by Utilizing Single-Cell Hi-C and RNA FISH Data.

Although multiplexed DNA fluorescence in situ hybridization (FISH) enables tracking the spatial localization of thousands of genomic loci using probes within individual cells, the high rates of undetected probes impede the depiction of 3D chromosome structures. Current data imputation methods neither utilize single-cell Hi-C data, which elucidate 3D genome architectures using sequencing nor leverage multimodal RNA FISH data that reflect cell-type information, limiting the effectiveness of these methods in complex tissues such as the mouse brain. To this end, a novel multiplexed DNA FISH imputation method named ImputeHiFI is proposed, which fully utilizes the complementary structural information from single-cell Hi-C data and the cell type signature from RNA FISH data to obtain a high-fidelity and complete spatial location of chromatin loci. ImputeHiFI enhances cell clustering, compartment identification, and cell subtype detection at the single-cell level in the mouse brain. ImputeHiFI improves the recognition of cell-type-specific loops in three high-resolution datasets. In short, ImputeHiFI is a powerful tool capable of imputing multiplexed DNA FISH data from various resolutions and imaging protocols, facilitating studies of 3D genome structures and functions.

Regulatory landscape enrichment analysis (RLEA): a computational toolkit for non-coding variant enrichment and cell type prioritization

BackgroundAs genomic studies continue to implicate non-coding sequences in disease, testing the roles of these variants requires insights into the cell type(s) in which they are likely to be mediating their effects. Prior methods for associating non-coding variants with cell types have involved approaches using linkage disequilibrium or ontological associations, incurring significant processing requirements. GaiaAssociation is a freely available, open-source software that enables thousands of genomic loci implicated in a phenotype to be tested for enrichment at regulatory loci of multiple cell types in minutes, permitting insights into the cell type(s) mediating the studied phenotype.ResultsIn this work, we present Regulatory Landscape Enrichment Analysis (RLEA) by GaiaAssociation and demonstrate its capability to test the enrichment of 12,133 variants across the cis-regulatory regions of 44 cell types. This analysis was completed in 134.0 ± 2.3 s, highlighting the efficient processing provided by GaiaAssociation. The intuitive interface requires only four inputs, offers a collection of customizable functions, and visualizes variant enrichment in cell-type regulatory regions through a heatmap matrix. GaiaAssociation is available on PyPi for download as a command line tool or Python package and the source code can also be installed from GitHub at https://github.com/GreallyLab/gaiaAssociation.ConclusionsGaiaAssociation is a novel package that provides an intuitive and efficient resource to understand the enrichment of non-coding variants across the cis-regulatory regions of different cells, empowering studies seeking to identify disease-mediating cell types.

DIPG-45. ELONGIN B REGULATES H3K27M CHROMATIN INCORPORATION TO PROMOTE MALIGNANCY IN DIFFUSE MIDLINE GLIOMA

Abstract Chromatin dysregulation is both a cause and consequence of aberrant RNA Polymerase II (Pol2) transcription in cancer. Some of the most striking examples of chromatin disruption in cancer include somatic mutations in genes encoding histone H3 (H3.3/H3.1K27M) observed in ~80% of diffuse midline glioma (DMG). H3K27M mutant histones act in a dominant manner to disrupt chromatin regulation and Pol2-dependent transcription in DMG. How H3K27M impinges on the Pol2 transcriptional machinery and whether aberrant Pol2 transcription can reciprocally induce chromatin dysregulation in DMG remain open questions. We conducted epigenome-focused CRISPR dropout screens and identified multiple regulators of Pol2 elongation, including Elongin B and C (ELOB/ELOC) as DMG genetic dependencies. Follow-up studies confirm that knockout (KO) of ELOB inhibits DMG proliferation in neurosphere culture and abrogates DMG tumorigenesis in orthotopic xenografts. Additional chromatin profiling studies reveal that ELOB binding sites are enriched in H3K27M mutant histones, suggesting that ELOB works together with H3K27M to promote malignant gene expression in DMG. Indeed, reciprocal knockout studies suggest that while ELOB chromatin binding patterns were not altered by H3K27M KO, ELOB KO significantly reduced H3K27M chromatin incorporation at thousands of genomic loci. Further RNAseq and PROseq analyses reveal that ELOB KO alters Pol2 transcription resulting in the up-regulation of genes associated with neuronal and muscle differentiation and the repression of glial stem cell genes. Together, these findings suggest that H3K27M chromatin binding patterns and aberrant gene regulation in DMG are, in part, dictated by ELOB and the Pol2 elongation machinery. Finally, analyses of patient datasets indicate that ELOB complex genes are overexpressed in high grade glioma and DMG and correlate with decreased patient survival. Collectively, these studies suggest that a functional interaction between ELOB and H3K27M mutant histones induces DMG-specific aberrations in chromatin and gene regulation to drive malignancy.

Loss of Ezh2 function remodels the DNA replication initiation landscape

In metazoan cells, DNA replication initiates from thousands of genomic loci scattered throughout the genome called DNA replication origins. Origins are strongly associated with euchromatin, particularly open genomic regions such as promoters and enhancers. However, over a third of transcriptionally silent genes are associated with DNA replication initiation. Most of these genes are bound and repressed by the Polycomb repressive complex-2 (PRC2) through the repressive H3K27me3 mark. This is the strongest overlap observed for a chromatin regulator with replication origin activity. Here, we asked whether Polycomb-mediated gene repression is functionally involved in recruiting DNA replication origins to transcriptionally silent genes. We show that the absence of EZH2, the catalytic subunit of PRC2, results in increased DNA replication initiation, specifically in the vicinity of EZH2 binding sites. The increase in DNA replication initiation does not correlate with transcriptional de-repression or the acquisition of activating histone marks but does correlate with loss of H3K27me3 from bivalent promoters.

Hacking hematopoiesis - emerging tools for examining variant effects.

Hematopoiesis is a continuous process of blood and immune cell production. It is orchestrated by thousands of gene products that respond to extracellular signals by guiding cell fate decisions to meet the needs of the organism. Although much of our knowledge of this process comes from work in model systems, we have learned a great deal from studies on human genetic variation. Considerable insight has emerged from studies on presumed monogenic blood disorders, which continue to provide key insights into the mechanisms critical for hematopoiesis. Furthermore, the emergence of large-scale biobanks and cohorts has uncovered thousands of genomic loci associated with blood cell traits and diseases. Some of these blood cell trait-associated loci act as modifiers of what were once thought to be monogenic blood diseases. However, most of these loci await functional validation. Here, we discuss the validation bottleneck and emerging methods to more effectively connect variant to function. In particular, we highlight recent innovations in genome editing, which have paved the path forward for high-throughput functional assessment of loci. Finally, we discuss existing barriers to progress, including challenges in manipulating the genomes of primary hematopoietic cells.

The H3K9me2-binding protein AGDP3 limits DNA methylation and transcriptional gene silencing in Arabidopsis.

DNA methylation, a conserved epigenetic mark, is critical for tuning temporal and spatial gene expression. The Arabidopsis thaliana DNA glycosylase/lyase REPRESSOR OF SILENCING 1 (ROS1) initiates active DNA demethylation and is required to prevent DNA hypermethylation at thousands of genomic loci. However, how ROS1 is recruited to specific loci is not well understood. Here, we report the discovery of Arabidopsis AGENET Domain Containing Protein 3 (AGDP3) as a cellular factor that is required to prevent gene silencing and DNA hypermethylation. AGDP3 binds to H3K9me2 marks in its target DNA via its AGD12 cassette. Analysis of the crystal structure of the AGD12 cassette of AGDP3 in complex with an H3K9me2 peptide revealed that dimethylated H3K9 and unmodified H3K4 are specifically anchored into two different surface pockets. A histidine residue located in the methyllysine binding aromatic cage provides AGDP3 with pH-dependent H3K9me2 binding capacity. Our results uncover a molecular mechanism for the regulation of DNA demethylation by the gene silencing mark H3K9me2.

A quantitative metric of pioneer activity reveals that HNF4A has stronger in vivo pioneer activity than FOXA1

BackgroundWe and others have suggested that pioneer activity — a transcription factor’s (TF’s) ability to bind and open inaccessible loci — is not a qualitative trait limited to a select class of pioneer TFs. We hypothesize that most TFs display pioneering activity that depends on the TF concentration and the motif content at their target loci.ResultsHere, we present a quantitative in vivo measure of pioneer activity that captures the relative difference in a TF’s ability to bind accessible versus inaccessible DNA. The metric is based on experiments that use CUT&Tag to measure the binding of doxycycline-inducible TFs. For each location across the genome, we determine the concentration of doxycycline required for a TF to reach half-maximal occupancy; lower concentrations reflect higher affinity. We propose that the relative difference in a TF’s affinity between ATAC-seq labeled accessible and inaccessible binding sites is a measure of its pioneer activity. We estimate binding affinities at tens of thousands of genomic loci for the endodermal TFs FOXA1 and HNF4A and show that HNF4A has stronger pioneer activity than FOXA1. We show that both FOXA1 and HNF4A display higher binding affinity at inaccessible sites with more copies of their respective motifs. The quantitative analysis of binding suggests different modes of binding for FOXA1, including an anti-cooperative mode of binding at certain accessible loci.ConclusionsOur results suggest that relative binding affinities are reasonable measures of pioneer activity and support the model wherein most TFs have some degree of context-dependent pioneer activity.

Hybridization decreases native cutthroat trout reproductive fitness.

Examining natural selection in wild populations is challenging, but crucial to understanding many ecological and evolutionary processes. Additionally, in hybridizing populations, natural selection may be an important determinant of the eventual outcome of hybridization. We characterized several components of relative fitness in hybridizing populations of Yellowstone cutthroat trout and rainbow trout in an effort to better understand the prolonged persistence of both parental species despite predictions of extirpation. Thousands of genomic loci enabled precise quantification of hybrid status in adult and subsequent juvenile generations; a subset of those data also identified parent-offspring relationships. We used linear models and simulations to assess the effects of ancestry on reproductive output and mate choice decisions. We found a relatively low number of late-stage (F3+) hybrids and an excess of F2 juveniles relative to the adult generation in one location, which suggests the presence of hybrid breakdown decreasing the fitness of F2+ hybrids later in life. Assessments of reproductive output showed that Yellowstone cutthroat trout are more likely to successfully reproduce and produce slightly more offspring than their rainbow trout and hybrid counterparts. Mate choice appeared to be largely random, though we did find statistical support for slight female preference for males of similar ancestry. Together, these results show that native Yellowstone cutthroat trout are able to outperform rainbow trout in terms of reproduction and suggest that management action to exclude rainbow trout from spawning locations may bolster the now-rare Yellowstone cutthroat trout.

EzQTL: A Web Platform for Interactive Visualization and Colocalization of QTLs and GWAS Loci

Genome-wide association studies (GWAS) have identified thousands of genomic loci associated with complex diseases and traits, including cancer. The vast majority of common trait-associated variants identified via GWAS fall in non-coding regions of the genome, posing a challenge in elucidating the causal variants, genes, and mechanisms involved. Expression quantitative trait locus (eQTL) and other molecular QTL studies have been valuable resources in identifying candidate causal genes from GWAS loci through statistical colocalization methods. While QTL colocalization is becoming a standard analysis in post-GWAS investigation, an easy web tool for users to perform formal colocalization analyses with either user-provided or public GWAS and eQTL datasets has been lacking. Here, we present ezQTL, a web-based bioinformatic application to interactively visualize and analyze genetic association data such as GWAS loci and molecular QTLs under different linkage disequilibrium (LD) patterns (1000 Genomes Project, UK Biobank, or user-provided data). This application allows users to perform data quality control for variants matched between different datasets, LD visualization, and two-trait colocalization analyses using two state-of-the-art methodologies (eCAVIAR and HyPrColoc), including batch processing. ezQTL is a free and publicly available cross-platform web tool, which can be accessed online at https://analysistools.cancer.gov/ezqtl.

Vitamin D and Its Target Genes.

The vitamin D metabolite 1α,25-dihydroxyvitamin D3 is the natural, high-affinity ligand of the transcription factor vitamin D receptor (VDR). In many tissues and cell types, VDR binds in a ligand-dependent fashion to thousands of genomic loci and modulates, via local chromatin changes, the expression of hundreds of primary target genes. Thus, the epigenome and transcriptome of VDR-expressing cells is directly affected by vitamin D. Vitamin D target genes encode for proteins with a large variety of physiological functions, ranging from the control of calcium homeostasis, innate and adaptive immunity, to cellular differentiation. This review will discuss VDR’s binding to genomic DNA, as well as its genome-wide locations and interaction with partner proteins, in the context of chromatin. This information will be integrated into a model of vitamin D signaling, explaining the regulation of vitamin D target genes.

Phylogenomics, introgression, and demographic history of South American true toads (Rhinella).

The effects of genetic introgression on species boundaries and how they affect species' integrity and persistence over evolutionary time have received increased attention. The increasing availability of genomic data has revealed contrasting patterns of gene flow across genomic regions, which impose challenges to inferences of evolutionary relationships and of patterns of genetic admixture across lineages. By characterizing patterns of variation across thousands of genomic loci in a widespread complex of true toads (Rhinella), we assess the true extent of genetic introgression across species thought to hybridize to extreme degrees based on natural history observations and multilocus analyses. Comprehensive geographic sampling of five large-ranged Neotropical taxa revealed multiple distinct evolutionary lineages that span large geographic areas and, at times, distinct biomes. The inferred major clades and genetic clusters largely correspond to currently recognized taxa; however, we also found evidence of cryptic diversity within taxa. While previous phylogenetic studies revealed extensive mitonuclear discordance, our genetic clustering analyses uncovered several admixed individuals within major genetic groups. Accordingly, historical demographic analyses supported that the evolutionary history of these toads involved cross-taxon gene flow both at ancient and recent times. Lastly, ABBA-BABA tests revealed widespread allele sharing across species boundaries, a pattern that can be confidently attributed to genetic introgression as opposed to incomplete lineage sorting. These results confirm previous assertions that the evolutionary history of Rhinella was characterized by various levels of hybridization even across environmentally heterogeneous regions, posing exciting questions about what factors prevent complete fusion of diverging yet highly interdependent evolutionary lineages.

Coordinated glucocorticoid receptor and MAFB action induces tolerogenesis and epigenome remodeling in dendritic cells.

Glucocorticoids (GCs) exert potent anti-inflammatory effects in immune cells through the glucocorticoid receptor (GR). Dendritic cells (DCs), central actors for coordinating immune responses, acquire tolerogenic properties in response to GCs. Tolerogenic DCs (tolDCs) have emerged as a potential treatment for various inflammatory diseases. To date, the underlying cell type-specific regulatory mechanisms orchestrating GC-mediated acquisition of immunosuppressive properties remain poorly understood. In this study, we investigated the transcriptomic and epigenomic remodeling associated with differentiation to DCs in the presence of GCs. Our analysis demonstrates a major role of MAFB in this process, in synergy with GR. GR and MAFB both interact with methylcytosine dioxygenase TET2 and bind to genomic loci that undergo specific demethylation in tolDCs. We also show that the role of MAFB is more extensive, binding to thousands of genomic loci in tolDCs. Finally, MAFB knockdown erases the tolerogenic properties of tolDCs and reverts the specific DNA demethylation and gene upregulation. The preeminent role of MAFB is also demonstrated in vivo for myeloid cells from synovium in rheumatoid arthritis following GC treatment. Our results imply that, once directly activated by GR, MAFB plays a critical role in orchestrating the epigenomic and transcriptomic remodeling that define the tolerogenic phenotype.

High-Throughput Nucleotide Resolution Predictions of Assay Limitations Increase the Reliability and Concordance of Clinical Tests.

The ability of next-generation sequencing (NGS) assays to interrogate thousands of genomic loci has revolutionized genetic testing. However, translation to the clinic is impeded by false-negative results that pose a risk to patients. In response, regulatory bodies are calling for reliability measures to be reported alongside NGS results. Existing methods to estimate reliability do not account for sample- and position-specific variability, which can be significant. Here, we report an approach that computes reliability metrics for every genomic position and sample interrogated by an NGS assay. Our approach predicts the limit of detection (LOD), the lowest reliably detectable variant fraction, by taking technical factors into account. We initially explored how LOD is affected by input material amount, library conversion rate, sequencing coverage, and sequencing error rate. This revealed that LOD depends heavily on genomic context and sample properties. Using these insights, we developed a computational approach to predict LOD on the basis of a biophysical model of the NGS workflow. We focused on targeted assays for cell-free DNA, but, in principle, this approach applies to any NGS assay. We validated our approach by showing that it accurately predicts LOD and distinguishes reliable from unreliable results when screening 580 lung cancer samples for actionable mutations. Compared with a standard variant calling workflow, our approach avoided most false negatives and improved interassay concordance from 94% to 99%. Our approach, which we name LAVA (LOD-aware variant analysis), reports the LOD for every position and sample interrogated by an NGS assay. This enables reliable results to be identified and improves the transparency and safety of genetic tests.

A Shift in Paradigms: Spatial Genomics Approaches to Reveal Single-Cell Principles of Genome Organization.

The genome tridimensional (3D) organization and its role towards the regulation of key cell processes such as transcription is currently a main question in biology. Interphase chromosomes are spatially segregated into “territories,” epigenetically-defined large domains of chromatin that interact to form “compartments” with common transcriptional status, and insulator-flanked domains called “topologically associating domains” (TADs). Moreover, chromatin organizes around nuclear structures such as lamina, speckles, or the nucleolus to acquire a higher-order genome organization. Due to recent technological advances, the different hierarchies are being solved. Particularly, advances in microscopy technologies are shedding light on the genome structure at multiple levels. Intriguingly, more and more reports point to high variability and stochasticity at the single-cell level. However, the functional consequences of such variability in genome conformation are still unsolved. Here, I will discuss the implication of the cell-to-cell heterogeneity at the different scales in the context of newly developed imaging approaches, particularly multiplexed Fluorescence in situ hybridization methods that enabled “chromatin tracing.” Extensions of these methods are now combining spatial information of dozens to thousands of genomic loci with the localization of nuclear features such as the nucleolus, nuclear speckles, or even histone modifications, creating the fast-moving field of “spatial genomics.” As our view of genome organization shifts the focus from ensemble to single-cell, new insights to fundamental questions begin to emerge.

Comparison of seven single cell whole genome amplification commercial kits using targeted sequencing

Advances in whole genome amplification (WGA) techniques enable understanding of the genomic sequence at a single cell level. Demand for single cell dedicated WGA kits (scWGA) has led to the development of several commercial kit. To this point, no robust comparison of all available kits was performed. Here, we benchmark an economical assay, comparing all commercially available scWGA kits. Our comparison is based on targeted sequencing of thousands of genomic loci, including highly mutable regions, from a large cohort of human single cells. Using this approach we have demonstrated the superiority of Ampli1 in genome coverage and of RepliG in reduced error rate. In summary, we show that no single kit is optimal across all categories, highlighting the need for a dedicated kit selection in accordance with experimental requirements.

In vivo locus-specific editing of the neuroepigenome.

Studies over the past several decades have identified numerous epigenetic mechanisms associated with pathological states in psychiatric and neurological disease. Until recently, studies investigating chromatin-regulatory proteins, using overexpression or knockdown approaches, did not establish causal roles for epigenetic modifications at specific genes because these techniques typically affect hundreds or thousands of genomic loci. In this Review, we describe recent efforts in using locus-specific neuroepigenome editing in vivo to, for the first time, define causal relationships between a single chromatin modification at a specific gene in a defined cell population and downstream measures at the molecular, cellular, circuit and behavioural levels. We briefly introduce three epigenome-editing platforms: zinc-finger proteins, transcriptional activator-like effectors and clustered regularly interspaced short palindromic repeats (CRISPR). We then explore the development of in vivo neuroepigenome-editing tools and their applications to resolve epigenetic contributions to the pathophysiology of brain diseases. We also discuss technical considerations for in vivo neuroepigenome-editing experiments and ongoing innovations in the field, including new tools to investigate chromatin marks, manipulate chromatin topology and induce epigenetic modifications at multiple genes in the same cell. Lastly, we explore the potential clinical applications of in vivo neuroepigenome editing for treating brain pathology.

Abundant expression of maternal siRNAs is a conserved feature of seed development

Small RNAs are abundant in plant reproductive tissues, especially 24-nucleotide (nt) small interfering RNAs (siRNAs). Most 24-nt siRNAs are dependent on RNA Pol IV and RNA-DEPENDENT RNA POLYMERASE 2 (RDR2) and establish DNA methylation at thousands of genomic loci in a process called RNA-directed DNA methylation (RdDM). In Brassica rapa, RdDM is required in the maternal sporophyte for successful seed development. Here, we demonstrate that a small number of siRNA loci account for over 90% of siRNA expression during B. rapa seed development. These loci exhibit unique characteristics with regard to their copy number and association with genomic features, but they resemble canonical 24-nt siRNA loci in their dependence on RNA Pol IV/RDR2 and role in RdDM. These loci are expressed in ovules before fertilization and in the seed coat, embryo, and endosperm following fertilization. We observed a similar pattern of 24-nt siRNA expression in diverse angiosperms despite rapid sequence evolution at siren loci. In the endosperm, siren siRNAs show a marked maternal bias, and siren expression in maternal sporophytic tissues is required for siren siRNA accumulation. Together, these results demonstrate that seed development occurs under the influence of abundant maternal siRNAs that might be transported to, and function in, filial tissues.

Glucocorticoid Receptor Orchestrates Interactions Between Glucocorticoid ‐Responsive and ‐Unresponsive Enhancers to Modulate Gene Expression

RATIONALEGlucocorticoids (GCs) are stress hormones that act on the GC receptor (GR; NR3C1) to elicit effects, including repression of inflammatory gene expression. The mechanisms by which GR modulate gene expression are not completely understood. GR binds to thousands of genomic loci, but only regulates the expression of hundreds of genes. Genome‐wide chromatin interaction studies revealed that the proximal promotor of various GC‐regulated genes interacts with multiple distal enhancer elements that harbour GC‐dependent GR binding. The significance of such interactions and the contribution of each locus to the achieved gene expression are yet to be investigated. The current work provides a focused perspective on genomic events that lead to the induction of KLF9 by GCs, a gene expression feature that is highly conserved in multiple cell types.RESULTSGC‐mediated induction of KLF9 was confirmed in the airway epithelial cell lines (A549 and BEAS‐2B), as well as in primary culture of various airway structural cells. This induction was abolished in A549 cells either by pretreating cells with the GR antagonist, Org34517, or by knocking down GR with siRNA. Chromatin immunoprecipitation of GR following GC treatment in A549, BEAS‐2B, and primary human bronchial epithelial (HBE) cells showed GR binding to sites 5.9, 6.7, 25, and 65 kb upstream of the KLF9 gene. Global nascent transcript analysis of BEAS‐2B cells following GC treatment showed marked induction of enhancer RNA (eRNA) from each of these sites. This observation was validated using qPCR in A549 and primary HBE cells. GR binding regions were then cloned upstream to a luciferase reporter to test their enhancer activity due to GCs in isolation from their genomic context. There was a clear distinction between proximal regions (i.e. 5.9 and 6.7 kb from KLF9 promotor) and the distal regions (25 and 65 kb) in terms of basal activity and response to GCs. Proximal regions showed very strong basal reporter activity that was not affected by GC treatment. Distal regions, however, showed minimal basal activity but produced dose‐dependent increase in reporter activity by GCs, which was abolished by mutating the GC response elements (GRE) from these regions. Kinetics of long‐range genomic interactions identified a significant interaction between the distal and the proximal regions in a slightly delayed fashion (4h onward), this was temporally associated with enhanced binding of EP300 (histone acetylase associated with gene activation) to the proximal region.CONCLUSIONSOur data suggests that each GR binding site upstream to KLF9 responds differently to GCs, but all sites may collectively contribute to the basal transcriptional activity, induction, and maintenance of KLF9 expression. Co‐localization of other transcriptional regulators may explain the distinction between different sites. Our data are essential in establishing a model to explain GR genomic binding behaviours and how they correlate with gene expression outcomes, thus providing insights on how to modulate GC actions to achieve clinical benefits.Support or Funding InformationFunding: Canadian Institutes of Health Research (CIHR) and Natural Sciences and Engineering Research Council of Canada (NSERC)

High-Throughput Screens of PAM-Flexible Cas9 Variants for Gene Knockout and Transcriptional Modulation

SUMMARYA key limitation of the widely used CRISPR enzyme S. pyogenes Cas9 is the strict requirement of an NGG protospacer-adjacent motif (PAM) at the target site. This constraint can be limiting for genome editing applications that require precise Cas9 positioning. Recently, two Cas9 variants with a relaxed PAM requirement (NG) have been developed (xCas9 and Cas9-NG), but their activity has been measured at only a small number of endogenous sites. Here, we devise a high-throughput Cas9 pooled competition screen to compare the performance of Cas9 variants at thousands of genomic loci for gene knockout, transcriptional activation, and inhibition. We show that PAM flexibility comes at a substantial cost of decreased DNA targeting and cleavage. Of the PAM-flexible variants, we find that Cas9-NG outperforms xCas9 regardless of genome engineering modality or PAM. Finally, we combine xCas9 mutations with those of Cas9-NG, creating a stronger transcriptional modulator than existing PAM-flexible Cas9 variants.

FOXA1 mutations alter pioneering activity, differentiation and prostate cancer phenotypes.

Mutations in the FOXA1 transcription factor define a unique subset of prostate cancers but the functional consequences of these mutations and whether they confer gain or loss of function is unknown1-9. By annotating the FOXA1 mutation landscape from 3086 human prostate cancers, we define two hotspots in the forkhead domain: Wing2 (~50% of all mutations) and R219 (~5%), a highly conserved DNA contact residue. Clinically, Wing2 mutations are seen in adenocarcinomas at all stages, whereas R219 mutations are enriched in metastatic tumors with neuroendocrine histology. Interrogation of the biologic properties of FOXA1WT and 14 FOXA1 mutants revealed gain-of-function in mouse prostate organoid proliferation assays. 12 of these mutants, as well as FOXA1WT, promoted an exaggerated pro-luminal differentiation program whereas two different R219 mutants blocked luminal differentiation and activate a mesenchymal and neuroendocrine transcriptional program. ATAC-seq of FOXA1WT and representative Wing2 and R219 mutants revealed dramatic, mutant-specific changes in open chromatin at thousands of genomic loci, together with novel sites of FOXA1 binding and associated increases in gene expression. Of note, peaks in R219 mutant expressing cells lack the canonical core FOXA1 binding motifs (GTAAAC/T) but are enriched for a related, non-canonical motif (GTAAAG/A), which is preferentially activated by R219 mutant FOXA1 in reporter assays. Thus, FOXA1 mutations alter its normal pioneering function through perturbation of normal luminal epithelial differentiation programs, providing further support to the role of lineage plasticity in cancer progression.