Human TSC2 Mutant Cells Exhibit Aberrations in Early Neurodevelopment Accompanied by Changes in the DNA Methylome.

Tuberous Sclerosis Complex (TSC) is a debilitating developmental disorder characterized by a variety of clinical manifestations. While benign tumors in the heart, lungs, kidney, and brain are all hallmarks of the disease, the most severe symptoms of TSC are often neurological, including seizures, autism, psychiatric disorders, and intellectual disabilities. TSC is caused by loss of function mutations in the TSC1 or TSC2 genes and consequent dysregulation of signaling via mechanistic Target of Rapamycin Complex 1 (mTORC1). While TSC neurological phenotypes are well-documented, it is not yet known how early in neural development TSC1/2-mutant cells diverge from the typical developmental trajectory. Another outstanding question is the contribution of homozygous-mutant cells to disease phenotypes and whether such phenotypes are also seen in the heterozygous-mutant populations that comprise the vast majority of cells in patients. Using TSC patient-derived isogenic induced pluripotent stem cells (iPSCs) with defined genetic changes, we observed aberrant early neurodevelopment in vitro, including misexpression of key proteins associated with lineage commitment and premature electrical activity. These alterations in differentiation were coincident with hundreds of differentially methylated DNA regions, including loci associated with key genes in neurodevelopment. Collectively, these data suggest that mutation or loss of TSC2 affects gene regulation and expression at earlier timepoints than previously appreciated, with implications for whether and how prenatal treatment should be pursued.

Genomic dissection and mutation-specific target discovery for breast cancer PIK3CA hotspot mutations.

BACKGROUND: Recent advancements in high-throughput genomics and targeted therapies have provided tremendous potential to identify and therapeutically target distinct mutations associated with cancers. However, to date the majority of targeted therapies are used to treat all functional mutations within the same gene, regardless of affected codon or phenotype. RESULTS: In this study, we developed a functional genomic analysis workflow with a unique isogenic cell line panel bearing two distinct hotspot PIK3CA mutations, E545K and H1047R, to accurately identify targetable differences between mutations within the same gene. We performed RNA-seq and ATAC-seq and identified distinct transcriptomic and epigenomic differences associated with each PIK3CA hotspot mutation. We used this data to curate a select CRISPR knock out screen to identify mutation-specific gene pathway vulnerabilities. These data revealed AREG as a E545K-preferential target that was further validated through in vitro analysis and publicly available patient databases. CONCLUSIONS: Using our multi-modal genomics framework, we discover distinct differences in genomic regulation between PIK3CA hotspot mutations, suggesting the PIK3CA mutations have different regulatory effects on the function and downstream signaling of the PI3K complex. Our results demonstrate the potential to rapidly uncover mutation specific molecular targets, specifically AREG and a proximal gene regulatory region, that may provide clinically relevant therapeutic targets. The methods outlined provide investigators with an integrative strategy to identify mutation-specific targets for the treatment of other oncogenic mutations in an isogenic system.

Human gene regulatory evolution is driven by the divergence of regulatory element function in both cis and trans.

Gene regulatory divergence between species can result from cis-acting local changes to regulatory element DNA sequences or global trans-acting changes to the regulatory environment. Understanding how these mechanisms drive regulatory evolution has been limited by challenges in identifying trans-acting changes. We present a comprehensive approach to directly identify cis- and trans-divergent regulatory elements between human and rhesus macaque lymphoblastoid cells using assay for transposase-accessible chromatin coupled to self-transcribing active regulatory region (ATAC-STARR) sequencing. In addition to thousands of cis changes, we discover an unexpected number (∼10,000) of trans changes and show that cis and trans elements exhibit distinct patterns of sequence divergence and function. We further identify differentially expressed transcription factors that underlie ∼37% of trans differences and trace how cis changes can produce cascades of trans changes. Overall, we find that most divergent elements (67%) experienced changes in both cis and trans, revealing a substantial role for trans divergence-alone and together with cis changes-in regulatory differences between species.

Endocrine islet ß-cell subtypes with differential function are derived from biochemically distinct embryonic endocrine islet progenitors that are regulated by maternal nutrients.

Endocrine islet b cells comprise heterogenous cell subsets. Yet when/how these subsets are produced and how stable they are remain unknown. Addressing these questions is important for preventing/curing diabetes, because lower numbers of b cells with better secretory function is a high risk of this disease. Using combinatorial cell lineage tracing, scRNA-seq, and DNA methylation analysis, we show here that embryonic islet progenitors with distinct gene expression and DNA methylation produce b-cell subtypes of different function and viability in adult mice. The subtype with better function is enriched for genes involved in vesicular production/trafficking, stress response, and Ca2+-secretion coupling, which further correspond to differential DNA methylation in putative enhancers of these genes. Maternal overnutrition, a major diabetes risk factor, reduces the proportion of endocrine progenitors of the b-cell subtype with better-function via deregulating DNA methyl transferase 3a. Intriguingly, the gene signature that defines mouse b-cell subtypes can reliably divide human cells into two sub-populations while the proportion of b cells with better-function is reduced in diabetic donors. The implication of these results is that modulating DNA methylation in islet progenitors using maternal food supplements can be explored to improve b-cell function in the prevention and therapy of diabetes.

Enhancer-promoter activation by the Kaposi sarcoma-associated herpesvirus episome maintenance protein LANA.

Higher-order genome structure influences the transcriptional regulation of cellular genes through the juxtaposition of regulatory elements, such as enhancers, close to promoters of target genes. While enhancer activation has emerged as an important facet of Kaposi sarcoma-associated herpesvirus (KSHV) biology, the mechanisms controlling enhancer-target gene expression remain obscure. Here, we discover that the KSHV genome tethering protein latency-associated nuclear antigen (LANA) potentiates enhancer-target gene expression in primary effusion lymphoma (PEL), a highly aggressive B cell lymphoma causally associated with KSHV. Genome-wide analyses demonstrate increased levels of enhancer RNA transcription as well as activating chromatin marks at LANA-bound enhancers. 3D genome conformation analyses identified genes critical for latency and tumorigenesis as targets of LANA-occupied enhancers, and LANA depletion results in their downregulation. These findings reveal a mechanism in enhancer-gene coordination and describe a role through which the main KSHV tethering protein regulates essential gene expression in PEL.

Genomic dissection and mutation-specific target discovery for breast cancer PIK3CA hotspot mutations.

Background: Recent advancements in high-throughput genomics and targeted therapies have provided tremendous potential to identify and therapeutically target distinct mutations associated with cancers. However, to date the majority of targeted therapies are used to treat all functional mutations within the same gene, regardless of affected codon or phenotype. Results: In this study, we developed a functional genomic analysis workflow with a unique isogenic cell line panel bearing two distinct hotspot PIK3CA mutations, E545K and H1047R, to accurately identify targetable differences between mutations within the same gene. We performed RNA-seq and ATAC-seq and identified distinct transcriptomic and epigenomic differences associated with each PIK3CA hotspot mutation. We used this data to curate a select CRISPR knock out screen to identify mutation-specific gene pathway vulnerabilities. These data revealed AREG as a E545K-preferential target that was further validated through in vitro analysis and publicly available patient databases. Conclusions: Using our multi-modal genomics framework, we discover distinct differences in genomic regulation between PIK3CA hotspot mutations, suggesting the PIK3CA mutations have different regulatory effects on the function and downstream signaling of the PI3K complex. Our results demonstrate the potential to rapidly uncover mutation specific molecular targets, specifically AREG and a proximal gene regulatory region, that may provide clinically relevant therapeutic targets. The methods outlined provide investigators with an integrative strategy to identify mutation-specific targets for the treatment of other oncogenic mutations in an isogenic system.

Polycomb repressor complex 2 suppresses interferon-responsive MHC-II expression in melanoma cells and is associated with anti-PD-1 resistance.

BACKGROUND: Despite the remarkable success of immunotherapy in treating melanoma, understanding of the underlying mechanisms of resistance remains limited. Emerging evidence suggests that upregulation of tumor-specific major histocompatibility complex-II (tsMHC-II) serves as a predictive marker for the response to anti-programmed death-1 (PD-1)/programmed death ligand 1 (PD-L1) therapy in various cancer types. The genetic and epigenetic pathways modulating tsMHC-II expression remain incompletely characterized. Here, we provide evidence that polycomb repressive complex 2 (PRC2)/EZH2 signaling and resulting H3K27 hypermethylation suppresses tsMHC-II. METHODS: RNA sequencing data from tumor biopsies from patients with cutaneous melanoma treated with or without anti-PD-1, targeted inhibition assays, and assays for transposase-accessible chromatin with sequencing were used to observe the relationship between EZH2 inhibition and interferon (IFN)-γ inducibility within the MHC-II pathway. RESULTS: We find that increased EZH2 pathway messenger RNA (mRNA) expression correlates with reduced mRNA expression of both presentation and T-cell genes. Notably, targeted inhibition assays revealed that inhibition of EZH2 influences the expression dynamics and inducibility of the MHC-II pathway following IFN-γ stimulation. Additionally, our analysis of patients with metastatic melanoma revealed a significant inverse association between PRC2-related gene expression and response to anti-PD-1 therapy. CONCLUSIONS: Collectively, our findings demonstrate that EZH2 inhibition leads to enhanced MHC-II expression potentially resulting from improved chromatin accessibility at CIITA, the master regulator of MHC-II. These insights shed light on the molecular mechanisms involved in tsMHC-II suppression and highlight the potential of targeting EZH2 as a therapeutic strategy to improve immunotherapy efficacy.