Loss of MGA repression mediated by an atypical polycomb complex promotes tumor progression and invasiveness

Genomic Characterization of High-Grade Serous Ovarian Carcinoma Reveals Distinct Somatic Features in Black Individuals

Black individuals experience worse survival after a diagnosis of high-grade serous ovarian carcinoma (HGSC) than White individuals and are underrepresented in ovarian cancer research. To date, the understanding of the molecular and genomic heterogeneity of HGSC is based primarily on the evaluation of tumors from White individuals. In the present study, we performed whole-exome sequencing on HGSC samples from 211 Black patients to identify significantly mutated genes and characterize mutational signatures, assessing their distributions by gene expression subtypes. The occurrence and frequency of somatic mutations and signatures by self-reported race were compared with historic data from The Cancer Genome Atlas (TCGA). Despite technical differences (e.g., formalin-fixed vs. fresh-frozen tissue), the distribution of mutations and their variant classifications for major HGSC genes were nearly identical across study populations. However, de novo significantly mutated gene analysis identified genes not previously reported in TCGA analysis, including the oncogene KRAS and the potential tumor suppressor OBSCN. The prevalence of the homologous recombination deficiency signature was higher among Black individuals with the immunoreactive gene expression subtype compared with the mesenchymal and proliferative subtypes. These findings were confirmed by comparing the data from Black patients with those from 123 White patients with identical tissue collection and processing. Overall, this study suggests that, although most features of HGSC tumor phenotypes are similar in Black and White populations, there may be clinically relevant differences. If validated, these phenotypes may be important for clinical decision-making and would have been missed by characterizing tumors from White individuals only. Significance: Elucidation of the somatic mutational landscape of high-grade serous ovarian carcinoma in Black individuals, who experience poor survival and are underrepresented in research, could inform patient prognosis and enable precision medicine opportunities.

MAX inactivation deregulates the MYC network and induces neuroendocrine neoplasia in multiple tissues

The MYC transcription factor requires MAX for DNA binding and widespread activation of gene expression in both normal and neoplastic cells. Inactivating mutations in MAX are associated with a subset of neuroendocrine cancers including pheochromocytoma, pituitary adenoma, and small cell lung cancer. Neither the extent nor the mechanisms of MAX tumor suppression are well understood. Deleting Max across multiple mouse neuroendocrine tissues, we find that Max inactivation alone produces pituitary adenomas, while Max inactivation cooperates with Rb1/Trp53 loss to accelerate medullary thyroid C cell and pituitary adenoma development. In the thyroid tumor cell lines, MAX loss triggers a marked shift in genomic occupancy by other members of the MYC network (MNT, MLX, MondoA) supporting metabolism, survival, and proliferation of neoplastic neuroendocrine cells. Our work reveals MAX as a broad suppressor of neuroendocrine tumorigenesis through its ability to maintain a balance of genomic occupancies among the diverse transcription factors in the MYC network.

Hao-Fountain syndrome protein USP7 controls neuronal differentiation via BCOR-ncPRC1.1

Pathogenic variants in the ubiquitin-specific protease 7 (USP7) gene cause a neurodevelopmental disorder called Hao-Fountain syndrome. However, it remains unclear which of USP7's pleiotropic functions are relevant for neurodevelopment. Here, we present a combination of quantitative proteomics, transcriptomics, and epigenomics to define the USP7 regulatory circuitry during neuronal differentiation. USP7 activity is required for the transcriptional programs that direct both the differentiation of embryonic stem cells into neural stem cells and the neuronal differentiation of SH-SY5Y neuroblastoma cells. USP7 controls the dosage of the Polycomb monubiquitylated histone H2A lysine 119 (H2AK119ub1) ubiquitin ligase complexes ncPRC1.1 and ncPRC1.6. Loss-of-function experiments revealed that BCOR-ncPRC1.1, but not ncPRC1.6, is a key effector of USP7 during neuronal differentiation. Indeed, BCOR-ncPRC1.1 mediates a major portion of USP7-dependent gene regulation during this process. Besides providing a detailed map of the USP7 regulome during neurodifferentiation, our results suggest that USP7- and ncPRC1.1-associated neurodevelopmental disorders involve dysregulation of a shared epigenetic network.

L3MBTL2 maintains MYCN-amplified neuroblastoma cell proliferation through silencing NRIP3 and BRME1 genes

Epigenetic alterations critically affect tumor development. Polycomb-group complexes constitute an evolutionarily conserved epigenetic machinery that regulates stem cell fate and development. They are implicated in tumorigenesis, primarily via histone modification. Polycomb repressive complex 1 (PRC1) complexes 1-6 (PRC1.1-6) mediate the ubiquitination of histone H2A on lysine 119 (H2AK119ub). Here, we studied the functional roles of a PRC1.6 molecule, L3MBTL2, in neuroblastoma (NB) cells. L3MBTL2-knockout and knockdown revealed that L3MBTL2 depletion suppressed NB cell proliferation via cell-cycle arrest and gamma-H2A.X upregulation. L3MBTL2-knockout profoundly suppressed xenograft tumor formation. Transcriptome analysis revealed suppressed cell-cycle-related and activated differentiation-related pathways. Break repair meiotic recombinase recruitment factor 1 (BRME1) and nuclear receptor interacting protein 3 (NRIP3) were notably de-repressed by L3MBTL2-knockout. The deletion of L3MBTL2 reduced enrichment of H2AK119ub and PCGF6 at transcriptional start site proximal regions of the targets. Add-back studies unveiled the importance of L3MBTL2-BRME1 and -NRIP3 axes for NB cell proliferation. We further manifested the association of MYCN with de-repression of NRIP3 in an L3MBTL2-deficient context. Therefore, this study first revealed the significance of L3MBTL2-mediated gene silencing in MYCN-amplified NB cells.

In vivo CRISPR screens identify Mga as an immunotherapy target in triple-negative breast cancer

Immune evasion is not only critical for tumor initiation and progression, but also determines the efficacy of immunotherapies. Through iterative in vivo CRISPR screens with seven syngeneic tumor models, we identified core and context-dependent immune evasion pathways across cancer types. This valuable high-confidence dataset is available for the further understanding of tumor intrinsic immunomodulators, which may lead to the discovery of effective anticancer therapeutic targets. With a focus on triple-negative breast cancer (TNBC), we found that Mga knock-out significantly enhances antitumor immunity and inhibits tumor growth. Transcriptomics and single-cell RNA sequencing analyses revealed that Mga influences various immune-related pathways in the tumor microenvironment. Our findings suggest that Mga may play a role in modulating the tumor immune landscape, though the precise mechanisms require further investigation. Interestingly, we observed that low MGA expression in breast cancer patients correlates with a favorable prognosis, particularly in those with active interferon-γ signaling. These observations provide insights into tumor immune escape mechanisms and suggest that further exploration of MGA's function could potentially lead to effective therapeutic strategies in TNBC.

MAX inactivation deregulates the MYC network and induces neuroendocrine neoplasia in multiple tissues

The MYC transcription factor requires MAX for DNA binding and widespread activation of gene expression in both normal and neoplastic cells. Surprisingly, inactivating mutations in MAX are associated with a subset of neuroendocrine cancers including pheochromocytoma, pituitary adenoma and small cell lung cancer. Neither the extent nor the mechanisms of MAX tumor suppression are well understood. Delet-ing Max across multiple mouse neuroendocrine tissues, we find Max inactivation alone produces pituitary adenomas while Max loss cooperates with Rb1/Trp53 loss to accelerate medullary thyroid C-cell and pituitary adenoma development. In the thyroid tumor cell lines, MAX loss triggers a striking shift in genomic occupancy by other members of the MYC network (MNT, MLX, MondoA) supporting metabolism, survival and proliferation of neoplastic neuroendocrine cells. Our work reveals MAX as a broad suppressor of neuroendocrine tumorigenesis through its ability to maintain a balance of genomic occupancies among the diverse transcription factors in the MYC network.

Development of an epigenetic clock resistant to changes in immune cell composition

Epigenetic clocks are age predictors that use machine-learning models trained on DNA CpG methylation values to predict chronological or biological age. Increases in predicted epigenetic age relative to chronological age (epigenetic age acceleration) are connected to aging-associated pathologies, and changes in epigenetic age are linked to canonical aging hallmarks. However, epigenetic clocks rely on training data from bulk tissues whose cellular composition changes with age. Here, we found that human naive CD8+ T cells, which decrease in frequency during aging, exhibit an epigenetic age 15-20 years younger than effector memory CD8+ T cells from the same individual. Importantly, homogenous naive T cells isolated from individuals of different ages show a progressive increase in epigenetic age, indicating that current epigenetic clocks measure two independent variables, aging and immune cell composition. To isolate the age-associated cell intrinsic changes, we created an epigenetic clock, the IntrinClock, that did not change among 10 immune cell types tested. IntrinClock shows a robust predicted epigenetic age increase in a model of replicative senescence in vitro and age reversal during OSKM-mediated reprogramming.

The genomic landscape of lung cancer in never-smokers from the Women's Health Initiative

Over 200,000 individuals are diagnosed with lung cancer in the United States every year, with a growing proportion of cases, especially lung adenocarcinoma, occurring in individuals who have never smoked. Women over the age of 50 comprise the largest affected demographic. To understand the genomic drivers of lung adenocarcinoma and therapeutic response in this population, we performed whole genome and/or whole exome sequencing on 73 matched lung tumor/normal pairs from postmenopausal women who participated in the Women's Health Initiative. Somatic copy number alterations showed little variation by smoking status, suggesting that aneuploidy may be a general characteristic of lung cancer regardless of smoke exposure. Similarly, clock-like and APOBEC mutation signatures were prevalent but did not differ in tumors from smokers and never-smokers. However, mutations in both EGFR and KRAS showed unique allelic differences determined by smoking status that are known to alter tumor response to targeted therapy. Mutations in the MYC-network member MGA were more prevalent in tumors from smokers. Fusion events in ALK, RET, and ROS1 were absent, likely due to age-related differences in fusion prevalence. Our work underscores the profound effect of smoking status, age, and sex on the tumor mutational landscape and identifies areas of unmet medical need.

Functional characterization of cooperating MGA mutations in RUNX1::RUNX1T1 acute myeloid leukemia

MGA (Max-gene associated) is a dual-specificity transcription factor that negatively regulates MYC-target genes to inhibit proliferation and promote differentiation. Loss-of-function mutations in MGA have been commonly identified in several hematological neoplasms, including acute myeloid leukemia (AML) with RUNX1::RUNX1T1, however, very little is known about the impact of these MGA alterations on normal hematopoiesis or disease progression. We show that representative MGA mutations identified in patient samples abolish protein-protein interactions and transcriptional activity. Using a series of human and mouse model systems, including a newly developed conditional knock-out mouse strain, we demonstrate that loss of MGA results in upregulation of MYC and E2F targets, cell cycle genes, mTOR signaling, and oxidative phosphorylation in normal hematopoietic cells, leading to enhanced proliferation. The loss of MGA induces an open chromatin state at promoters of genes involved in cell cycle and proliferation. RUNX1::RUNX1T1 expression in Mga-deficient murine hematopoietic cells leads to a more aggressive AML with a significantly shortened latency. These data show that MGA regulates multiple pro-proliferative pathways in hematopoietic cells and cooperates with the RUNX1::RUNX1T1 fusion oncoprotein to enhance leukemogenesis.

Mutational processes of tobacco smoking and APOBEC activity generate protein-truncating mutations in cancer genomes

Mutational signatures represent a genomic footprint of endogenous and exogenous mutational processes through tumor evolution. However, their functional impact on the proteome remains incompletely understood. We analyzed the protein-coding impact of single-base substitution (SBS) signatures in 12,341 cancer genomes from 18 cancer types. Stop-gain mutations (SGMs) (i.e., nonsense mutations) were strongly enriched in SBS signatures of tobacco smoking, APOBEC cytidine deaminases, and reactive oxygen species. These mutational processes alter specific trinucleotide contexts and thereby substitute serines and glutamic acids with stop codons. SGMs frequently affect cancer hallmark pathways and tumor suppressors such as TP53, FAT1, and APC. Tobacco-driven SGMs in lung cancer correlate with smoking history and highlight a preventable determinant of these harmful mutations. APOBEC-driven SGMs are enriched in YTCA motifs and associate with APOBEC3A expression. Our study exposes SGM expansion as a genetic mechanism by which endogenous and carcinogenic mutational processes directly contribute to protein loss of function, oncogenesis, and tumor heterogeneity.

Functional Characterization of Cooperating MGA Mutations in RUNX1::RUNX1T1 Acute Myeloid Leukemia

MGA (Max-gene associated) is a dual-specificity transcription factor that negatively regulates MYC-target genes to inhibit proliferation and promote differentiation. Loss-of-function mutations in MGA have been commonly identified in several hematological neoplasms, including acute myeloid leukemia (AML) with RUNX1::RUNX1T1, however, very little is known about the impact of these MGA alterations on normal hematopoiesis or disease progression. We show that representative MGA mutations identified in patient samples abolish protein-protein interactions and transcriptional activity. Using a series of human and mouse model systems, including a newly developed conditional knock-out mouse strain, we demonstrate that loss of MGA results in upregulation of MYC and E2F targets, cell cycle genes, mTOR signaling, and oxidative phosphorylation in normal hematopoietic cells, leading to enhanced proliferation. The loss of MGA induces an open chromatin state at promotors of genes involved in cell cycle and proliferation. RUNX1::RUNX1T1 expression in Mga-deficient murine hematopoietic cells leads to a more aggressive AML with a significantly shortened latency. These data show that MGA regulates multiple pro-proliferative pathways in hematopoietic cells and cooperates with the RUNX1::RUNX1 T1 fusion oncoprotein to enhance leukemogenesis.

The evolution of two transmissible cancers in Tasmanian devils

Tasmanian devils have spawned two transmissible cancer lineages, named devil facial tumor 1 (DFT1) and devil facial tumor 2 (DFT2). We investigated the genetic diversity and evolution of these clones by analyzing 78 DFT1 and 41 DFT2 genomes relative to a newly assembled, chromosome-level reference. Time-resolved phylogenetic trees reveal that DFT1 first emerged in 1986 (1982 to 1989) and DFT2 in 2011 (2009 to 2012). Subclone analysis documents transmission of heterogeneous cell populations. DFT2 has faster mutation rates than DFT1 across all variant classes, including substitutions, indels, rearrangements, transposable element insertions, and copy number alterations, and we identify a hypermutated DFT1 lineage with defective DNA mismatch repair. Several loci show plausible evidence of positive selection in DFT1 or DFT2, including loss of chromosome Y and inactivation of MGA, but none are common to both cancers. This study reveals the parallel long-term evolution of two transmissible cancers inhabiting a common niche in Tasmanian devils.

In Vivo Modeling of CLL Transformation to Richter Syndrome Reveals Convergent Evolutionary Paths and Therapeutic Vulnerabilities

Transformation to aggressive disease histologies generates formidable clinical challenges across cancers, but biological insights remain few. We modeled the genetic heterogeneity of chronic lymphocytic leukemia (CLL) through multiplexed in vivo CRISPR-Cas9 B-cell editing of recurrent CLL loss-of-function drivers in mice and recapitulated the process of transformation from indolent CLL into large cell lymphoma [i.e., Richter syndrome (RS)]. Evolutionary trajectories of 64 mice carrying diverse combinatorial gene assortments revealed coselection of mutations in Trp53, Mga, and Chd2 and the dual impact of clonal Mga/Chd2 mutations on E2F/MYC and interferon signaling dysregulation. Comparative human and murine RS analyses demonstrated tonic PI3K signaling as a key feature of transformed disease, with constitutive activation of the AKT and S6 kinases, downmodulation of the PTEN phosphatase, and convergent activation of MYC/PI3K transcriptional programs underlying enhanced sensitivity to MYC/mTOR/PI3K inhibition. This robust experimental system presents a unique framework to study lymphoid biology and therapy.

SIGNIFICANCE

Mouse models reflective of the genetic complexity and heterogeneity of human tumors remain few, including those able to recapitulate transformation to aggressive disease histologies. Herein, we model CLL transformation into RS through multiplexed in vivo gene editing, providing key insight into the pathophysiology and therapeutic vulnerabilities of transformed disease. This article is highlighted in the In This Issue feature, p. 101.

Identification of a novel target of SETD1A histone methyltransferase and the clinical significance in pancreatic cancer

Although histone H3K4 methyltransferase SETD1A is overexpressed in various cancer types, the molecular mechanism underlying its overexpression and its target genes in pancreatic ductal adenocarcinoma (PDAC) remain unclarified. We conducted immunohistochemical staining for SETD1A in 105 human PDAC specimens to assess the relationship between SETD1A overexpression and clinicopathological features. The function and target genes of SETD1A were investigated using human pancreatic cancer cell lines. SETD1A expression was upregulated in 51.4% of patients with PDAC and was an independent prognostic factor associated with shorter disease-free survival after resection (p < 0.05). Knockdown and overexpression of SETD1A showed that SETD1A plays a crucial role in increasing the proliferation and motility of PDAC cells. SETD1A overexpression increased tumorigenicity. RNA sequencing of SETD1A-knockdown cells revealed downregulation of RUVBL1, an oncogenic protein ATP-dependent DNA helicase gene. ChIP analysis revealed that SETD1A binds to the RUVBL1 promoter region, resulting in increased H3K4me3 levels. Knockdown of RUVBL1 showed inhibition of cell proliferation, migration, and invasion of PDAC cells, which are similar biological effects to SETD1A knockdown. High expression of both SETD1A and RUVBL1 was an independent prognostic factor not only for disease-free survival but also for overall survival (p < 0.05). In conclusion, we identified RUVBL1 as a novel downstream target gene of the SETD1A-H3K4me3 pathway. Co-expression of SETD1A and RUVBL1 is an important factor for predicting the prognosis of patients with PDAC.

Changes of Mutations and Copy-Number and Enhanced Cell Migration during Breast Tumorigenesis

Although cancer stem cells (CSCs) play a major role in tumorigenesis and metastasis, the role of genetic alterations in invasiveness of CSCs is still unclear. Tumor microenvironment signals, such as extracellular matrix (ECM) composition, significantly influence cell behaviors. Unfortunately, these signals are often lost in in vitro cell culture. This study determines putative CSC populations, examines genetic changes during tumorigenesis of human breast epithelial stem cells, and investigates single-cell migration properties on ECM-mimetic platforms. Whole exome sequencing data indicate that tumorigenic cells have a higher somatic mutation burden than non-tumorigenic cells, and that mutations exclusive to tumorigenic cells exhibit higher predictive deleterious scores. Tumorigenic cells exhibit distinct somatic copy number variations (CNVs) including gain of duplications in chromosomes 5 and 8. ECM-mimetic topography selectively enhances migration speed of tumorigenic cells, but not of non-tumorigenic cells, and results in a wide distribution of tumorigenic single-cell migration speeds, suggesting heterogeneity in cellular sensing of contact guidance cues. This study identifies mutations and CNVs acquired during breast tumorigenesis, which can be associated with enhanced migration of breast tumorigenic cells, and demonstrates that a nanotopographically-defined platform can be applied to recapitulate an ECM structure for investigating cellular migration in the simulated tumor microenvironment.

Cannabinoid receptor 2 plays a pro-tumorigenic role in non-small cell lung cancer by limiting anti-tumor activity of CD8+ T and NK cells

Cannabinoid (CB) receptors (CB1 and CB2) are expressed on cancer cells and their expression influences carcinogenesis in various tumor entities. Cells of the tumor microenvironment (TME) also express CB receptors, however, their role in tumor development is still unclear. We, therefore, investigated the role of TME-derived CB1 and CB2 receptors in a model of non-small cell lung cancer (NSCLC). Leukocytes in the TME of mouse and human NSCLC express CB receptors, with CB2 showing higher expression than CB1. In the tumor model, using CB1- (CB1 -/-) and CB2-knockout (CB2 -/-) mice, only deficiency of CB2, but not of CB1, resulted in reduction of tumor burden vs. wild type (WT) littermates. This was accompanied by increased accumulation and tumoricidal activity of CD8+ T and natural killer cells, as well as increased expression of programmed death-1 (PD-1) and its ligand on lymphoid and myeloid cells, respectively. CB2 -/- mice responded significantly better to anti-PD-1 therapy than WT mice. The treatment further increased infiltration of cytotoxic lymphocytes into the TME of CB2 -/- mice. Our findings demonstrate that TME-derived CB2 dictates the immune cell recruitment into tumors and the responsiveness to anti-PD-1 therapy in a model of NSCLC. CB2 could serve as an adjuvant target for immunotherapy.

Somatic mutation but not aneuploidy differentiates lung cancer in never-smokers and smokers

Lung cancer in never-smokers disproportionately affects older women. To understand the mutational landscape of this cohort, we performed detailed genome characterization of 73 lung adenocarcinomas from participants of the Women’s Health Initiative (WHI). We find enrichment of EGFR mutations in never-/light-smokers and KRAS mutations in heavy smokers as expected, but we also show that the specific variants of these genes differ by smoking status, with important therapeutic implications. Mutational signature analysis revealed signatures of clock, APOBEC, and DNA repair deficiency in never-/light-smokers; however, the mutational load of these signatures did not differ significantly from those found in smokers. Last, tumors from both smokers and never-/light-smokers shared copy number subtypes, with no significant differences in aneuploidy. Thus, the genomic landscape of lung cancer in never-/light-smokers and smokers is predominantly differentiated by somatic mutations and not copy number alterations.

Genomic Characterization of Prostatic Basal Cell Carcinoma

Basal cell carcinoma (BCC) of the prostate is a rare tumor. Compared with the more common acinar adenocarcinoma (AAC) of the prostate, BCCs show features of basal cell differentiation and are thought to be biologically distinct from AAC. The spectrum of molecular alterations of BCC has not been comprehensively described, and genomic studies are lacking. Herein, whole genome sequencing was performed on archival formalin-fixed, paraffin-embedded specimens of two cases with BCC. Prostatic BCCs were characterized by an overall low copy number and mutational burden. Recurrent copy number loss of chromosome 16 was observed. In addition, putative driver gene alterations in KIT, DENND3, PTPRU, MGA, and CYLD were identified. Mechanistically, depletion of the CYLD protein resulted in increased proliferation of prostatic basal cells in vitro. Collectively, these studies show that prostatic BCC displays distinct genomic alterations from AAC and highlight a potential role for loss of chromosome 16 in the pathogenesis of this rare tumor type.

USP7 regulates the ncPRC1 Polycomb axis to stimulate genomic H2AK119ub1 deposition uncoupled from H3K27me3

Ubiquitin-specific protease 7 (USP7) has been implicated in cancer progression and neurodevelopment. However, its molecular targets remain poorly characterized. We combined quantitative proteomics, transcriptomics, and epigenomics to define the core USP7 network. Our multi-omics analysis reveals USP7 as a control hub that links genome regulation, tumor suppression, and histone H2A ubiquitylation (H2AK119ub1) by noncanonical Polycomb-repressive complexes (ncPRC1s). USP7 strongly stabilizes ncPRC1.6 and, to a lesser extent, ncPRC1.1. Moreover, USP7 represses expression of AUTS2, which suppresses H2A ubiquitylation by ncPRC1.3/5. Collectively, these USP7 activities promote the genomic deposition of H2AK119ub1 by ncPRC1, especially at transcriptionally repressed loci. Notably, USP7-dependent changes in H2AK119ub1 levels are uncoupled from H3K27me3. Even complete loss of the PRC1 catalytic core and H2AK119ub1 has only a limited effect on H3K27me3. Besides defining the USP7 regulome, our results reveal that H2AK119ub1 dosage is largely disconnected from H3K27me3.

Atlas of regulated target genes of transcription factors (ART-TF) in human ES cells

Background

Transcription factors (TFs) play central roles in maintaining “stemness” of embryonic stem (ES) cells and their differentiation into several hundreds of adult cell types. The regulatory competence of TFs is routinely assessed by detecting target genes to which they bind. However, these data do not indicate which target genes are activated, repressed, or not affected by the change of TF abundance. There is a lack of large-scale studies that compare the genome binding of TFs with the expression change of target genes after manipulation of each TF.

Results

In this paper we associated human TFs with their target genes by two criteria: binding to genes, evaluated from published ChIP-seq data (n = 1868); and change of target gene expression shortly after induction of each TF in human ES cells. Lists of direction- and strength-specific regulated target genes are generated for 311 TFs (out of 351 TFs tested) with expected proportion of false positives less than or equal to 0.30, including 63 new TFs not present in four existing databases of target genes. Our lists of direction-specific targets for 152 TFs (80.0%) are larger that in the TRRUST database. In average, 30.9% of genes that respond greater than or equal to twofold to the induction of TFs are regulated targets. Regulated target genes indicate that the majority of TFs are either strong activators or strong repressors, whereas sets of genes that responded greater than or equal to twofold to the induction of TFs did not show strong asymmetry in the direction of expression change. The majority of human TFs (82.1%) regulated their target genes primarily via binding to enhancers. Repression of target genes is more often mediated by promoter-binding than activation of target genes. Enhancer-promoter loops are more abundant among strong activator and repressor TFs.

Conclusions

We developed an atlas of regulated targets of TFs (ART-TF) in human ES cells by combining data on TF binding with data on gene expression change after manipulation of individual TFs. Sets of regulated gene targets were identified with a controlled rate of false positives. This approach contributes to the understanding of biological functions of TFs and organization of gene regulatory networks. This atlas should be a valuable resource for ES cell-based regenerative medicine studies.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12859-022-04924-3.

Polycomb group ring finger protein 6 suppresses Myc-induced lymphomagenesis

Max dimerizes with Mga to form the repressive complex PRC1.6; another PRC1.6 subunit, Pcgf6, suppresses Myc-induced lymphomagenesis but, unexpectedly, does so in a Mga- and PRC1.6-independent manner.

Max is an obligate dimerization partner for the Myc transcription factors and for several repressors, such as Mnt, Mxd1-4, and Mga, collectively thought to antagonize Myc function in transcription and oncogenesis. Mga, in particular, is part of the variant Polycomb group repressive complex PRC1.6. Here, we show that ablation of the distinct PRC1.6 subunit Pcgf6–but not Mga–accelerates Myc-induced lymphomagenesis in Eµ-myc transgenic mice. Unexpectedly, however, Pcgf6 loss shows no significant impact on transcriptional profiles, in neither pre-tumoral B-cells, nor lymphomas. Altogether, these data unravel an unforeseen, Mga- and PRC1.6-independent tumor suppressor activity of Pcgf6.

Mechanisms of Polycomb group protein function in cancer

Cancer arises from a multitude of disorders resulting in loss of differentiation and a stem cell-like phenotype characterized by uncontrolled growth. Polycomb Group (PcG) proteins are members of multiprotein complexes that are highly conserved throughout evolution. Historically, they have been described as essential for maintaining epigenetic cellular memory by locking homeotic genes in a transcriptionally repressed state. What was initially thought to be a function restricted to a few target genes, subsequently turned out to be of much broader relevance, since the main role of PcG complexes is to ensure a dynamically choregraphed spatio-temporal regulation of their numerous target genes during development. Their ability to modify chromatin landscapes and refine the expression of master genes controlling major switches in cellular decisions under physiological conditions is often misregulated in tumors. Surprisingly, their functional implication in the initiation and progression of cancer may be either dependent on Polycomb complexes, or specific for a subunit that acts independently of other PcG members. In this review, we describe how misregulated Polycomb proteins play a pleiotropic role in cancer by altering a broad spectrum of biological processes such as the proliferation-differentiation balance, metabolism and the immune response, all of which are crucial in tumor progression. We also illustrate how interfering with PcG functions can provide a powerful strategy to counter tumor progression.

Normal and Neoplastic Growth Suppression by the Extended Myc Network

Among the first discovered and most prominent cellular oncogenes is MYC, which encodes a bHLH-ZIP transcription factor (Myc) that both activates and suppresses numerous genes involved in proliferation, energy production, metabolism and translation. Myc belongs to a small group of bHLH-ZIP transcriptional regulators (the Myc Network) that includes its obligate heterodimerization partner Max and six "Mxd proteins" (Mxd1-4, Mnt and Mga), each of which heterodimerizes with Max and largely opposes Myc's functions. More recently, a second group of bHLH-ZIP proteins (the Mlx Network) has emerged that bears many parallels with the Myc Network. It is comprised of the Myc-like factors ChREBP and MondoA, which, in association with the Max-like member Mlx, regulate smaller and more functionally restricted repertoires of target genes, some of which are shared with Myc. Opposing ChREBP and MondoA are heterodimers comprised of Mlx and Mxd1, Mxd4 and Mnt, which also structurally and operationally link the two Networks. We discuss here the functions of these "Extended Myc Network" members, with particular emphasis on their roles in suppressing normal and neoplastic growth. These roles are complex due to the temporal- and tissue-restricted expression of Extended Myc Network proteins in normal cells, their regulation of both common and unique target genes and, in some cases, their functional redundancy.

Polycomb group proteins in cancer: multifaceted functions and strategies for modulation

Polycomb repressive complexes (PRCs) are a heterogenous collection of dozens, if not hundreds, of protein complexes composed of various combinations of subunits. PRCs are transcriptional repressors important for cell-type specificity during development, and as such, are commonly mis-regulated in cancer. PRCs are broadly characterized as PRC1 with histone ubiquitin ligase activity, or PRC2 with histone methyltransferase activity; however, the mechanism by which individual PRCs, particularly the highly diverse set of PRC1s, alter gene expression has not always been clear. Here we review the current understanding of how PRCs act, both individually and together, to establish and maintain gene repression, the biochemical contribution of individual PRC subunits, the mis-regulation of PRC function in different cancers, and the current strategies for modulating PRC activity. Increased mechanistic understanding of PRC function, as well as cancer-specific roles for individual PRC subunits, will uncover better targets and strategies for cancer therapies.

The Multiple Faces of MNT and Its Role as a MYC Modulator

MNT is a crucial modulator of MYC, controls several cellular functions, and is activated in most human cancers. It is the largest, most divergent, and most ubiquitously expressed protein of the MXD family. MNT was first described as a MYC antagonist and tumor suppressor. Indeed, 10% of human tumors present deletions of one MNT allele. However, some reports show that MNT functions in cooperation with MYC by maintaining cell proliferation, promoting tumor cell survival, and supporting MYC-driven tumorigenesis in cellular and animal models. Although MAX was originally considered MNT's obligate partner, our recent findings demonstrate that MNT also works independently. MNT forms homodimers and interacts with proteins both outside and inside of the proximal MYC network. These complexes are involved in a wide array of cellular processes, from transcriptional repression via SIN3 to the modulation of metabolism through MLX as well as immunity and apoptosis via REL. In this review, we discuss the present knowledge of MNT with a special focus on its interactome, which sheds light on the complex and essential role of MNT in cell biology.