Mammalian CDC14 phosphatases control exit from stemness in pluripotent cells.

Maintenance of stemness is tightly linked to cell cycle regulation through protein phosphorylation by cyclin-dependent kinases (CDKs). However, how this process is reversed during differentiation is unknown. We report here that exit from stemness and differentiation of pluripotent cells along the neural lineage are controlled by CDC14, a CDK-counteracting phosphatase whose function in mammals remains obscure. Lack of the two CDC14 family members, CDC14A and CDC14B, results in deficient development of the neural system in the mouse and impairs neural differentiation from embryonic stem cells (ESCs). Mechanistically, CDC14 directly dephosphorylates specific proline-directed Ser/Thr residues of undifferentiated embryonic transcription Factor 1 (UTF1) during the exit from stemness, triggering its proteasome-dependent degradation. Multiomic single-cell analysis of transcription and chromatin accessibility in differentiating ESCs suggests that increased UTF1 levels in the absence of CDC14 prevent the proper firing of bivalent promoters required for differentiation. CDC14 phosphatases are dispensable for mitotic exit, suggesting that CDC14 phosphatases have evolved to control stemness rather than cell cycle exit and establish the CDK-CDC14 axis as a critical molecular switch for linking cell cycle regulation and self-renewal.

Transient exposure to miR-203 enhances the differentiation capacity of established pluripotent stem cells.

Full differentiation potential along with self-renewal capacity is a major property of pluripotent stem cells (PSCs). However, the differentiation capacity frequently decreases during expansion of PSCs in vitro. We show here that transient exposure to a single microRNA, expressed at early stages during normal development, improves the differentiation capacity of already-established murine and human PSCs. Short exposure to miR-203 in PSCs (miPSCs) induces a transient expression of 2C markers that later results in expanded differentiation potency to multiple lineages, as well as improved efficiency in tetraploid complementation and human-mouse interspecies chimerism assays. Mechanistically, these effects are at least partially mediated by direct repression of de novo DNA methyltransferases Dnmt3a and Dnmt3b, leading to transient and reversible erasure of DNA methylation. These data support the use of transient exposure to miR-203 as a versatile method to reset the epigenetic memory in PSCs, and improve their effectiveness in regenerative medicine.

Pseudoalignment tools as an efficient alternative to detect repeated transposable elements in scRNAseq data.

MOTIVATION: Transposable elements (TE) have played a major role in configuring the structures of mammalian genomes through evolution. In normal conditions, the expression of these elements is repressed by different epigenetic regulation mechanisms such as DNA methylation, histone modification and regulation by small RNAs. TE re-activation is associated with stemness potential acquisition, regulation of innate immunity and disease, such as cancer. However, the vast majority of current knowledge in the field is based on bulk expression studies, and very little is known on cell-type- or state-specific expression of TE-derived transcripts. Therefore, cost-efficient single-cell-resolution TE expression analytical approaches are needed. RESULTS: We have implemented an analytical approach based on pseudoalignment to consensus sequences to incorporate TE expression information to scRNAseq data. AVAILABILITY AND IMPLEMENTATION: All the data and code implemented are available as Supplementary data and in: https://github.com/jmzvillarreal/kallisto_TE_scRNAseq. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

3D chromatin connectivity underlies replication origin efficiency in mouse embryonic stem cells.

In mammalian cells, chromosomal replication starts at thousands of origins at which replisomes are assembled. Replicative stress triggers additional initiation events from 'dormant' origins whose genomic distribution and regulation are not well understood. In this study, we have analyzed origin activity in mouse embryonic stem cells in the absence or presence of mild replicative stress induced by aphidicolin, a DNA polymerase inhibitor, or by deregulation of origin licensing factor CDC6. In both cases, we observe that the majority of stress-responsive origins are also active in a small fraction of the cell population in a normal S phase, and stress increases their frequency of activation. In a search for the molecular determinants of origin efficiency, we compared the genetic and epigenetic features of origins displaying different levels of activation, and integrated their genomic positions in three-dimensional chromatin interaction networks derived from high-depth Hi-C and promoter-capture Hi-C data. We report that origin efficiency is directly proportional to the proximity to transcriptional start sites and to the number of contacts established between origin-containing chromatin fragments, supporting the organization of origins in higher-level DNA replication factories.

miR-203 drives breast cancer cell differentiation.

A hallmark of many malignant tumors is dedifferentiated (immature) cells bearing slight or no resemblance to the normal cells from which the cancer originated. Tumor dedifferentiated cells exhibit a higher capacity to survive to chemo and radiotherapies and have the ability to incite tumor relapse. Inducing cancer cell differentiation would abolish their self-renewal and invasive capacity and could be combined with the current standard of care, especially in poorly differentiated and aggressive tumors (with worst prognosis). However, differentiation therapy is still in its early stages and the intrinsic complexity of solid tumor heterogeneity demands innovative approaches in order to be efficiently translated into the clinic. We demonstrate here that microRNA 203, a potent driver of differentiation in pluripotent stem cells (ESCs and iPSCs), promotes the differentiation of mammary gland tumor cells. Combining mouse in vivo approaches and both mouse and human-derived tridimensional organoid cultures, we report that miR-203 influences the self-renewal capacity, plasticity and differentiation potential of breast cancer cells and prevents tumor cell growth in vivo. Our work sheds light on differentiation-based antitumor therapies and offers miR-203 as a promising tool for directly confronting the tumor-maintaining and regeneration capability of cancer cells.

Lymphatic vessels interact dynamically with the hair follicle stem cell niche during skin regeneration in vivo.

Lymphatic vessels are essential for skin fluid homeostasis and immune cell trafficking. Whether the lymphatic vasculature is associated with hair follicle regeneration is, however, unknown. Here, using steady and live imaging approaches in mouse skin, we show that lymphatic vessels distribute to the anterior permanent region of individual hair follicles, starting from development through all cycle stages and interconnecting neighboring follicles at the bulge level, in a stem cell-dependent manner. Lymphatic vessels further connect hair follicles in triads and dynamically flow across the skin. At the onset of the physiological stem cell activation, or upon pharmacological or genetic induction of hair follicle growth, lymphatic vessels transiently expand their caliber suggesting an increased tissue drainage capacity. Interestingly, the physiological caliber increase is associated with a distinct gene expression correlated with lymphatic vessel reorganization. Using mouse genetics, we show that lymphatic vessel depletion blocks hair follicle growth. Our findings point toward the lymphatic vasculature being important for hair follicle development, cycling, and organization, and define lymphatic vessels as stem cell niche components, coordinating connections at tissue-level, thus provide insight into their functional contribution to skin regeneration.

Chromatin regulation by Histone H4 acetylation at Lysine 16 during cell death and differentiation in the myeloid compartment.

Histone H4 acetylation at Lysine 16 (H4K16ac) is a key epigenetic mark involved in gene regulation, DNA repair and chromatin remodeling, and though it is known to be essential for embryonic development, its role during adult life is still poorly understood. Here we show that this lysine is massively hyperacetylated in peripheral neutrophils. Genome-wide mapping of H4K16ac in terminally differentiated blood cells, along with functional experiments, supported a role for this histone post-translational modification in the regulation of cell differentiation and apoptosis in the hematopoietic system. Furthermore, in neutrophils, H4K16ac was enriched at specific DNA repeats. These DNA regions presented an accessible chromatin conformation and were associated with the cleavage sites that generate the 50 kb DNA fragments during the first stages of programmed cell death. Our results thus suggest that H4K16ac plays a dual role in myeloid cells as it not only regulates differentiation and apoptosis, but it also exhibits a non-canonical structural role in poising chromatin for cleavage at an early stage of neutrophil cell death.

Hsa-miR-139-5p is a prognostic thyroid cancer marker involved in HNRNPF-mediated alternative splicing.

It is critical to identify biomarkers and functional networks associated with aggressive thyroid cancer to anticipate disease progression and facilitate personalized patient management. We performed miRNome sequencing of 46 thyroid tumors enriched with advanced disease patients with a median follow-up of 96 months. MiRNome profiles correlated with tumor-specific histopathological and molecular features, such as stromal cell infiltration and tumor driver mutation. Differential expression analysis revealed a consistent hsa-miR-139-5p downexpression in primary carcinomas from patients with recurrent/metastatic disease compared to disease-free patients, sustained in paired local metastases and validated in publicly available thyroid cancer series. Exogenous expression of hsa-miR-139-5p significantly reduced migration and proliferation of anaplastic thyroid cancer cells. Proteomic analysis indicated RICTOR, SMAD2/3 and HNRNPF as putative hsa-miR-139-5p targets in our cell system. Abundance of HNRNPF mRNA, encoding an alternative splicing factor involved in cryptic exon inclusion/exclusion, inversely correlated with hsa-miR-139-5p expression in human tumors. RNA sequencing analysis revealed 174 splicing events differentially regulated upon HNRNPF repression in our cell system, affecting genes involved in RTK/RAS/MAPK and PI3K/AKT/MTOR signaling cascades among others. These results point at the hsa-miR-139-5p/HNRNPF axis as a novel regulatory mechanism associated with the modulation of major thyroid cancer signaling pathways and tumor virulence.

RNA-Seq reveals the existence of a CDKN1C-E2F1-TP53 axis that is altered in human T-cell lymphoblastic lymphomas.

BACKGROUND: Precursor T-cell lymphoblastic lymphomas (T-LBL) are rare aggressive hematological malignancies that mainly develop in children. As in other cancers, the loss of cell cycle control plays a prominent role in the pathogenesis in these malignancies that is primarily attributed to loss of CDKN2A (encoding protein p16INK4A). However, the impact of the deregulation of other genes such as CDKN1C, E2F1, and TP53 remains to be clarified. Interestingly, experiments in mouse models have proven that conditional T-cell specific deletion of Cdkn1c gene may induce a differentiation block at the DN3 to DN4 transition, and that the loss of this gene in the absence of Tp53 led to aggressive thymic lymphomas. RESULTS: In this manuscript, we demonstrated that the simultaneous deregulation of CDKN1C, E2F1, and TP53 genes by epigenetic mechanisms and/or the deregulation of specific microRNAs, together with additional impairing of TP53 function by the expression of dominant-negative isoforms are common features in primary human T-LBLs. CONCLUSIONS: Previous experimental work in mice revealed that T-cell specific deletion of Cdkn1c accelerates lymphomagenesis in the absence of Tp53. If, as expected, the consequences of the deregulation of the CDKN1C-E2F1-TP53 axis were the same as those experimentally demonstrated in mouse models, the disruption of this axis might be useful to predict tumor aggressiveness, and to provide the basis towards the development of potential therapeutic strategiesin human T-LBL.

Renal fibroblasts are involved in fibrogenic changes in kidney fibrosis associated with dysfunctional telomeres.

Tubulointerstitial fibrosis associated with chronic kidney disease (CKD) represents a global health care problem. We previously reported that short and dysfunctional telomeres lead to interstitial renal fibrosis; however, the cell-of-origin of kidney fibrosis associated with telomere dysfunction is currently unknown. We induced telomere dysfunction by deleting the Trf1 gene encoding a telomere-binding factor specifically in renal fibroblasts in both short-term and long-term life-long experiments in mice to identify the role of fibroblasts in renal fibrosis. Short-term Trf1 deletion in renal fibroblasts was not sufficient to trigger kidney fibrosis but was sufficient to induce inflammatory responses, ECM deposition, cell cycle arrest, fibrogenesis, and vascular rarefaction. However, long-term persistent deletion of Trf1 in fibroblasts resulted in kidney fibrosis accompanied by an elevated urinary albumin-to-creatinine ratio (uACR) and a decrease in mouse survival. These cellular responses lead to the macrophage-to-myofibroblast transition (MMT), endothelial-to-mesenchymal transition (EndMT), and partial epithelial-to-mesenchymal transition (EMT), ultimately causing kidney fibrosis at the humane endpoint (HEP) when the deletion of Trf1 in fibroblasts is maintained throughout the lifespan of mice. Our findings contribute to a better understanding of the role of dysfunctional telomeres in the onset of the profibrotic alterations that lead to kidney fibrosis.

Hepatic nutrient and hormone signaling to mTORC1 instructs the postnatal metabolic zonation of the liver.

The metabolic functions of the liver are spatially organized in a phenomenon called zonation, linked to the differential exposure of portal and central hepatocytes to nutrient-rich blood. The mTORC1 signaling pathway controls cellular metabolism in response to nutrients and insulin fluctuations. Here we show that simultaneous genetic activation of nutrient and hormone signaling to mTORC1 in hepatocytes results in impaired establishment of postnatal metabolic and zonal identity of hepatocytes. Mutant hepatocytes fail to upregulate postnatally the expression of Frizzled receptors 1 and 8, and show reduced Wnt/ß-catenin activation. This defect, alongside diminished paracrine Wnt2 ligand expression by endothelial cells, underlies impaired postnatal maturation. Impaired zonation is recapitulated in a model of constant supply of nutrients by parenteral nutrition to piglets. Our work shows the role of hepatocyte sensing of fluctuations in nutrients and hormones for triggering a latent metabolic zonation program.

A mild increase in nutrient signaling to mTORC1 in mice leads to parenchymal damage, myeloid inflammation and shortened lifespan.

The mechanistic target of rapamycin complex 1 controls cellular anabolism in response to growth factor signaling and to nutrient sufficiency signaled through the Rag GTPases. Inhibition of mTOR reproducibly extends longevity across eukaryotes. Here we report that mice that endogenously express active mutant variants of RagC exhibit multiple features of parenchymal damage that include senescence, expression of inflammatory molecules, increased myeloid inflammation with extensive features of inflammaging and a ~30% reduction in lifespan. Through bone marrow transplantation experiments, we show that myeloid cells are abnormally activated by signals emanating from dysfunctional RagC-mutant parenchyma, causing neutrophil extravasation that inflicts additional inflammatory damage. Therapeutic suppression of myeloid inflammation in aged RagC-mutant mice attenuates parenchymal damage and extends survival. Together, our findings link mildly increased nutrient signaling to limited lifespan in mammals, and support a two-component process of parenchymal damage and myeloid inflammation that together precipitate a time-dependent organ deterioration that limits longevity.

TIMP1 mediates astrocyte-dependent local immunosuppression in brain metastasis acting on infiltrating CD8+ T cells.

Immunotherapies against brain metastases have shown clinical benefits when applied to asymptomatic patients, but they are largely ineffective in symptomatic cases for unknown reasons. Here we dissect the heterogeneity in metastasis-associated astrocytes using scRNAseq and report a population that blocks the antitumoral activity of infiltrating T cells. This pro-tumoral activity is mediated by the secretion of TIMP1 from a cluster of pSTAT3+ astrocytes that acts on CD63+ CD8+ T cells to modulate their function. Using genetic and pharmacologic approaches in mouse and human brain metastasis models, we demonstrate that combining immune checkpoint blockade antibodies with the inhibition of astrocyte-mediated local immunosuppression may benefit patients with symptomatic brain metastases. We further reveal that the presence of TIMP1 in liquid biopsies provides a biomarker to select patients for this combined immunotherapy. Overall, our findings demonstrate an unexpected immunomodulatory role for astrocytes in brain metastases with clinical implications.

Luminal Rank loss decreases cell fitness leading to basal cell bipotency in parous mammary glands.

Rank signaling pathway regulates mammary gland homeostasis and epithelial cell differentiation. Although Rank receptor is expressed by basal cells and luminal progenitors, its role in each individual cell lineage remains unclear. By combining temporal/lineage specific Rank genetic deletion with lineage tracing techniques, we found that loss of luminal Rank reduces the luminal progenitor pool and leads to aberrant alveolar-like differentiation with high protein translation capacity in virgin mammary glands. These Rank-deleted luminal cells are unable to expand during the first pregnancy, leading to lactation failure and impairment of protein synthesis potential in the parous stage. The unfit parous Rank-deleted luminal cells in the alveoli are progressively replaced by Rank-proficient cells early during the second pregnancy, thereby restoring lactation. Transcriptomic analysis and functional assays point to the awakening of basal bipotency after pregnancy by the induction of Rank/NF-κB signaling in basal parous cell to restore lactation and tissue homeostasis.

Endoglin, a Novel Biomarker and Therapeutical Target to Prevent Malignant Peripheral Nerve Sheath Tumor Growth and Metastasis.

PURPOSE: Malignant peripheral nerve sheath tumors (MPNST) are highly aggressive soft-tissue sarcomas that lack effective treatments, underscoring the urgent need to uncover novel mediators of MPNST pathogenesis that may serve as potential therapeutic targets. Tumor angiogenesis is considered a critical event in MPNST transformation and progression. Here, we have investigated whether endoglin (ENG), a TGFß coreceptor with a crucial role in angiogenesis, could be a novel therapeutic target in MPNSTs. EXPERIMENTAL DESIGN: ENG expression was evaluated in human peripheral nerve sheath tumor tissues and plasma samples. Effects of tumor cell-specific ENG expression on gene expression, signaling pathway activation and in vivo MPNST growth and metastasis, were investigated. The efficacy of ENG targeting in monotherapy or in combination with MEK inhibition was analyzed in xenograft models. RESULTS: ENG expression was found to be upregulated in both human MPNST tumor tissues and plasma-circulating small extracellular vesicles. We demonstrated that ENG modulates Smad1/5 and MAPK/ERK pathway activation and pro-angiogenic and pro-metastatic gene expression in MPNST cells and plays an active role in tumor growth and metastasis in vivo. Targeting with ENG-neutralizing antibodies (TRC105/M1043) decreased MPNST growth and metastasis in xenograft models by reducing tumor cell proliferation and angiogenesis. Moreover, combination of anti-ENG therapy with MEK inhibition effectively reduced tumor cell growth and angiogenesis. CONCLUSIONS: Our data unveil a tumor-promoting function of ENG in MPNSTs and support the use of this protein as a novel biomarker and a promising therapeutic target for this disease.

Folliculin-interacting protein FNIP2 impacts on overweight and obesity through a polymorphism in a conserved 3' untranslated region.

BACKGROUND: Overweight and obesity are defined by an anomalous or excessive fat accumulation that may compromise health. To find single-nucleotide polymorphisms (SNPs) influencing metabolic phenotypes associated with the obesity state, we analyze multiple anthropometric and clinical parameters in a cohort of 790 healthy volunteers and study potential associations with 48 manually curated SNPs, in metabolic genes functionally associated with the mechanistic target of rapamycin (mTOR) pathway. RESULTS: We identify and validate rs2291007 within a conserved region in the 3'UTR of folliculin-interacting protein FNIP2 that correlates with multiple leanness parameters. The T-to-C variant represents the major allele in Europeans and disrupts an ancestral target sequence of the miRNA miR-181b-5p, thus resulting in increased FNIP2 mRNA levels in cancer cell lines and in peripheral blood from carriers of the C allele. Because the miRNA binding site is conserved across vertebrates, we engineered the T-to-C substitution in the endogenous Fnip2 allele in mice. Primary cells derived from Fnip2 C/C mice show increased mRNA stability, and more importantly, Fnip2 C/C mice replicate the decreased adiposity and increased leanness observed in human volunteers. Finally, expression levels of FNIP2 in both human samples and mice negatively associate with leanness parameters, and moreover, are the most important contributor in a multifactorial model of body mass index prediction. CONCLUSIONS: We propose that rs2291007 influences human leanness through an evolutionarily conserved modulation of FNIP2 mRNA levels.

A clinically compatible drug-screening platform based on organotypic cultures identifies vulnerabilities to prevent and treat brain metastasis.

We report a medium-throughput drug-screening platform (METPlatform) based on organotypic cultures that allows to evaluate inhibitors against metastases growing in situ. By applying this approach to the unmet clinical need of brain metastasis, we identified several vulnerabilities. Among them, a blood-brain barrier permeable HSP90 inhibitor showed high potency against mouse and human brain metastases at clinically relevant stages of the disease, including a novel model of local relapse after neurosurgery. Furthermore, in situ proteomic analysis applied to metastases treated with the chaperone inhibitor uncovered a novel molecular program in brain metastasis, which includes biomarkers of poor prognosis and actionable mechanisms of resistance. Our work validates METPlatform as a potent resource for metastasis research integrating drug-screening and unbiased omic approaches that is compatible with human samples. Thus, this clinically relevant strategy is aimed to personalize the management of metastatic disease in the brain and elsewhere.

Stratification of radiosensitive brain metastases based on an actionable S100A9/RAGE resistance mechanism.

Whole-brain radiotherapy (WBRT) is the treatment backbone for many patients with brain metastasis; however, its efficacy in preventing disease progression and the associated toxicity have questioned the clinical impact of this approach and emphasized the need for alternative treatments. Given the limited therapeutic options available for these patients and the poor understanding of the molecular mechanisms underlying the resistance of metastatic lesions to WBRT, we sought to uncover actionable targets and biomarkers that could help to refine patient selection. Through an unbiased analysis of experimental in vivo models of brain metastasis resistant to WBRT, we identified activation of the S100A9-RAGE-NF-κB-JunB pathway in brain metastases as a potential mediator of resistance in this organ. Targeting this pathway genetically or pharmacologically was sufficient to revert the WBRT resistance and increase therapeutic benefits in vivo at lower doses of radiation. In patients with primary melanoma, lung or breast adenocarcinoma developing brain metastasis, endogenous S100A9 levels in brain lesions correlated with clinical response to WBRT and underscored the potential of S100A9 levels in the blood as a noninvasive biomarker. Collectively, we provide a molecular framework to personalize WBRT and improve its efficacy through combination with a radiosensitizer that balances therapeutic benefit and toxicity.

Short and dysfunctional telomeres sensitize the kidneys to develop fibrosis.

Accumulation of short telomeres is a hallmark of aging. Mutations in telomerase or telomere-binding proteins lead to telomere shortening or dysfunction and are at the origin of human pathologies known as 'telomere syndromes', which are characterized by loss of the regenerative capacity of tissues and fibrotic pathologies. Here, we generated two mouse models of kidney fibrosis, either by combining telomerase deficiency to induce telomere shortening and a low dose of folic acid, or by conditionally deleting Trf1, a component of the shelterin telomere protective complex, from the kidneys. We find that short telomeres sensitize the kidneys to develop fibrosis in response to folic acid and exacerbate the epithelial-to-mesenchymal transition (EMT) program. Trf1 deletion in kidneys led to fibrosis and EMT activation. Our findings suggest that telomere shortening or dysfunction may contribute to pathological, age-associated renal fibrosis by influencing the EMT program.

Compi: a framework for portable and reproducible pipelines.

Compi is an application framework to develop end-user, pipeline-based applications with a primary emphasis on: (i) user interface generation, by automatically generating a command-line interface based on the pipeline specific parameter definitions; (ii) application packaging, with compi-dk, which is a version-control-friendly tool to package the pipeline application and its dependencies into a Docker image; and (iii) application distribution provided through a public repository of Compi pipelines, named Compi Hub, which allows users to discover, browse and reuse them easily. By addressing these three aspects, Compi goes beyond traditional workflow engines, having been specially designed for researchers who want to take advantage of common workflow engine features (such as automatic job scheduling or logging, among others) while keeping the simplicity and readability of shell scripts without the need to learn a new programming language. Here we discuss the design of various pipelines developed with Compi to describe its main functionalities, as well as to highlight the similarities and differences with similar tools that are available. An open-source distribution under the Apache 2.0 License is available from GitHub (available at https://github.com/sing-group/compi). Documentation and installers are available from https://www.sing-group.org/compi. A specific repository for Compi pipelines is available from Compi Hub (available at https://www.sing-group.org/compihub.