Perivascular niche cells sense thrombocytopenia and activate hematopoietic stem cells in an IL-1 dependent manner.

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.

High-depth sequencing characterization of viral dynamics across tissues in fatal COVID-19 reveals compartmentalized infection.

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.

Natural history of Ebola virus disease in rhesus monkeys shows viral variant emergence dynamics and tissue-specific host responses.

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.

Analysis of 6.4 million SARS-CoV-2 genomes identifies mutations associated with fitness.

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 PyR0, 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 PyR0 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. PyR0 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.

Analysis of 6.4 million SARS-CoV-2 genomes identifies mutations associated with fitness.

Repeated emergence of SARS-CoV-2 variants with increased fitness necessitates rapid detection and characterization of new lineages. To address this need, we 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 non-spike mutations within the nucleocapsid and nonstructural proteins. PyR 0 forecasts growth of new lineages from their mutational profile, identifies viral lineages of concern as they emerge, and prioritizes mutations of biological and public health concern for functional characterization. ONE SENTENCE SUMMARY: A Bayesian hierarchical model of all SARS-CoV-2 viral genomes predicts lineage fitness and identifies associated mutations.

Transmission from vaccinated individuals in a large SARS-CoV-2 Delta variant outbreak.

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.

External signals regulate continuous transcriptional states in hematopoietic stem cells.

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.

Most organs in the human body are maintained by a type of immature cells known as adult stem cells, which ensure a constant supply of new, mature cells. Adult stem cells monitor their environment through external signalling molecules and replace damaged cells as needed. Stem cell therapy takes advantage of the regenerative ability of immature stem cells and can be helpful for conditions such as blood diseases, autoimmune diseases, neurodegeneration and cancer. For example, hematopoietic stem-cell transplantation is a treatment for some types of cancer and blood disorders, in which stem cells are harvested from the blood or bone marrow and reintroduced into the body, where they can develop into all types of blood cells, including white blood cells, red blood cells and platelets. Hematopoietic stem-cell transplants have been in use for over 30 years, but they remain a highly risky procedure. One of the challenges is that outcomes can vary between patients and many of the factors that can influence the 'regenerative' potential of hematopoietic stem cells, such as external signalling molecules, are not well understood. To fill this gap, Fast et al. analysed which genes are turned on and off in hematopoietic stem cells in response to several external signalling molecules. To do so, three signalling pathways in mice were altered by injecting them with different chemicals. After two hours, the hematopoietic stem cells were purified and the gene expression for each cell was analysed. This revealed that the types of genes and the strength at which they were affected by each chemical was unique. Moreover, hematopoietic stem cells responded rapidly to external signals, with substantial differences in gene expression between individual groups of cells. Contrary to more specialised cells, the external signalling genes in some hematopoietic stem cells were already activated without being injected with external signalling molecules. This suggest that low levels of external signalling molecules released from their microenvironment may prepare stem cells to better respond to future stress or injuries. These results help to better understand stem cells and to evaluate how the signalling state of hematopoietic stem cells affects regeneration, and ultimately improve hematopoietic stem cell transplantation for patients.

Evidence of transmission from fully vaccinated individuals in a large outbreak of the SARS-CoV-2 Delta variant in Provincetown, Massachusetts.

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.

Imprinted Gene Expression and Function of the Dopa Decarboxylase Gene in the Developing Heart.

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.

COVID-19 tissue atlases reveal SARS-CoV-2 pathology and cellular targets.

COVID-19, which is caused by SARS-CoV-2, can result in acute respiratory distress syndrome and multiple organ failure1-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.

Single cell transcriptomics of primate sensory neurons identifies cell types associated with chronic pain.

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.

Distinct evolutionary paths in chronic lymphocytic leukemia during resistance to the graft-versus-leukemia effect.

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.

Transcription of intragenic CpG islands influences spatiotemporal host gene pre-mRNA processing.

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.

Myocardial differentiation is dependent upon endocardial signaling during early cardiogenesis in vitro.

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 Nfatc1 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.

Unravelling Intratumoral Heterogeneity through High-Sensitivity Single-Cell Mutational Analysis and Parallel RNA Sequencing.

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.

Linking transcriptional and genetic tumor heterogeneity through allele analysis of single-cell RNA-seq data.

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.

Ezh2 and Runx1 Mutations Collaborate to Initiate Lympho-Myeloid Leukemia in Early Thymic Progenitors.

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.

Niche-mediated depletion of the normal hematopoietic stem cell reservoir by Flt3-ITD-induced myeloproliferation.

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.

Single-cell transcriptomics uncovers distinct molecular signatures of stem cells in chronic myeloid leukemia.

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.