Single-cell landscape of the ecosystem in early-relapse hepatocellular carcinoma.

Hepatocellular carcinoma (HCC) has high relapse and low 5-year survival rates. Single-cell profiling in relapsed HCC may aid in the design of effective anticancer therapies, including immunotherapies. We profiled the transcriptomes of ∼17,000 cells from 18 primary or early-relapse HCC cases. Early-relapse tumors have reduced levels of regulatory T cells, increased dendritic cells (DCs), and increased infiltrated CD8+ T cells, compared with primary tumors, in two independent cohorts. Remarkably, CD8+ T cells in recurrent tumors overexpressed KLRB1 (CD161) and displayed an innate-like low cytotoxic state, with low clonal expansion, unlike the classical exhausted state observed in primary HCC. The enrichment of these cells was associated with a worse prognosis. Differential gene expression and interaction analyses revealed potential immune evasion mechanisms in recurrent tumor cells that dampen DC antigen presentation and recruit innate-like CD8+ T cells. Our comprehensive picture of the HCC ecosystem provides deeper insights into immune evasion mechanisms associated with tumor relapse.

Senescence atlas reveals an aged-like inflamed niche that blunts muscle regeneration.

Tissue regeneration requires coordination between resident stem cells and local niche cells1,2. Here we identify that senescent cells are integral components of the skeletal muscle regenerative niche that repress regeneration at all stages of life. The technical limitation of senescent-cell scarcity3 was overcome by combining single-cell transcriptomics and a senescent-cell enrichment sorting protocol. We identified and isolated different senescent cell types from damaged muscles of young and old mice. Deeper transcriptome, chromatin and pathway analyses revealed conservation of cell identity traits as well as two universal senescence hallmarks (inflammation and fibrosis) across cell type, regeneration time and ageing. Senescent cells create an aged-like inflamed niche that mirrors inflammation associated with ageing (inflammageing4) and arrests stem cell proliferation and regeneration. Reducing the burden of senescent cells, or reducing their inflammatory secretome through CD36 neutralization, accelerates regeneration in young and old mice. By contrast, transplantation of senescent cells delays regeneration. Our results provide a technique for isolating in vivo senescent cells, define a senescence blueprint for muscle, and uncover unproductive functional interactions between senescent cells and stem cells in regenerative niches that can be overcome. As senescent cells also accumulate in human muscles, our findings open potential paths for improving muscle repair throughout life.

Cell transcriptomic atlas of the non-human primate Macaca fascicularis.

Studying tissue composition and function in non-human primates (NHPs) is crucial to understand the nature of our own species. Here we present a large-scale cell transcriptomic atlas that encompasses over 1 million cells from 45 tissues of the adult NHP Macaca fascicularis. This dataset provides a vast annotated resource to study a species phylogenetically close to humans. To demonstrate the utility of the atlas, we have reconstructed the cell-cell interaction networks that drive Wnt signalling across the body, mapped the distribution of receptors and co-receptors for viruses causing human infectious diseases, and intersected our data with human genetic disease orthologues to establish potential clinical associations. Our M. fascicularis cell atlas constitutes an essential reference for future studies in humans and NHPs.

MYB insufficiency disrupts proteostasis in hematopoietic stem cells, leading to age-related neoplasia.

MYB plays a key role in gene regulation throughout the hematopoietic hierarchy and is critical for the maintenance of normal hematopoietic stem cells (HSC). Acquired genetic dysregulation of MYB is involved in the etiology of a number of leukemias, although inherited noncoding variants of the MYB gene are a susceptibility factor for many hematological conditions, including myeloproliferative neoplasms (MPN). The mechanisms that connect variations in MYB levels to disease predisposition, especially concerning age dependency in disease initiation, are completely unknown. Here, we describe a model of Myb insufficiency in mice that leads to MPN, myelodysplasia, and leukemia in later life, mirroring the age profile of equivalent human diseases. We show that this age dependency is intrinsic to HSC, involving a combination of an initial defective cellular state resulting from small effects on the expression of multiple genes and a progressive accumulation of further subtle changes. Similar to previous studies showing the importance of proteostasis in HSC maintenance, we observed altered proteasomal activity and elevated proliferation indicators, followed by elevated ribosome activity in young Myb-insufficient mice. We propose that these alterations combine to cause an imbalance in proteostasis, potentially creating a cellular milieu favoring disease initiation.

NR1H3 (LXRα) is associated with pro-inflammatory macrophages, predicts survival and suggests potential therapeutic rationales in diffuse large b-cell lymphoma.

The role of macrophages (Mo) and their prognostic impact in diffuse large B-cell lymphomas (DLBCL) remain controversial. By regulating the lipid metabolism, Liver-X-Receptors (LXRs) control Mo polarization/inflammatory response, and their pharmacological modulation is under clinical investigation to treat human cancers, including lymphomas. Herein, we surveyed the role of LXRs in DLBCL for prognostic purposes. Comparing bulk tumors with purified malignant and normal B-cells, we found an intriguing association of NR1H3, encoding for the LXR-α isoform, with the tumor microenvironment (TME). CIBERSORTx-based purification on large DLBCL datasets revealed a high expression of the receptor transcript in M1-like pro-inflammatory Mo. By determining an expression cut-off of NR1H3, we used digital measurement to validate its prognostic capacity on two large independent on-trial and real-world cohorts. Independently of classical prognosticators, NR1H3high patients displayed longer survival compared with NR1H3low cases and a high-resolution Mo GEP dissection suggested a remarkable transcriptional divergence between subgroups. Overall, our findings indicate NR1H3 as a Mo-related biomarker identifying patients at higher risk and prompt future preclinical studies investigating its mouldability for therapeutic purposes.

Expression characteristics and potential function of non-coding RNA in mouse cortical cells.

Non-coding RNAs (ncRNAs) play essential regulatory functions in various physiological and pathological processes in the brain. To systematically characterize the ncRNA profile in cortical cells, we downloaded single-cell SMART-Seq v4 data of mouse cerebral cortex. Our results revealed that the ncRNAs alone are sufficient to define the identity of most cortical cell types. We identified 1,600 ncRNAs that exhibited cell type specificity, even yielding to distinguish microglia from perivascular macrophages with ncRNA. Moreover, we characterized cortical layer and region specific ncRNAs, in line with the results by spatial transcriptome (ST) data. By constructing a co-expression network of ncRNAs and protein-coding genes, we predicted the function of ncRNAs. By integrating with genome-wide association studies data, we established associations between cell type-specific ncRNAs and traits related to neurological disorders. Collectively, our study identified differentially expressed ncRNAs at multiple levels and provided the valuable resource to explore the functions and dysfunctions of ncRNAs in cortical cells.

A spatiotemporal atlas of mouse liver homeostasis and regeneration.

The mechanism by which mammalian liver cell responses are coordinated during tissue homeostasis and perturbation is poorly understood, representing a major obstacle in our understanding of many diseases. This knowledge gap is caused by the difficulty involved with studying multiple cell types in different states and locations, particularly when these are transient. We have combined Stereo-seq (spatiotemporal enhanced resolution omics-sequencing) with single-cell transcriptomic profiling of 473,290 cells to generate a high-definition spatiotemporal atlas of mouse liver homeostasis and regeneration at the whole-lobe scale. Our integrative study dissects in detail the molecular gradients controlling liver cell function, systematically defining how gene networks are dynamically modulated through intercellular communication to promote regeneration. Among other important regulators, we identified the transcriptional cofactor TBL1XR1 as a rheostat linking inflammation to Wnt/ß-catenin signaling for facilitating hepatocyte proliferation. Our data and analytical pipelines lay the foundation for future high-definition tissue-scale atlases of organ physiology and malfunction.

Cell atlas of CCl 4-induced progressive liver fibrosis reveals stage-specific responses.

Chronic liver injury leads to progressive liver fibrosis and ultimately cirrhosis, a major cause of morbidity and mortality worldwide. However, there are currently no effective anti-fibrotic therapies available, especially for late-stage patients, which is partly attributed to the major knowledge gap regarding liver cell heterogeneity and cell-specific responses in different fibrosis stages. To reveal the multicellular networks regulating mammalian liver fibrosis from mild to severe phenotypes, we generated a single-nucleus transcriptomic atlas encompassing 49 919 nuclei corresponding to all main liver cell types at different stages of murine carbon tetrachloride (CCl 4)-induced progressive liver fibrosis. Integrative analysis distinguished the sequential responses to injury of hepatocytes, hepatic stellate cells and endothelial cells. Moreover, we reconstructed cell-cell interactions and gene regulatory networks implicated in these processes. These integrative analyses uncovered previously overlooked aspects of hepatocyte proliferation exhaustion and disrupted pericentral metabolic functions, dysfunction for clearance by apoptosis of activated hepatic stellate cells, accumulation of pro-fibrotic signals, and the switch from an anti-angiogenic to a pro-angiogenic program during CCl 4-induced progressive liver fibrosis. Our dataset thus constitutes a useful resource for understanding the molecular basis of progressive liver fibrosis using a relevant animal model.

Spatially resolved gene regulatory and disease-related vulnerability map of the adult Macaque cortex.

Single cell approaches have increased our knowledge about the cell type composition of the non-human primate (NHP), but a detailed characterization of area-specific regulatory features remains outstanding. We generated single-cell transcriptomic and chromatin accessibility (single-cell ATAC) data of 358,237 cells from prefrontal cortex (PFC), primary motor cortex (M1) and primary visual cortex (V1) of adult female cynomolgus monkey brain, and integrated this dataset with Stereo-seq (spatial enhanced resolution omics-sequencing) of the corresponding cortical areas to assign topographic information to molecular states. We identified area-specific chromatin accessible sites and their targeted genes, including the cell type-specific transcriptional regulatory network associated with excitatory neurons heterogeneity. We reveal calcium ion transport and axon guidance genes related to specialized functions of PFC and M1, identified the similarities and differences between adult macaque and human oligodendrocyte trajectories, and mapped the genetic variants and gene perturbations of human diseases to NHP cortical cells. This resource establishes a transcriptomic and chromatin accessibility combinatory regulatory landscape at a single-cell and spatially resolved resolution in NHP cortex.

Global Profiling of the Lysine Crotonylome in Different Pluripotent States.

Pluripotent stem cells (PSCs) can be expanded in vitro in different culture conditions, resulting in a spectrum of cell states with distinct properties. Understanding how PSCs transition from one state to another, ultimately leading to lineage-specific differentiation, is important for developmental biology and regenerative medicine. Although there is significant information regarding gene expression changes controlling these transitions, less is known about post-translational modifications of proteins. Protein crotonylation is a newly discovered post-translational modification where lysine residues are modified with a crotonyl group. Here, we employed affinity purification of crotonylated peptides and liquid chromatography-tandem mass spectrometry (LC-MS/MS) to systematically profile protein crotonylation in mouse PSCs in different states including ground, metastable, and primed states, as well as metastable PSCs undergoing early pluripotency exit. We successfully identified 3628 high-confidence crotonylated sites in 1426 proteins. These crotonylated proteins are enriched for factors involved in functions/processes related to pluripotency such as RNA biogenesis, central carbon metabolism, and proteasome function. Moreover, we found that increasing the cellular levels of crotonyl-coenzyme A (crotonyl-CoA) through crotonic acid treatment promotes proteasome activity in metastable PSCs and delays their differentiation, consistent with previous observations showing that enhanced proteasome activity helps to sustain pluripotency. Our atlas of protein crotonylation will be valuable for further studies of pluripotency regulation and may also provide insights into the role of metabolism in other cell fate transitions.

Capture of the newly transcribed RNA interactome using click chemistry.

Application of synthetic nucleoside analogues to capture newly transcribed RNAs has unveiled key features of RNA metabolism. Whether this approach could be adapted to isolate the RNA-bound proteome (RNA interactome) was, however, unexplored. We have developed a new method (capture of the newly transcribed RNA interactome using click chemistry, or RICK) for the systematic identification of RNA-binding proteins based on the incorporation of 5-ethynyluridine into newly transcribed RNAs followed by UV cross-linking and click chemistry-mediated biotinylation. The RNA-protein adducts are then isolated by affinity capture using streptavidin-coated beads. Through high-throughput RNA sequencing and mass spectrometry, the RNAs and proteins can be elucidated globally. A typical RICK experimental procedure takes only 1 d, excluding the steps of cell preparation, 5-ethynyluridine labeling, validation (silver staining, western blotting, quantitative reverse-transcription PCR (qRT-PCR) or RNA sequencing (RNA-seq)) and proteomics. Major advantages of RICK are the capture of RNA-binding proteins interacting with any type of RNA and, particularly, the ability to discern between newly transcribed and steady-state RNAs through controlled labeling. Thanks to its versatility, RICK will facilitate the characterization of the total and newly transcribed RNA interactome in different cell types and conditions.

Single-Nucleus Chromatin Accessibility Landscape Reveals Diversity in Regulatory Regions Across Distinct Adult Rat Cortex.

Rats have been widely used as an experimental organism in psychological, pharmacological, and behavioral studies by modeling human diseases such as neurological disorders. It is critical to identify and characterize cell fate determinants and their regulatory mechanisms in single-cell resolutions across rat brain regions. Here, we applied droplet-based single-nucleus assay for transposase-accessible chromatin using sequencing (snATAC-seq) to systematically profile the single-cell chromatin accessibility across four dissected brain areas in adult Sprague-Dawley (SD) rats with a total of 59,023 single nuclei and identified 16 distinct cell types. Interestingly, we found that different cortex regions exhibit diversity in both cellular compositions and gene regulatory regions. Several cell-type-specific transcription factors (TFs), including SPI1, KLF4, KLF6, and NEUROD2, have been shown to play important roles during the pathogenesis of various neurological diseases, such as Alzheimer's disease (AD), astrocytic gliomas, autism spectrum disorder (ASD), and intellectual disabilities. Therefore, our single-nucleus atlas of rat cortex could serve as an invaluable resource for dissecting the regulatory mechanisms underlying diverse cortex cell fates and further revealing the regulatory networks of neuropathogenesis.

PHC1 maintains pluripotency by organizing genome-wide chromatin interactions of the Nanog locus.

Polycomb group (PcG) proteins maintain cell identity by repressing gene expression during development. Surprisingly, emerging studies have recently reported that a number of PcG proteins directly activate gene expression during cell fate determination process. However, the mechanisms by which they direct gene activation in pluripotency remain poorly understood. Here, we show that Phc1, a subunit of canonical polycomb repressive complex 1 (cPRC1), can exert its function in pluripotency maintenance via a PRC1-independent activation of Nanog. Ablation of Phc1 reduces the expression of Nanog and overexpression of Nanog partially rescues impaired pluripotency caused by Phc1 depletion. We find that Phc1 interacts with Nanog and activates Nanog transcription by stabilizing the genome-wide chromatin interactions of the Nanog locus. This adds to the already known canonical function of PRC1 in pluripotency maintenance via a PRC1-dependent repression of differentiation genes. Overall, our study reveals a function of Phc1 to activate Nanog transcription through regulating chromatin architecture and proposes a paradigm for PcG proteins to maintain pluripotency.

Role of Long Non-coding RNAs in Reprogramming to Induced Pluripotency.

The generation of induced pluripotent stem cells through somatic cell reprogramming requires a global reorganization of cellular functions. This reorganization occurs in a multi-phased manner and involves a gradual revision of both the epigenome and transcriptome. Recent studies have shown that the large-scale transcriptional changes observed during reprogramming also apply to long non-coding RNAs (lncRNAs), a type of traditionally neglected RNA species that are increasingly viewed as critical regulators of cellular function. Deeper understanding of lncRNAs in reprogramming may not only help to improve this process but also have implications for studying cell plasticity in other contexts, such as development, aging, and cancer. In this review, we summarize the current progress made in profiling and analyzing the role of lncRNAs in various phases of somatic cell reprogramming, with emphasis on the re-establishment of the pluripotency gene network and X chromosome reactivation.

Oxidised metabolites of the omega-6 fatty acid linoleic acid activate dFOXO.

Obesity-induced inflammation, or meta-inflammation, plays key roles in metabolic syndrome and is a significant risk factor in diabetes and cardiovascular disease. To investigate causal links between obesity, meta-inflammation, and insulin signaling we established a Drosophila model to determine how elevated dietary fat and changes in the levels and balance of saturated fatty acids (SFAs) and polyunsaturated fatty acids (PUFAs) influence inflammation. We observe negligible effect of saturated fatty acid on inflammation but marked enhancement or suppression by omega-6 and omega-3 PUFAs, respectively. Using combined lipidomic and genetic analysis, we show omega-6 PUFA enhances meta-inflammation by producing linoleic acid-derived lipid mediator 9-hydroxy-octadecadienoic acid (9-HODE). Transcriptome analysis reveals 9-HODE functions by regulating FOXO family transcription factors. We show 9-HODE activates JNK, triggering FOXO nuclear localisation and chromatin binding. FOXO TFs are important transducers of the insulin signaling pathway that are normally down-regulated by insulin. By activating FOXO, 9-HODE could antagonise insulin signaling providing a molecular conduit linking changes in dietary fatty acid balance, meta-inflammation, and insulin resistance.

Generation of an induced pluripotent stem cell line (GIBHi004-A) from a Parkinson's disease patient with mutant DJ-1/PARK7 (p.L10P).

Mutations occurring in the gene body of PARK7 (encoding DJ-1/PARK7) cause autosomal recessive early-onset parkinsonism (AREP). These mutations produce a loss of function and have been reported to lead to dopaminergic neuron degeneration in the substantia nigra. However, the underlying mechanisms are largely unknown. Here, we report the generation of a patient-derived induced pluripotent stem cell (iPSC) line carrying mutant DJ-1 (p.L10P). This cell line is a valuable tool for in vitro Parkinson's disease (PD) modeling and mechanistic studies.

JMJD3 acts in tandem with KLF4 to facilitate reprogramming to pluripotency.

The interplay between the Yamanaka factors (OCT4, SOX2, KLF4 and c-MYC) and transcriptional/epigenetic co-regulators in somatic cell reprogramming is incompletely understood. Here, we demonstrate that the histone H3 lysine 27 trimethylation (H3K27me3) demethylase JMJD3 plays conflicting roles in mouse reprogramming. On one side, JMJD3 induces the pro-senescence factor Ink4a and degrades the pluripotency regulator PHF20 in a reprogramming factor-independent manner. On the other side, JMJD3 is specifically recruited by KLF4 to reduce H3K27me3 at both enhancers and promoters of epithelial and pluripotency genes. JMJD3 also promotes enhancer-promoter looping through the cohesin loading factor NIPBL and ultimately transcriptional elongation. This competition of forces can be shifted towards improved reprogramming by using early passage fibroblasts or boosting JMJD3's catalytic activity with vitamin C. Our work, thus, establishes a multifaceted role for JMJD3, placing it as a key partner of KLF4 and a scaffold that assists chromatin interactions and activates gene transcription.

CEBPA-mutated leukemia is sensitive to genetic and pharmacological targeting of the MLL1 complex.

The gene encoding the transcription factor C/EBPα is mutated in 10-15% of acute myeloid leukemia (AML) patients. N-terminal CEBPA mutations cause ablation of full-length C/EBPα without affecting the expression of a shorter oncogenic isoform, termed p30. The mechanistic basis of p30-induced leukemogenesis is incompletely understood. Here, we demonstrate that the MLL1 histone-methyltransferase complex represents a critical actionable vulnerability in CEBPA-mutated AML. Oncogenic C/EBPα p30 and MLL1 show global co-localization on chromatin and p30 exhibits robust physical interaction with the MLL1 complex. CRISPR/Cas9-mediated mutagenesis of MLL1 results in proliferation arrest and myeloid differentiation in C/EBPα p30-expressing cells. In line, CEBPA-mutated hematopoietic progenitor cells are hypersensitive to pharmacological targeting of the MLL1 complex. Inhibitor treatment impairs proliferation and restores myeloid differentiation potential in mouse and human AML cells with CEBPA mutations. Finally, we identify the transcription factor GATA2 as a direct critical target of the p30-MLL1 interaction. Altogether, we show that C/EBPα p30 requires the MLL1 complex to regulate oncogenic gene expression and that CEBPA-mutated AML is hypersensitive to perturbation of the MLL1 complex. These findings identify the MLL1 complex as a potential therapeutic target in AML with CEBPA mutations.

Dependence on Myb expression is attenuated in myeloid leukaemia with N-terminal CEBPA mutations.

Mutations at the N- or C-terminus of C/EBPα are frequent in acute myeloid leukaemia (AML) with normal karyotype. Here, we investigate the role of the transcription factor Myb in AMLs driven by different combinations of CEBPA mutations. Using knockdown of Myb in murine cell lines modelling the spectrum of CEBPA mutations, we show that the effect of reduced Myb depends on the mutational status of the two Cebpa alleles. Importantly, Myb knockdown fails to override the block in myeloid differentiation in cells with biallelic N-terminal C/EBPα mutations, demonstrating for the first time that the dependency on Myb is much lower in AML with this mutational profile. By comparing gene expression following Myb knockdown and chromatin immunoprecipitation sequencing data for the binding of C/EBPα isoforms, we provide evidence for a functional cooperation between C/EBPα and Myb in the maintenance of AML. This co-dependency breaks down when both alleles of CEBPA harbour N-terminal mutations, as a subset of C/EBPα-regulated genes only bind the short p30 C/EBPα isoform and, unlike other C/EBPα-regulated genes, do so without a requirement for Myb.