Summary The molecular complexity of the bone marrow (BM) microenvironment and its response to stress are incompletely understood, despite its key role in the regulation of hematopoiesis. Here we map the transcriptional landscape of BM vascular, perivascular, and osteoblast niche populations at single-cell resolution at both homeostasis and under stress hematopoiesis. This analysis revealed a previously unappreciated level of cellular heterogeneity within the BM niche, identified novel cellular subsets, and resolved cellular sources of pro-hematopoietic growth factors, chemokines, and membrane-bound ligands. Under conditions of stress, our studies revealed a significant transcriptional remodeling of these niche elements, including an adipocytic skewing of the perivascular cells. Among the stress-induced changes, we observed that vascular Notch ligand delta-like ligands (Dll1,4) were downregulated. In the absence of vascular Dll4, hematopoietic stem cells (HSC) prematurely induced a myeloid transcriptional program. These findings refine our understanding of the cellular architecture of the BM niche, reveal a dynamic and heterogeneous molecular landscape that is highly sensitive to stress, and illustrate the utility of single cell transcriptomic data in systematically evaluating the regulation of hematopoiesis by discrete niche populations.

Effective public response to a pandemic relies upon accurate measurement of the extent and dynamics of an outbreak. Viral genome sequencing has emerged as a powerful approach to link seemingly unrelated cases, and large-scale sequencing surveillance can inform on critical epidemiological parameters. Here, we report the analysis of 864 SARS-CoV-2 sequences from cases in the New York City metropolitan area during the COVID-19 outbreak in Spring 2020. The majority of cases had no recent travel history or known exposure, and genetically linked cases were spread throughout the region. Comparison to global viral sequences showed that early transmission was most linked to cases from Europe. Our data are consistent with numerous seeds from multiple sources and a prolonged period of unrecognized community spreading. This work highlights the complementary role of genomic surveillance in addition to traditional epidemiological indicators.

Splicing alterations are common in diseases such as cancer, where mutations in splicing factor genes are frequently responsible for aberrant splicing. Here we present an alternative mechanism for splicing regulation in T-cell acute lymphoblastic leukemia (T-ALL) that involves posttranslational stabilization of the splicing machinery via deubiquitination. We demonstrate there are extensive exon skipping changes in disease, affecting proteasomal subunits, cell-cycle regulators, and the RNA machinery. We present that the serine/arginine-rich splicing factors (SRSF), controlling exon skipping, are critical for leukemia cell survival. The ubiquitin-specific peptidase 7 (USP7) regulates SRSF6 protein levels via active deubiquitination, and USP7 inhibition alters the exon skipping pattern and blocks T-ALL growth. The splicing inhibitor H3B-8800 affects splicing of proteasomal transcripts and proteasome activity and acts synergistically with proteasome inhibitors in inhibiting T-ALL growth. Our study provides the proof-of-principle for regulation of splicing factors via deubiquitination and suggests new therapeutic modalities in T-ALL. SignifiCAnCE: Our study provides a new proof-of-principle for posttranslational regulation of splicing factors independently of mutations in aggressive T-cell leukemia. It further suggests a new drug combination of splicing and proteasomal inhibitors, a concept that might apply to other diseases with or without mutations affecting the splicing machinery.

New mutations conferring resistance to SARS-CoV-2 therapeutics have important clinical implications. We describe the first cases of an independently acquired V792I RNA-dependent RNA polymerase mutation developing in renal transplant recipients after remdesivir exposure. Our work underscores the need for augmented efforts to identify concerning mutations and address their clinical implications.

During the first phase of the COVID-19 epidemic, New York City rapidly became the epicenter of the pandemic in the United States. While molecular phylogenetic analyses have previously highlighted multiple introductions and a period of cryptic community transmission within New York City, little is known about the circulation of SARS-CoV-2 within and among its boroughs. We here perform phylogeographic investigations to gain insights into the circulation of viral lineages during the first months of the New York City outbreak. Our analyses describe the dispersal dynamics of viral lineages at the state and city levels, illustrating that peripheral samples likely correspond to distinct dispersal events originating from the main metropolitan city areas. In line with the high prevalence recorded in this area, our results highlight the relatively important role of the borough of Queens as a transmission hub associated with higher local circulation and dispersal of viral lineages toward the surrounding boroughs.

Effective public response to a pandemic relies upon accurate measurement of the extent and dynamics of an outbreak. Viral genome sequencing has emerged as a powerful approach to link seemingly unrelated cases, and large-scale sequencing surveillance can inform on critical epidemiological parameters. Here, we report the analysis of 236 SARS-CoV2 sequences from cases in the New York City metropolitan area during the initial stages of the 2020 COVID-19 outbreak. The majority of cases throughout the region had no recent travel history or known exposure, and genetically linked cases were spread throughout the region. Comparison to global viral sequences showed that the majority were most related to cases from Europe. Our data are consistent with numerous seed transmissions from multiple sources and a prolonged period of unrecognized community spreading. This work highlights the complementary role of real-time genomic surveillance in addition to traditional epidemiological indicators.

COVID-19 is a respiratory illness caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) and declared by the World Health Organization a global public health emergency. Among the severe outbreaks across South America, Uruguay has become known for curtailing SARS-CoV-2 exceptionally well. To understand the SARS-CoV-2 introductions, local transmissions, and associations with genomic and clinical parameters in Uruguay, we sequenced the viral genomes of 44 outpatients and inpatients in a private healthcare system in its capital, Montevideo, from March to May 2020. We performed a phylogeographic analysis using sequences from our cohort and other studies that indicate a minimum of 23 independent introductions into Uruguay, resulting in five major transmission clusters. Our data suggest that most introductions resulting in chains of transmission originate from other South American countries, with the earliest seeding of the virus in late February 2020, weeks before the borders were closed to all non-citizens and a partial lockdown implemented. Genetic analyses suggest a dominance of S and G clades (G, GH, GR) that make up >90% of the viral strains in our study. In our cohort, lethal outcome of SARS-CoV-2 infection significantly correlated with arterial hypertension, kidney failure, and ICU admission (FDR < 0.01), but not with any mutation in a structural or non-structural protein, such as the spike D614G mutation. Our study contributes genetic, phylodynamic, and clinical correlation data about the exceptionally well-curbed SARS-CoV-2 outbreak in Uruguay, which furthers the understanding of disease patterns and regional aspects of the pandemic in Latin America.

The efficacy of COVID-19 mRNA vaccines is high, but breakthrough infections still occur. We compared the SARS-CoV-2 genomes of 76 breakthrough cases after full vaccination with BNT162b2 (Pfizer/BioNTech), mRNA-1273 (Moderna), or JNJ-78436735 (Janssen) to unvaccinated controls (February-April 2021) in metropolitan New York, including their phylogenetic relationship, distribution of variants, and full spike mutation profiles. Their median age was 48 years; seven required hospitalization and one died. Most breakthrough infections (57/76) occurred with B.1.1.7 (Alpha) or B.1.526 (Iota). Among the 7 hospitalized cases, 4 were infected with B.1.1.7, including 1 death. Both unmatched and matched statistical analyses considering age, sex, vaccine type, and study month as covariates supported the null hypothesis of equal variant distributions between vaccinated and unvaccinated in chi-squared and McNemar tests (p>0.1) highlighting a high vaccine efficacy against B.1.1.7 and B.1.526. There was no clear association among breakthroughs between type of vaccine received and variant. In the vaccinated group, spike mutations in the N-terminal domain and receptor-binding domain that have been associated with immune evasion were overrepresented. The evolving dynamic of SARS-CoV-2 variants requires broad genomic analyses of breakthrough infections to provide real-life information on immune escape mediated by circulating variants and their spike mutations.