Search

BackgroundGlioblastoma (GBM) evades the immune system by creating an immune-suppressive tumor microenvironment (TME), where GBM-associated myeloid cells are geared towards tumor-supportive roles. However, it is unclear whether recruited myeloid cells are phenotypically and functionally identical. Here, we aim to understand the TME heterogeneity in GBM patients recapitulated in patient-derived orthotopic xenografts (PDOXs) and systematically characterize myeloid cell type identities at the molecular and functional level.MethodsWe applied single-cell RNA-sequencing and spatial transcriptomics, multicolor flow cytometry, immunohistochemistry and functional assays to examine the heterogeneity of the TME in GBM. Various GBM PDOXs representing different tumor phenotypes were analyzed and compared to the patient tumors, normal brain and mouse GL261 glioma model.ResultsPDOX models recapitulate the major components of the TME detected in human GBM, where tumor cells reciprocally interact with host cells to create a GBM-specific TME. We detect the most prominent transcriptomic adaptations in myeloid cells, which are largely of microglial origin. We reveal intra-tumoral heterogeneity of microglia and identify diverse phenotypic states across distinct GBM landscapes and tumor niches. GBM-educated microglia acquire dendritic cell-like features, displaying increased migration and phagocytosis. We further find novel microglial states expressing astrocytic and endothelial markers. Lastly, we show that temozolomide (TMZ) treatment leads to transcriptomic plasticity of both GBM tumor cells and adjacent TME components.ConclusionsOur data provide insight into the phenotypic adaptation of the heterogeneous TME instructed by GBM. We uncover that GBM-educated microglia are represented by various concomitant states, both in patients and recapitulated in PDOXs, displaying different pro- or anti-tumoral properties that are modulated by anti-neoplastic treatments, such as TMZ.KEY POINTSGBM-educated tumor microenvironment is faithfully recapitulated and modulated in PDOX modelsMicroglia represent an essential myeloid cell population in the GBM microenvironmentGBM-educated microglia acquire heterogeneous transcriptomic states across distinct tumor nichesGBM-educated microglia subsets display phagocytic and dendritic cell-like gene expression programs, which are modulated upon TMZ treatmentIMPORTANCE OF THE STUDYThis manuscript addresses tumor-immune interactions in GBM, focusing on the molecular changes of the myeloid compartment. We find that myeloid cells, the most abundant immune cell population in brain tumors, undergo the most prominent transcriptional adaptation in the TME. Resident microglia represent the main myeloid cell population in the cellular tumor, while peripheral-derived myeloid cells appear to infiltrate the brain at sites of blood-brain barrier disruption. We identify reactive dendritic cell-like gene expression programs associated with enhanced phagocytic and antigen-presentation features in GBM-educated microglia subsets that might be harnessed for novel immunotherapeutic approaches. Overall, PDOX models faithfully recapitulate the major components of the GBM-educated TME and allow assessment of phenotypic changes in the GBM ecosystem upon treatment.

Human stem cell-derived organoids have great potential for modelling physiological and pathological processes. They recapitulate in vitro the organization and function of a respective organ or part of an organ. Human midbrain organoids (hMOs) have been described to contain midbrain-specific dopaminergic neurons that release the neurotransmitter dopamine. However, the human midbrain contains also additional neuronal cell types, which are functionally interacting with each other. Here, we analysed hMOs at high-resolution by means of single-cell RNA sequencing (scRNA-seq), imaging and electrophysiology to unravel cell heterogeneity. Our findings demonstrate that hMOs show essential neuronal functional properties as spontaneous electrophysiological activity of different neuronal subtypes, including dopaminergic, GABAergic, glutamatergic and serotonergic neurons. Recapitulating these in vivo features makes hMOs an excellent tool for in vitro disease phenotyping and drug discovery. Electronic supplementary material The online version of this article (10.1007/s00441-020-03249-y) contains supplementary material, which is available to authorized users.

Objective: To pinpoint the earliest cellular defects underlying seizure onset (epileptogenic period) during perinatal brain development in a new zebrafish model of Dravet syndrome (DS) and to investigate potential disease-modifying activity of the 5HT2 receptor agonist fenfluramine. Methods: We used CRISPR/Cas9 mutagenesis to introduce a missense mutation, designed to perturb ion transport function in all channel isoforms, into scn1lab, the zebrafish orthologue of SCN1A (encoding voltage-gated sodium channel alpha subunit 1). We performed behavioral analysis and electroencephalographic recordings to measure convulsions and epileptiform discharges, followed by single-cell RNA-Seq, morphometric analysis of transgenic reporter-labeled γ-aminobutyric acidergic (GABAergic) neurons, and pharmacological profiling of mutant larvae. Results: Homozygous mutant (scn1labmut/mut ) larvae displayed spontaneous seizures with interictal, preictal, and ictal discharges (mean = 7.5 per 20-minute recording; P < .0001; one-way analysis of variance). Drop-Seq analysis revealed a 2:1 shift in the ratio of glutamatergic to GABAergic neurons in scn1labmut/mut larval brains versus wild type (WT), with dynamic changes in neuronal, glial, and progenitor cell populations. To explore disease pathophysiology further, we quantified dendritic arborization in GABAergic neurons and observed a 40% reduction in arbor number compared to WT (P < .001; n = 15 mutant, n = 16 WT). We postulate that the significant reduction in inhibitory arbors causes an inhibitory to excitatory neurotransmitter imbalance that contributes to seizures and enhanced electrical brain activity in scn1labmut/mut larvae (high-frequency range), with subsequent GABAergic neuronal loss and astrogliosis. Chronic fenfluramine administration completely restored dendritic arbor numbers to normal in scn1labmut/mut larvae, whereas similar treatment with the benzodiazepine diazepam attenuated seizures, but was ineffective in restoring neuronal cytoarchitecture. BrdU labeling revealed cell overproliferation in scn1labmut/mut larval brains that were rescued by fenfluramine but not diazepam. Significance: Our findings provide novel insights into early mechanisms of DS pathogenesis, describe dynamic cell population changes in the scn1labmut/mut brain, and present first-time evidence for potential disease modification by fenfluramine. Keywords: Dravet syndrome; epileptogenesis; fenfluramine; sodium channel; zebrafish. publishedVersion

Stress hormones are believed to skew the CD4 T-cell differentiation towards a Th2 response via a T-cell-extrinsic mechanism. Using isolated primary human naïve and memory CD4 T cells, here we show that both adrenergic- and glucocorticoid-mediated stress signalling pathways play a CD4 naïve T-cell-intrinsic role in regulating the Th1/Th2 differentiation balance. Both stress hormones reduced the Th1 programme and cytokine production by inhibiting mTORC1 signalling via two parallel mechanisms. Stress hormone signalling inhibited mTORC1 in naïve CD4 T cells (1) by affecting the PI3K/AKT pathway and (2) by regulating the expression of the circadian rhythm gene, period circadian regulator 1 (PER1). Both stress hormones induced the expression of PER1, which inhibited mTORC1 signalling, thus reducing Th1 differentiation. This previously unrecognized cell-autonomous mechanism connects stress hormone signalling with CD4 T-cell differentiation via mTORC1 and a specific circadian clock gene, namely PER1.

Increasing evidence suggests that neurodevelopmental alterations might contribute to increase the susceptibility to develop neurodegenerative diseases. We investigate the occurrence of developmental abnormalities in dopaminergic neurons in a model of Parkinson’s disease (PD). We monitor the differentiation of human patient-specific neuroepithelial stem cells (NESCs) into dopaminergic neurons. Using high-throughput image analyses and single-cell RNA sequencing, we observe that the PD-associated LRRK2-G2019S mutation alters the initial phase of neuronal differentiation by accelerating cell-cycle exit with a concomitant increase in cell death. We identify the NESC-specific core regulatory circuit and a molecular mechanism underlying the observed phenotypes. The expression of NR2F1, a key transcription factor involved in neurogenesis, decreases in LRRK2-G2019S NESCs, neurons, and midbrain organoids compared to controls. We also observe accelerated dopaminergic differentiation in vivo in NR2F1-deficient mouse embryos. This suggests a pathogenic mechanism involving the LRRK2-G2019S mutation, where the dynamics of dopaminergic differentiation are modified via NR2F1. publishedVersion

Human midbrain organoid models represent a new tool to study the underlaying etiology of Parkinson’s disease in physiological conditions.

While immunopathology has been widely studied in severe COVID-19 patients, immunoprotective factors in non-hospitalized patients have remained largely elusive. We systematically analyzed 484 peripheral immune cell signatures, various serological parameters and TCR repertoire in a longitudinal cohort of 63 mild and 15 hospitalized patients versus 14 asymptomatic and 26 control individuals. Within three days following PCR diagnosis, we observed coordinated responses of CD4 and CD8 T cells, various antigen presenting cells and antibody-secreting cells in mild, but not hospitalized COVID-19 patients. This early-stage SARS-CoV-2-specific response was predominantly characterized by substantially expanded clonotypes of CD4 and less of CD8 T cells. The early-stage responses of T cells and dendritic cells were highly predictive for later seroconversion and protective antibody levels after three weeks in mild non-hospitalized, but not in hospitalized patients. Our systemic analysis provides the first full picture and early-stage trajectory of highly coordinated immune responses in mild COVID-19 patients.

BackgroundCell types in ventral midbrain are involved in diseases with variable genetic susceptibility, such as Parkinson’s disease and schizophrenia. Many genetic variants affect regulatory regions and alter gene expression in a cell-type-specific manner depending on the chromatin structure and accessibility.ResultsWe report 20,658 single-nuclei chromatin accessibility profiles of ventral midbrain from two genetically and phenotypically distinct mouse strains. We distinguish ten cell types based on chromatin profiles and analysis of accessible regions controlling cell identity genes highlights cell-type-specific key transcription factors. Regulatory variation segregating the mouse strains manifests more on transcriptome than chromatin level. However, cell-type-level data reveals changes not captured at tissue level. To discover the scope and cell-type specificity ofcis-acting variation in midbrain gene expression, we identify putative regulatory variants and show them to be enriched at differentially expressed loci. Finally, we find TCF7L2 to mediatetrans-acting variation selectively in midbrain neurons.ConclusionsOur data set provides an extensive resource to study gene regulation in mesencephalon and provides insights into control of cell identity in the midbrain and identifies cell-type-specific regulatory variation possibly underlying phenotypic and behavioural differences between mouse strains.

Human brain organoid models that recapitulate the physiology and complexity of the human brain have a great potential for in vitro disease modeling, in particular for neurodegenerative diseases, such as Parkinson disease. In the present study, we compare single-cell RNA-sequencing data of human midbrain organoids to the developing human embryonic midbrain. We demonstrate that the in vitro model is comparable to its in vivo equivalents in terms of developmental path and cellular composition. Moreover, we investigate the potential of midbrain organoids for modeling early developmental changes in Parkinson disease. Therefore, we compare the single-cell RNA-sequencing data of healthy-individual-derived midbrain organoids to their isogenic LRRK2-p.Gly2019Ser-mutant counterparts. We show that the LRRK2 p.Gly2019Ser variant alters neurodevelopment, resulting in an untimely and incomplete differentiation with reduced cellular variability. Finally, we present four candidate genes, APP, DNAJC6, GATA3, and PTN, that might contribute to the LRRK2-p.Gly2019Ser-associated transcriptome changes that occur during early neurodevelopment.

While immunopathology has been widely studied in severe COVID-19 patients, immunoprotective factors in non-hospitalized patients have remained largely elusive. We systematically analyzed 484 peripheral immune cell signatures, various serological parameters and TCR repertoire in a longitudinal cohort of 63 mild and 15 hospitalized patients versus 14 asymptomatic and 26 control individuals. Within three days following PCR diagnosis, we observed coordinated responses of CD4 and CD8 T cells, various antigen presenting cells and antibody-secreting cells in mild, but not hospitalized COVID-19 patients. This early-stage SARS-CoV-2-specific response was predominantly characterized by substantially expanded clonotypes of CD4 and less of CD8 T cells. The early-stage responses of T cells and dendritic cells were highly predictive for later seroconversion and protective antibody levels after three weeks in mild non-hospitalized, but not in hospitalized patients. Our systemic analysis reveals the first full picture and early-stage trajectory of highly-coordinated immune responses in mild COVID-19 patients. Funding: The Predi-COVID study is supported by the Luxembourg National Research Fund (FNR) (Predi-COVID, 14716273) and the Andre Losch Fondation. We also highly appreciate the expert support of the IBBL processing and biorepository teams. F.Q.H. was partially supported by FNR CORE programme grant (CORE/14/BM/8231540/GeDES), FNR AFR- RIKEN bilateral programme (TregBAR, 11228353, F.Q.H. and M.O.) and PRIDE programme grants (PRIDE/11012546/NEXTIMMUNE and PRIDE/10907093/CRITICS). C.H. was partially supported by the FNR fast-track call COVID-19/2020-1/14703957/COV-Immun. Declaration of Interests: The authors declare that they have no conflict of interest. Ethics Approval Statement: All collections were performed with approval from relevant ethic organizations. Informed consent was obtained from each participant prior to collection. The blood sampling was performed by nurses from Clinical and Epidemiological Investigation Centre (CIEC) of LIH.

Parkinson’s disease (PD) is the second most prevalent neurodegenerative disorder, characterized by the loss of dopaminergic neurons (mDA) in the midbrain. The heterogenous pathology and complex underlying mechanisms are only partly understood and there is no treatment able to reverse PD progression. Here, we targeted the disease mechanisms by focusing on the ILE368ASN mutation within the PINK1 (PARK6) gene and systematically characterized midbrain dopaminergic neurons obtained from human induced pluripotent stem cells (iPSCs). Single-cell RNA sequencing (RNAseq) and pairwise analysis of gene expression identified genes consistently differentially expressed during the mDA neuron differentiation process. Subsequent network analysis revealed that these genes form a core network, which interacts with all known 19 protein-coding Parkinson’s disease-associated genes and includes ubiquitination, mitochondrial, protein processing, RNA metabolism, and secretory pathways as important subnetworks. Our findings indicate a unified network underlying PD pathology and offers new interpretation of the phenotypic heterogeneity of PD.

Cell types in ventral midbrain are involved in diseases with variable genetic susceptibility such as Parkinson’s disease and schizophrenia. Many genetic variants affect regulatory regions and alter gene expression. We report 20 658 single nuclei chromatin accessibility profiles of ventral midbrain from two genetically and phenotypically distinct mouse strains. We distinguish ten cell types based on chromatin profiles and analysis of accessible regions controlling cell identity genes highlights cell type-specific key transcription factors. Regulatory variation segregating the mouse strains manifests more on transcriptome than chromatin level. However, cell type-level data reveals changes not captured at tissue level. To discover the scope and cell-type specificity ofcis-acting variation in midbrain gene expression, we identify putative regulatory variants and show them to be enriched at differentially expressed loci. Finally, we find TCF7L2 to mediatetrans-acting variation selectively in midbrain neurons. Our dataset provides an extensive resource to study gene regulation in mesencephalon.

Patient-derived cancer models are essential tools for studying tumor biology and preclinical interventions. Here, we show that glioma patient-derived orthotopic xenografts (PDOXs) enable long-term propagation of patient tumors and represent clinically relevant patient avatars. We created a large collection of PDOXs from primary and recurrent gliomas with and without mutations in IDH1, which retained histopathological, genetic, epigenetic and transcriptomic features of patient tumors with no mouse-specific clonal evolution. Longitudinal PDOX models recapitulate the limited genetic evolution of gliomas observed in patient tumors following treatment. PDOX-derived standardized tumor organoid cultures enabled assessment of drug responses, which were validated in mice. PDOXs showed clinically relevant responses to Temozolomide and to targeted treatments such as EGFR and CDK4/6 inhibitors in (epi)genetically defined groups, according to MGMT promoter and EGFR/CDK status respectively. Dianhydrogalactitol, a bifunctional alkylating agent, showed promising potential against glioblastoma. Our study underlines the clinical relevance of glioma PDOX models for translational research and personalized treatment studies.

While immunopathology has been widely studied in patients with severe COVID-19, immune responses in non-hospitalized patients have remained largely elusive. We systematically analyze 484 peripheral cellular or soluble immune features in a longitudinal cohort of 63 mild and 15 hospitalized patients versus 14 asymptomatic and 26 household controls. We observe a transient increase of IP10/CXCL10 and interferon-β levels, coordinated responses of dominant SARS-CoV-2-specific CD4 and fewer CD8 T cells, and various antigen-presenting and antibody-secreting cells in mild patients within 3 days of PCR diagnosis. The frequency of key innate immune cells and their functional marker expression are impaired in hospitalized patients at day 1 of inclusion. T cell and dendritic cell responses at day 1 are highly predictive for SARS-CoV-2-specific antibody responses after 3 weeks in mild but not hospitalized patients. Our systematic analysis reveals a combinatorial picture and trajectory of various arms of the highly coordinated early-stage immune responses in mild COVID-19 patients.

Patient-based cancer models are essential tools for studying tumor biology and for the assessment of drug responses in a translational context. We report the establishment a large cohort of unique organoids and patient-derived orthotopic xenografts (PDOX) of various glioma subtypes, including gliomas with mutations in IDH1, and paired longitudinal PDOX from primary and recurrent tumors of the same patient. We show that glioma PDOXs enable long-term propagation of patient tumors and represent clinically relevant patient avatars that retain histopathological, genetic, epigenetic, and transcriptomic features of parental tumors. We find no evidence of mouse-specific clonal evolution in glioma PDOXs. Our cohort captures individual molecular genotypes for precision medicine including mutations in IDH1, ATRX, TP53, MDM2/4, amplification of EGFR, PDGFRA, MET, CDK4/6, MDM2/4, and deletion of CDKN2A/B, PTCH, and PTEN. Matched longitudinal PDOX recapitulate the limited genetic evolution of gliomas observed in patients following treatment. At the histological level, we observe increased vascularization in the rat host as compared to mice. PDOX-derived standardized glioma organoids are amenable to high-throughput drug screens that can be validated in mice. We show clinically relevant responses to temozolomide (TMZ) and to targeted treatments, such as EGFR and CDK4/6 inhibitors in (epi)genetically defined subgroups, according to MGMT promoter and EGFR/CDK status, respectively. Dianhydrogalactitol (VAL-083), a promising bifunctional alkylating agent in the current clinical trial, displayed high therapeutic efficacy, and was able to overcome TMZ resistance in glioblastoma. Our work underscores the clinical relevance of glioma organoids and PDOX models for translational research and personalized treatment studies and represents a unique publicly available resource for precision oncology.

Cell types in ventral midbrain are involved in diseases with variable genetic susceptibility such as Parkinson's disease and schizophrenia. Many genetic variants affect regulatory regions and alter gene expression. We report 20 658 single nuclei chromatin accessibility profiles of ventral midbrain from two genetically and phenotypically distinct mouse strains. We distinguish ten cell types based on chromatin profiles and analysis of accessible regions controlling cell identity genes highlights cell type-specific key transcription factors. Regulatory variation segregating the mouse strains manifests more on transcriptome than chromatin level. However, cell type-level data reveals changes not captured at tissue level. To discover the scope and cell-type specificity of cis-acting variation in midbrain gene expression, we identify putative regulatory variants and show them to be enriched at differentially expressed loci. Finally, we find TCF7L2 to mediate trans-acting variation selectively in midbrain neurons. Our dataset provides an extensive resource to study gene regulation in mesencephalon.

Summary Ca2+ is a central second messenger in eukaryotic cells that regulates many cellular processes. Recently, we have indicated that typical Ca2+ signals are not purely oscillatory as widely assumed, but exhibit stochastic spiking with cell type and pathway specific characteristics. Here, we present the Calcium Signaling Analyzer (CaSiAn), an open source software tool that allows for quantifying these signal characteristics including individual spike properties and time course statistics in a semi-automated manner. CaSiAn provides an intuitive graphical user interface allowing experimentalists to easily process a large amount of Ca2+ signals, interactively tune peak detection, revise statistical measures and access the quantified signal properties as excel or text files. Availability and implementation CaSiAn is implemented in Java and available on Github (https://github.com/mmahsa/CaSiAn) as well as on the project page (http://r3lab.uni.lu/web/casa). Supplementary information Supplementary data are available at Bioinformatics online.

Human stem cell-derived organoids have great potential for modelling physiological and pathological processes. They recapitulatein vitrothe organisation and function of a respective organ or part of an organ. Human midbrain organoids (hMOs) have been described to contain midbrain-specific dopaminergic neurons that release the neurotransmitter dopamine. However, the human midbrain contains also additional neuronal cell types, which are functionally interacting with each other. Here, we analysed hMOs at high-resolution by means of single-cell RNA-sequencing (scRNA-seq), imaging and electrophysiology to unravel cell heterogeneity. Our findings demonstrate that hMOs show essential neuronal functional properties as spontaneous electrophysiological activity of different neuronal subtypes, including dopaminergic, GABAergic, and glutamatergic neurons. Recapitulating thesein vivofeatures makes hMOs an excellent tool forin vitrodisease phenotyping and drug discovery.

Serine catabolism via the folate cycle provides formate that is essential for nucleotide synthesis in proliferating cells. In addition to this canonical function to support biomass production in anabolic cells, we have recently demonstrated in vitro and in vivo that formate production in cancer cells is often in excess of the anabolic demand. This excess formate production is characterized by formate overflow and thus, net formate excretion into the tumor microenvironment. Interestingly, we observe increased rates of formate overflow upon different chemical perturbations that induce growth arrest. Thus, stressed cancer cells that encounter growth restriction such as upon chemotherapy, are often characterized by increased formate release rates. We demonstrated that such high formate levels in the extracellular space promote invasion of glioblastoma cells. Using ex vivo brain slice cultures and an orthotopic brain tumor model, we demonstrate that silencing MTHFD1L, the essential enzyme to enable formate overflow, results in decreased invasiveness of the tumor. Embarking from this observation, we investigated the underlying mechanism and now provide evidence that the formate-dependent increase of cell motility is mediated by an activation of Ca2+ signaling. Activation of Ca2+ signaling triggers integrin and matrix metallopeptidase (MMP) responses enabling the invasion process. Targeting either the Ca2+ response or MMP release can suppress the formate dependent increase in invasion. Finally, we tested the effect of formate also in context of breast cancer where we were able to recapitulate our observation of increased invasiveness and, in this case, formate also promoted the metastatic potential. We conclude that excreted formate might serve as a cellular stress signal that represents a promotive trigger to support tumor escape mechanisms.

Intratumoral heterogeneity is one of the major barriers in the treatment of many tumor types. In glioblastoma (GBM) – the most aggressive brain tumor - survival of patients following standard therapy has stagnated to an average of 14 months, largely due to the inherent heterogeneity and phenotypic plasticity. Patient-derived orthotopic xenografts (PDOX) are important tools for preclinical cancer studies. We have shown that they recapitulate histological and molecular features of patient tumors, thus producing more clinically relevant study outcomes when compared to other preclinical models. To understand the impact of treatment on intratumoral heterogeneity and interactions of tumor cells with the microenvironment, we applied single-cell RNA-seq in GBM PDOX models generated from naive and treatment exposed patient tumors. Here we show that PDOX models recapitulate all the major cell types and transcriptional programs reported in GBM patient samples. We demonstrate that GBM tumor cells occur in distinct phenotypic states dictated by microenvironmental cues and treatment-induced pressure. Non-neoplastic cells present within xenografts undergo transcriptomic adaptation similar to the GBM microenvironment in patients. Analysis of treated PDOXs and longitudinal PDOX models derived from patients prior and after treatment allows us to investigate short-term and stable long-term transcriptomic adaptations. We aim to reveal key molecular regulators enhancing tumor cell plasticity and resistance, that could be targeted in combinatory treatments. Our GBM PDOX models provide a platform for high-quality pre-clinical research and an improved clinical relevance in the development and testing of novel therapeutics.