Mediator kinase inhibition impedes transcriptional plasticity and prevents resistance to ERK/MAPK-targeted therapy in KRAS-mutant cancers.

Acquired resistance remains a major challenge for therapies targeting oncogene activated pathways. KRAS is the most frequently mutated oncogene in human cancers, yet strategies targeting its downstream signaling kinases have failed to produce durable treatment responses. Here, we developed multiple models of acquired resistance to dual-mechanism ERK/MAPK inhibitors across KRAS-mutant pancreatic, colorectal, and lung cancers, and then probed the long-term events enabling survival against this class of drugs. These studies revealed that resistance emerges secondary to large-scale transcriptional adaptations that are diverse and cell line-specific. Transcriptional reprogramming extends beyond the well-established early response, and instead represents a dynamic, evolved process that is refined to attain a stably resistant phenotype. Mechanistic and translational studies reveal that resistance to dual-mechanism ERK/MAPK inhibition is broadly susceptible to manipulation of the epigenetic machinery, and that Mediator kinase, in particular, can be co-targeted at a bottleneck point to prevent diverse, cell line-specific resistance programs.

Transcriptional and epigenetic regulators of human CD8+ T cell function identified through orthogonal CRISPR screens.

Clinical response to adoptive T cell therapies is associated with the transcriptional and epigenetic state of the cell product. Thus, discovery of regulators of T cell gene networks and their corresponding phenotypes has potential to improve T cell therapies. Here we developed pooled, epigenetic CRISPR screening approaches to systematically profile the effects of activating or repressing 120 transcriptional and epigenetic regulators on human CD8+ T cell state. We found that BATF3 overexpression promoted specific features of memory T cells and attenuated gene programs associated with cytotoxicity, regulatory T cell function, and exhaustion. Upon chronic antigen stimulation, BATF3 overexpression countered phenotypic and epigenetic signatures of T cell exhaustion. Moreover, BATF3 enhanced the potency of CAR T cells in both in vitro and in vivo tumor models and programmed a transcriptional profile that correlates with positive clinical response to adoptive T cell therapy. Finally, we performed CRISPR knockout screens that defined cofactors and downstream mediators of the BATF3 gene network.

Cooperative regulation of coupled oncoprotein synthesis and stability in triple-negative breast cancer by EGFR and CDK12/13.

Evidence has long suggested that epidermal growth factor receptor (EGFR) may play a prominent role in triple-negative breast cancer (TNBC) pathogenesis, but clinical trials of EGFR inhibitors have yielded disappointing results. Using a candidate drug screen, we identified that inhibition of cyclin-dependent kinases 12 and 13 (CDK12/13) dramatically sensitizes diverse models of TNBC to EGFR blockade. This combination therapy drives cell death through the 4E-BP1-dependent suppression of the translation and translation-linked turnover of driver oncoproteins, including MYC. A genome-wide CRISPR/Cas9 screen identified the CCR4-NOT complex as a major determinant of sensitivity to the combination therapy whose loss renders 4E-BP1 unresponsive to drug-induced dephosphorylation, thereby rescuing MYC translational suppression and promoting MYC stability. The central roles of CCR4-NOT and 4E-BP1 in response to the combination therapy were further underscored by the observation of CNOT1 loss and rescue of 4E-BP1 phosphorylation in TNBC cells that naturally evolved therapy resistance. Thus, pharmacological inhibition of CDK12/13 reveals a long-proposed EGFR dependence in TNBC that functions through the cooperative regulation of translation-coupled oncoprotein stability.

Orthogonal CRISPR screens to identify transcriptional and epigenetic regulators of human CD8 T cell function.

The clinical response to adoptive T cell therapies is strongly associated with transcriptional and epigenetic state. Thus, technologies to discover regulators of T cell gene networks and their corresponding phenotypes have great potential to improve the efficacy of T cell therapies. We developed pooled CRISPR screening approaches with compact epigenome editors to systematically profile the effects of activation and repression of 120 transcription factors and epigenetic modifiers on human CD8+ T cell state. These screens nominated known and novel regulators of T cell phenotypes with BATF3 emerging as a high confidence gene in both screens. We found that BATF3 overexpression promoted specific features of memory T cells such as increased IL7R expression and glycolytic capacity, while attenuating gene programs associated with cytotoxicity, regulatory T cell function, and T cell exhaustion. In the context of chronic antigen stimulation, BATF3 overexpression countered phenotypic and epigenetic signatures of T cell exhaustion. CAR T cells overexpressing BATF3 significantly outperformed control CAR T cells in both in vitro and in vivo tumor models. Moreover, we found that BATF3 programmed a transcriptional profile that correlated with positive clinical response to adoptive T cell therapy. Finally, we performed CRISPR knockout screens with and without BATF3 overexpression to define co-factors and downstream factors of BATF3, as well as other therapeutic targets. These screens pointed to a model where BATF3 interacts with JUNB and IRF4 to regulate gene expression and illuminated several other novel targets for further investigation.

ZBTB7A promotes virus-host homeostasis during human coronavirus 229E infection.

The cellular fate after infection with human coronaviruses (HCoVs) is typically death. Previous data suggest, however, that the transcriptional state of an individual cell may sometimes allow additional outcomes of infection. Here, to probe the range of interactions a permissive cell type can have with a HCoV, we perform a CRISPR activation screen with HCoV-229E. The screen identified the transcription factor ZBTB7A, which strongly promotes cell survival after infection. Rather than suppressing viral infection, ZBTB7A upregulation allows the virus to induce a persistent infection and homeostatic state with the cell. We also find that control of oxidative stress is a primary driver of cellular survival during HCoV-229E infection. These data illustrate that, in addition to the nature of the infecting virus and the type of cell that it encounters, the cellular gene expression profile prior to infection can affect the eventual fate.

Epigenome editing of the CFTR-locus for treatment of cystic fibrosis.

BACKGROUND: Mechanisms governing the diversity of CFTR gene expression throughout the body are complex. Multiple intronic and distal regulatory elements are responsible for regulating differential CFTR expression across tissues. METHODS: Drawing on published data, 18 high-priority genomic regions were identified and interrogated for CFTR-enhancer function using CRISPR/dCas9-based epigenome editing tools. Each region was evaluated by dCas9p300 and dCas9KRAB for its ability to enhance or repress CFTR expression, respectively. RESULTS: Multiple genomic regions were tested for enhancer activity using CRISPR/dCas9 epigenome editing. dCas9p300 mediates a significant increase in CFTR mRNA levels when targeted to the promoter and a region 44 kb upstream of the transcriptional start site in a CFTR-low expressing cell line. Multiple gRNAs targeting the promoter induced a robust increase in CFTR protein levels. In contrast, dCas9KRAB-mediated repression is much more robust with 10 of the 18 evaluated genomic regions inducing CFTR protein knockdown. To evaluate the therapeutic efficacy of modulating CFTR gene regulation, dCas9p300 was used to induce elevated levels of CFTR from the endogenous locus in ΔF508/ΔF508 human bronchial epithelial cells. Ussing chamber studies demonstrated a synergistic increase in ion transport in response to CRISPR-induced expression of ΔF508 CFTR mRNA along with VX809 treatment. CONCLUSIONS: CRISPR/dCas9-based epigenome-editing provides a previously unexplored tool for interrogating CFTR enhancer function. Here, we demonstrate that therapeutic interventions that increase the expression of CFTR may improve the efficacy of CFTR modulators. A better understanding CFTR regulatory mechanisms could uncover novel therapeutic interventions for the development of cystic fibrosis therapies.

c-Maf regulates the plasticity of group 3 innate lymphoid cells by restraining the type 1 program.

CCR6- group 3 innate lymphoid cells (ILC3s) are mediators of intestinal immunity and barrier function that possess the capacity to acquire type 1 effector features and fully convert into ILC1s. The molecular mechanisms governing such plasticity are undefined. Here, we identified c-Maf as an essential regulator of ILC3 homeostasis and plasticity that limits physiological ILC1 conversion. Phenotypic analysis of effector status in Maf-deficient CCR6- ILC3s, coupled with evaluation of global changes in transcriptome, chromatin accessibility, and transcription factor motif enrichment, revealed that c-Maf enforces ILC3 identity. c-Maf promoted ILC3 accessibility and supported RORγt activity and expression of type 3 effector genes. Conversely, c-Maf antagonized type 1 programming, largely through restraint of T-bet expression and function. Mapping of the dynamic changes in chromatin landscape accompanying CCR6- ILC3 development and ILC1 conversion solidified c-Maf as a gatekeeper of type 1 regulatory transformation and a controller of ILC3 fate.

Glucocorticoid receptor recruits to enhancers and drives activation by motif-directed binding.

Glucocorticoids are potent steroid hormones that regulate immunity and metabolism by activating the transcription factor (TF) activity of glucocorticoid receptor (GR). Previous models have proposed that DNA binding motifs and sites of chromatin accessibility predetermine GR binding and activity. However, there are vast excesses of both features relative to the number of GR binding sites. Thus, these features alone are unlikely to account for the specificity of GR binding and activity. To identify genomic and epigenetic contributions to GR binding specificity and the downstream changes resultant from GR binding, we performed hundreds of genome-wide measurements of TF binding, epigenetic state, and gene expression across a 12-h time course of glucocorticoid exposure. We found that glucocorticoid treatment induces GR to bind to nearly all pre-established enhancers within minutes. However, GR binds to only a small fraction of the set of accessible sites that lack enhancer marks. Once GR is bound to enhancers, a combination of enhancer motif composition and interactions between enhancers then determines the strength and persistence of GR binding, which consequently correlates with dramatic shifts in enhancer activation. Over the course of several hours, highly coordinated changes in TF binding and histone modification occupancy occur specifically within enhancers, and these changes correlate with changes in the expression of nearby genes. Following GR binding, changes in the binding of other TFs precede changes in chromatin accessibility, suggesting that other TFs are also sensitive to genomic features beyond that of accessibility.

Clustering gene expression time series data using an infinite Gaussian process mixture model.

Transcriptome-wide time series expression profiling is used to characterize the cellular response to environmental perturbations. The first step to analyzing transcriptional response data is often to cluster genes with similar responses. Here, we present a nonparametric model-based method, Dirichlet process Gaussian process mixture model (DPGP), which jointly models data clusters with a Dirichlet process and temporal dependencies with Gaussian processes. We demonstrate the accuracy of DPGP in comparison to state-of-the-art approaches using hundreds of simulated data sets. To further test our method, we apply DPGP to published microarray data from a microbial model organism exposed to stress and to novel RNA-seq data from a human cell line exposed to the glucocorticoid dexamethasone. We validate our clusters by examining local transcription factor binding and histone modifications. Our results demonstrate that jointly modeling cluster number and temporal dependencies can reveal shared regulatory mechanisms. DPGP software is freely available online at https://github.com/PrincetonUniversity/DP_GP_cluster.

Coordinated regulatory variation associated with gestational hyperglycaemia regulates expression of the novel hexokinase HKDC1.

Maternal glucose levels during pregnancy impact the developing fetus, affecting metabolic health both early and later on in life. Both genetic and environmental factors influence maternal metabolism, but little is known about the genetic mechanisms that alter glucose metabolism during pregnancy. Here, we report that haplotypes previously associated with gestational hyperglycaemia in the third trimester disrupt regulatory element activity and reduce expression of the nearby HKDC1 gene. We further find that experimentally reducing or increasing HKDC1 expression reduces or increases hexokinase activity, respectively, in multiple cellular models; in addition, purified HKDC1 protein has hexokinase activity in vitro. Together, these results suggest a novel mechanism of gestational glucose regulation in which the effects of genetic variants in multiple regulatory elements alter glucose homeostasis by coordinately reducing expression of the novel hexokinase HKDC1.

Distinct properties of cell-type-specific and shared transcription factor binding sites.

Most human transcription factors bind a small subset of potential genomic sites and often use different subsets in different cell types. To identify mechanisms that govern cell-type-specific transcription factor binding, we used an integrative approach to study estrogen receptor α (ER). We found that ER exhibits two distinct modes of binding. Shared sites, bound in multiple cell types, are characterized by high-affinity estrogen response elements (EREs), inaccessible chromatin, and a lack of DNA methylation, while cell-specific sites are characterized by a lack of EREs, co-occurrence with other transcription factors, and cell-type-specific chromatin accessibility and DNA methylation. These observations enabled accurate quantitative models of ER binding that suggest tethering of ER to one-third of cell-specific sites. The distinct properties of cell-specific binding were also observed with glucocorticoid receptor and for ER in primary mouse tissues, representing an elegant genomic encoding scheme for generating cell-type-specific gene regulation.

Evidence for multiple roles for grainyhead-like 2 in the establishment and maintenance of human mucociliary airway epithelium.[corrected].

Most of the airways of the human lung are lined by an epithelium made up of ciliated and secretory luminal cells and undifferentiated basal progenitor cells. The integrity of this epithelium and its ability to act as a selective barrier are critical for normal lung function. In other epithelia, there is evidence that transcription factors of the evolutionarily conserved grainyheadlike (GRHL) family play key roles in coordinating multiple cellular processes required for epithelial morphogenesis, differentiation, remodeling, and repair. However, only a few target genes have been identified, and little is known about GRHL function in the adult lung. Here we focus on the role of GRHL2 in primary human bronchial epithelial cells, both as undifferentiated progenitors and as they differentiate in air-liquid interface culture into an organized mucociliary epithelium with transepithelial resistance. Using a dominant-negative protein or shRNA to inhibit GRHL2, we follow changes in epithelial phenotype and gene transcription using RNA sequencing or microarray analysis. We identify several hundreds of genes that are directly or indirectly regulated by GRHL2 in both undifferentiated cells and air-liquid interface cultures. Using ChIP sequencing to map sites of GRHL2 binding in the basal cells, we identify 7,687 potential primary targets and confirm that GRHL2 binding is strongly enriched near GRHL2-regulated genes. Taken together, the results support the hypothesis that GRHL2 plays a key role in regulating many physiological functions of human airway epithelium, including those involving cell morphogenesis, adhesion, and motility.

The chromosome 3q25 genomic region is associated with measures of adiposity in newborns in a multi-ethnic genome-wide association study.

Newborns characterized as large and small for gestational age are at risk for increased mortality and morbidity during the first year of life as well as for obesity and dysglycemia as children and adults. The intrauterine environment and fetal genes contribute to the fetal size at birth. To define the genetic architecture underlying the newborn size, we performed a genome-wide association study (GWAS) in 4281 newborns in four ethnic groups from the Hyperglycemia and Adverse Pregnancy Outcome Study. We tested for association with newborn anthropometric traits (birth length, head circumference, birth weight, percent fat mass and sum of skinfolds) and newborn metabolic traits (cord glucose and C-peptide) under three models. Model 1 adjusted for field center, ancestry, neonatal gender, gestational age at delivery, parity, maternal age at oral glucose tolerance test (OGTT); Model 2 adjusted for Model 1 covariates, maternal body mass index (BMI) at OGTT, maternal height at OGTT, maternal mean arterial pressure at OGTT, maternal smoking and drinking; Model 3 adjusted for Model 2 covariates, maternal glucose and C-peptide at OGTT. Strong evidence for association was observed with measures of newborn adiposity (sum of skinfolds model 3 Z-score 7.356, P = 1.90×10⁻¹³, and to a lesser degree fat mass and birth weight) and a region on Chr3q25.31 mapping between CCNL and LEKR1. These findings were replicated in an independent cohort of 2296 newborns. This region has previously been shown to be associated with birth weight in Europeans. The current study suggests that association of this locus with birth weight is secondary to an effect on fat as opposed to lean body mass.

Dynamic DNA methylation across diverse human cell lines and tissues.

As studies of DNA methylation increase in scope, it has become evident that methylation has a complex relationship with gene expression, plays an important role in defining cell types, and is disrupted in many diseases. We describe large-scale single-base resolution DNA methylation profiling on a diverse collection of 82 human cell lines and tissues using reduced representation bisulfite sequencing (RRBS). Analysis integrating RNA-seq and ChIP-seq data illuminates the functional role of this dynamic mark. Loci that are hypermethylated across cancer types are enriched for sites bound by NANOG in embryonic stem cells, which supports and expands the model of a stem/progenitor cell signature in cancer. CpGs that are hypomethylated across cancer types are concentrated in megabase-scale domains that occur near the telomeres and centromeres of chromosomes, are depleted of genes, and are enriched for cancer-specific EZH2 binding and H3K27me3 (repressive chromatin). In noncancer samples, there are cell-type specific methylation signatures preserved in primary cell lines and tissues as well as methylation differences induced by cell culture. The relationship between methylation and expression is context-dependent, and we find that CpG-rich enhancers bound by EP300 in the bodies of expressed genes are unmethylated despite the dense gene-body methylation surrounding them. Non-CpG cytosine methylation occurs in human somatic tissue, is particularly prevalent in brain tissue, and is reproducible across many individuals. This study provides an atlas of DNA methylation across diverse and well-characterized samples and enables new discoveries about DNA methylation and its role in gene regulation and disease.

Genistein and bisphenol A exposure cause estrogen receptor 1 to bind thousands of sites in a cell type-specific manner.

Endogenous estrogens that are synthesized in the body impact gene regulation by activating estrogen receptors in diverse cell types. Exogenous compounds that have estrogenic properties can also be found circulating in the blood in both children and adults. The genome-wide impact of these environmental estrogens on gene regulation is unclear. To obtain an integrated view of gene regulation in response to environmental and endogenous estrogens on a genome-wide scale, we performed ChIP-seq to identify estrogen receptor 1 (ESR1; previously estrogen receptor α) binding sites, and RNA-seq in endometrial cancer cells exposed to bisphenol A (BPA; found in plastics), genistein (GEN; found in soybean), or 17ß-estradiol (E2; an endogenous estrogen). GEN and BPA treatment induces thousands of ESR1 binding sites and >50 gene expression changes, representing a subset of E2-induced gene regulation changes. Genes affected by E2 were highly enriched for ribosome-associated proteins; however, GEN and BPA failed to regulate most ribosome-associated proteins and instead enriched for transporters of carboxylic acids. Treatment-dependent changes in gene expression were associated with treatment-dependent ESR1 binding sites, with the exception that many genes up-regulated by E2 harbored a BPA-induced ESR1 binding site but failed to show any expression change after BPA treatment. GEN and BPA exhibited a similar relationship to E2 in the breast cancer line T-47D, where cell type specificity played a much larger role than treatment specificity. Overall, both environmental estrogens clearly regulate gene expression through ESR1 on a genome-wide scale, although with lower potency resulting in less ESR1 binding sites and less gene expression changes compared to the endogenous estrogen, E2.

The hypersensitive glucocorticoid response specifically regulates period 1 and expression of circadian genes.

Glucocorticoids regulate gene expression by binding and activating the glucocorticoid receptor (GR). While ligand affinity determines the global sensitivity of the response, additional proteins act on the genome to tune sensitivity of some genes. However, the genomic extent and specificity of dose-specific glucocorticoid responses are unknown. We show that dose-specific glucocorticoid responses are extraordinarily specific at the genomic scale, able to distinctly express a single gene, the circadian rhythm gene for Period 1 (PER1), at concentrations consistent with the nighttime nadir of human cortisol. We mapped the PER1 response to a single GR binding site. The specific GR binding sequence did not impact sensitivity, and we instead attributed the response to a combination of additional transcription factors and chromatin accessibility acting in the same locus. The PER1 hypersensitive response element is conserved in the mouse, where we found similar upregulation of Per1 in pituitary cells. Targeted and transient overexpression of PER1 led to regulation of additional circadian rhythm genes hours later, suggesting that hypersensitive expression of PER1 impacts circadian gene expression. These findings show that hypersensitive GR binding occurs throughout the genome, drives targeted gene expression, and may be important to endocrine mediation of peripheral circadian rhythms.

Analysis of DNA methylation in a three-generation family reveals widespread genetic influence on epigenetic regulation.

The methylation of cytosines in CpG dinucleotides is essential for cellular differentiation and the progression of many cancers, and it plays an important role in gametic imprinting. To assess variation and inheritance of genome-wide patterns of DNA methylation simultaneously in humans, we applied reduced representation bisulfite sequencing (RRBS) to somatic DNA from six members of a three-generation family. We observed that 8.1% of heterozygous SNPs are associated with differential methylation in cis, which provides a robust signature for Mendelian transmission and relatedness. The vast majority of differential methylation between homologous chromosomes (>92%) occurs on a particular haplotype as opposed to being associated with the gender of the parent of origin, indicating that genotype affects DNA methylation of far more loci than does gametic imprinting. We found that 75% of genotype-dependent differential methylation events in the family are also seen in unrelated individuals and that overall genotype can explain 80% of the variation in DNA methylation. These events are under-represented in CpG islands, enriched in intergenic regions, and located in regions of low evolutionary conservation. Even though they are generally not in functionally constrained regions, 22% (twice as many as expected by chance) of genes harboring genotype-dependent DNA methylation exhibited allele-specific gene expression as measured by RNA-seq of a lymphoblastoid cell line, indicating that some of these events are associated with gene expression differences. Overall, our results demonstrate that the influence of genotype on patterns of DNA methylation is widespread in the genome and greatly exceeds the influence of imprinting on genome-wide methylation patterns.

Agent-based modeling of the context dependency in T cell recognition.

Antigen recognition by T cells is a key event in the adaptive immune response. T cells scan the surface of antigen-presenting cells (APCs) or target cells for specific peptides bound to MHC molecules. In the physiological setting, a typical APC presents tens of thousands of diverse endogenous self-derived peptides complexed to MHC (pMHC complexes). When 'foreign' peptides are presented, they constitute a small fraction of the total surface peptide repertoire. As T cells seem to be capable of discerning minute amounts of 'foreign' peptides among a complex background of self-peptides, endogenous peptides are generally assumed to play no role in recognition. However, recent results suggest that these background peptides may alter the sensitivity of T cells to foreign peptides. Current experimental limitations preclude analysis of peptide mixtures approaching physiological complexity, making it difficult to further address the role of complex background peptides. In this paper, we present a computational model to test how complex, varied peptide populations on an APC could potentially modulate a T cell's ability to detect the presence of small numbers of agonist peptides among a diverse population. We use the model to investigate the notion that under physiological conditions, T cell recognition of foreign peptides is context dependent, that is, T cells process signals gathered from all pMHC interactions, not just from a few agonist peptides while ignoring all others.