Genomic context sensitizes regulatory elements to genetic disruption.

Genomic context critically modulates regulatory function but is difficult to manipulate systematically. The murine insulin-like growth factor 2 (Igf2)/H19 locus is a paradigmatic model of enhancer selectivity, whereby CTCF occupancy at an imprinting control region directs downstream enhancers to activate either H19 or Igf2. We used synthetic regulatory genomics to repeatedly replace the native locus with 157-kb payloads, and we systematically dissected its architecture. Enhancer deletion and ectopic delivery revealed previously uncharacterized long-range regulatory dependencies at the native locus. Exchanging the H19 enhancer cluster with the Sox2 locus control region (LCR) showed that the H19 enhancers relied on their native surroundings while the Sox2 LCR functioned autonomously. Analysis of regulatory DNA actuation across cell types revealed that these enhancer clusters typify broader classes of context sensitivity genome wide. These results show that unexpected dependencies influence even well-studied loci, and our approach permits large-scale manipulation of complete loci to investigate the relationship between regulatory architecture and function.

Synthetic regulatory genomics uncovers enhancer context dependence at the Sox2 locus.

Sox2 expression in mouse embryonic stem cells (mESCs) depends on a distal cluster of DNase I hypersensitive sites (DHSs), but their individual contributions and degree of interdependence remain a mystery. We analyzed the endogenous Sox2 locus using Big-IN to scarlessly integrate large DNA payloads incorporating deletions, rearrangements, and inversions affecting single or multiple DHSs, as well as surgical alterations to transcription factor (TF) recognition sequences. Multiple mESC clones were derived for each payload, sequence-verified, and analyzed for Sox2 expression. We found that two DHSs comprising a handful of key TF recognition sequences were each sufficient for long-range activation of Sox2 expression. By contrast, three nearby DHSs were entirely context dependent, showing no activity alone but dramatically augmenting the activity of the autonomous DHSs. Our results highlight the role of context in modulating genomic regulatory element function, and our synthetic regulatory genomics approach provides a roadmap for the dissection of other genomic loci.

Mouse genome rewriting and tailoring of three important disease loci.

Genetically engineered mouse models (GEMMs) help us to understand human pathologies and develop new therapies, yet faithfully recapitulating human diseases in mice is challenging. Advances in genomics have highlighted the importance of non-coding regulatory genome sequences, which control spatiotemporal gene expression patterns and splicing in many human diseases1,2. Including regulatory extensive genomic regions, which requires large-scale genome engineering, should enhance the quality of disease modelling. Existing methods set limits on the size and efficiency of DNA delivery, hampering the routine creation of highly informative models that we call genomically rewritten and tailored GEMMs (GREAT-GEMMs). Here we describe 'mammalian switching antibiotic resistance markers progressively for integration' (mSwAP-In), a method for efficient genome rewriting in mouse embryonic stem cells. We demonstrate the use of mSwAP-In for iterative genome rewriting of up to 115 kb of a tailored Trp53 locus, as well as for humanization of mice using 116 kb and 180 kb human ACE2 loci. The ACE2 model recapitulated human ACE2 expression patterns and splicing, and notably, presented milder symptoms when challenged with SARS-CoV-2 compared with the existing K18-hACE2 model, thus representing a more human-like model of infection. Finally, we demonstrated serial genome writing by humanizing mouse Tmprss2 biallelically in the ACE2 GREAT-GEMM, highlighting the versatility of mSwAP-In in genome writing.

A versatile platform for locus-scale genome rewriting and verification.

Routine rewriting of loci associated with human traits and diseases would facilitate their functional analysis. However, existing DNA integration approaches are limited in terms of scalability and portability across genomic loci and cellular contexts. We describe Big-IN, a versatile platform for targeted integration of large DNAs into mammalian cells. CRISPR/Cas9-mediated targeting of a landing pad enables subsequent recombinase-mediated delivery of variant payloads and efficient positive/negative selection for correct clones in mammalian stem cells. We demonstrate integration of constructs up to 143 kb, and an approach for one-step scarless delivery. We developed a staged pipeline combining PCR genotyping and targeted capture sequencing for economical and comprehensive verification of engineered stem cells. Our approach should enable combinatorial interrogation of genomic functional elements and systematic locus-scale analysis of genome function.

Sequencing identifies multiple early introductions of SARS-CoV-2 to the New York City region.

Effective public response to a pandemic relies upon accurate measurement of the extent and dynamics of an outbreak. Viral genome sequencing has emerged as a powerful approach to link seemingly unrelated cases, and large-scale sequencing surveillance can inform on critical epidemiological parameters. Here, we report the analysis of 864 SARS-CoV-2 sequences from cases in the New York City metropolitan area during the COVID-19 outbreak in spring 2020. The majority of cases had no recent travel history or known exposure, and genetically linked cases were spread throughout the region. Comparison to global viral sequences showed that early transmission was most linked to cases from Europe. Our data are consistent with numerous seeds from multiple sources and a prolonged period of unrecognized community spreading. This work highlights the complementary role of genomic surveillance in addition to traditional epidemiological indicators.

Chromatin activation as a unifying principle underlying pathogenic mechanisms in multiple myeloma.

Multiple myeloma (MM) is a plasma cell neoplasm associated with a broad variety of genetic lesions. In spite of this genetic heterogeneity, MMs share a characteristic malignant phenotype whose underlying molecular basis remains poorly characterized. In the present study, we examined plasma cells from MM using a multi-epigenomics approach and demonstrated that, when compared to normal B cells, malignant plasma cells showed an extensive activation of regulatory elements, in part affecting coregulated adjacent genes. Among target genes up-regulated by this process, we found members of the NOTCH, NF-kB, MTOR signaling, and TP53 signaling pathways. Other activated genes included sets involved in osteoblast differentiation and response to oxidative stress, all of which have been shown to be associated with the MM phenotype and clinical behavior. We functionally characterized MM-specific active distant enhancers controlling the expression of thioredoxin (TXN), a major regulator of cellular redox status and, in addition, identified PRDM5 as a novel essential gene for MM. Collectively, our data indicate that aberrant chromatin activation is a unifying feature underlying the malignant plasma cell phenotype.

Epigenomic profiling of myelofibrosis reveals widespread DNA methylation changes in enhancer elements and ZFP36L1 as a potential tumor suppressor gene that is epigenetically regulated.

In this study we interrogated the DNA methylome of myelofibrosis patients using high-density DNA methylation arrays. We detected 35,215 differentially methylated CpG, corresponding to 10,253 genes, between myelofibrosis patients and healthy controls. These changes were present both in primary and secondary myelofibrosis, which showed no differences between them. Remarkably, most differentially methylated CpG were located outside gene promoter regions and showed significant association with enhancer regions. This aberrant enhancer hypermethylation was negatively correlated with the expression of 27 genes in the myelofibrosis cohort. Of these, we focused on the ZFP36L1 gene and validated its decreased expression and enhancer DNA hypermethylation in an independent cohort of patients and myeloid cell-lines. In vitro reporter assay and 5'-azacitidine treatment confirmed the functional relevance of hyper-methylation of ZFP36L1 enhancer. Furthermore, in vitro rescue of ZFP36L1 expression had an impact on cell proliferation and induced apoptosis in SET-2 cell line indicating a possible role of ZFP36L1 as a tumor suppressor gene in myelofibrosis. Collectively, we describe the DNA methylation profile of myelofibrosis, identifying extensive changes in enhancer elements and revealing ZFP36L1 as a novel candidate tumor suppressor gene.

Molecular characterization of autophagic and apoptotic signaling induced by sorafenib in liver cancer cells.

Sorafenib is the unique accepted molecular targeted drug for the treatment of patients in advanced stage of hepatocellular carcinoma. The current study evaluated cell signaling regulation of endoplasmic reticulum (ER) stress, c-Jun-N-terminal kinase (JNK), Akt, and 5'AMP-activated protein kinase (AMPK) leading to autophagy and apoptosis induced by sorafenib. Sorafenib induced early (3-12 hr) ER stress characterized by an increase of Ser51 P-eIF2α/eIF2α, C/EBP homologous protein (CHOP), IRE1α, and sXBP1, but a decrease of activating transcription factor 6 expression, overall temporally associated with the increase of Thr183,Tyr185 P-JNK1/2/JNK1/2, Thr172 P-AMPKα, Ser413 P-Foxo3a, Thr308 P-AKt/AKt and Thr32 P-Foxo3a/Foxo3a ratios, and reduction of Ser2481 P-mammalian target of rapamycin (mTOR)/mTOR and protein translation. This pattern was related to a transient increase of tBid, Bim EL , Beclin-1, Bcl-xL, Bcl-2, autophagy markers, and reduction of myeloid cell leukemia-1 (Mcl-1) expression. The progressive increase of CHOP expression, and reduction of Thr308 P-AKt/AKt and Ser473 P-AKt/AKt ratios were associated with the reduction of autophagic flux and an additional upregulation of Bim EL expression and caspase-3 activity (24 hr). Small interfering-RNA (si-RNA) assays showed that Bim, but not Bak and Bax, was involved in the induction of caspase-3 in sorafenib-treated HepG2 cells. Sorafenib increased autophagic and apoptotic markers in tumor-derived xenograft model. In conclusion, the early sorafenib-induced ER stress and regulation of JNK and AMPK-dependent signaling were related to the induction of survival autophagic process. The sustained drug treatment induced a progressive increase of ER stress and PERK-CHOP-dependent rise of Bim EL , which was associated with the shift from autophagy to apoptosis. The kinetic of Bim EL expression profile might also be related to the tight balance between AKt- and AMPK-related signaling leading to Foxo3a-dependent BIM EL upregulation.

Melatonin-induced increase in sensitivity of human hepatocellular carcinoma cells to sorafenib is associated with reactive oxygen species production and mitophagy.

Effects of sorafenib in hepatocellular carcinoma (HCC) are frequently transient due to tumor-acquired resistance, a phenotype that could be targeted by other molecules to reduce this adaptive response. Because melatonin is known to exert antitumor effects in HCC cells, this study investigated whether and how melatonin reduces resistance to sorafenib. Susceptibility to sorafenib (10 nmol/L to 50 µmol/L) in the presence of melatonin (1 and 2 mmol/L) was assessed in HCC cell lines HepG2, HuH7, and Hep3B. Cell viability was reduced by sorafenib from 1 µmol/L in HepG2 or HuH7 cells, and 2.5 µmol/L in Hep3B cells. Co-administration of melatonin and sorafenib exhibited a synergistic cytotoxic effect on HepG2 and HuH7 cells, while Hep3B cells displayed susceptibility to doses of sorafenib that had no effect when administrated alone. Co-administration of 2.5 µmol/L sorafenib and 1 mmol/L melatonin induced apoptosis in Hep3B cells, increasing PARP hydrolysis and BAX expression. We also observed an early colocalization of mitochondria with lysosomes, correlating with the expression of mitophagy markers PINK1 and Parkin and a reduction of mitofusin-2 and mtDNA compared with sorafenib administration alone. Moreover, increased reactive oxygen species production and mitochondrial membrane depolarization were elicited by drug combination, suggesting their contribution to mitophagy induction. Interestingly, Parkin silencing by siRNA to impair mitophagy significantly reduced cell killing, PARP cleavage, and BAX expression. These results demonstrate that the pro-oxidant capacity of melatonin and its impact on mitochondria stability and turnover via mitophagy increase sensitivity to the cytotoxic effect of sorafenib.

Melatonin and endoplasmic reticulum stress: relation to autophagy and apoptosis.

Endoplasmic reticulum (ER) is a dynamic organelle that participates in a number of cellular functions by controlling lipid metabolism, calcium stores, and proteostasis. Under stressful situations, the ER environment is compromised, and protein maturation is impaired; this causes misfolded proteins to accumulate and a characteristic stress response named unfolded protein response (UPR). UPR protects cells from stress and contributes to cellular homeostasis re-establishment; however, during prolonged ER stress, UPR activation promotes cell death. ER stressors can modulate autophagy which in turn, depending of the situation, induces cell survival or death. Interactions of different autophagy- and apoptosis-related proteins and also common signaling pathways have been found, suggesting an interplay between these cellular processes, although their dynamic features are still unknown. A number of pathologies including metabolic, neurodegenerative and cardiovascular diseases, cancer, inflammation, and viral infections are associated with ER stress, leading to a growing interest in targeting components of the UPR as a therapeutic strategy. Melatonin has a variety of antioxidant, anti-inflammatory, and antitumor effects. As such, it modulates apoptosis and autophagy in cancer cells, neurodegeneration and the development of liver diseases as well as other pathologies. Here, we review the effects of melatonin on the main ER stress mechanisms, focusing on its ability to regulate the autophagic and apoptotic processes. As the number of studies that have analyzed ER stress modulation by this indole remains limited, further research is necessary for a better understanding of the crosstalk between ER stress, autophagy, and apoptosis and to clearly delineate the mechanisms by which melatonin modulates these responses.

Ceramide metabolism regulates autophagy and apoptotic cell death induced by melatonin in liver cancer cells.

Autophagy is a process that maintains homeostasis during stress, although it also contributes to cell death under specific contexts. Ceramides have emerged as important effectors in the regulation of autophagy, mediating the crosstalk with apoptosis. Melatonin induces apoptosis of cancer cells; however, its role in autophagy and ceramide metabolism has yet to be clearly elucidated. This study was aimed to evaluate the effect of melatonin administration on autophagy and ceramide metabolism and its possible link with melatonin-induced apoptotic cell death in hepatocarcinoma (HCC) cells. Melatonin (2 mm) transiently induced autophagy in HepG2 cells through JNK phosphorylation, characterized by increased Beclin-1 expression, p62 degradation, and LC3II and LAMP-2 colocalization, which translated in decreased cell viability. Moreover, ATG5 silencing sensitized HepG2 cells to melatonin-induced apoptosis, suggesting a dual role of autophagy in cell death. Melatonin enhanced ceramide levels through both de novo synthesis and acid sphingomyelinase (ASMase) stimulation. Serine palmitoyltransferase (SPT) inhibition with myriocin prevented melatonin-induced autophagy and ASMase inhibition with imipramine-impaired autophagy flux. However, ASMase inhibition partially protected HepG2 cells against melatonin, while SPT inhibition significantly enhanced cell death. Findings suggest a crosstalk between SPT-mediated ceramide generation and autophagy in protecting against melatonin, while specific ASMase-induced ceramide production participates in melatonin-mediated cell death. Thus, dual blocking of SPT and autophagy emerges as a potential strategy to potentiate the apoptotic effects of melatonin in liver cancer cells.

Inhibition of matrix metalloproteinase-9 and nuclear factor kappa B contribute to melatonin prevention of motility and invasiveness in HepG2 liver cancer cells.

Hepatocellular carcinoma (HCC) is one of the most lethal human cancers worldwide because of its high incidence and its metastatic potential. Extracellular matrix degradation by matrix metalloproteinases (MMPs) has been connected with cancer cell invasion, and it has been suggested that inhibition of MMPs by synthetic and natural inhibitors may be of great importance in the HCC therapies. Melatonin, the main product of the pineal gland, exerts antiproliferative, proapoptotic, and antiangiogenic properties in HepG2 human hepatocellular cells, and exhibits anti-invasive and antimetastatic activities by suppressing the enzymatic activity of MMP-9 in different tumor types. However, the underlying mechanism of anti-invasive activity in HCC models has not been fully elucidated. Here, we demonstrate that 1 mm melatonin dosage reduced in IL-1ß-induced HepG2 cells MMP-9 gelatinase activity and inhibited cell invasion and motility through downregulation of MMP-9 gene expression and upregulation of the MMP-9-specific inhibitor tissue inhibitor of metalloproteinases (TIMP)-1. No significant changes were observed in the expression and activity of MMP-2, the other proteinase implicated in matrix collagen degradation, and its tissue inhibitor, TIMP-2. Also, melatonin significantly suppressed IL-1ß-induced nuclear factor-kappaB (NF-κB) translocation and transcriptional activity. In summary, we demonstrate that melatonin modulates motility and invasiveness of HepG2 cell in vitro through a molecular mechanism that involves TIMP-1 upregulation and attenuation of MMP-9 expression and activity via NF-κB signal pathway inhibition.

Genomic context sensitizes regulatory elements to genetic disruption.

Enhancer function is frequently investigated piecemeal using truncated reporter assays or single deletion analysis. Thus it remains unclear to what extent enhancer function at native loci relies on surrounding genomic context. Using the Big-IN technology for targeted integration of large DNAs, we analyzed the regulatory architecture of the murine Igf2/H19 locus, a paradigmatic model of enhancer selectivity. We assembled payloads containing a 157-kb functional Igf2/H19 locus and engineered mutations to genetically direct CTCF occupancy at the imprinting control region (ICR) that switches the target gene of the H19 enhancer cluster. Contrasting activity of payloads delivered at the endogenous Igf2/H19 locus or ectopically at Hprt revealed that the Igf2/H19 locus includes additional, previously unknown long-range regulatory elements. Exchanging components of the Igf2/H19 locus with the well-studied Sox2 locus showed that the H19 enhancer cluster functioned poorly out of context, and required its native surroundings to activate Sox2 expression. Conversely, the Sox2 locus control region (LCR) could activate both Igf2 and H19 outside its native context, but its activity was only partially modulated by CTCF occupancy at the ICR. Analysis of regulatory DNA actuation across different cell types revealed that, while the H19 enhancers are tightly coordinated within their native locus, the Sox2 LCR acts more independently. We show that these enhancer clusters typify broader classes of loci genome-wide. Our results show that unexpected dependencies may influence even the most studied functional elements, and our synthetic regulatory genomics approach permits large-scale manipulation of complete loci to investigate the relationship between locus architecture and function.

An automated network-based tool to search for metabolic vulnerabilities in cancer.

The development of computational tools for the systematic prediction of metabolic vulnerabilities of cancer cells constitutes a central question in systems biology. Here, we present gmctool, a freely accessible online tool that allows us to accomplish this task in a simple, efficient and intuitive environment. gmctool exploits the concept of genetic Minimal Cut Sets (gMCSs), a theoretical approach to synthetic lethality based on genome-scale metabolic networks, including a unique database of synthetic lethals computed from Human1, the most recent metabolic reconstruction of human cells. gmctool introduces qualitative and quantitative improvements over our previously developed algorithms to predict, visualize and analyze metabolic vulnerabilities in cancer, demonstrating a superior performance than competing algorithms. A detailed illustration of gmctool is presented for multiple myeloma (MM), an incurable hematological malignancy. We provide in vitro experimental evidence for the essentiality of CTPS1 (CTPS synthase) and UAP1 (UDP-N-Acetylglucosamine Pyrophosphorylase 1) in specific MM patient subgroups.

Preclinical models for prediction of immunotherapy outcomes and immune evasion mechanisms in genetically heterogeneous multiple myeloma.

The historical lack of preclinical models reflecting the genetic heterogeneity of multiple myeloma (MM) hampers the advance of therapeutic discoveries. To circumvent this limitation, we screened mice engineered to carry eight MM lesions (NF-κB, KRAS, MYC, TP53, BCL2, cyclin D1, MMSET/NSD2 and c-MAF) combinatorially activated in B lymphocytes following T cell-driven immunization. Fifteen genetically diverse models developed bone marrow (BM) tumors fulfilling MM pathogenesis. Integrative analyses of ∼500 mice and ∼1,000 patients revealed a common MAPK-MYC genetic pathway that accelerated time to progression from precursor states across genetically heterogeneous MM. MYC-dependent time to progression conditioned immune evasion mechanisms that remodeled the BM microenvironment differently. Rapid MYC-driven progressors exhibited a high number of activated/exhausted CD8+ T cells with reduced immunosuppressive regulatory T (Treg) cells, while late MYC acquisition in slow progressors was associated with lower CD8+ T cell infiltration and more abundant Treg cells. Single-cell transcriptomics and functional assays defined a high ratio of CD8+ T cells versus Treg cells as a predictor of response to immune checkpoint blockade (ICB). In clinical series, high CD8+ T/Treg cell ratios underlie early progression in untreated smoldering MM, and correlated with early relapse in newly diagnosed patients with MM under Len/Dex therapy. In ICB-refractory MM models, increasing CD8+ T cell cytotoxicity or depleting Treg cells reversed immunotherapy resistance and yielded prolonged MM control. Our experimental models enable the correlation of MM genetic and immunological traits with preclinical therapy responses, which may inform the next-generation immunotherapy trials.

The transcription factor DDIT3 is a potential driver of dyserythropoiesis in myelodysplastic syndromes.

Myelodysplastic syndromes (MDS) are hematopoietic stem cell (HSC) malignancies characterized by ineffective hematopoiesis, with increased incidence in older individuals. Here we analyze the transcriptome of human HSCs purified from young and older healthy adults, as well as MDS patients, identifying transcriptional alterations following different patterns of expression. While aging-associated lesions seem to predispose HSCs to myeloid transformation, disease-specific alterations may trigger MDS development. Among MDS-specific lesions, we detect the upregulation of the transcription factor DNA Damage Inducible Transcript 3 (DDIT3). Overexpression of DDIT3 in human healthy HSCs induces an MDS-like transcriptional state, and dyserythropoiesis, an effect associated with a failure in the activation of transcriptional programs required for normal erythroid differentiation. Moreover, DDIT3 knockdown in CD34+ cells from MDS patients with anemia is able to restore erythropoiesis. These results identify DDIT3 as a driver of dyserythropoiesis, and a potential therapeutic target to restore the inefficient erythroid differentiation characterizing MDS patients.

Tissue context determines the penetrance of regulatory DNA variation.

Functional assessment of disease-associated sequence variation at non-coding regulatory elements is complicated by their high degree of context sensitivity to both the local chromatin and nuclear environments. Allelic profiling of DNA accessibility across individuals has shown that only a select minority of sequence variation affects transcription factor (TF) occupancy, yet low sequence diversity in human populations means that no experimental assessment is available for the majority of disease-associated variants. Here we describe high-resolution in vivo maps of allelic DNA accessibility in liver, kidney, lung and B cells from 5 increasingly diverged strains of F1 hybrid mice. The high density of heterozygous sites in these hybrids enables precise quantification of effect size and cell-type specificity for hundreds of thousands of variants throughout the mouse genome. We show that chromatin-altering variants delineate characteristic sensitivity profiles for hundreds of TF motifs. We develop a compendium of TF-specific sensitivity profiles accounting for genomic context effects. Finally, we link maps of allelic accessibility to allelic transcript levels in the same samples. This work provides a foundation for quantitative prediction of cell-type specific effects of non-coding variation on TF activity, which will facilitate both fine-mapping and systems-level analyses of common disease-associated variation in human genomes.

Design and Synthesis of Novel Epigenetic Inhibitors Targeting Histone Deacetylases, DNA Methyltransferase 1, and Lysine Methyltransferase G9a with In Vivo Efficacy in Multiple Myeloma.

Concomitant inhibition of key epigenetic pathways involved in silencing tumor suppressor genes has been recognized as a promising strategy for cancer therapy. Herein, we report a first-in-class series of quinoline-based analogues that simultaneously inhibit histone deacetylases (from a low nanomolar range) and DNA methyltransferase-1 (from a mid-nanomolar range, IC50 < 200 nM). Additionally, lysine methyltransferase G9a inhibitory activity is achieved (from a low nanomolar range) by introduction of a key lysine mimic group at the 7-position of the quinoline ring. The corresponding epigenetic functional cellular responses are observed: histone-3 acetylation, DNA hypomethylation, and decreased histone-3 methylation at lysine-9. These chemical probes, multitarget epigenetic inhibitors, were validated against the multiple myeloma cell line MM1.S, demonstrating promising in vitro activity of 12a (CM-444) with GI50 of 32 nM, an adequate therapeutic window (>1 log unit), and a suitable pharmacokinetic profile. In vivo, 12a achieved significant antitumor efficacy in a xenograft mouse model of human multiple myeloma.