Integrative metabolomics-genomics analysis identifies key networks in a stem cell-based model of schizophrenia.

Schizophrenia (SCZ) is a neuropsychiatric disorder, caused by a combination of genetic and environmental factors. The etiology behind the disorder remains elusive although it is hypothesized to be associated with the aberrant response to neurotransmitters, such as dopamine and glutamate. Therefore, investigating the link between dysregulated metabolites and distorted neurodevelopment holds promise to offer valuable insights into the underlying mechanism of this complex disorder. In this study, we aimed to explore a presumed correlation between the transcriptome and the metabolome in a SCZ model based on patient-derived induced pluripotent stem cells (iPSCs). For this, iPSCs were differentiated towards cortical neurons and samples were collected longitudinally at various developmental stages, reflecting neuroepithelial-like cells, radial glia, young and mature neurons. The samples were analyzed by both RNA-sequencing and targeted metabolomics and the two modalities were used to construct integrative networks in silico. This multi-omics analysis revealed significant perturbations in the polyamine and gamma-aminobutyric acid (GABA) biosynthetic pathways during rosette maturation in SCZ lines. We particularly observed the downregulation of the glutamate decarboxylase encoding genes GAD1 and GAD2, as well as their protein product GAD65/67 and their biochemical product GABA in SCZ samples. Inhibition of ornithine decarboxylase resulted in further decrease of GABA levels suggesting a compensatory activation of the ornithine/putrescine pathway as an alternative route for GABA production. These findings indicate an imbalance of cortical excitatory/inhibitory dynamics occurring during early neurodevelopmental stages in SCZ. Our study supports the hypothesis of disruption of inhibitory circuits to be causative for SCZ and establishes a novel in silico approach that enables for integrative correlation of metabolic and transcriptomic data of psychiatric disease models.

Spatially resolved multi-omics deciphers bidirectional tumor-host interdependence in glioblastoma.

Glioblastomas are malignant tumors of the central nervous system hallmarked by subclonal diversity and dynamic adaptation amid developmental hierarchies. The source of dynamic reorganization within the spatial context of these tumors remains elusive. Here, we characterized glioblastomas by spatially resolved transcriptomics, metabolomics, and proteomics. By deciphering regionally shared transcriptional programs across patients, we infer that glioblastoma is organized by spatial segregation of lineage states and adapts to inflammatory and/or metabolic stimuli, reminiscent of the reactive transformation in mature astrocytes. Integration of metabolic imaging and imaging mass cytometry uncovered locoregional tumor-host interdependence, resulting in spatially exclusive adaptive transcriptional programs. Inferring copy-number alterations emphasizes a spatially cohesive organization of subclones associated with reactive transcriptional programs, confirming that environmental stress gives rise to selection pressure. A model of glioblastoma stem cells implanted into human and rodent neocortical tissue mimicking various environments confirmed that transcriptional states originate from dynamic adaptation to various environments.

Wiskott-Aldrich syndrome protein forms nuclear condensates and regulates alternative splicing.

The diverse functions of WASP, the deficiency of which causes Wiskott-Aldrich syndrome (WAS), remain poorly defined. We generated three isogenic WAS models using patient induced pluripotent stem cells and genome editing. These models recapitulated WAS phenotypes and revealed that WASP deficiency causes an upregulation of numerous RNA splicing factors and widespread altered splicing. Loss of WASP binding to splicing factor gene promoters frequently leads to aberrant epigenetic activation. WASP interacts with dozens of nuclear speckle constituents and constrains SRSF2 mobility. Using an optogenetic system, we showed that WASP forms phase-separated condensates that encompasses SRSF2, nascent RNA and active Pol II. The role of WASP in gene body condensates is corroborated by ChIPseq and RIPseq. Together our data reveal that WASP is a nexus regulator of RNA splicing that controls the transcription of splicing factors epigenetically and the dynamics of the splicing machinery through liquid-liquid phase separation.

Exhaustion of CD39-Expressing CD8+ T Cells in Crohn's Disease Is Linked to Clinical Outcome.

BACKGROUND & AIMS: Exhaustion of CD8 T cells has been suggested to inform different clinical outcomes in Crohn's disease, but detailed analyses are lacking. This study aimed to identify the role of exhaustion on a single-cell level and identify relevant CD8 T cell populations in Crohn's disease. METHODS: Blood and intestinal tissue from 58 patients with Crohn's disease (active disease or remission) were assessed for CD8 T cell expression of exhaustion markers and their cytokine profile by highly multiplexed flow and mass cytometry. Key disease-associated subsets were sorted and analyzed by RNA sequencing. CD39 inhibition assays were performed in vitro. RESULTS: Activated CD39+ and CD39+PD-1+ CD8 T cell subsets expressing multiple exhaustion markers were enriched at low frequency in active Crohn's disease. Their cytokine production capacity was inversely linked to the Harvey-Bradshaw Index. Subset-level protein and transcriptome profiling revealed co-existence of effector and exhaustion programs in CD39+ and CD39+ PD-1+CD8 T cells, with CD39+ cells likely originating from the intestine. CD39 enzymatic activity controlled T cell cytokine production. Importantly, transcriptional exhaustion signatures were enriched in remission in CD39-expressing subsets with up-regulation of TOX. Subset-level transcriptomics revealed a CD39-related gene module that is associated with the clinical course. CONCLUSIONS: These data showed a role for the exhaustion of peripheral CD39-expressing CD8 T cell subsets in Crohn's disease. Their low frequency illustrated the utility of single-cell cytometry methods for identification of relevant immune populations. Importantly, the link of their exhaustion status to the clinical activity and their specific gene signatures have implications for exhaustion-based personalized medicine approaches.

Dynamic transcriptome analysis reveals signatures of paradoxical effect of vemurafenib on human dermal fibroblasts.

BACKGROUND: Vemurafenib (PLX4032) is one of the most frequently used treatments for late-stage melanoma patients with the BRAFV600E mutation; however, acquired resistance to the drug poses as a major challenge. It remains to be determined whether off-target effects of vemurafenib on normal stroma components could reshape the tumor microenvironment in a way that contributes to cancer progression and drug resistance. METHODS: By using temporally-resolved RNA- and ATAC-seq, we studied the early molecular changes induced by vemurafenib in human dermal fibroblast (HDF), a main stromal component in melanoma and other tumors with high prevalence of BRAFV600 mutations. RESULTS: Transcriptomics analyses revealed a stepwise up-regulation of proliferation signatures, together with a down-regulation of autophagy and proteolytic processes. The gene expression changes in HDF strongly correlated in an inverse way with those in BRAFV600E mutant malignant melanoma (MaMel) cell lines, consistent with the observation of a paradoxical effect of vemurafenib, leading to hyperphosphorylation of MEK1/2 and ERK1/2. The transcriptional changes in HDF were not strongly determined by alterations in chromatin accessibility; rather, an already permissive chromatin landscape seemed to facilitate the early accessibility to MAPK/ERK-regulated transcription factor binding sites. Combinatorial treatment with the MEK inhibitor trametinib did not preclude the paradoxical activation of MAPK/ERK signaling in HDF. When administered together, vemurafenib partially compensated for the reduction of cell viability and proliferation induced by trametinib. These paradoxical changes were restrained by using the third generation BRAF inhibitor PLX8394, a so-called paradox breaker compound. However, the advantageous effects on HDF during combination therapies were also lost. CONCLUSIONS: Vemurafenib induces paradoxical changes in HDF, enabled by a permissive chromatin landscape. These changes might provide an advantage during combination therapies, by compensating for the toxicity induced in stromal cells by less specific MAPK/ERK inhibitors. Our results highlight the relevance of evaluating the effects of the drugs on non-transformed stromal components, carefully considering the implications of their administration either as mono- or combination therapies. Video Abstract.

Clinical evolution, genetic landscape and trajectories of clonal hematopoiesis in SAMD9/SAMD9L syndromes.

Germline SAMD9 and SAMD9L mutations (SAMD9/9Lmut) predispose to myelodysplastic syndromes (MDS) with propensity for somatic rescue. In this study, we investigated a clinically annotated pediatric MDS cohort (n = 669) to define the prevalence, genetic landscape, phenotype, therapy outcome and clonal architecture of SAMD9/9L syndromes. In consecutively diagnosed MDS, germline SAMD9/9Lmut accounted for 8% and were mutually exclusive with GATA2 mutations present in 7% of the cohort. Among SAMD9/9Lmut cases, refractory cytopenia was the most prevalent MDS subtype (90%); acquired monosomy 7 was present in 38%; constitutional abnormalities were noted in 57%; and immune dysfunction was present in 28%. The clinical outcome was independent of germline mutations. In total, 67 patients had 58 distinct germline SAMD9/9Lmut clustering to protein middle regions. Despite inconclusive in silico prediction, 94% of SAMD9/9Lmut suppressed HEK293 cell growth, and mutations expressed in CD34+ cells induced overt cell death. Furthermore, we found that 61% of SAMD9/9Lmut patients underwent somatic genetic rescue (SGR) resulting in clonal hematopoiesis, of which 95% was maladaptive (monosomy 7 ± cancer mutations), and 51% had adaptive nature (revertant UPD7q, somatic SAMD9/9Lmut). Finally, bone marrow single-cell DNA sequencing revealed multiple competing SGR events in individual patients. Our findings demonstrate that SGR is common in SAMD9/9Lmut MDS and exemplify the exceptional plasticity of hematopoiesis in children.

Negative correlation of single-cell PAX3:FOXO1 expression with tumorigenicity in rhabdomyosarcoma.

Rhabdomyosarcomas (RMS) are phenotypically and functionally heterogeneous. Both primary human RMS cultures and low-passage Myf6Cre,Pax3:Foxo1,p53 mouse RMS cell lines, which express the fusion oncoprotein Pax3:Foxo1 and lack the tumor suppressor Tp53 (Myf6Cre,Pax3:Foxo1,p53), exhibit marked heterogeneity in PAX3:FOXO1 (P3F) expression at the single cell level. In mouse RMS cells, P3F expression is directed by the Pax3 promoter and coupled to eYFP YFPlow/P3Flow mouse RMS cells included 87% G0/G1 cells and reorganized their actin cytoskeleton to produce a cellular phenotype characterized by more efficient adhesion and migration. This translated into higher tumor-propagating cell frequencies of YFPlow/P3Flow compared with YFPhigh/P3Fhigh cells. Both YFPlow/P3Flow and YFPhigh/P3Fhigh cells gave rise to mixed clones in vitro, consistent with fluctuations in P3F expression over time. Exposure to the anti-tropomyosin compound TR100 disrupted the cytoskeleton and reversed enhanced migration and adhesion of YFPlow/P3Flow RMS cells. Heterogeneous expression of PAX3:FOXO1 at the single cell level may provide a critical advantage during tumor progression.

Age-dependent instability of mature neuronal fate in induced neurons from Alzheimer's patients.

Sporadic Alzheimer's disease (AD) exclusively affects elderly people. Using direct conversion of AD patient fibroblasts into induced neurons (iNs), we generated an age-equivalent neuronal model. AD patient-derived iNs exhibit strong neuronal transcriptome signatures characterized by downregulation of mature neuronal properties and upregulation of immature and progenitor-like signaling pathways. Mapping iNs to longitudinal neuronal differentiation trajectory data demonstrated that AD iNs reflect a hypo-mature neuronal identity characterized by markers of stress, cell cycle, and de-differentiation. Epigenetic landscape profiling revealed an underlying aberrant neuronal state that shares similarities with malignant transformation and age-dependent epigenetic erosion. To probe for the involvement of aging, we generated rejuvenated iPSC-derived neurons that showed no significant disease-related transcriptome signatures, a feature that is consistent with epigenetic clock and brain ontogenesis mapping, which indicate that fibroblast-derived iNs more closely reflect old adult brain stages. Our findings identify AD-related neuronal changes as age-dependent cellular programs that impair neuronal identity.

Single-Cell RNA-Seq Reveals that CD9 Is a Negative Marker of Glucose-Responsive Pancreatic ß-like Cells Derived from Human Pluripotent Stem Cells.

To date, it remains unclear if there are specific cell-surface markers for purifying glucose-responsive pancreatic ß-like cells derived from human pluripotent stem cells (hPSCs). In searching for this, we generated an efficient protocol for differentiating ß-like cells from human embryonic stem cells. We performed single-cell RNA sequencing and found that CD9 is a negative cell-surface marker of ß-like cells, as most INS+ cells are CD9-. We purified ß-like cells for spontaneous formation of islet-like organoids against CD9, and found significantly more NKX6.1+MAFA+C-PEPTIDE+ ß-like cells in the CD9- than in the CD9+ population. CD9- cells also demonstrate better glucose responsiveness than CD9+ cells. In humans, we observe more CD9+C-PEPTIDE+ ß cells in the fetal than in the adult cadaveric islets and more Ki67+ proliferating cells among CD9+ fetal ß cells. Taken together, our experiments show that CD9 is a cell-surface marker for negative enrichment of glucose-responsive ß-like cells differentiated from hPSCs.

Single-cell transcriptomics uncover distinct innate and adaptive cell subsets during tissue homeostasis and regeneration.

Recently, immune cell-mediated tissue repair and regeneration has been an emerging paradigm of regenerative medicine. Immune cells form an essential part of the wound as induction of inflammation is a necessary step to elicit tissue healing. Rapid progress in transcriptomic analyses by high-throughput next-generation sequencing has been developed to study gene regulatory network and establish molecular signatures of immune cells that could potentially predict their functional roles in tissue repair and regeneration. However, the identification of cellular heterogeneity especially on the rare cell subsets has been limited in transcriptomic analyses of bulk cell populations. Therefore, genome-wide, single-cell RNA sequencing (scRNA-Seq) has offered an unprecedented approach to unravel cellular diversity and to study novel immune cell populations involved in tissue repair and regeneration through unsupervised sampling of individual cells without the need to rely on prior knowledge about cell-specific markers. The analysis of gene expression patterns at a single-cell resolution also holds promises to uncover the mechanisms and therefore the development of therapeutic strategy promoting immunoregenerative medicine. In this review, we will discuss how scRNA-Seq facilitates the characterization of immune cells, including macrophages, innate lymphoid cells and T and B lymphocytes, discovery of immune cell heterogeneity, identification of novel subsets, and tracking of developmental trajectories of distinct immune cells during tissue homeostasis, repair, and regeneration.

Premature Activation of Immune Transcription Programs in Autoimmune-Predisposed Mouse Embryonic Stem Cells and Blastocysts.

Autoimmune diabetes is a complex multifactorial disease with genetic and environmental factors playing pivotal roles. While many genes associated with the risk of diabetes have been identified to date, the mechanisms by which external triggers contribute to the genetic predisposition remain unclear. Here, we derived embryonic stem (ES) cell lines from diabetes-prone non-obese diabetic (NOD) and healthy C57BL/6 (B6) mice. While overall pluripotency markers were indistinguishable between newly derived NOD and B6 ES cells, we discovered several differentially expressed genes that normally are not expressed in ES cells. Several genes that reside in previously identified insulin-dependent diabetics (Idd) genomic regions were up-regulated in NOD ES cells. Gene set enrichment analysis showed that different groups of genes associated with immune functions are differentially expressed in NOD. Transcriptomic analysis of NOD blastocysts validated several differentially overexpressed Idd genes compared to B6. Genome-wide mapping of active histone modifications using ChIP-Seq supports active expression as the promoters and enhancers of activated genes are also marked by active histone modifications. We have also found that NOD ES cells secrete more inflammatory cytokines. Our data suggest that the known genetic predisposition of NOD to autoimmune diabetes leads to epigenetic instability of several Idd regions.

Prediabetes Induced by a Single Autoimmune B Cell Clone.

While B cells play a significant role in the onset of type-1 diabetes (T1D), little is know about their role in those early stages. Thus, to gain new insights into the role of B cells in T1D, we converted a physiological early pancreas-infiltrating B cell into a novel BCR mouse model using Somatic Cell Nuclear Transfer (SCNT). Strikingly, SCNT-derived B1411 model displayed neither developmental block nor anergy. Instead, B1411 underwent spontaneous germinal center reactions. Without T cell help, B1411-Rag1-/- was capable of forming peri-/intra-pancreatic lymph nodes, and undergoing class-switching. RNA-Seq analysis identified 93 differentially expressed genes in B1411 compared to WT B cells, including Irf7, Usp18, and Mda5 that had been linked to a potential viral etiology of T1D. We also found various members of the oligoadenylate synthase (OAS) family to be enriched in B1411, such as Oas1, which had recently also been linked to T1D. Strikingly, when challenged with glucose B1411-Rag1-/- mice displayed impaired glucose tolerance.

Extraction of Distinct Neuronal Cell Types from within a Genetically Continuous Population.

Single-cell transcriptomics of neocortical neurons have revealed more than 100 clusters corresponding to putative cell types. For inhibitory and subcortical projection neurons (SCPNs), there is a strong concordance between clusters and anatomical descriptions of cell types. In contrast, cortico-cortical projection neurons (CCPNs) separate into surprisingly few transcriptomic clusters, despite their diverse anatomical projection types. We used projection-dependent single-cell transcriptomic analyses and monosynaptic rabies tracing to compare mouse primary visual cortex CCPNs projecting to different higher visual areas. We find that layer 2/3 CCPNs with different anatomical projections differ systematically in their gene expressions, despite forming only a single genetic cluster. Furthermore, these neurons receive feedback selectively from the same areas to which they project. These findings demonstrate that gene-expression analysis in isolation is insufficient to identify neuron types and have important implications for understanding the functional role of cortical feedback circuits.

MicroRNA-146a regulates immune-related adverse events caused by immune checkpoint inhibitors.

Immune checkpoint inhibitor (ICI) therapy has shown a significant benefit in the treatment of a variety of cancer entities. However, immune-related adverse events (irAEs) occur frequently and can lead to ICI treatment termination. MicroRNA-146a (miR-146a) has regulatory functions in immune cells. We observed that mice lacking miR-146a developed markedly more severe irAEs compared with WT mice in several irAE target organs in 2 different murine models. miR-146a-/- mice exhibited increased T cell activation and effector function upon ICI treatment. Moreover, neutrophil numbers in the spleen and the inflamed intestine were highly increased in ICI-treated miR-146a-/- mice. Therapeutic administration of a miR-146a mimic reduced irAE severity. To validate our preclinical findings in patients, we analyzed the effect of a SNP in the MIR146A gene on irAE severity in 167 patients treated with ICIs. We found that the SNP rs2910164 leading to reduced miR-146a expression was associated with an increased risk of developing severe irAEs, reduced progression-free survival, and increased neutrophil counts both at baseline and during ICI therapy. In conclusion, we characterized miR-146a as a molecular target for preventing ICI-mediated autoimmune dysregulation. Furthermore, we identified the MIR146A SNP rs2910164 as a biomarker to predict severe irAE development in ICI-treated patients.

Chemical modulation of transcriptionally enriched signaling pathways to optimize the conversion of fibroblasts into neurons.

Direct conversion of human somatic fibroblasts into induced neurons (iNs) allows for the generation of functional neurons while bypassing any stem cell intermediary stages. Although iN technology has an enormous potential for modeling age-related diseases, as well as therapeutic approaches, the technology faces limitations due to variable conversion efficiencies and a lack of thorough understanding of the signaling pathways directing iN conversion. Here, we introduce a new all-in-one inducible lentiviral system that simplifies fibroblast transgenesis for the two pioneer transcription factors, Ngn2 and Ascl1, and markedly improves iN yields. Further, our timeline RNA-Seq data across the course of conversion has identified signaling pathways that become transcriptionally enriched during iN conversion. Small molecular modulators were identified for four signaling pathways that reliably increase the yield of iNs. Taken together, these advances provide an improved toolkit for iN technology and new insight into the mechanisms influencing direct iN conversion.

Pathological priming causes developmental gene network heterochronicity in autistic subject-derived neurons.

Autism spectrum disorder (ASD) is thought to emerge during early cortical development. However, the exact developmental stages and associated molecular networks that prime disease propensity are elusive. To profile early neurodevelopmental alterations in ASD with macrocephaly, we monitored subject-derived induced pluripotent stem cells (iPSCs) throughout the recapitulation of cortical development. Our analysis revealed ASD-associated changes in the maturational sequence of early neuron development, involving temporal dysregulation of specific gene networks and morphological growth acceleration. The observed changes tracked back to a pathologically primed stage in neural stem cells (NSCs), reflected by altered chromatin accessibility. Concerted over-representation of network factors in control NSCs was sufficient to trigger ASD-like features, and circumventing the NSC stage by direct conversion of ASD iPSCs into induced neurons abolished ASD-associated phenotypes. Our findings identify heterochronic dynamics of a gene network that, while established earlier in development, contributes to subsequent neurodevelopmental aberrations in ASD.

Single-cell transcriptomics reveal that PD-1 mediates immune tolerance by regulating proliferation of regulatory T cells.

BACKGROUND: We have previously reported an antigen-specific protocol to induce transplant tolerance and linked suppression to human embryonic stem cell (hESC)-derived tissues in immunocompetent mice through coreceptor and costimulation blockade. However, the exact mechanisms of acquired immune tolerance in this model have remained unclear. METHODS: We utilize the NOD.Foxp3hCD2 reporter mouse line and an ablative anti-hCD2 antibody to ask if CD4+FOXP3+ regulatory T cells (Treg) are required for coreceptor and costimulation blockade-induced immune tolerance. We also perform genome-wide single-cell RNA-sequencing to interrogate Treg during immune rejection and tolerance and to indicate possible mechanisms involved in sustaining Treg function. RESULTS: We show that Treg are indispensable for tolerance induced by coreceptor and costimulation blockade as depletion of which with an anti-hCD2 antibody resulted in rejection of hESC-derived pancreatic islets in NOD.Foxp3hCD2 mice. Single-cell transcriptomic profiling of 12,964 intragraft CD4+ T cells derived from rejecting and tolerated grafts reveals that Treg are heterogeneous and functionally distinct in the two outcomes of transplant rejection and tolerance. Treg appear to mainly promote chemotactic and ubiquitin-dependent protein catabolism during transplant rejection while seeming to harness proliferative and immunosuppressive function during tolerance. We also demonstrate that this form of acquired transplant tolerance is associated with increased proliferation and PD-1 expression by Treg. Blocking PD-1 signaling with a neutralizing anti-PD-1 antibody leads to reduced Treg proliferation and graft rejection. CONCLUSIONS: Our results suggest that short-term coreceptor and costimulation blockade mediates immune tolerance to hESC-derived pancreatic islets by promoting Treg proliferation through engagement of PD-1. Our findings could give new insights into clinical development of hESC-derived pancreatic tissues, combined with immunotherapies that expand intragraft Treg, as a potentially sustainable alternative treatment for T1D.

Regulatory T Cells Promote Apelin-Mediated Sprouting Angiogenesis in Type 2 Diabetes.

The role of CD4+ T cells in the ischemic tissues of T2D patients remains unclear. Here, we report that T2D patients' vascular density was negatively correlated with the number of infiltrating CD4+ T cells after ischemic injury. Th1 was the predominant subset, and Th1-derived IFN-γ and TNF-α directly impaired human angiogenesis. We then blocked CD4+ T cell infiltration into the ischemic tissues of both Leprdb/db and diet-induced obese T2D mice. Genome-wide RNA sequencing shows an increased proliferative and angiogenic capability of diabetic ECs in ischemic tissues. Moreover, wire myography shows enhanced EC function and laser Doppler imaging reveals improved post-ischemic blood reperfusion. Mechanistically, functional revascularization after CD4 coreceptor blockade was mediated by Tregs. Genetic lineage tracing via Cdh5-CreER and Apln-CreER and coculture assays further illustrate that Tregs increased vascular density and induced de novo sprouting angiogenesis in a paracrine manner. Taken together, our results reveal that Th1 impaired while Tregs promoted functional post-ischemic revascularization in obesity and diabetes.