Epigenomics-Based Identification of Major Cell Identity Regulators within Heterogeneous Cell Populations

Overexpression of ZNF488 supports pancreatic cancer cell proliferation and tumorigenesis through inhibition of ferroptosis via regulating SCD1-mediated unsaturated fatty acid metabolism

BACKGROUND

Pancreatic cancer is a malignancy with high mortality. Once diagnosed, effective treatment strategies are limited and the five-year survival is extremely poor. Recent studies have shown that zinc finger proteins play important roles in tumorigenesis, including pancreatic cancer. However, it remains unknown on the clinical significance, function and underlying mechanisms of zinc finger protein 488 (ZNF488) during the development of pancreatic cancer.

METHODS

The clinical relevance of ZNF488 and stearoyl-CoA desaturase 1 (SCD1) was examined by analyzing the data from The Cancer Genome Atlas (TCGA) and immunohistochemical staining of the tissue microarray. Gain-of-function and loss-of-function experiments were performed by transfecting the cells with overexpressing lentivirus and siRNAs or shRNA lentivirus, respectively. The function of ZNF488 in pancreatic cancer was assessed by CCK8, colony formation, EdU staining, PI/Annexin V staining and xenografted tumorigenesis. Chip-qPCR assay was conducted to examine the interaction between ZNF488 and the promoter sequence of SCD1. Transcription activity was measured by dual luciferase reporter assay. mRNA and protein expression was detected by qRT-PCR and immunoblotting experiment, respectively. Fatty acid was quantified by gas chromatography mass spectrometry.

RESULTS

ZNF488 was overexpressed in pancreatic cancer samples compared with normal tissues. High expression of ZNF488 predicted the poor prognosis of the patients. In vitro, ZNF488 upregulation contributed to the EuU cooperation, proliferation and colony formation of MIAPaCa-2 and PANC-1 cells. Based on PI/Annexin V and trypan blue staining results, we showed that ZNF488 suppressed the ferroptosis and apoptosis of pancreatic cancer cells. Mechanistically, ZNF488 directly interacted with the promoter sequence of SCD1 gene and promoted its transcription activity, which resulted in enhanced palmitoleic and oleic acid production, as well as the peroxidation of fatty acid. In vivo, ZNF488 overexpression promoted the xenograted tumorigenesis of PANC-1, which was reversed by SCD1 knockdown. Importantly, combination of erastin and SCD1 inhibitors A939572 completely blunted the growth of ZNF488 overexpressed MIAPaCa-2 and PANC-1 cells. Usage of A939572 or erastin recovered the sensitivity of pancreatic cancer cells to the treatment of gemcitabine. Lastly, we found a positive correlation between ZNF488 and SCD1 in pancreatic cancer patients based on TCGA and immunohistochemical staining results.

CONCLUSION

Overexpression of ZNF488 suppresses the ferroptosis and apoptosis to support the growth and tumorigenesis of pancreatic cancer through augmentation of SCD1-mediated unsaturated fatty acid metabolism. Combination of SCD1 inhibitors, ferroptosis inducers or gemcitabine could be applied for the treatment of pancreatic cancer with overexpression of ZNF488.

Low RNA stability signifies increased post-transcriptional regulation of cell identity genes

Cell identity genes are distinct from other genes with respect to the epigenetic mechanisms to activate their transcription, e.g. by super-enhancers and broad H3K4me3 domains. However, it remains unclear whether their post-transcriptional regulation is also unique. We performed a systematic analysis of transcriptome-wide RNA stability in nine cell types and found that unstable transcripts were enriched in cell identity-related pathways while stable transcripts were enriched in housekeeping pathways. Joint analyses of RNA stability and chromatin state revealed significant enrichment of super-enhancers and broad H3K4me3 domains at the gene loci of unstable transcripts. Intriguingly, the RNA m6A methyltransferase, METTL3, preferentially binds to chromatin at super-enhancers, broad H3K4me3 domains and their associated genes. METTL3 binding intensity is positively correlated with RNA m6A methylation and negatively correlated with RNA stability of cell identity genes, probably due to co-transcriptional m6A modifications promoting RNA decay. Nanopore direct RNA-sequencing showed that METTL3 knockdown has a stronger effect on RNA m6A and mRNA stability for cell identity genes. Our data suggest a run-and-brake model, where cell identity genes undergo both frequent transcription and fast RNA decay to achieve precise regulation of RNA expression.

Epigenetic dosage identifies two major and functionally distinct β cell subtypes

The mechanisms that specify and stabilize cell subtypes remain poorly understood. Here, we identify two major subtypes of pancreatic β cells based on histone mark heterogeneity (βHI and βLO). βHI cells exhibit ∼4-fold higher levels of H3K27me3, distinct chromatin organization and compaction, and a specific transcriptional pattern. βHI and βLO cells also differ in size, morphology, cytosolic and nuclear ultrastructure, epigenomes, cell surface marker expression, and function, and can be FACS separated into CD24+ and CD24- fractions. Functionally, βHI cells have increased mitochondrial mass, activity, and insulin secretion in vivo and ex vivo. Partial loss of function indicates that H3K27me3 dosage regulates βHI/βLO ratio in vivo, suggesting that control of β cell subtype identity and ratio is at least partially uncoupled. Both subtypes are conserved in humans, with βHI cells enriched in humans with type 2 diabetes. Thus, epigenetic dosage is a novel regulator of cell subtype specification and identifies two functionally distinct β cell subtypes.

Identification of limb-specific Lmx1b auto-regulatory modules with Nail-patella syndrome pathogenicity

LMX1B haploinsufficiency causes Nail-patella syndrome (NPS; MIM 161200), characterized by nail dysplasia, absent/hypoplastic patellae, chronic kidney disease, and glaucoma. Accordingly in mice, Lmx1b has been shown to play crucial roles in the development of the limb, kidney and eye. Although one functional allele of Lmx1b appears adequate for development, Lmx1b null mice display ventral-ventral distal limbs with abnormal kidney, eye and cerebellar development, more disruptive, but fully concordant with NPS. In Lmx1b functional knockouts (KOs), Lmx1b transcription in the limb is decreased nearly 6-fold, indicating autoregulation. Herein, we report on two conserved Lmx1b-associated cis-regulatory modules (LARM1 and LARM2) that are bound by Lmx1b, amplify Lmx1b expression with unique spatial modularity in the limb, and are necessary for Lmx1b-mediated limb dorsalization. These enhancers, being conserved across vertebrates (including coelacanth, but not other fish species), and required for normal locomotion, provide a unique opportunity to study the role of dorsalization in the fin to limb transition. We also report on two NPS patient families with normal LMX1B coding sequence, but with loss-of-function variations in the LARM1/2 region, stressing the role of regulatory modules in disease pathogenesis.

Nail-patella syndrome (NPS) is characterized by nail dysplasia, absent/hypoplastic patellae, chronic kidney disease, and glaucoma and can be caused by haploinsufficiency of LMX1B; however, not all patients harbor pathogenic LMX1B mutations. Here the authors show that loss-of-function variations in upstream enhancer sequences are responsible for a limb specific form of human NPS.

The chromatin, topological and regulatory properties of pluripotency-associated poised enhancers are conserved in vivo

Poised enhancers (PEs) represent a genetically distinct set of distal regulatory elements that control the expression of major developmental genes. Before becoming activated in differentiating cells, PEs are already bookmarked in pluripotent cells with unique chromatin and topological features that could contribute to their privileged regulatory properties. However, since PEs were originally characterized in embryonic stem cells (ESC), it is currently unknown whether PEs are functionally conserved in vivo. Here, we show that the chromatin and 3D structural features of PEs are conserved among mouse pluripotent cells both in vitro and in vivo. We also uncovered that the interactions between PEs and their target genes are globally controlled by the combined action of Polycomb, Trithorax and architectural proteins. Moreover, distal regulatory sequences located close to developmental genes and displaying the typical genetic (i.e. CpG islands) and chromatin (i.e. high accessibility and H3K27me3 levels) features of PEs are commonly found across vertebrates. These putative PEs show high sequence conservation within specific vertebrate clades, with only a few being evolutionary conserved across all vertebrates. Lastly, by genetically disrupting PEs in mouse and chicken embryos, we demonstrate that these regulatory elements play essential roles during the induction of major developmental genes in vivo.

Overarching control of autophagy and DNA damage response by CHD6 revealed by modeling a rare human pathology

Members of the chromodomain-helicase-DNA binding (CHD) protein family are chromatin remodelers implicated in human pathologies, with CHD6 being one of its least studied members. We discovered a de novo CHD6 missense mutation in a patient clinically presenting the rare Hallermann-Streiff syndrome (HSS). We used genome editing to generate isogenic iPSC lines and model HSS in relevant cell types. By combining genomics with functional in vivo and in vitro assays, we show that CHD6 binds a cohort of autophagy and stress response genes across cell types. The HSS mutation affects CHD6 protein folding and impairs its ability to recruit co-remodelers in response to DNA damage or autophagy stimulation. This leads to accumulation of DNA damage burden and senescence-like phenotypes. We therefore uncovered a molecular mechanism explaining HSS onset via chromatin control of autophagic flux and genotoxic stress surveillance.

Conserved Epigenetic Regulatory Logic Infers Genes Governing Cell Identity

Determining genes that orchestrate cell differentiation in development and disease remains a fundamental goal of cell biology. This study establishes a genome-wide metric based on the gene-repressive trimethylation of histone H3 at lysine 27 (H3K27me3) across hundreds of diverse cell types to identify genetic regulators of cell differentiation. We introduce a computational method, TRIAGE, which uses discordance between gene-repressive tendency and expression to identify genetic drivers of cell identity. We apply TRIAGE to millions of genome-wide single-cell transcriptomes, diverse omics platforms, and eukaryotic cells and tissue types. Using a wide range of data, we validate the performance of TRIAGE in identifying cell-type-specific regulatory factors across diverse species including human, mouse, boar, bird, fish, and tunicate. Using CRISPR gene editing, we use TRIAGE to experimentally validate RNF220 as a regulator of Ciona cardiopharyngeal development and SIX3 as required for differentiation of endoderm in human pluripotent stem cells. A record of this paper’s transparent peer review process is included in the Supplemental Information.

Perturbing genes controlling cell decisions have major implications in development or disease. However, identifying key regulatory genes from the thousands expressed in a cell is challenging. TRIAGE is a computational method that distills patterns of epigenetic repression across diverse cell types to infer regulatory genes using input gene expression data from any cell type. Demonstrating its utility, we combine single-cell RNA-seq and TRIAGE to identify and experimentally confirm novel regulators of heart development in evolutionarily distant species.

GPR17 is an essential regulator for the temporal adaptation of sonic hedgehog signalling in neural tube development

Dorsal-ventral pattern formation of the neural tube is regulated by temporal and spatial activities of extracellular signalling molecules. Sonic hedgehog (Shh) assigns ventral neural subtypes via activation of the Gli transcription factors. Shh activity in the neural progenitor cells changes dynamically during differentiation, but the mechanisms regulating this dynamicity are not fully understood. Here, we show that temporal change of intracellular cAMP levels confers the temporal Shh signal, and the purinergic G-protein-coupled receptor GPR17 plays an essential role in this regulation. GPR17 is highly expressed in the ventral progenitor regions of the neural tube and acts as a negative regulator of the Shh signal in chick embryos. Although the activation of the GPR17-related signal inhibits ventral identity, perturbation of Gpr17 expression leads to aberrant expansion of ventral neural domains. Notably, perturbation of Gpr17 expression partially inhibits the negative feedback of Gli activity. Moreover, GPR17 increases cAMP activity, suggesting that it exerts its function by inhibiting the processing of Gli3 protein. GPR17 also negatively regulates Shh signalling in neural cells differentiated from mouse embryonic stem cells, suggesting that GPR17 function is conserved among different organisms. Our results demonstrate that GPR17 is a novel negative regulator of Shh signalling in a wide range of cellular contexts.

Development of glial restricted human neural stem cells for oligodendrocyte differentiation in vitro and in vivo

In this study, we have developed highly expandable neural stem cells (NSCs) from HESCs and iPSCs that artificially express the oligodendrocyte (OL) specific transcription factor gene Zfp488. This is enough to restrict them to an exclusive oligodendrocyte progenitor cell (OPC) fate during differentiation in vitro and in vivo. During CNS development, Zfp488 is induced during the early stages of OL generation, and then again during terminal differentiation of OLs. Interestingly, the human ortholog Znf488, crucial for OL development in human, has been recently identified to function as a dorsoventral pattering regulator in the ventral spinal cord for the generation of P1, P2/pMN, and P2 neural progenitor domains. Forced expression of Zfp488 gene in human NSCs led to the robust generation of OLs and suppression of neuronal and astrocyte fate in vitro and in vivo. Zfp488 expressing NSC derived oligodendrocytes are functional and can myelinate rat dorsal root ganglion neurons in vitro, and form myelin in Shiverer mice brain in vivo. After transplantation near a site of demyelination, Zfp488 expressing hNSCs migrated to the lesion and differentiated into premyelinating OLs. A certain fraction also homed in the subventricular zone (SVZ). Zfp488-ZsGreen1-hNSC derived OLs formed compact myelin in Shiverer mice brain seen under the electron microscope. Transplanted human neural stem cells (NSC) that have the potential to differentiate into functional oligodendrocytes in response to remyelinating signals can be a powerful therapeutic intervention for disorders where oligodendrocyte (OL) replacement is beneficial.

Corrupted coordination of epigenetic modifications leads to diverging chromatin states and transcriptional heterogeneity in CLL

Cancer evolution is fueled by epigenetic as well as genetic diversity. In chronic lymphocytic leukemia (CLL), intra-tumoral DNA methylation (DNAme) heterogeneity empowers evolution. Here, to comprehensively study the epigenetic dimension of cancer evolution, we integrate DNAme analysis with histone modification mapping and single cell analyses of RNA expression and DNAme in 22 primary CLL and 13 healthy donor B lymphocyte samples. Our data reveal corrupted coherence across different layers of the CLL epigenome. This manifests in decreased mutual information across epigenetic modifications and gene expression attributed to cell-to-cell heterogeneity. Disrupted epigenetic-transcriptional coordination in CLL is also reflected in the dysregulation of the transcriptional output as a function of the combinatorial chromatin states, including incomplete Polycomb-mediated gene silencing. Notably, we observe unexpected co-mapping of typically mutually exclusive activating and repressing histone modifications, suggestive of intra-tumoral epigenetic diversity. Thus, CLL epigenetic diversification leads to decreased coordination across layers of epigenetic information, likely reflecting an admixture of cells with diverging cellular identities.

In chronic lymphocytic leukemia (CLL), evolution is driven by transcriptional and epigenetic heterogeneity. Here, the authors integrate epigenomic analyses to show how intra-tumoral epigenetic diversity results in divergent chromatin states in CLL cells, increasing cell-to-cell transcriptional heterogeneity.

Pathological ASXL1 Mutations and Protein Variants Impair Neural Crest Development

The neural crest (NC) gives rise to a multitude of fetal tissues, and its misregulation is implicated in congenital malformations. Here, we investigated molecular mechanisms pertaining to NC-related symptoms in Bohring-Opitz syndrome (BOS), a developmental disorder linked to mutations in the Polycomb group factor Additional sex combs-like 1 (ASXL1). Genetically edited human pluripotent stem cell lines that were differentiated to NC progenitors and then xenotransplanted into chicken embryos demonstrated an impairment of NC delamination and emigration. Molecular analysis showed that ASXL1 mutations correlated with reduced activation of the transcription factor ZIC1 and the NC gene regulatory network. These findings were supported by differentiation experiments using BOS patient-derived induced pluripotent stem cell lines. Expression of truncated ASXL1 isoforms (amino acids 1-900) recapitulated the NC phenotypes in vitro and in ovo, raising the possibility that truncated ASXL1 variants contribute to BOS pathology. Collectively, we expand the understanding of truncated ASXL1 in BOS and in the human NC.

Modeling the Pathological Long-Range Regulatory Effects of Human Structural Variation with Patient-Specific hiPSCs

The pathological consequences of structural variants disrupting 3D genome organization can be difficult to elucidate in vivo due to differences in gene dosage sensitivity between mice and humans. This is illustrated by branchiooculofacial syndrome (BOFS), a rare congenital disorder caused by heterozygous mutations within TFAP2A, a neural crest regulator for which humans, but not mice, are haploinsufficient. Here, we present a BOFS patient carrying a heterozygous inversion with one breakpoint located within a topologically associating domain (TAD) containing enhancers essential for TFAP2A expression in human neural crest cells (hNCCs). Using patient-specific hiPSCs, we show that, although the inversion shuffles the TFAP2A hNCC enhancers with novel genes within the same TAD, this does not result in enhancer adoption. Instead, the inversion disconnects one TFAP2A allele from its cognate enhancers, leading to monoallelic and haploinsufficient TFAP2A expression in patient hNCCs. Our work illustrates the power of hiPSC differentiation to unveil long-range pathomechanisms.

Chromatin Immunoprecipitation (ChIP) Protocol for Low-abundance Embryonic Samples

Chromatin immunoprecipitation (ChIP) is a widely-used technique for mapping the localization of post-translationally modified histones, histone variants, transcription factors, or chromatin-modifying enzymes at a given locus or on a genome-wide scale. The combination of ChIP assays with next-generation sequencing (i.e., ChIP-Seq) is a powerful approach to globally uncover gene regulatory networks and to improve the functional annotation of genomes, especially of non-coding regulatory sequences. ChIP protocols normally require large amounts of cellular material, thus precluding the applicability of this method to investigating rare cell types or small tissue biopsies. In order to make the ChIP assay compatible with the amount of biological material that can typically be obtained in vivo during early vertebrate embryogenesis, we describe here a simplified ChIP protocol in which the number of steps required to complete the assay were reduced to minimize sample loss. This ChIP protocol has been successfully used to investigate different histone modifications in various embryonic chicken and adult mouse tissues using low to medium cell numbers (5 x 10⁴ - 5 x 10⁵ cells). Importantly, this protocol is compatible with ChIP-seq technology using standard library preparation methods, thus providing global epigenomic maps in highly relevant embryonic tissues.