Transcription factor antagonism regulates heterogeneity in embryonic stem cell states.

Gene expression heterogeneity underlies cell states and contributes to developmental robustness. While heterogeneity can arise from stochastic transcriptional processes, the extent to which it is regulated is unclear. Here, we characterize the regulatory program underlying heterogeneity in murine embryonic stem cell (mESC) states. We identify differentially active and transcribed enhancers (DATEs) across states. DATEs regulate differentially expressed genes and are distinguished by co-binding of transcription factors Klf4 and Zfp281. In contrast to other factors that interact in a positive feedback network stabilizing mESC cell-type identity, Klf4 and Zfp281 drive opposing transcriptional and chromatin programs. Abrogation of factor binding to DATEs dampens variation in gene expression, and factor loss alters kinetics of switching between states. These results show antagonism between factors at enhancers results in gene expression heterogeneity and formation of cell states, with implications for the generation of diverse cell types during development.

MITF deficiency accelerates GNAQ-driven uveal melanoma.

Cutaneous melanoma (CM) and uveal melanoma (UM) both originate from the melanocytic lineage but are primarily driven by distinct oncogenic drivers, BRAF/NRAS or GNAQ/GNA11, respectively. The melanocytic master transcriptional regulator, MITF, is essential for both CM development and maintenance, but its role in UM is largely unexplored. Here, we use zebrafish models to dissect the key UM oncogenic signaling events and establish the role of MITF in UM tumors. Using a melanocytic lineage expression system, we showed that patient-derived mutations of GNAQ (GNAQQ209L) or its upstream CYSLTR2 receptor (CYSLTR2L129Q) both drive UM when combined with a cooperating mutation, tp53M214K/M214K. The tumor-initiating potential of the major GNAQ/11 effector pathways, YAP, and phospholipase C-ß (PLCß)­ERK was also investigated in this system and thus showed that while activated YAP (YAPAA) induced UM with high potency, the patient-derived PLCß4 mutation (PLCB4D630Y) very rarely yielded UM tumors in the tp53M214K/M214K context. Remarkably, mitfa deficiency was profoundly UM promoting, dramatically accelerating the onset and progression of tumors induced by Tg(mitfa:GNAQQ209L);tp53M214K/M214K or Tg(mitfa:CYSLTR2L129Q);tp53M214K/M214K. Moreover, mitfa loss was sufficient to cooperate with GNAQQ209L to drive tp53­wild type UM development and allowed Tg(mitfa:PLCB4D630Y);tp53M214K/M214K melanocyte lineage cells to readily form tumors. Notably, all of the mitfa−/− UM tumors, including those arising in Tg(mitfa:PLCB4D630Y);tp53M214K/M214K;mitfa−/− zebrafish, displayed nuclear YAP while lacking hyperactive ERK indicative of PLCß signaling. Collectively, these data show that YAP signaling is the major mediator of UM and that MITF acts as a bona fide tumor suppressor in UM in direct opposition to its essential role in CM.

EMT and primary ciliogenesis: For better or worse in sickness and in health.

Epithelial-mesenchymal transition (EMT) and primary ciliogenesis are two cell-biological programs that are essential for development of multicellular organisms and whose abnormal regulation results in many diseases (i.e., developmental anomalies and cancers). Emerging studies suggest an intricate interplay between these two processes. Here, we discuss physiological and pathological contexts in which their interconnections promote normal development or disease progression. We describe underlying molecular mechanisms of the interplay and EMT/ciliary signaling axes that influence EMT-related processes (i.e., stemness, motility and invasion). Understanding the molecular and cellular mechanisms of the relationship between EMT and primary ciliogenesis may provide new insights in the etiology of diseases related to EMT and cilia dysfunction.

Acquired resistance to PRMT5 inhibition induces concomitant collateral sensitivity to paclitaxel.

Epigenetic regulators play key roles in cancer and are increasingly being targeted for treatment. However, for many, little is known about mechanisms of resistance to the inhibition of these regulators. We have generated a model of resistance to inhibitors of protein arginine methyltransferase 5 (PRMT5). This study was conducted in KrasG12D;Tp53-null lung adenocarcinoma (LUAD) cell lines. Resistance to PRMT5 inhibitors (PRMT5i) arose rapidly, and barcoding experiments showed that this resulted from a drug-induced transcriptional state switch, not selection of a preexisting population. This resistant state is both stable and conserved across variants arising from distinct LUAD lines. Moreover, it brought with it vulnerabilities to other chemotherapeutics, especially the taxane paclitaxel. This paclitaxel sensitivity depended on the presence of stathmin 2 (STMN2), a microtubule regulator that is specifically expressed in the resistant state. Remarkably, STMN2 was also essential for resistance to PRMT5 inhibition. Thus, a single gene is required for both acquisition of resistance to PRMT5i and collateral sensitivity to paclitaxel in our LUAD cells. Accordingly, the combination of PRMT5i and paclitaxel yielded potent and synergistic killing of the murine LUAD cells. Importantly, the synergy between PRMT5i and paclitaxel also extended to human cancer cell lines. Finally, analysis of The Cancer Genome Atlas patient data showed that high STMN2 levels correlate with complete regression of tumors in response to taxane treatment. Collectively, this study reveals a recurring mechanism of PRMT5i resistance in LUAD and identifies collateral sensitivities that have potential clinical relevance.

BMI1 induces an invasive signature in melanoma that promotes metastasis and chemoresistance.

Melanoma can switch between proliferative and invasive states, which have identifying gene expression signatures that correlate with good and poor prognosis, respectively. However, the mechanisms controlling these signatures are poorly understood. In this study, we identify BMI1 as a key determinant of melanoma metastasis by which its overexpression enhanced and its deletion impaired dissemination. Remarkably, in this tumor type, BMI1 had no effect on proliferation or primary tumor growth but enhanced every step of the metastatic cascade. Consistent with the broad spectrum of effects, BMI1 activated widespread gene expression changes, which are characteristic of melanoma progression and also chemoresistance. Accordingly, we showed that up-regulation or down-regulation of BMI1 induced resistance or sensitivity to BRAF inhibitor treatment and that induction of noncanonical Wnt by BMI1 is required for this resistance. Finally, we showed that our BMI1-induced gene signature encompasses all of the hallmarks of the previously described melanoma invasive signature. Moreover, our signature is predictive of poor prognosis in human melanoma and is able to identify primary tumors that are likely to become metastatic. These data yield key insights into melanoma biology and establish BMI1 as a compelling drug target whose inhibition would suppress both metastasis and chemoresistance of melanoma.

Targeting the cyclin-dependent kinase 5 in metastatic melanoma.

The cyclin-dependent kinase 5 (CDK5), originally described as a neuronal-specific kinase, is also frequently activated in human cancers. Using conditional CDK5 knockout mice and a mouse model of highly metastatic melanoma, we found that CDK5 is dispensable for the growth of primary tumors. However, we observed that ablation of CDK5 completely abrogated the metastasis, revealing that CDK5 is essential for the metastatic spread. In mouse and human melanoma cells CDK5 promotes cell invasiveness by directly phosphorylating an intermediate filament protein, vimentin, thereby inhibiting assembly of vimentin filaments. Chemical inhibition of CDK5 blocks the metastatic spread of patient-derived melanomas in patient-derived xenograft (PDX) mouse models. Hence, inhibition of CDK5 might represent a very potent therapeutic strategy to impede the metastatic dissemination of malignant cells.

Proteomic analysis of pRb loss highlights a signature of decreased mitochondrial oxidative phosphorylation.

The retinoblastoma tumor suppressor (pRb) protein associates with chromatin and regulates gene expression. Numerous studies have identified Rb-dependent RNA signatures, but the proteomic effects of Rb loss are largely unexplored. We acutely ablated Rb in adult mice and conducted a quantitative analysis of RNA and proteomic changes in the colon and lungs, where Rb(KO) was sufficient or insufficient to induce ectopic proliferation, respectively. As expected, Rb(KO) caused similar increases in classic pRb/E2F-regulated transcripts in both tissues, but, unexpectedly, their protein products increased only in the colon, consistent with its increased proliferative index. Thus, these protein changes induced by Rb loss are coupled with proliferation but uncoupled from transcription. The proteomic changes in common between Rb(KO) tissues showed a striking decrease in proteins with mitochondrial functions. Accordingly, RB1 inactivation in human cells decreased both mitochondrial mass and oxidative phosphorylation (OXPHOS) function. RB(KO) cells showed decreased mitochondrial respiratory capacity and the accumulation of hypopolarized mitochondria. Additionally, RB/Rb loss altered mitochondrial pyruvate oxidation from (13)C-glucose through the TCA cycle in mouse tissues and cultured cells. Consequently, RB(KO) cells have an enhanced sensitivity to mitochondrial stress conditions. In summary, proteomic analyses provide a new perspective on Rb/RB1 mutation, highlighting the importance of pRb for mitochondrial function and suggesting vulnerabilities for treatment.

The retinoblastoma protein induces apoptosis directly at the mitochondria.

The retinoblastoma protein gene RB-1 is mutated in one-third of human tumors. Its protein product, pRB (retinoblastoma protein), functions as a transcriptional coregulator in many fundamental cellular processes. Here, we report a nonnuclear role for pRB in apoptosis induction via pRB's direct participation in mitochondrial apoptosis. We uncovered this activity by finding that pRB potentiated TNFα-induced apoptosis even when translation was blocked. This proapoptotic function was highly BAX-dependent, suggesting a role in mitochondrial apoptosis, and accordingly, a fraction of endogenous pRB constitutively associated with mitochondria. Remarkably, we found that recombinant pRB was sufficient to trigger the BAX-dependent permeabilization of mitochondria or liposomes in vitro. Moreover, pRB interacted with BAX in vivo and could directly bind and conformationally activate BAX in vitro. Finally, by targeting pRB specifically to mitochondria, we generated a mutant that lacked pRB's classic nuclear roles. This mito-tagged pRB retained the ability to promote apoptosis in response to TNFα and also additional apoptotic stimuli. Most importantly, induced expression of mito-tagged pRB in Rb(-/-);p53(-/-) tumors was sufficient to block further tumor development. Together, these data establish a nontranscriptional role for pRB in direct activation of BAX and mitochondrial apoptosis in response to diverse stimuli, which is profoundly tumor-suppressive.

EMT programs promote basal mammary stem cell and tumor-initiating cell stemness by inducing primary ciliogenesis and Hedgehog signaling.

Tissue regeneration relies on adult stem cells (SCs) that possess the ability to self-renew and produce differentiating progeny. In an analogous manner, the development of certain carcinomas depends on a small subset of tumor cells, called "tumor-initiating cells" (TICs), with SC-like properties. Mammary SCs (MaSCs) reside in the basal compartment of the mammary epithelium, and their neoplastic counterparts, mammary TICs (MaTICs), are thought to serve as the TICs for the claudin-low subtype of breast cancer. MaSCs and MaTICs both use epithelial-mesenchymal transition (EMT) programs to acquire SC properties, but the mechanism(s) connecting EMT programs to stemness remain unclear. Here we show that this depends on primary cilia, which are nonmotile, cell-surface structures that serve as platforms for receiving cues and enable activation of various signaling pathways. We show that MaSC and MaTIC EMT programs induce primary cilia formation and Hedgehog (Hh) signaling, which has previously been implicated in both MaSC and MaTIC function. Moreover, ablation of these primary cilia is sufficient to repress Hh signaling, the stemness of MaSCs, and the tumor-forming potential of MaTICs. Together, our findings establish primary ciliogenesis and consequent Hh signaling as a key mechanism by which MaSC and MaTIC EMT programs promote stemness and thereby support mammary tissue outgrowth and tumors of basal origin.

The Rb tumor suppressor regulates epithelial cell migration and polarity.

Altered cell polarity and migration are hallmarks of cancer and metastases. Here we show that inactivation of the retinoblastoma gene (Rb) tumor suppressor causes defects in tissue closure that reflect the inability of Rb null epithelial cells to efficiently migrate and polarize. These defects occur independently of pRB's anti-proliferative role and instead correlate with upregulation of RhoA signaling and mislocalization of apical-basal polarity proteins. Notably, concomitant inactivation of tp53 specifically overrides the motility defect, and not the aberrant polarity, thereby uncovering previously unappreciated mechanisms by which Rb and tp53 mutations cooperate to promote cancer development and metastases.

A homozygous deleterious CDK10 mutation in a patient with agenesis of corpus callosum, retinopathy, and deafness.

The primary cilium is a key organelle in numerous physiological and developmental processes. Genetic defects in the formation of this non-motile structure, in its maintenance and function, underlie a wide array of ciliopathies in human, including craniofacial, brain and heart malformations, and retinal and hearing defects. We used exome sequencing to study the molecular basis of disease in an 11-year-old female patient who suffered from growth retardation, global developmental delay with absent speech acquisition, agenesis of corpus callosum and paucity of white matter, sensorineural deafness, retinitis pigmentosa, vertebral anomalies, patent ductus arteriosus, and facial dysmorphism reminiscent of STAR syndrome, a suspected ciliopathy. A homozygous variant, c.870_871del, was identified in the CDK10 gene, predicted to cause a frameshift, p.Trp291Alafs*18, in the cyclin-dependent kinase 10 protein. CDK10 mRNAs were detected in patient cells and do not seem to undergo non-sense mediated decay. CDK10 is the binding partner of Cyclin M (CycM) and CDK10/CycM protein kinase regulates ciliogenesis and primary cilium elongation. Notably, CycM gene is mutated in patients with STAR syndrome. Following incubation, the patient cells appeared less elongated and more densely populated than the control cells suggesting that the CDK10 mutation affects the cytoskeleton. Upon starvation and staining with acetylated-tubulin, γ-tubulin, and Arl13b, the patient cells exhibited fewer and shorter cilia than control cells. These findings underscore the importance of CDK10 for the regulation of ciliogenesis. CDK10 defect is likely associated with a new form of ciliopathy phenotype; additional patients may further validate this association.

A vertebrate gene, ticrr, is an essential checkpoint and replication regulator.

Eukaryotes have numerous checkpoint pathways to protect genome fidelity during normal cell division and in response to DNA damage. Through a screen for G2/M checkpoint regulators in zebrafish, we identified ticrr (for TopBP1-interacting, checkpoint, and replication regulator), a previously uncharacterized gene that is required to prevent mitotic entry after treatment with ionizing radiation. Ticrr deficiency is embryonic-lethal in the absence of exogenous DNA damage because it is essential for normal cell cycle progression. Specifically, the loss of ticrr impairs DNA replication and disrupts the S/M checkpoint, leading to premature mitotic entry and mitotic catastrophe. We show that the human TICRR ortholog associates with TopBP1, a known checkpoint protein and a core component of the DNA replication preinitiation complex (pre-IC), and that the TICRR-TopBP1 interaction is stable without chromatin and requires BRCT motifs essential for TopBP1's replication and checkpoint functions. Most importantly, we find that ticrr deficiency disrupts chromatin binding of pre-IC, but not prereplication complex, components. Taken together, our data show that TICRR acts in association with TopBP1 and plays an essential role in pre-IC formation. It remains to be determined whether Ticrr represents the vertebrate ortholog of the yeast pre-IC component Sld3, or a hitherto unknown metazoan replication and checkpoint regulator.

Rb regulates fate choice and lineage commitment in vivo.

Mutation of the retinoblastoma gene (RB1) tumour suppressor occurs in one-third of all human tumours and is particularly associated with retinoblastoma and osteosarcoma. Numerous functions have been ascribed to the product of the human RB1 gene, the retinoblastoma protein (pRb). The best known is pRb's ability to promote cell-cycle exit through inhibition of the E2F transcription factors and the transcriptional repression of genes encoding cell-cycle regulators. In addition, pRb has been shown in vitro to regulate several transcription factors that are master differentiation inducers. Depending on the differentiation factor and cellular context, pRb can either suppress or promote their transcriptional activity. For example, pRb binds to Runx2 and potentiates its ability to promote osteogenic differentiation in vitro. In contrast, pRb acts with E2F to suppress peroxisome proliferator-activated receptor gamma subunit (PPAR-gamma), the master activator of adipogenesis. Because osteoblasts and adipocytes can both arise from mesenchymal stem cells, these observations suggest that pRb might play a role in the choice between these two fates. However, so far, there is no evidence for this in vivo. Here we use mouse models to address this hypothesis in mesenchymal tissue development and tumorigenesis. Our data show that Rb status plays a key role in establishing fate choice between bone and brown adipose tissue in vivo.

Comparative oncogenomic analysis of copy number alterations in human and zebrafish tumors enables cancer driver discovery.

The identification of cancer drivers is a major goal of current cancer research. Finding driver genes within large chromosomal events is especially challenging because such alterations encompass many genes. Previously, we demonstrated that zebrafish malignant peripheral nerve sheath tumors (MPNSTs) are highly aneuploid, much like human tumors. In this study, we examined 147 zebrafish MPNSTs by massively parallel sequencing and identified both large and focal copy number alterations (CNAs). Given the low degree of conserved synteny between fish and mammals, we reasoned that comparative analyses of CNAs from fish versus human MPNSTs would enable elimination of a large proportion of passenger mutations, especially on large CNAs. We established a list of orthologous genes between human and zebrafish, which includes approximately two-thirds of human protein-coding genes. For the subset of these genes found in human MPNST CNAs, only one quarter of their orthologues were co-gained or co-lost in zebrafish, dramatically narrowing the list of candidate cancer drivers for both focal and large CNAs. We conclude that zebrafish-human comparative analysis represents a powerful, and broadly applicable, tool to enrich for evolutionarily conserved cancer drivers.

p53-deficient cells rely on ATM- and ATR-mediated checkpoint signaling through the p38MAPK/MK2 pathway for survival after DNA damage.

In response to DNA damage, eukaryotic cells activate ATM-Chk2 and/or ATR-Chk1 to arrest the cell cycle and initiate DNA repair. We show that, in the absence of p53, cells depend on a third cell-cycle checkpoint pathway involving p38MAPK/MK2 for cell-cycle arrest and survival after DNA damage. MK2 depletion in p53-deficient cells, but not in p53 wild-type cells, caused abrogation of the Cdc25A-mediated S phase checkpoint after cisplatin exposure and loss of the Cdc25B-mediated G2/M checkpoint following doxorubicin treatment, resulting in mitotic catastrophe and pronounced regression of murine tumors in vivo. We show that the Chk1 inhibitor UCN-01 also potently inhibits MK2, suggesting that its clinical efficacy results from the simultaneous disruption of two critical checkpoint pathways in p53-defective cells.

Life and death decisions by the E2F transcription factors.

The E2F transcription factors are critical regulators of genes required for appropriate progression through the cell cycle, and in special circumstances they can also promote the expression of another class of genes that function in the apoptotic program. Since E2Fs can initiate both cell proliferation and cell death, it is not surprising that the pro-apoptotic capacity of these proteins is subject to complex regulation. Recent study has expanded our knowledge of the factors influencing E2F-induced apoptosis as well as downstream targets of E2F in this process.

Highly aneuploid zebrafish malignant peripheral nerve sheath tumors have genetic alterations similar to human cancers.

Aneuploidy is a hallmark of human cancers, but most mouse cancer models lack the extensive aneuploidy seen in many human tumors. The zebrafish is becoming an increasingly popular model for studying cancer. Here we report that malignant peripheral nerve sheath tumors (MPNSTs) that arise in zebrafish as a result of mutations in either ribosomal protein (rp) genes or in p53 are highly aneuploid. Karyotyping reveals that these tumors frequently harbor near-triploid numbers of chromosomes, and they vary in chromosome number from cell to cell within a single tumor. Using array comparative genomic hybridization, we found that, as in human cancers, certain fish chromosomes are preferentially overrepresented, whereas others are underrepresented in many MPNSTs. In addition, we obtained evidence for recurrent subchromosomal amplifications and deletions that may contain genes involved in cancer initiation or progression. These focal amplifications encompassed several genes whose amplification is observed in human tumors, including met, cyclinD2, slc45a3, and cdk6. One focal amplification included fgf6a. Increasing fgf signaling via a mutation that overexpresses fgf8 accelerated the onset of MPNSTs in fish bearing a mutation in p53, suggesting that fgf6a itself may be a driver of MPNSTs. Our results suggest that the zebrafish is a useful model in which to study aneuploidy in human cancer and in which to identify candidate genes that may act as drivers in fish and potentially also in human tumors.

BMI1 is required for melanocyte stem cell maintenance and hair pigmentation.

The epigenetic repressor BMI1 plays an integral role in promoting the self-renewal and proliferation of many adult stem cell populations, and also tumor types, primarily through silencing the Cdkn2a locus, which encodes the tumor suppressors p16Ink4a and p19Arf . However, in cutaneous melanoma, BMI1 drives epithelial-mesenchymal transition programs, and thus metastasis, while having little impact on proliferation or primary tumor growth. This raised questions about the requirement and role for BMI1 in melanocyte stem cell (McSC) biology. Here, we demonstrate that murine melanocyte-specific Bmi1 deletion causes premature hair greying and gradual loss of melanocyte lineage cells. Depilation enhances this hair greying defect, accelerating depletion of McSCs in early hair cycles, suggesting that BMI1 acts to protect McSCs against stress. RNA-seq of McSCs, harvested before onset of detectable phenotypic defects, revealed that Bmi1 deletion derepresses p16Ink4a and p19Arf , as observed in many other stem cell contexts. Additionally, BMI1 loss downregulated the glutathione S-transferase enzymes, Gsta1 and Gsta2, which can suppress oxidative stress. Accordingly, treatment with the antioxidant N-acetyl cysteine (NAC) partially rescued melanocyte expansion. Together, our data establish a critical function for BMI1 in McSC maintenance that reflects a partial role for suppression of oxidative stress, and likely transcriptional repression of Cdkn2a.

E2F4 loss suppresses tumorigenesis in Rb mutant mice.

The E2F transcription factors mediate the activation or repression of key cell cycle regulatory genes under the control of the retinoblastoma protein (pRB) tumor suppressor and its relatives, p107 and p130. Here we investigate how E2F4, the major "repressive" E2F, contributes to pRB's tumor-suppressive properties. Remarkably, E2F4 loss suppresses the development of both pituitary and thyroid tumors in Rb(+/-) mice. Importantly, E2F4 loss also suppresses the inappropriate gene expression and proliferation of pRB-deficient cells. Biochemical analyses suggest that this tumor suppression occurs via a novel mechanism: E2F4 loss allows p107 and p130 to regulate the pRB-specific, activator E2Fs. We also detect these novel E2F complexes in pRB-deficient cells, suggesting that they play a significant role in the regulation of tumorigenesis in vivo.

Mutation of p107 exacerbates the consequences of Rb loss in embryonic tissues and causes cardiac and blood vessel defects.

The retinoblastoma tumor-suppressor protein, pRb, is a member of the pocket protein family that includes p107 and p130. These proteins have well-defined roles in regulating entry into and exit from the cell cycle and also have cell cycle-independent roles in facilitating differentiation. Here we investigate the overlap between pocket protein's function during embryonic development by using conditional mutant alleles to generate Rb;p107 double-mutant embryos (DKOs) that develop in the absence of placental defects. These DKOs die between e13.5 and e14.5, much earlier than either the conditional Rb or the germline p107 single mutants, which survive to birth or are largely viable, respectively. Analyses of the e13.5 DKOs shows that p107 mutation exacerbates the phenotypes resulting from pRb loss in the central nervous system and lens, but not in the peripheral nervous system. In addition, these embryos exhibit novel phenotypes, including increased proliferation of blood vessel endothelial cells, and heart defects, including double-outlet right ventricle (DORV). The DORV is caused, at least in part, by a defect in blood vessel endothelial cells and/or heart mesenchymal cells. These findings demonstrate novel, overlapping functions for pRb and p107 in numerous murine tissues.