Search

Identifying initial triggering events in neurodegenerative disorders is critical to developing preventive therapies. In Huntington's disease (HD), hyperdopaminergia-probably triggered by the dysfunction of the most affected neurons, indirect pathway spiny projection neurons (iSPNs)-is believed to induce hyperkinesia, an early stage HD symptom. However, how this change arises and contributes to HD pathogenesis is unclear. Here, we demonstrate that genetic disruption of iSPNs function by Ntrk2/Trkb deletion in mice results in increased striatal dopamine and midbrain dopaminergic neurons, preceding hyperkinetic dysfunction. Transcriptomic analysis of iSPNs at the pre-symptomatic stage showed de-regulation of metabolic pathways, including upregulation of Gsto2, encoding glutathione S-transferase omega-2 (GSTO2). Selectively reducing Gsto2 in iSPNs in vivo effectively prevented dopaminergic dysfunction and halted the onset and progression of hyperkinetic symptoms. This study uncovers a functional link between altered iSPN BDNF-TrkB signalling, glutathione-ascorbate metabolism and hyperdopaminergic state, underscoring the vital role of GSTO2 in maintaining dopamine balance., Competing Interests: Competing interests: The authors declare no competing interests., (© 2024. The Author(s).)

Ebola virus (EBOV) causes Ebola virus disease (EVD), marked by severe hemorrhagic fever; however, the mechanisms underlying the disease remain unclear. To assess the molecular basis of EVD across time, we performed RNA sequencing on 17 tissues from a natural history study of 21 rhesus monkeys, developing new methods to characterize host-pathogen dynamics. We identified alterations in host gene expression with previously unknown tissue-specific changes, including downregulation of genes related to tissue connectivity. EBOV was widely disseminated throughout the body; using a new, broadly applicable deconvolution method, we found that viral load correlated with increased monocyte presence. Patterns of viral variation between tissues differentiated primary infections from compartmentalized infections, and several variants impacted viral fitness in a EBOV/Kikwit minigenome system, suggesting that functionally significant variants can emerge during early infection. This comprehensive portrait of host-pathogen dynamics in EVD illuminates new features of pathogenesis and establishes resources to study other emerging pathogens., Competing Interests: P.C.S. is a co-founder and shareholder of Sherlock Biosciences and Delve Bio, a board member and shareholder of Danaher Corporation, and has filed IP related to genomic sequencing and diagnostic technologies. A.K.S. reports compensation for consulting and/or scientific advisory board (SAB) membership from Merck, Honeycomb Biotechnologies, Cellarity, Repertoire Immune Medicines, Ochre Bio, Third Rock Ventures, Hovione, Relation Therapeutics, FL82, FL86, Empress Therapeutics, IntrECate Biotherapeutics, Senda Biosciences, and Dahlia Biosciences unrelated to this work. F.Z. is a scientific advisor and cofounder of Editas Medicine, Beam Therapeutics, Pairwise Plants, Arbor Biotechnologies, and Aera Therapeutics. F.Z. is a scientific advisor for Octant., (© 2023 The Authors.)

Hematopoietic stem cells (HSCs) residing in specialized niches in the bone marrow are responsible for the balanced output of multiple short-lived blood cell lineages in steady-state and in response to different challenges. However, feedback mechanisms by which HSCs, through their niches, sense acute losses of specific blood cell lineages remain to be established. While all HSCs replenish platelets, previous studies have shown that a large fraction of HSCs are molecularly primed for the megakaryocyte-platelet lineage and are rapidly recruited into proliferation upon platelet depletion. Platelets normally turnover in an activation-dependent manner, herein mimicked by antibodies inducing platelet activation and depletion. Antibody-mediated platelet activation upregulates expression of Interleukin-1 (IL-1) in platelets, and in bone marrow extracellular fluid in vivo. Genetic experiments demonstrate that rather than IL-1 directly activating HSCs, activation of bone marrow Lepr + perivascular niche cells expressing IL-1 receptor is critical for the optimal activation of quiescent HSCs upon platelet activation and depletion. These findings identify a feedback mechanism by which activation-induced depletion of a mature blood cell lineage leads to a niche-dependent activation of HSCs to reinstate its homeostasis., (© 2023. The Author(s).)

SARS-CoV-2 distribution and circulation dynamics are not well understood due to challenges in assessing genomic data from tissue samples. We develop experimental and computational workflows for high-depth viral sequencing and high-resolution genomic analyses from formalin-fixed, paraffin-embedded tissues and apply them to 120 specimens from six subjects with fatal COVID-19. To varying degrees, viral RNA is present in extrapulmonary tissues from all subjects. The majority of the 180 viral variants identified within subjects are unique to individual tissue samples. We find more high-frequency (>10%) minor variants in subjects with a longer disease course, with one subject harboring ten such variants, exclusively in extrapulmonary tissues. One tissue-specific high-frequency variant was a nonsynonymous mutation in the furin-cleavage site of the spike protein. Our findings suggest adaptation and/or compartmentalized infection, illuminating the basis of extrapulmonary COVID-19 symptoms and potential for viral reservoirs, and have broad utility for investigating human pathogens., (© 2023. The Author(s).)

Repeated emergence of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) variants with increased fitness underscores the value of rapid detection and characterization of new lineages. We have developed PyR 0 , a hierarchical Bayesian multinomial logistic regression model that infers relative prevalence of all viral lineages across geographic regions, detects lineages increasing in prevalence, and identifies mutations relevant to fitness. Applying PyR 0 to all publicly available SARS-CoV-2 genomes, we identify numerous substitutions that increase fitness, including previously identified spike mutations and many nonspike mutations within the nucleocapsid and nonstructural proteins. PyR 0 forecasts growth of new lineages from their mutational profile, ranks the fitness of lineages as new sequences become available, and prioritizes mutations of biological and public health concern for functional characterization.

An outbreak of over 1,000 COVID-19 cases in Provincetown, Massachusetts (MA), in July 2021-the first large outbreak mostly in vaccinated individuals in the US-prompted a comprehensive public health response, motivating changes to national masking recommendations and raising questions about infection and transmission among vaccinated individuals. To address these questions, we combined viral genomic and epidemiological data from 467 individuals, including 40% of outbreak-associated cases. The Delta variant accounted for 99% of cases in this dataset; it was introduced from at least 40 sources, but 83% of cases derived from a single source, likely through transmission across multiple settings over a short time rather than a single event. Genomic and epidemiological data supported multiple transmissions of Delta from and between fully vaccinated individuals. However, despite its magnitude, the outbreak had limited onward impact in MA and the US overall, likely due to high vaccination rates and a robust public health response., Competing Interests: Declaration of interests P.C.S. is a co-founder of, shareholder in, and scientific advisor to Sherlock Biosciences, as well as a board member of and shareholder in Danaher Corporation. J.E.L. has received consulting fees from Sherlock Biosciences. A.P. is a venture partner at Google Ventures. W.P.H. is a member of the scientific advisory board of Biobot. Other authors report no competing interests., (Copyright © 2021 The Authors. Published by Elsevier Inc. All rights reserved.)

Hematopoietic stem cells (HSCs) must ensure adequate blood cell production following distinct external stressors. A comprehensive understanding of in vivo heterogeneity and specificity of HSC responses to external stimuli is currently lacking. We performed single-cell RNA sequencing (scRNA-Seq) on functionally validated mouse HSCs and LSK (Lin-, c-Kit+, Sca1+) progenitors after in vivo pharmacological perturbation of niche signals interferon, granulocyte colony-stimulating factor (G-CSF), and prostaglandin. We identified six HSC states that are characterized by enrichment but not exclusive expression of marker genes. External signals induced rapid transitions between HSC states but transcriptional response varied both between external stimulants and within the HSC population for a given perturbation. In contrast to LSK progenitors, HSCs were characterized by a greater link between molecular signatures at baseline and in response to external stressors. Chromatin analysis of unperturbed HSCs and LSKs by scATAC-Seq suggested some HSC-specific, cell intrinsic predispositions to niche signals. We compiled a comprehensive resource of HSC- and LSK progenitor-specific chromatin and transcriptional features that represent determinants of signal receptiveness and regenerative potential during stress hematopoiesis., Competing Interests: EF, AS, MM, ER, SY, YZ, JG, NB, NB, FC No competing interests declared, DS is a director and equity holder of Agios Pharmaceuticals, Magenta Therapeutics, Editas Medicines, ClearCreekBio, and Life-VaultBio; a founder of Fate Therapeutics and Magenta Therapeutics; and a consultant to FOG Pharma and VCanBio, LZ is founder and stockholder of Fate, Inc, Scholar Rock, Camp4 therapeutics and a scientific advisor for Stemgent, (© 2021, Fast et al.)

Bone metastases are devastating complications of cancer. They are particularly common in prostate cancer (PCa), represent incurable disease, and are refractory to immunotherapy. We seek to define distinct features of the bone marrow (BM) microenvironment by analyzing single cells from bone metastatic prostate tumors, involved BM, uninvolved BM, and BM from cancer-free, orthopedic patients, and healthy individuals. Metastatic PCa is associated with multifaceted immune distortion, specifically exhaustion of distinct T cell subsets, appearance of macrophages with states specific to PCa bone metastases. The chemokine CCL20 is notably overexpressed by myeloid cells, as is its cognate CCR6 receptor on T cells. Disruption of the CCL20-CCR6 axis in mice with syngeneic PCa bone metastases restores T cell reactivity and significantly prolongs animal survival. Comparative high-resolution analysis of PCa bone metastases shows a targeted approach for relieving local immunosuppression for therapeutic effect., Competing Interests: Declaration of interests P.V.K. serves on the Scientific Advisory Board to Celsius Therapeutics Inc. and Biomage Inc. D.T.S. is a director and shareholder for Agios Therapeutics and Editas Medicines; a founder, director, shareholder, and scientific advisory board member for Magenta Therapeutics and LifeVault Bio, a shareholder and founder of Fate Therapeutics, and a director, founder, and shareholder for Clear Creek Bio, a consultant for FOG Pharma and VCanBio, and a recipient of sponsored research funding from Novartis. D.B.S. is a founder, consultant, and shareholder for Clear Creek Bio., (Copyright © 2021 Elsevier Inc. All rights reserved.)

Multiple summer events, including large indoor gatherings, in Provincetown, Massachusetts (MA), in July 2021 contributed to an outbreak of over one thousand COVID-19 cases among residents and visitors. Most cases were fully vaccinated, many of whom were also symptomatic, prompting a comprehensive public health response, motivating changes to national masking recommendations, and raising questions about infection and transmission among vaccinated individuals. To characterize the outbreak and the viral population underlying it, we combined genomic and epidemiological data from 467 individuals, including 40% of known outbreak-associated cases. The Delta variant accounted for 99% of sequenced outbreak-associated cases. Phylogenetic analysis suggests over 40 sources of Delta in the dataset, with one responsible for a single cluster containing 83% of outbreak-associated genomes. This cluster was likely not the result of extensive spread at a single site, but rather transmission from a common source across multiple settings over a short time. Genomic and epidemiological data combined provide strong support for 25 transmission events from, including many between, fully vaccinated individuals; genomic data alone provides evidence for an additional 64. Together, genomic epidemiology provides a high-resolution picture of the Provincetown outbreak, revealing multiple cases of transmission of Delta from fully vaccinated individuals. However, despite its magnitude, the outbreak was restricted in its onward impact in MA and the US, likely due to high vaccination rates and a robust public health response.

COVID-19, which is caused by SARS-CoV-2, can result in acute respiratory distress syndrome and multiple organ failure 1-4 , but little is known about its pathophysiology. Here we generated single-cell atlases of 24 lung, 16 kidney, 16 liver and 19 heart autopsy tissue samples and spatial atlases of 14 lung samples from donors who died of COVID-19. Integrated computational analysis uncovered substantial remodelling in the lung epithelial, immune and stromal compartments, with evidence of multiple paths of failed tissue regeneration, including defective alveolar type 2 differentiation and expansion of fibroblasts and putative TP63 + intrapulmonary basal-like progenitor cells. Viral RNAs were enriched in mononuclear phagocytic and endothelial lung cells, which induced specific host programs. Spatial analysis in lung distinguished inflammatory host responses in lung regions with and without viral RNA. Analysis of the other tissue atlases showed transcriptional alterations in multiple cell types in heart tissue from donors with COVID-19, and mapped cell types and genes implicated with disease severity based on COVID-19 genome-wide association studies. Our foundational dataset elucidates the biological effect of severe SARS-CoV-2 infection across the body, a key step towards new treatments.

Dopa decarboxylase (DDC) synthesizes serotonin in the developing mouse heart where it is encoded by Ddc_exon1a , a tissue-specific paternally expressed imprinted gene. Ddc_exon1a shares an imprinting control region (ICR) with the imprinted, maternally expressed (outside of the central nervous system) Grb10 gene on mouse chromosome 11, but little else is known about the tissue-specific imprinted expression of Ddc_exon1a . Fluorescent immunostaining localizes DDC to the developing myocardium in the pre-natal mouse heart, in a region susceptible to abnormal development and implicated in congenital heart defects in human. Ddc_exon1a and Grb10 are not co-expressed in heart nor in brain where Grb10 is also paternally expressed, despite sharing an ICR, indicating they are mechanistically linked by their shared ICR but not by Grb10 gene expression. Evidence from a Ddc_exon1a gene knockout mouse model suggests that it mediates the growth of the developing myocardium and a thinning of the myocardium is observed in a small number of mutant mice examined, with changes in gene expression detected by microarray analysis. Comparative studies in the human developing heart reveal a paternal expression bias with polymorphic imprinting patterns between individual human hearts at DDC_EXON1a , a finding consistent with other imprinted genes in human., Competing Interests: The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (Copyright © 2021 Prickett, Montibus, Barkas, Amante, Franco, Cowley, Puszyk, Shannon, Irving, Madon-Simon, Ward, Schulz, Baldwin and Oakey.)

Distinct types of dorsal root ganglion sensory neurons may have unique contributions to chronic pain. Identification of primate sensory neuron types is critical for understanding the cellular origin and heritability of chronic pain. However, molecular insights into the primate sensory neurons are missing. Here we classify non-human primate dorsal root ganglion sensory neurons based on their transcriptome and map human pain heritability to neuronal types. First, we identified cell correlates between two major datasets for mouse sensory neuron types. Machine learning exposes an overall cross-species conservation of somatosensory neurons between primate and mouse, although with differences at individual gene level, highlighting the importance of primate data for clinical translation. We map genomic loci associated with chronic pain in human onto primate sensory neuron types to identify the cellular origin of chronic pain. Genome-wide associations for chronic pain converge on two different neuronal types distributed between pain disorders that display different genetic susceptibilities, suggesting both unique and shared mechanisms between different pain conditions.

The SARS-CoV-2 pandemic has caused over 1 million deaths globally, mostly due to acute lung injury and acute respiratory distress syndrome, or direct complications resulting in multiple-organ failures. Little is known about the host tissue immune and cellular responses associated with COVID-19 infection, symptoms, and lethality. To address this, we collected tissues from 11 organs during the clinical autopsy of 17 individuals who succumbed to COVID-19, resulting in a tissue bank of approximately 420 specimens. We generated comprehensive cellular maps capturing COVID-19 biology related to patients' demise through single-cell and single-nucleus RNA-Seq of lung, kidney, liver and heart tissues, and further contextualized our findings through spatial RNA profiling of distinct lung regions. We developed a computational framework that incorporates removal of ambient RNA and automated cell type annotation to facilitate comparison with other healthy and diseased tissue atlases. In the lung, we uncovered significantly altered transcriptional programs within the epithelial, immune, and stromal compartments and cell intrinsic changes in multiple cell types relative to lung tissue from healthy controls. We observed evidence of: alveolar type 2 (AT2) differentiation replacing depleted alveolar type 1 (AT1) lung epithelial cells, as previously seen in fibrosis; a concomitant increase in myofibroblasts reflective of defective tissue repair; and, putative TP63 + intrapulmonary basal-like progenitor (IPBLP) cells, similar to cells identified in H1N1 influenza, that may serve as an emergency cellular reserve for severely damaged alveoli. Together, these findings suggest the activation and failure of multiple avenues for regeneration of the epithelium in these terminal lungs. SARS-CoV-2 RNA reads were enriched in lung mononuclear phagocytic cells and endothelial cells, and these cells expressed distinct host response transcriptional programs. We corroborated the compositional and transcriptional changes in lung tissue through spatial analysis of RNA profiles in situ and distinguished unique tissue host responses between regions with and without viral RNA, and in COVID-19 donor tissues relative to healthy lung. Finally, we analyzed genetic regions implicated in COVID-19 GWAS with transcriptomic data to implicate specific cell types and genes associated with disease severity. Overall, our COVID-19 cell atlas is a foundational dataset to better understand the biological impact of SARS-CoV-2 infection across the human body and empowers the identification of new therapeutic interventions and prevention strategies., Competing Interests: Competing Interests P.D., R.F., E.M.M., M.R., E.H.R., L.P., T.He., J.R., J.B., and S.W. are employees and stockholders at Nanostring Technologies Inc. D.Z., is a former employee and stockholder at NanoString Technologies. N.H., holds equity in BioNTech and Related Sciences. T.H.is an employee and stockholder of Prime Medicine as of Oct. 13, 2020. G.H. is an employee of Genentech as of Nov 16, 2020. R.N. is a founder, shareholder, and member of the board at Rhinostics Inc. A.R. is a co-founder and equity holder of Celsius Therapeutics, an equity holder in Immunitas, and was an SAB member of ThermoFisher Scientific, Syros Pharmaceuticals, Neogene Therapeutics and Asimov until July 31, 2020. From August 1, 2020, A.R. is an employee of Genentech. From October 19, 2020, O.R.-R is an employee of Genentech. P.C.S is a co-founder and shareholder of Sherlock Biosciences, and a Board member and shareholder of Danaher Corporation. A.K.S. reports compensation for consulting and/or SAB membership from Honeycomb Biotechnologies, Cellarity, Repertoire Immune Medicines, Ochre Bio, and Dahlia Biosciences. Z.G.J. reports grant support from Gilead Science, Pfizer, compensation for consulting from Olix Pharmaceuticals. Y.V.P. reports grant support from Enanta Pharmaceuticals, CymaBay Therapeutics, Morphic Therapeutic; consulting and/or SAB in Ambys Medicines, Morphic Therapeutics, Enveda Therapeutics, BridgeBio Pharma, as well as being an Editor – American Journal of Physiology-Gastrointestinal and Liver Physiology. GS reports consultant service in Alnylam Pharmaceuticals, Merck, Generon, Glympse Bio, Inc., Mayday Foundation, Novartis Pharmaceuticals, Quest Diagnostics, Surrozen, Terra Firma, Zomagen Bioscience, Pandion Therapeutics, Inc. Durect Corporation; royalty from UpToDate Inc., and Editor service in Hepatology Communications. P.R.T. receives consulting fees from Cellarity Inc., and Surrozen Inc., for work not related to this manuscript.

Leukemic relapse remains a major barrier to successful allogeneic hematopoietic stem cell transplantation (allo-HSCT) for aggressive hematologic malignancies. The basis for relapse of advanced lymphoid malignancies remains incompletely understood and may involve escape from the graft-versus-leukemia (GvL) effect. We hypothesized that for patients with chronic lymphocytic leukemia (CLL) treated with allo-HSCT, leukemic cell-intrinsic features influence transplant outcomes by directing the evolutionary trajectories of CLL cells. Integrated genetic, transcriptomic, and epigenetic analyses of CLL cells from 10 patients revealed that the clinical kinetics of post-HSCT relapse are shaped by distinct molecular dynamics. Early relapses after allo-HSCT exhibited notable genetic stability; single CLL cell transcriptional analysis demonstrated a cellular heterogeneity that was static over time. In contrast, CLL cells relapsing late after allo-HSCT displayed notable genetic evolution and evidence of neoantigen depletion, consistent with marked single-cell transcriptional shifts that were unique to each patient. We observed a greater rate of epigenetic change for late relapses not seen in early relapses or relapses after chemotherapy alone, suggesting that the selection pressures of the GvL bottleneck are unlike those imposed by chemotherapy. No selective advantage for human leukocyte antigen (HLA) loss was observed, even when present in pretransplant subpopulations. Gain of stem cell modules was a common signature associated with leukemia relapse regardless of posttransplant relapse kinetics. These data elucidate the biological pathways that underlie GvL resistance and posttransplant relapse., (Copyright © 2020 The Authors, some rights reserved; exclusive licensee American Association for the Advancement of Science. No claim to original U.S. Government Works.)

Alternative splicing (AS) and alternative polyadenylation (APA) generate diverse transcripts in mammalian genomes during development and differentiation. Epigenetic marks such as trimethylation of histone H3 lysine 36 (H3K36me3) and DNA methylation play a role in generating transcriptome diversity. Intragenic CpG islands (iCGIs) and their corresponding host genes exhibit dynamic epigenetic and gene expression patterns during development and between different tissues. We hypothesise that iCGI-associated H3K36me3, DNA methylation and transcription can influence host gene AS and/or APA. We investigate H3K36me3 and find that this histone mark is not a major regulator of AS or APA in our model system. Genomewide, we identify over 4000 host genes that harbour an iCGI in the mammalian genome, including both previously annotated and novel iCGI/host gene pairs. The transcriptional activity of these iCGIs is tissue- and developmental stage-specific and, for the first time, we demonstrate that the premature termination of host gene transcripts upstream of iCGIs is closely correlated with the level of iCGI transcription in a DNA-methylation independent manner. These studies suggest that iCGI transcription, rather than H3K36me3 or DNA methylation, interfere with host gene transcription and pre-mRNA processing genomewide and contributes to the spatiotemporal diversification of both the transcriptome and proteome., (© The Author(s) 2020. Published by Oxford University Press on behalf of Nucleic Acids Research.)

Single-cell RNA sequencing is often applied in study designs that include multiple individuals, conditions or tissues. To identify recurrent cell subpopulations in such heterogeneous collections, we developed Conos, an approach that relies on multiple plausible inter-sample mappings to construct a global graph connecting all measured cells. The graph enables identification of recurrent cell clusters and propagation of information between datasets in multi-sample or atlas-scale collections.

The endocardium interacts with the myocardium to promote proliferation and morphogenesis during the later stages of heart development. However, the role of the endocardium in early cardiac ontogeny remains under-explored. Given the shared origin, subsequent juxtaposition, and essential cell-cell interactions of endocardial and myocardial cells throughout heart development, we hypothesized that paracrine signaling from the endocardium to the myocardium is crucial for initiating early differentiation of myocardial cells. To test this, we generated an in vitro , endocardial-specific ablation model using the diphtheria toxin receptor under the regulatory elements of the Nfat c1 genomic locus ( NFATc1-DTR ). Early treatment of NFATc1-DTR mouse embryoid bodies with diphtheria toxin efficiently ablated endocardial cells, which significantly attenuated the percentage of beating EBs in culture and expression of early and late myocardial differentiation markers. The addition of Bmp2 during endocardial ablation partially rescued myocyte differentiation, maturation and function. Therefore, we conclude that early stages of myocardial differentiation rely on endocardial paracrine signaling mediated in part by Bmp2. Our findings provide novel insight into early endocardial-myocardial interactions that can be explored to promote early myocardial development and growth., Competing Interests: Competing interestsThe authors declare no competing or financial interests., (© 2019. Published by The Company of Biologists Ltd.)

Single-cell RNA sequencing (scRNA-seq) has emerged as a powerful tool for resolving transcriptional heterogeneity. However, its application to studying cancerous tissues is currently hampered by the lack of coverage across key mutation hotspots in the vast majority of cells; this lack of coverage prevents the correlation of genetic and transcriptional readouts from the same single cell. To overcome this, we developed TARGET-seq, a method for the high-sensitivity detection of multiple mutations within single cells from both genomic and coding DNA, in parallel with unbiased whole-transcriptome analysis. Applying TARGET-seq to 4,559 single cells, we demonstrate how this technique uniquely resolves transcriptional and genetic tumor heterogeneity in myeloproliferative neoplasms (MPN) stem and progenitor cells, providing insights into deregulated pathways of mutant and non-mutant cells. TARGET-seq is a powerful tool for resolving the molecular signatures of genetically distinct subclones of cancer cells., (Copyright © 2019 The Author(s). Published by Elsevier Inc. All rights reserved.)

Characterization of intratumoral heterogeneity is critical to cancer therapy, as the presence of phenotypically diverse cell populations commonly fuels relapse and resistance to treatment. Although genetic variation is a well-studied source of intratumoral heterogeneity, the functional impact of most genetic alterations remains unclear. Even less understood is the relative importance of other factors influencing heterogeneity, such as epigenetic state or tumor microenvironment. To investigate the relationship between genetic and transcriptional heterogeneity in a context of cancer progression, we devised a computational approach called HoneyBADGER to identify copy number variation and loss of heterozygosity in individual cells from single-cell RNA-sequencing data. By integrating allele and normalized expression information, HoneyBADGER is able to identify and infer the presence of subclone-specific alterations in individual cells and reconstruct the underlying subclonal architecture. By examining several tumor types, we show that HoneyBADGER is effective at identifying deletions, amplifications, and copy-neutral loss-of-heterozygosity events and is capable of robustly identifying subclonal focal alterations as small as 10 megabases. We further apply HoneyBADGER to analyze single cells from a progressive multiple myeloma patient to identify major genetic subclones that exhibit distinct transcriptional signatures relevant to cancer progression. Other prominent transcriptional subpopulations within these tumors did not line up with the genetic subclonal structure and were likely driven by alternative, nonclonal mechanisms. These results highlight the need for integrative analysis to understand the molecular and phenotypic heterogeneity in cancer., (© 2018 Fan et al.; Published by Cold Spring Harbor Laboratory Press.)

Lympho-myeloid restricted early thymic progenitors (ETPs) are postulated to be the cell of origin for ETP leukemias, a therapy-resistant leukemia associated with frequent co-occurrence of EZH2 and RUNX1 inactivating mutations, and constitutively activating signaling pathway mutations. In a mouse model, we demonstrate that Ezh2 and Runx1 inactivation targeted to early lymphoid progenitors causes a marked expansion of pre-leukemic ETPs, showing transcriptional signatures characteristic of ETP leukemia. Addition of a RAS-signaling pathway mutation (Flt3-ITD) results in an aggressive leukemia co-expressing myeloid and lymphoid genes, which can be established and propagated in vivo by the expanded ETPs. Both mouse and human ETP leukemias show sensitivity to BET inhibition in vitro and in vivo, which reverses aberrant gene expression induced by Ezh2 inactivation., (Copyright © 2018 Elsevier Inc. All rights reserved.)

Although previous studies suggested that the expression of FMS-like tyrosine kinase 3 (Flt3) initiates downstream of mouse hematopoietic stem cells (HSCs), FLT3 internal tandem duplications ( FLT3 ITDs) have recently been suggested to intrinsically suppress HSCs. Herein, single-cell interrogation found Flt3 mRNA expression to be absent in the large majority of phenotypic HSCs, with a strong negative correlation between Flt3 and HSC-associated gene expression. Flt3-ITD knock-in mice showed reduced numbers of phenotypic HSCs, with an even more severe loss of long-term repopulating HSCs, likely reflecting the presence of non-HSCs within the phenotypic HSC compartment. Competitive transplantation experiments established that Flt3-ITD compromises HSCs through an extrinsically mediated mechanism of disrupting HSC-supporting bone marrow stromal cells, with reduced numbers of endothelial and mesenchymal stromal cells showing increased inflammation-associated gene expression. Tumor necrosis factor (TNF), a cell-extrinsic potent negative regulator of HSCs, was overexpressed in bone marrow niche cells from FLT3-ITD mice, and anti-TNF treatment partially rescued the HSC phenotype. These findings, which establish that Flt3-ITD-driven myeloproliferation results in cell-extrinsic suppression of the normal HSC reservoir, are of relevance for several aspects of acute myeloid leukemia biology., (© 2017 Mead et al.)

Recent advances in single-cell transcriptomics are ideally placed to unravel intratumoral heterogeneity and selective resistance of cancer stem cell (SC) subpopulations to molecularly targeted cancer therapies. However, current single-cell RNA-sequencing approaches lack the sensitivity required to reliably detect somatic mutations. We developed a method that combines high-sensitivity mutation detection with whole-transcriptome analysis of the same single cell. We applied this technique to analyze more than 2,000 SCs from patients with chronic myeloid leukemia (CML) throughout the disease course, revealing heterogeneity of CML-SCs, including the identification of a subgroup of CML-SCs with a distinct molecular signature that selectively persisted during prolonged therapy. Analysis of nonleukemic SCs from patients with CML also provided new insights into cell-extrinsic disruption of hematopoiesis in CML associated with clinical outcome. Furthermore, we used this single-cell approach to identify a blast-crisis-specific SC population, which was also present in a subclone of CML-SCs during the chronic phase in a patient who subsequently developed blast crisis. This approach, which might be broadly applied to any malignancy, illustrates how single-cell analysis can identify subpopulations of therapy-resistant SCs that are not apparent through cell-population analysis.

Background: Recent advances in single-cell techniques have provided the opportunity to finely dissect cellular heterogeneity within populations previously defined by "bulk" assays and to uncover rare cell types. In human hematopoiesis, megakaryocytes and erythroid cells differentiate from a shared precursor, the megakaryocyte-erythroid progenitor (MEP), which remains poorly defined., Results: To clarify the cellular pathway in erythro-megakaryocyte differentiation, we correlate the surface immunophenotype, transcriptional profile, and differentiation potential of individual MEP cells. Highly purified, single MEP cells were analyzed using index fluorescence-activated cell sorting and parallel targeted transcriptional profiling of the same cells was performed using a specifically designed panel of genes. Differentiation potential was tested in novel, single-cell differentiation assays. Our results demonstrate that immunophenotypic MEP comprise three distinct subpopulations: "Pre-MEP," enriched for erythroid/megakaryocyte progenitors but with residual myeloid differentiation capacity; "E-MEP," strongly biased towards erythroid differentiation; and "MK-MEP," a previously undescribed, rare population of cells that are bipotent but primarily generate megakaryocytic progeny. Therefore, conventionally defined MEP are a mixed population, as a minority give rise to mixed-lineage colonies while the majority of cells are transcriptionally primed to generate exclusively single-lineage output., Conclusions: Our study clarifies the cellular hierarchy in human megakaryocyte/erythroid lineage commitment and highlights the importance of using a combination of single-cell approaches to dissect cellular heterogeneity and identify rare cell types within a population. We present a novel immunophenotyping strategy that enables the prospective identification of specific intermediate progenitor populations in erythro-megakaryopoiesis, allowing for in-depth study of disorders including inherited cytopenias, myeloproliferative disorders, and erythromegakaryocytic leukemias.

DNA binding factors are essential for regulating gene expression. CTCF and cohesin are DNA binding factors with central roles in chromatin organization and gene expression. We determined the sites of CTCF and cohesin binding to DNA in mouse brain, genome wide and in an allele-specific manner with high read-depth ChIP-seq. By comparing our results with existing data for mouse liver and embryonic stem (ES) cells, we investigated the tissue specificity of CTCF binding sites. ES cells have fewer unique CTCF binding sites occupied than liver and brain, consistent with a ground-state pattern of CTCF binding that is elaborated during differentiation. CTCF binding sites without the canonical consensus motif were highly tissue specific. In brain, a third of CTCF and cohesin binding sites coincide, consistent with the potential for many interactions between cohesin and CTCF but also many instances of independent action. In the context of genomic imprinting, CTCF and/or cohesin bind to a majority but not all differentially methylated regions, with preferential binding to the unmethylated parental allele. Whether the parental allele-specific methylation was established in the parental germlines or post-fertilization in the embryo is not a determinant in CTCF or cohesin binding. These findings link CTCF and cohesin with the control regions of a subset of imprinted genes, supporting the notion that imprinting control is mechanistically diverse.