Identity Noise and Adipogenic Traits Characterize Dermal Fibroblast Aging.

During aging, stromal functions are thought to be impaired, but little is known whether this stems from changes of fibroblasts. Using population- and single-cell transcriptomics, as well as long-term lineage tracing, we studied whether murine dermal fibroblasts are altered during physiological aging under different dietary regimes that affect longevity. We show that the identity of old fibroblasts becomes undefined, with the fibroblast states present in young skin no longer clearly demarcated. In addition, old fibroblasts not only reduce the expression of genes involved in the formation of the extracellular matrix, but also gain adipogenic traits, paradoxically becoming more similar to neonatal pro-adipogenic fibroblasts. These alterations are sensitive to systemic metabolic changes: long-term caloric restriction reversibly prevents them, whereas a high-fat diet potentiates them. Our results therefore highlight loss of cell identity and the acquisition of adipogenic traits as a mechanism underlying cellular aging, which is influenced by systemic metabolism.

Roastgsa: a comparison of rotation-based scores for gene set enrichment analysis.

BACKGROUND: Gene-wise differential expression is usually the first major step in the statistical analysis of high-throughput data obtained from techniques such as microarrays or RNA-sequencing. The analysis at gene level is often complemented by interrogating the data in a broader biological context that considers as unit of measure groups of genes that may have a common function or biological trait. Among the vast number of publications about gene set analysis (GSA), the rotation test for gene set analysis, also referred to as roast, is a general sample randomization approach that maintains the integrity of the intra-gene set correlation structure in defining the null distribution of the test. RESULTS: We present roastgsa, an R package that contains several enrichment score functions that feed the roast algorithm for hypothesis testing. These implemented methods are evaluated using both simulated and benchmarking data in microarray and RNA-seq datasets. We find that computationally intensive measures based on Kolmogorov-Smirnov (KS) statistics fail to improve the rates of simpler measures of GSA like mean and maxmean scores. We also show the importance of accounting for the gene linear dependence structure of the testing set, which is linked to the loss of effective signature size. Complete graphical representation of the results, including an approximation for the effective signature size, can be obtained as part of the roastgsa output. CONCLUSIONS: We encourage the usage of the absmean (non-directional), mean (directional) and maxmean (directional) scores for roast GSA analysis as these are simple measures of enrichment that have presented dominant results in all provided analyses in comparison to the more complex KS measures.

Differential expression of human tRNA genes drives the abundance of tRNA-derived fragments.

The human genome encodes hundreds of transfer RNA (tRNA) genes but their individual contribution to the tRNA pool is not fully understood. Deep sequencing of tRNA transcripts (tRNA-Seq) can estimate tRNA abundance at single gene resolution, but tRNA structures and posttranscriptional modifications impair these analyses. Here we present a bioinformatics strategy to investigate differential tRNA gene expression and use it to compare tRNA-Seq datasets from cultured human cells and human brain. We find that sequencing caveats affect quantitation of only a subset of human tRNA genes. Unexpectedly, we detect several cases where the differences in tRNA expression among samples do not involve variations at the level of isoacceptor tRNA sets (tRNAs charged with the same amino acid but using different anticodons), but rather among tRNA genes within the same isodecoder set (tRNAs having the same anticodon sequence). Because isodecoder tRNAs are functionally equal in terms of genetic translation, their differential expression may be related to noncanonical tRNA functions. We show that several instances of differential tRNA gene expression result in changes in the abundance of tRNA-derived fragments (tRFs) but not of mature tRNAs. Examples of differentially expressed tRFs include PIWI-associated RNAs, tRFs present in tissue samples but not in cells cultured in vitro, and somatic tissue-specific tRFs. Our data support that differential expression of tRNA genes regulate noncanonical tRNA functions performed by tRFs.

Designing alternative splicing RNA-seq studies. Beyond generic guidelines.

MOTIVATION: Designing an RNA-seq study depends critically on its specific goals, technology and underlying biology, which renders general guidelines inadequate. We propose a Bayesian framework to customize experiments so that goals can be attained and resources are not wasted, with a focus on alternative splicing. RESULTS: We studied how read length, sequencing depth, library preparation and the number of replicates affects cost-effectiveness of single-sample and group comparison studies. Optimal settings varied strongly according to the target organism or tissue (potential 50-500% cost cuts) and, interestingly, short reads outperformed long reads for standard analyses. Our framework learns key characteristics for study design from the data, and predicts if and how to continue experimentation. These predictions matched several follow-up experimental datasets that were used for validation. We provide default pipelines, but the framework can be combined with other data analysis methods and can help assess their relative merits. AVAILABILITY AND IMPLEMENTATION: casper package at www.bioconductor.org/packages/release/bioc/html/casper.html, Supplementary Manual by typing casperDesign() at the R prompt. CONTACT: rosselldavid@gmail.com SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Iro/IRX transcription factors negatively regulate Dpp/TGF-ß pathway activity during intestinal tumorigenesis.

Activating mutations in Wnt and EGFR/Ras signaling pathways are common in colorectal cancer (CRC). Remarkably, clonal co-activation of these pathways in the adult Drosophila midgut induces "tumor-like" overgrowths. Here, we show that, in these clones and in CRC cell lines, Dpp/TGF-ß acts as a tumor suppressor. Moreover, we discover that the Iroquois/IRX-family-protein Mirror downregulates the transcription of core components of the Dpp pathway, reducing its tumor suppressor activity. We also show that this genetic interaction is conserved in human CRC cells, where the Iro/IRX proteins IRX3 and IRX5 diminish the response to TGF-ß. IRX3 and IRX5 are upregulated in human adenomas, and their levels correlate inversely with the gene expression signature of response to TGF-ß. In addition, Irx5 expression confers a growth advantage in the presence of TGF-ß, but is selected against in its absence. Together, our results identify a set of Iro/IRX proteins as conserved negative regulators of Dpp/TGF-ß activity. We propose that during the characteristic adenoma-to-carcinoma transition of human CRC, the activity of IRX proteins could reduce the sensitivity to the cytostatic effect of TGF-ß, conferring a growth advantage to tumor cells prior to the acquisition of mutations in TGF-ß pathway components.

Structure and mechanisms of transport of human Asc1/CD98hc amino acid transporter.

Recent cryoEM studies elucidated details of the structural basis for the substrate selectivity and translocation of heteromeric amino acid transporters. However, Asc1/CD98hc is the only neutral heteromeric amino acid transporter that can function through facilitated diffusion, and the only one that efficiently transports glycine and D-serine, and thus has a regulatory role in the central nervous system. Here we use cryoEM, ligand-binding simulations, mutagenesis, transport assays, and molecular dynamics to define human Asc1/CD98hc determinants for substrate specificity and gain insights into the mechanisms that govern substrate translocation by exchange and facilitated diffusion. The cryoEM structure of Asc1/CD98hc is determined at 3.4-3.8 Å resolution, revealing an inward-facing semi-occluded conformation. We find that Ser 246 and Tyr 333 are essential for Asc1/CD98hc substrate selectivity and for the exchange and facilitated diffusion modes of transport. Taken together, these results reveal the structural bases for ligand binding and transport features specific to human Asc1.

Membrane to cortex attachment determines different mechanical phenotypes in LGR5+ and LGR5- colorectal cancer cells.

Colorectal cancer (CRC) tumors are composed of heterogeneous and plastic cell populations, including a pool of cancer stem cells that express LGR5. Whether these distinct cell populations display different mechanical properties, and how these properties might contribute to metastasis is poorly understood. Using CRC patient derived organoids (PDOs), we find that compared to LGR5- cells, LGR5+ cancer stem cells are stiffer, adhere better to the extracellular matrix (ECM), move slower both as single cells and clusters, display higher nuclear YAP, show a higher survival rate in response to mechanical confinement, and form larger transendothelial gaps. These differences are largely explained by the downregulation of the membrane to cortex attachment proteins Ezrin/Radixin/Moesin (ERMs) in the LGR5+ cells. By analyzing single cell RNA-sequencing (scRNA-seq) expression patterns from a patient cohort, we show that this downregulation is a robust signature of colorectal tumors. Our results show that LGR5- cells display a mechanically dynamic phenotype suitable for dissemination from the primary tumor whereas LGR5+ cells display a mechanically stable and resilient phenotype suitable for extravasation and metastatic growth.

The efficacy of chemotherapy is limited by intratumoral senescent cells expressing PD-L2.

Chemotherapy often generates intratumoral senescent cancer cells that strongly modify the tumor microenvironment, favoring immunosuppression and tumor growth. We discovered, through an unbiased proteomics screen, that the immune checkpoint inhibitor programmed cell death 1 ligand 2 (PD-L2) is highly upregulated upon induction of senescence in different types of cancer cells. PD-L2 is not required for cells to undergo senescence, but it is critical for senescent cells to evade the immune system and persist intratumorally. Indeed, after chemotherapy, PD-L2-deficient senescent cancer cells are rapidly eliminated and tumors do not produce the senescence-associated chemokines CXCL1 and CXCL2. Accordingly, PD-L2-deficient pancreatic tumors fail to recruit myeloid-derived suppressor cells and undergo regression driven by CD8 T cells after chemotherapy. Finally, antibody-mediated blockade of PD-L2 strongly synergizes with chemotherapy causing remission of mammary tumors in mice. The combination of chemotherapy with anti-PD-L2 provides a therapeutic strategy that exploits vulnerabilities arising from therapy-induced senescence.

Vitamin B12 is a limiting factor for induced cellular plasticity and tissue repair.

Transient reprogramming by the expression of OCT4, SOX2, KLF4 and MYC (OSKM) is a therapeutic strategy for tissue regeneration and rejuvenation, but little is known about its metabolic requirements. Here we show that OSKM reprogramming in mice causes a global depletion of vitamin B12 and molecular hallmarks of methionine starvation. Supplementation with vitamin B12 increases the efficiency of reprogramming both in mice and in cultured cells, the latter indicating a cell-intrinsic effect. We show that the epigenetic mark H3K36me3, which prevents illegitimate initiation of transcription outside promoters (cryptic transcription), is sensitive to vitamin B12 levels, providing evidence for a link between B12 levels, H3K36 methylation, transcriptional fidelity and efficient reprogramming. Vitamin B12 supplementation also accelerates tissue repair in a model of ulcerative colitis. We conclude that vitamin B12, through its key role in one-carbon metabolism and epigenetic dynamics, improves the efficiency of in vivo reprogramming and tissue repair.

Iron accumulation drives fibrosis, senescence and the senescence-associated secretory phenotype.

Fibrogenesis is part of a normal protective response to tissue injury that can become irreversible and progressive, leading to fatal diseases. Senescent cells are a main driver of fibrotic diseases through their secretome, known as senescence-associated secretory phenotype (SASP). Here, we report that cellular senescence, and multiple types of fibrotic diseases in mice and humans are characterized by the accumulation of iron. We show that vascular and hemolytic injuries are efficient in triggering iron accumulation, which in turn can cause senescence and promote fibrosis. Notably, we find that senescent cells persistently accumulate iron, even when the surge of extracellular iron has subdued. Indeed, under normal conditions of extracellular iron, cells exposed to different types of senescence-inducing insults accumulate abundant ferritin-bound iron, mostly within lysosomes, and present high levels of labile iron, which fuels the generation of reactive oxygen species and the SASP. Finally, we demonstrate that detection of iron by magnetic resonance imaging might allow non-invasive assessment of fibrotic burden in the kidneys of mice and in patients with renal fibrosis. Our findings suggest that iron accumulation plays a central role in senescence and fibrosis, even when the initiating events may be independent of iron, and identify iron metabolism as a potential therapeutic target for senescence-associated diseases.

Multi-omic rejuvenation of naturally aged tissues by a single cycle of transient reprogramming.

The expression of the pluripotency factors OCT4, SOX2, KLF4, and MYC (OSKM) can convert somatic differentiated cells into pluripotent stem cells in a process known as reprogramming. Notably, partial and reversible reprogramming does not change cell identity but can reverse markers of aging in cells, improve the capacity of aged mice to repair tissue injuries, and extend longevity in progeroid mice. However, little is known about the mechanisms involved. Here, we have studied changes in the DNA methylome, transcriptome, and metabolome in naturally aged mice subject to a single period of transient OSKM expression. We found that this is sufficient to reverse DNA methylation changes that occur upon aging in the pancreas, liver, spleen, and blood. Similarly, we observed reversion of transcriptional changes, especially regarding biological processes known to change during aging. Finally, some serum metabolites and biomarkers altered with aging were also restored to young levels upon transient reprogramming. These observations indicate that a single period of OSKM expression can drive epigenetic, transcriptomic, and metabolomic changes toward a younger configuration in multiple tissues and in the serum.

Mex3a marks drug-tolerant persister colorectal cancer cells that mediate relapse after chemotherapy.

Colorectal cancer (CRC) patient-derived organoids predict responses to chemotherapy. Here we used them to investigate relapse after treatment. Patient-derived organoids expand from highly proliferative LGR5+ tumor cells; however, we discovered that lack of optimal growth conditions specifies a latent LGR5+ cell state. This cell population expressed the gene MEX3A, is chemoresistant and regenerated the organoid culture after treatment. In CRC mouse models, Mex3a+ cells contributed marginally to metastatic outgrowth; however, after chemotherapy, Mex3a+ cells produced large cell clones that regenerated the disease. Lineage-tracing analysis showed that persister Mex3a+ cells downregulate the WNT/stem cell gene program immediately after chemotherapy and adopt a transient state reminiscent to that of YAP+ fetal intestinal progenitors. In contrast, Mex3a-deficient cells differentiated toward a goblet cell-like phenotype and were unable to resist chemotherapy. Our findings reveal that adaptation of cancer stem cells to suboptimal niche environments protects them from chemotherapy and identify a candidate cell of origin of relapse after treatment in CRC.

pTINCR microprotein promotes epithelial differentiation and suppresses tumor growth through CDC42 SUMOylation and activation.

The human transcriptome contains thousands of small open reading frames (sORFs) that encode microproteins whose functions remain largely unexplored. Here, we show that TINCR lncRNA encodes pTINCR, an evolutionary conserved ubiquitin-like protein (UBL) expressed in many epithelia and upregulated upon differentiation and under cellular stress. By gain- and loss-of-function studies, we demonstrate that pTINCR is a key inducer of epithelial differentiation in vitro and in vivo. Interestingly, low expression of TINCR associates with worse prognosis in several epithelial cancers, and pTINCR overexpression reduces malignancy in patient-derived xenografts. At the molecular level, pTINCR binds to SUMO through its SUMO interacting motif (SIM) and to CDC42, a Rho-GTPase critical for actin cytoskeleton remodeling and epithelial differentiation. Moreover, pTINCR increases CDC42 SUMOylation and promotes its activation, triggering a pro-differentiation cascade. Our findings suggest that the microproteome is a source of new regulators of cell identity relevant for cancer.

Functional patient-derived organoid screenings identify MCLA-158 as a therapeutic EGFR × LGR5 bispecific antibody with efficacy in epithelial tumors.

Patient-derived organoids (PDOs) recapitulate tumor architecture, contain cancer stem cells and have predictive value supporting personalized medicine. Here we describe a large-scale functional screen of dual-targeting bispecific antibodies (bAbs) on a heterogeneous colorectal cancer PDO biobank and paired healthy colonic mucosa samples. More than 500 therapeutic bAbs generated against Wingless-related integration site (WNT) and receptor tyrosine kinase (RTK) targets were functionally evaluated by high-content imaging to capture the complexity of PDO responses. Our drug discovery strategy resulted in the generation of MCLA-158, a bAb that specifically triggers epidermal growth factor receptor degradation in leucine-rich repeat-containing G-protein-coupled receptor 5-positive (LGR5+) cancer stem cells but shows minimal toxicity toward healthy LGR5+ colon stem cells. MCLA-158 exhibits therapeutic properties such as growth inhibition of KRAS-mutant colorectal cancers, blockade of metastasis initiation and suppression of tumor outgrowth in preclinical models for several epithelial cancer types.

A differentiation-based phylogeny of cancer subtypes.

Histopathological classification of human tumors relies in part on the degree of differentiation of the tumor sample. To date, there is no objective systematic method to categorize tumor subtypes by maturation. In this paper, we introduce a novel computational algorithm to rank tumor subtypes according to the dissimilarity of their gene expression from that of stem cells and fully differentiated tissue, and thereby construct a phylogenetic tree of cancer. We validate our methodology with expression data of leukemia, breast cancer and liposarcoma subtypes and then apply it to a broader group of sarcomas. This ranking of tumor subtypes resulting from the application of our methodology allows the identification of genes correlated with differentiation and may help to identify novel therapeutic targets. Our algorithm represents the first phylogeny-based tool to analyze the differentiation status of human tumors.

Stability of Imprinting and Differentiation Capacity in Naïve Human Cells Induced by Chemical Inhibition of CDK8 and CDK19.

Pluripotent stem cells can be stabilized in vitro at different developmental states by the use of specific chemicals and soluble factors. The naïve and primed states are the best characterized pluripotency states. Naïve pluripotent stem cells (PSCs) correspond to the early pre-implantation blastocyst and, in mice, constitute the optimal starting state for subsequent developmental applications. However, the stabilization of human naïve PSCs remains challenging because, after short-term culture, most current methods result in karyotypic abnormalities, aberrant DNA methylation patterns, loss of imprinting and severely compromised developmental potency. We have recently developed a novel method to induce and stabilize naïve human PSCs that consists in the simple addition of a chemical inhibitor for the closely related CDK8 and CDK19 kinases (CDK8/19i). Long-term cultured CDK8/19i-naïve human PSCs preserve their normal karyotype and do not show widespread DNA demethylation. Here, we investigate the long-term stability of allele-specific methylation at imprinted loci and the differentiation potency of CDK8/19i-naïve human PSCs. We report that long-term cultured CDK8/19i-naïve human PSCs retain the imprinting profile of their parental primed cells, and imprints are further retained upon differentiation in the context of teratoma formation. We have also tested the capacity of long-term cultured CDK8/19i-naïve human PSCs to differentiate into primordial germ cell (PGC)-like cells (PGCLCs) and trophoblast stem cells (TSCs), two cell types that are accessible from the naïve state. Interestingly, long-term cultured CDK8/19i-naïve human PSCs differentiated into PGCLCs with a similar efficiency to their primed counterparts. Also, long-term cultured CDK8/19i-naïve human PSCs were able to differentiate into TSCs, a transition that was not possible for primed PSCs. We conclude that inhibition of CDK8/19 stabilizes human PSCs in a functional naïve state that preserves imprinting and potency over long-term culture.

Global hyperactivation of enhancers stabilizes human and mouse naive pluripotency through inhibition of CDK8/19 Mediator kinases.

Pluripotent stem cells (PSCs) transition between cell states in vitro, reflecting developmental changes in the early embryo. PSCs can be stabilized in the naive state by blocking extracellular differentiation stimuli, particularly FGF-MEK signalling. Here, we report that multiple features of the naive state in human and mouse PSCs can be recapitulated without affecting FGF-MEK signalling or global DNA methylation. Mechanistically, chemical inhibition of CDK8 and CDK19 (hereafter CDK8/19) kinases removes their ability to repress the Mediator complex at enhancers. CDK8/19 inhibition therefore increases Mediator-driven recruitment of RNA polymerase II (RNA Pol II) to promoters and enhancers. This efficiently stabilizes the naive transcriptional program and confers resistance to enhancer perturbation by BRD4 inhibition. Moreover, naive pluripotency during embryonic development coincides with a reduction in CDK8/19. We conclude that global hyperactivation of enhancers drives naive pluripotency, and this can be achieved in vitro by inhibiting CDK8/19 kinase activity. These principles may apply to other contexts of cellular plasticity.

CD98hc (SLC3A2) sustains amino acid and nucleotide availability for cell cycle progression.

CD98 heavy chain (CD98hc) forms heteromeric amino acid (AA) transporters by interacting with different light chains. Cancer cells overexpress CD98hc-transporters in order to meet their increased nutritional and antioxidant demands, since they provide branched-chain AA (BCAA) and aromatic AA (AAA) availability while protecting cells from oxidative stress. Here we show that BCAA and AAA shortage phenocopies the inhibition of mTORC1 signalling, protein synthesis and cell proliferation caused by CD98hc ablation. Furthermore, our data indicate that CD98hc sustains glucose uptake and glycolysis, and, as a consequence, the pentose phosphate pathway (PPP). Thus, loss of CD98hc triggers a dramatic reduction in the nucleotide pool, which leads to replicative stress in these cells, as evidenced by the enhanced DNA Damage Response (DDR), S-phase delay and diminished rate of mitosis, all recovered by nucleoside supplementation. In addition, proper BCAA and AAA availability sustains the expression of the enzyme ribonucleotide reductase. In this regard, BCAA and AAA shortage results in decreased content of deoxynucleotides that triggers replicative stress, also recovered by nucleoside supplementation. On the basis of our findings, we conclude that CD98hc plays a central role in AA and glucose cellular nutrition, redox homeostasis and nucleotide availability, all key for cell proliferation.

Blimp-1 Mediates Tracheal Lumen Maturation in Drosophila melanogaster.

The specification of tissue identity during embryonic development requires precise spatio-temporal coordination of gene expression. Many transcription factors required for the development of organs have been identified and their expression patterns are known; however, the mechanisms through which they coordinate gene expression in time remain poorly understood. Here, we show that hormone-induced transcription factor Blimp-1 participates in the temporal coordination of tubulogenesis in Drosophila melanogaster by regulating the expression of many genes involved in tube maturation. In particular, we demonstrate that Blimp-1 regulates the expression of genes involved in chitin deposition and F-actin organization. We show that Blimp-1 is involved in the temporal control of lumen maturation by regulating the beginning of chitin deposition. We also report that Blimp-1 represses a variety of genes involved in tracheal maturation. Finally, we reveal that the kinase Btk29A serves as a link between Blimp-1 transcriptional repression and apical extracellular matrix organization.

Regulation of Mammary Luminal Cell Fate and Tumorigenesis by p38α.

Mammary stem and progenitor cells are essential for mammary gland homeostasis and are also candidates for cells of origin of mammary tumors. Here, we have investigated the function of the protein kinase p38α in the mammary gland using mice that delete this protein in the luminal epithelial cells. We show that p38α regulates the fate of luminal progenitor cells through modulation of the transcription factor RUNX1, an important controller of the estrogen receptor-positive cell lineage. We also provide evidence that the regulation of RUNX1 by p38α probably involves the kinase MSK1, which phosphorylates histone H3 at the RUNX1 promoter. Moreover, using a mouse model for breast cancer initiated by luminal cells, we show that p38α downregulation in mammary epithelial cells reduces tumor burden, which correlates with decreased numbers of tumor-initiating cells. Collectively, our results define a key role for p38α in luminal progenitor cell fate that affects mammary tumor formation.