Alternative platelet differentiation pathways initiated by nonhierarchically related hematopoietic stem cells.

Rare multipotent stem cells replenish millions of blood cells per second through a time-consuming process, passing through multiple stages of increasingly lineage-restricted progenitors. Although insults to the blood-forming system highlight the need for more rapid blood replenishment from stem cells, established models of hematopoiesis implicate only one mandatory differentiation pathway for each blood cell lineage. Here, we establish a nonhierarchical relationship between distinct stem cells that replenish all blood cell lineages and stem cells that replenish almost exclusively platelets, a lineage essential for hemostasis and with important roles in both the innate and adaptive immune systems. These distinct stem cells use cellularly, molecularly and functionally separate pathways for the replenishment of molecularly distinct megakaryocyte-restricted progenitors: a slower steady-state multipotent pathway and a fast-track emergency-activated platelet-restricted pathway. These findings provide a framework for enhancing platelet replenishment in settings in which slow recovery of platelets remains a major clinical challenge.

Profiling spatiotemporal gene expression of the developing human spinal cord and implications for ependymoma origin.

The spatiotemporal regulation of cell fate specification in the human developing spinal cord remains largely unknown. In this study, by performing integrated analysis of single-cell and spatial multi-omics data, we used 16 prenatal human samples to create a comprehensive developmental cell atlas of the spinal cord during post-conceptional weeks 5-12. This revealed how the cell fate commitment of neural progenitor cells and their spatial positioning are spatiotemporally regulated by specific gene sets. We identified unique events in human spinal cord development relative to rodents, including earlier quiescence of active neural stem cells, differential regulation of cell differentiation and distinct spatiotemporal genetic regulation of cell fate choices. In addition, by integrating our atlas with pediatric ependymomas data, we identified specific molecular signatures and lineage-specific genes of cancer stem cells during progression. Thus, we delineate spatiotemporal genetic regulation of human spinal cord development and leverage these data to gain disease insight.

A roadmap for the Human Developmental Cell Atlas.

The Human Developmental Cell Atlas (HDCA) initiative, which is part of the Human Cell Atlas, aims to create a comprehensive reference map of cells during development. This will be critical to understanding normal organogenesis, the effect of mutations, environmental factors and infectious agents on human development, congenital and childhood disorders, and the cellular basis of ageing, cancer and regenerative medicine. Here we outline the HDCA initiative and the challenges of mapping and modelling human development using state-of-the-art technologies to create a reference atlas across gestation. Similar to the Human Genome Project, the HDCA will integrate the output from a growing community of scientists who are mapping human development into a unified atlas. We describe the early milestones that have been achieved and the use of human stem-cell-derived cultures, organoids and animal models to inform the HDCA, especially for prenatal tissues that are hard to acquire. Finally, we provide a roadmap towards a complete atlas of human development.

Effects of In Vivo Gluten Challenge on PBMC Gene Expression Profiles in Diet Treated Celiac Disease.

The pathological mechanisms that lead to the onset and reactivation of celiac disease (CD) remain largely unknown. While gluten free diet (GFD) improves the intestinal damage and associated clinical symptoms in majority of cases, it falls short of providing full recovery. Additionally, late or misdiagnosis is also common as CD presents with a wide range of symptoms. Clear understanding of CD pathogenesis is thus critical to address both diagnostic and treatment concerns. We aimed to study the molecular impact of short gluten exposure in GFD treated CD patients, as well as identify biological pathways that remain altered constitutively in CD regardless of treatment. Using RNAseq profiling of PBMC samples collected from treated CD patients and gluten challenged patient and healthy controls, we explored the peripheral transcriptome in CD patients following a short gluten exposure. Short gluten exposure of just three days was enough to alter the genome-wide PBMC transcriptome of patients. Pathway analysis revealed gluten-induced upregulation of mainly immune response related pathways, both innate and adaptive, in CD patients. We evaluated the perturbation of biological pathways in sample-specific manner. Compared to gluten exposed healthy controls, pathways related to tight junction, olfactory transduction, metabolism of unsaturated fatty acids (such as arachidonic acid), metabolism of amino acids (such as cysteine and glutamate), and microbial infection were constitutively altered in CD patients regardless of treatment, while GFD treatment appears to mostly normalize immune response pathways to "healthy" state. Upstream regulator prediction analysis using differentially expressed genes identified constitutively activated regulators relatively proximal to previously reported CD associated loci, particularly SMARCA4 on 19p13.2 and CSF2 on 5q31. We also found constitutively upregulated genes in CD that are in CD associated genetic loci such as MEF2BNB-MEF2B (BORCS8-MEF2B) on 19p13.11 and CSTB on 21q22.3. RNAseq revealed strong effects of short oral gluten challenge on whole PBMC fraction and constitutively altered pathways in CD PBMC suggesting important factors other than gluten in CD pathogenesis.

LifeTime and improving European healthcare through cell-based interceptive medicine.

Here we describe the LifeTime Initiative, which aims to track, understand and target human cells during the onset and progression of complex diseases, and to analyse their response to therapy at single-cell resolution. This mission will be implemented through the development, integration and application of single-cell multi-omics and imaging, artificial intelligence and patient-derived experimental disease models during the progression from health to disease. The analysis of large molecular and clinical datasets will identify molecular mechanisms, create predictive computational models of disease progression, and reveal new drug targets and therapies. The timely detection and interception of disease embedded in an ethical and patient-centred vision will be achieved through interactions across academia, hospitals, patient associations, health data management systems and industry. The application of this strategy to key medical challenges in cancer, neurological and neuropsychiatric disorders, and infectious, chronic inflammatory and cardiovascular diseases at the single-cell level will usher in cell-based interceptive medicine in Europe over the next decade.

A cell fitness selection model for neuronal survival during development.

Developmental cell death plays an important role in the construction of functional neural circuits. In vertebrates, the canonical view proposes a selection of the surviving neurons through stochastic competition for target-derived neurotrophic signals, implying an equal potential for neurons to compete. Here we show an alternative cell fitness selection of neurons that is defined by a specific neuronal heterogeneity code. Proprioceptive sensory neurons that will undergo cell death and those that will survive exhibit different molecular signatures that are regulated by retinoic acid and transcription factors, and are independent of the target and neurotrophins. These molecular features are genetically encoded, representing two distinct subgroups of neurons with contrasted functional maturation states and survival outcome. Thus, in this model, a heterogeneous code of intrinsic cell fitness in neighboring neurons provides differential competitive advantage resulting in the selection of cells with higher capacity to survive and functionally integrate into neural networks.

High-throughput chromatin accessibility profiling at single-cell resolution.

Here we develop a high-throughput single-cell ATAC-seq (assay for transposition of accessible chromatin) method to measure physical access to DNA in whole cells. Our approach integrates fluorescence imaging and addressable reagent deposition across a massively parallel (5184) nano-well array, yielding a nearly 20-fold improvement in throughput (up to ~1800 cells/chip, 4-5 h on-chip processing time) and library preparation cost (~81¢ per cell) compared to prior microfluidic implementations. We apply this method to measure regulatory variation in peripheral blood mononuclear cells (PBMCs) and show robust, de novo clustering of single cells by hematopoietic cell type.

Classes and continua of hippocampal CA1 inhibitory neurons revealed by single-cell transcriptomics.

Understanding any brain circuit will require a categorization of its constituent neurons. In hippocampal area CA1, at least 23 classes of GABAergic neuron have been proposed to date. However, this list may be incomplete; additionally, it is unclear whether discrete classes are sufficient to describe the diversity of cortical inhibitory neurons or whether continuous modes of variability are also required. We studied the transcriptomes of 3,663 CA1 inhibitory cells, revealing 10 major GABAergic groups that divided into 49 fine-scale clusters. All previously described and several novel cell classes were identified, with three previously described classes unexpectedly found to be identical. A division into discrete classes, however, was not sufficient to describe the diversity of these cells, as continuous variation also occurred between and within classes. Latent factor analysis revealed that a single continuous variable could predict the expression levels of several genes, which correlated similarly with it across multiple cell types. Analysis of the genes correlating with this variable suggested it reflects a range from metabolically highly active faster-spiking cells that proximally target pyramidal cells to slower-spiking cells targeting distal dendrites or interneurons. These results elucidate the complexity of inhibitory neurons in one of the simplest cortical structures and show that characterizing these cells requires continuous modes of variation as well as discrete cell classes.

Exome sequencing of primary breast cancers with paired metastatic lesions reveals metastasis-enriched mutations in the A-kinase anchoring protein family (AKAPs).

BACKGROUND: Tumor heterogeneity in breast cancer tumors is today widely recognized. Most of the available knowledge in genetic variation however, relates to the primary tumor while metastatic lesions are much less studied. Many studies have revealed marked alterations of standard prognostic and predictive factors during tumor progression. Characterization of paired primary- and metastatic tissues should therefore be fundamental in order to understand mechanisms of tumor progression, clonal relationship to tumor evolution as well as the therapeutic aspects of systemic disease. METHODS: We performed full exome sequencing of primary breast cancers and their metastases in a cohort of ten patients and further confirmed our findings in an additional cohort of 20 patients with paired primary and metastatic tumors. Furthermore, we used gene expression from the metastatic lesions and a primary breast cancer data set to study the gene expression of the AKAP gene family. RESULTS: We report that somatic mutations in A-kinase anchoring proteins are enriched in metastatic lesions. The frequency of mutation in the AKAP gene family was 10% in the primary tumors and 40% in metastatic lesions. Several copy number variations, including deletions in regions containing AKAP genes were detected and showed consistent patterns in both investigated cohorts. In a second cohort containing 20 patients with paired primary and metastatic lesions, AKAP mutations showed an increasing variant allele frequency after multiple relapses. Furthermore, gene expression profiles from the metastatic lesions (n = 120) revealed differential expression patterns of AKAPs relative to the tumor PAM50 intrinsic subtype, which were most apparent in the basal-like subtype. This pattern was confirmed in primary tumors from TCGA (n = 522) and in a third independent cohort (n = 182). CONCLUSION: Several studies from primary cancers have reported individual AKAP genes to be associated with cancer risk and metastatic relapses as well as direct involvement in cellular invasion and migration processes. Our findings reveal an enrichment of mutations in AKAP genes in metastatic breast cancers and suggest the involvement of AKAPs in the metastatic process. In addition, we report an AKAP gene expression pattern that consistently follows the tumor intrinsic subtype, further suggesting AKAP family members as relevant players in breast cancer biology.

Comparative cell cycle transcriptomics reveals synchronization of developmental transcription factor networks in cancer cells.

The cell cycle coordinates core functions such as replication and cell division. However, cell-cycle-regulated transcription in the control of non-core functions, such as cell identity maintenance through specific transcription factors (TFs) and signalling pathways remains unclear. Here, we provide a resource consisting of mapped transcriptomes in unsynchronized HeLa and U2OS cancer cells sorted for cell cycle phase by Fucci reporter expression. We developed a novel algorithm for data analysis that enables efficient visualization and data comparisons and identified cell cycle synchronization of Notch signalling and TFs associated with development. Furthermore, the cell cycle synchronizes with the circadian clock, providing a possible link between developmental transcriptional networks and the cell cycle. In conclusion we find that cell cycle synchronized transcriptional patterns are temporally compartmentalized and more complex than previously anticipated, involving genes, which control cell identity and development.

Alternative TSSs are co-regulated in single cells in the mouse brain.

Alternative transcription start sites (TSSs) have been extensively studied genome-wide for many cell types and have been shown to be important during development and to regulate transcript abundance between cell types. Likewise, single-cell gene expression has been extensively studied for many cell types. However, how single cells use TSSs has not yet been examined. In particular, it is unknown whether alternative TSSs are independently expressed, or whether they are co-activated or even mutually exclusive in single cells. Here, we use a previously published single-cell RNA-seq dataset, comprising thousands of cells, to study alternative TSS usage. We find that alternative TSS usage is a regulated process, and the correlation between two TSSs expressed in single cells of the same cell type is surprisingly high. Our findings indicate that TSSs generally are regulated by common factors rather than being independently regulated or stochastically expressed.

Membrane-Depolarizing Channel Blockers Induce Selective Glioma Cell Death by Impairing Nutrient Transport and Unfolded Protein/Amino Acid Responses.

Glioma-initiating cells (GIC) are considered the underlying cause of recurrences of aggressive glioblastomas, replenishing the tumor population and undermining the efficacy of conventional chemotherapy. Here we report the discovery that inhibiting T-type voltage-gated Ca2+ and KCa channels can effectively induce selective cell death of GIC and increase host survival in an orthotopic mouse model of human glioma. At present, the precise cellular pathways affected by the drugs affecting these channels are unknown. However, using cell-based assays and integrated proteomics, phosphoproteomics, and transcriptomics analyses, we identified the downstream signaling events these drugs affect. Changes in plasma membrane depolarization and elevated intracellular Na+, which compromised Na+-dependent nutrient transport, were documented. Deficits in nutrient deficit acted in turn to trigger the unfolded protein response and the amino acid response, leading ultimately to nutrient starvation and GIC cell death. Our results suggest new therapeutic targets to attack aggressive gliomas. Cancer Res; 77(7); 1741-52. ©2017 AACR.

Visceral motor neuron diversity delineates a cellular basis for nipple- and pilo-erection muscle control.

Despite the variety of physiological and target-related functions, little is known regarding the cellular complexity in the sympathetic ganglion. We explored the heterogeneity of mouse stellate and thoracic ganglia and found an unexpected variety of cell types. We identified specialized populations of nipple- and pilo-erector muscle neurons. These neurons extended axonal projections and were born among other neurons during embryogenesis, but remained unspecialized until target organogenesis occurred postnatally. Target innervation and cell-type specification was coordinated by an intricate acquisition of unique combinations of growth factor receptors and the initiation of expression of concomitant ligands by the nascent erector muscles. Overall, our results provide compelling evidence for a highly sophisticated organization of the sympathetic nervous system into discrete outflow channels that project to well-defined target tissues and offer mechanistic insight into how diversity and connectivity are established during development.

Visualization and analysis of gene expression in tissue sections by spatial transcriptomics.

Analysis of the pattern of proteins or messengerRNAs (mRNAs) in histological tissue sections is a cornerstone in biomedical research and diagnostics. This typically involves the visualization of a few proteins or expressed genes at a time. We have devised a strategy, which we call "spatial transcriptomics," that allows visualization and quantitative analysis of the transcriptome with spatial resolution in individual tissue sections. By positioning histological sections on arrayed reverse transcription primers with unique positional barcodes, we demonstrate high-quality RNA-sequencing data with maintained two-dimensional positional information from the mouse brain and human breast cancer. Spatial transcriptomics provides quantitative gene expression data and visualization of the distribution of mRNAs within tissue sections and enables novel types of bioinformatics analyses, valuable in research and diagnostics.

A PBX1 transcriptional network controls dopaminergic neuron development and is impaired in Parkinson's disease.

Pre-B-cell leukemia homeobox (PBX) transcription factors are known to regulate organogenesis, but their molecular targets and function in midbrain dopaminergic neurons (mDAn) as well as their role in neurodegenerative diseases are unknown. Here, we show that PBX1 controls a novel transcriptional network required for mDAn specification and survival, which is sufficient to generate mDAn from human stem cells. Mechanistically, PBX1 plays a dual role in transcription by directly repressing or activating genes, such as Onecut2 to inhibit lateral fates during embryogenesis, Pitx3 to promote mDAn development, and Nfe2l1 to protect from oxidative stress. Notably, PBX1 and NFE2L1 levels are severely reduced in dopaminergic neurons of the substantia nigra of Parkinson's disease (PD) patients and decreased NFE2L1 levels increases damage by oxidative stress in human midbrain cells. Thus, our results reveal novel roles for PBX1 and its transcriptional network in mDAn development and PD, opening the door for new therapeutic interventions.

Origin, fate and dynamics of macrophages at central nervous system interfaces.

Perivascular, subdural meningeal and choroid plexus macrophages are non-parenchymal macrophages that mediate immune responses at brain boundaries. Although the origin of parenchymal microglia has recently been elucidated, much less is known about the precursors, the underlying transcriptional program and the dynamics of the other macrophages in the central nervous system (CNS). It was assumed that they have a high turnover from blood-borne monocytes. However, using parabiosis and fate-mapping approaches in mice, we found that CNS macrophages arose from hematopoietic precursors during embryonic development and established stable populations, with the notable exception of choroid plexus macrophages, which had dual origins and a shorter life span. The generation of CNS macrophages relied on the transcription factor PU.1, whereas the MYB, BATF3 and NR4A1 transcription factors were not required.

Rewired metabolism in drug-resistant leukemia cells: a metabolic switch hallmarked by reduced dependence on exogenous glutamine.

Cancer cells that escape induction therapy are a major cause of relapse. Understanding metabolic alterations associated with drug resistance opens up unexplored opportunities for the development of new therapeutic strategies. Here, we applied a broad spectrum of technologies including RNA sequencing, global untargeted metabolomics, and stable isotope labeling mass spectrometry to identify metabolic changes in P-glycoprotein overexpressing T-cell acute lymphoblastic leukemia (ALL) cells, which escaped a therapeutically relevant daunorubicin treatment. We show that compared with sensitive ALL cells, resistant leukemia cells possess a fundamentally rewired central metabolism characterized by reduced dependence on glutamine despite a lack of expression of glutamate-ammonia ligase (GLUL), a higher demand for glucose and an altered rate of fatty acid ß-oxidation, accompanied by a decreased pantothenic acid uptake capacity. We experimentally validate our findings by selectively targeting components of this metabolic switch, using approved drugs and starvation approaches followed by cell viability analyses in both the ALL cells and in an acute myeloid leukemia (AML) sensitive/resistant cell line pair. We demonstrate how comparative metabolomics and RNA expression profiling of drug-sensitive and -resistant cells expose targetable metabolic changes and potential resistance markers. Our results show that drug resistance is associated with significant metabolic costs in cancer cells, which could be exploited using new therapeutic strategies.

Amplification-free sequencing of cell-free DNA for prenatal non-invasive diagnosis of chromosomal aberrations.

Cell-free DNA has been used for fetal rhesus factor and sex determination, fetal aneuploidy screening, cancer diagnostics and monitoring, and other applications. However current methods of using cell free DNA require amplification, which leads to allelic dropout and bias especially when starting with small amounts of DNA. Here we describe an amplification-free method for sequencing of cell-free DNA, even from low levels of starting material. We evaluated this method in the context of prenatal diagnosis of fetal aneuploidy and compared it with a PCR-based library preparation method as well as a recently described method using unique molecular identifiers (UMI). All methods performed well, however coverage was increased by the amplification-free method and GC-induced bias was reduced by both the amplification-free method and the UMI method. Future diagnostic applications including whole genome sequencing of cell-free DNA will benefit from amplification-free sequencing.

Selective calcium sensitivity in immature glioma cancer stem cells.

Tumor-initiating cells are a subpopulation in aggressive cancers that exhibit traits shared with stem cells, including the ability to self-renew and differentiate, commonly referred to as stemness. In addition, such cells are resistant to chemo- and radiation therapy posing a therapeutic challenge. To uncover stemness-associated functions in glioma-initiating cells (GICs), transcriptome profiles were compared to neural stem cells (NSCs) and gene ontology analysis identified an enrichment of Ca2+ signaling genes in NSCs and the more stem-like (NSC-proximal) GICs. Functional analysis in a set of different GIC lines regarding sensitivity to disturbed homeostasis using A23187 and Thapsigargin, revealed that NSC-proximal GICs were more sensitive, corroborating the transcriptome data. Furthermore, Ca2+ drug sensitivity was reduced in GICs after differentiation, with most potent effect in the NSC-proximal GIC, supporting a stemness-associated Ca2+ sensitivity. NSCs and the NSC-proximal GIC line expressed a larger number of ion channels permeable to potassium, sodium and Ca2+. Conversely, a higher number of and higher expression levels of Ca2+ binding genes that may buffer Ca2+, were expressed in NSC-distal GICs. In particular, expression of the AMPA glutamate receptor subunit GRIA1, was found to associate with Ca2+ sensitive NSC-proximal GICs, and decreased as GICs differentiated along with reduced Ca2+ drug sensitivity. The correlation between high expression of Ca2+ channels (such as GRIA1) and sensitivity to Ca2+ drugs was confirmed in an additional nine novel GIC lines. Calcium drug sensitivity also correlated with expression of the NSC markers nestin (NES) and FABP7 (BLBP, brain lipid-binding protein) in this extended analysis. In summary, NSC-associated NES+/FABP7+/GRIA1+ GICs were selectively sensitive to disturbances in Ca2+ homeostasis, providing a potential target mechanism for eradication of an immature population of malignant cells.

Myelodysplastic syndromes are propagated by rare and distinct human cancer stem cells in vivo.

Evidence for distinct human cancer stem cells (CSCs) remains contentious and the degree to which different cancer cells contribute to propagating malignancies in patients remains unexplored. In low- to intermediate-risk myelodysplastic syndromes (MDS), we establish the existence of rare multipotent MDS stem cells (MDS-SCs), and their hierarchical relationship to lineage-restricted MDS progenitors. All identified somatically acquired genetic lesions were backtracked to distinct MDS-SCs, establishing their distinct MDS-propagating function in vivo. In isolated del(5q)-MDS, acquisition of del(5q) preceded diverse recurrent driver mutations. Sequential analysis in del(5q)-MDS revealed genetic evolution in MDS-SCs and MDS-progenitors prior to leukemic transformation. These findings provide definitive evidence for rare human MDS-SCs in vivo, with extensive implications for the targeting of the cells required and sufficient for MDS-propagation.