Variant-specific introduction and dispersal dynamics of SARS-CoV-2 in New York City - from Alpha to Omicron.

Since the latter part of 2020, SARS-CoV-2 evolution has been characterised by the emergence of viral variants associated with distinct biological characteristics. While the main research focus has centred on the ability of new variants to increase in frequency and impact the effective reproductive number of the virus, less attention has been placed on their relative ability to establish transmission chains and to spread through a geographic area. Here, we describe a phylogeographic approach to estimate and compare the introduction and dispersal dynamics of the main SARS-CoV-2 variants - Alpha, Iota, Delta, and Omicron - that circulated in the New York City area between 2020 and 2022. Notably, our results indicate that Delta had a lower ability to establish sustained transmission chains in the NYC area and that Omicron (BA.1) was the variant fastest to disseminate across the study area. The analytical approach presented here complements non-spatially-explicit analytical approaches that seek a better understanding of the epidemiological differences that exist among successive SARS-CoV-2 variants of concern.

Author Correction: The bone marrow microenvironment at single-cell resolution.

An Amendment to this paper has been published and can be accessed via a link at the top of the paper.

The bone marrow microenvironment at single-cell resolution.

The bone marrow microenvironment has a key role in regulating haematopoiesis, but its molecular complexity and response to stress are incompletely understood. Here we map the transcriptional landscape of mouse bone marrow vascular, perivascular and osteoblast cell populations at single-cell resolution, both at homeostasis and under conditions of stress-induced haematopoiesis. This analysis revealed previously unappreciated levels of cellular heterogeneity within the bone marrow niche and resolved cellular sources of pro-haematopoietic growth factors, chemokines and membrane-bound ligands. Our studies demonstrate a considerable transcriptional remodelling of niche elements under stress conditions, including an adipocytic skewing of perivascular cells. Among the stress-induced changes, we observed that vascular Notch delta-like ligands (encoded by Dll1 and Dll4) were downregulated. In the absence of vascular Dll4, haematopoietic stem cells prematurely induced a myeloid transcriptional program. These findings refine our understanding of the cellular architecture of the bone marrow niche, reveal a dynamic and heterogeneous molecular landscape that is highly sensitive to stress and illustrate the utility of single-cell transcriptomic data in evaluating the regulation of haematopoiesis by discrete niche populations.

Remdesivir Resistance in Transplant Recipients With Persistent Coronavirus Disease 2019.

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.

De novo assembly and annotation of the singing mouse genome.

BACKGROUND: Developing genomic resources for a diverse range of species is an important step towards understanding the mechanisms underlying complex traits. Specifically, organisms that exhibit unique and accessible phenotypes-of-interest allow researchers to address questions that may be ill-suited to traditional model organisms. We sequenced the genome and transcriptome of Alston's singing mouse (Scotinomys teguina), an emerging model for social cognition and vocal communication. In addition to producing advertisement songs used for mate attraction and male-male competition, these rodents are diurnal, live at high-altitudes, and are obligate insectivores, providing opportunities to explore diverse physiological, ecological, and evolutionary questions. RESULTS: Using PromethION, Illumina, and PacBio sequencing, we produced an annotated genome and transcriptome, which were validated using gene expression and functional enrichment analyses. To assess the usefulness of our assemblies, we performed single nuclei sequencing on cells of the orofacial motor cortex, a brain region implicated in song coordination, identifying 12 cell types. CONCLUSIONS: These resources will provide the opportunity to identify the molecular basis of complex traits in singing mice as well as to contribute data that can be used for large-scale comparative analyses.

Sequencing identifies multiple early introductions of SARS-CoV-2 to the New York City region.

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.

Dispersal dynamics of SARS-CoV-2 lineages during the first epidemic wave in New York City.

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.

Amplification Artifact in SARS-CoV-2 Omicron Sequences Carrying P681R Mutation, New York, USA.

Of 379 severe acute respiratory syndrome coronavirus 2 samples collected in New York, USA, we detected 86 Omicron variant sequences containing Delta variant mutation P681R. Probable explanations were co-infection with 2 viruses or contamination/amplification artifact. Repeated library preparation with fewer cycles showed the P681R calls were artifactual. Unusual mutations should be interpreted with caution.

Serial single-cell profiling analysis of metastatic TNBC during Nab-paclitaxel and pembrolizumab treatment.

PURPOSE: Immunotherapy has recently been shown to improve outcomes for advanced PD-L1-positive triple-negative breast cancer (TNBC) in the Impassion130 trial, leading to FDA approval of the first immune checkpoint inhibitor in combination with taxane chemotherapy. To further develop predictive biomarkers and improve therapeutic efficacy of the combination, interrogation of the tumor immune microenvironment before therapy as well as during each component of treatment is crucial. Here we use single-cell RNA sequencing (scRNA-seq) on tumor biopsies to assess immune cell changes from two patients with advanced TNBC treated in a prospective trial at predefined serial time points, before treatment, on taxane chemotherapy and on chemo-immunotherapy. METHODS: Both patients (one responder and one progressor) received the trial therapy, in cycle 1 nab-paclitaxel given as single agent, in cycle 2 nab-paclitaxel in combination with pembrolizumab. Tumor core biopsies were obtained at baseline, 3 weeks (after cycle 1, chemotherapy alone) and 6 weeks (after cycle 2, chemo-immunotherapy). Single-cell RNA sequencing (scRNA-seq) of both cancer cells and infiltrating immune cells isolated were performed from fresh tumor core biopsy specimens by 10 × chromium sequencing. RESULTS: ScRNA-seq analysis showed significant baseline heterogeneity of tumor-infiltrating immune cell populations between the two patients as well as modulation of the tumor microenvironment by chemotherapy and immunotherapy. In the responding patient there was a population of PD-1high-expressing T cells which significantly decreased after nab-paclitaxel plus pembrolizumab treatment as well as a presence of tissue-resident memory T cells (TRM). In contrast, tumors from the patient with rapid disease progression showed a prevalent and persistent myeloid compartment. CONCLUSIONS: Our study provides a deep cellular analysis of on-treatment changes during chemo-immunotherapy for advanced TNBC, demonstrating not only feasibility of single-cell analyses on serial tumor biopsies but also the heterogeneity of TNBC and differences in on-treatment changes in responder versus progressor.

How well do you know your mutation? Complex effects of genetic background on expressivity, complementation, and ordering of allelic effects.

For a given gene, different mutations influence organismal phenotypes to varying degrees. However, the expressivity of these variants not only depends on the DNA lesion associated with the mutation, but also on factors including the genetic background and rearing environment. The degree to which these factors influence related alleles, genes, or pathways similarly, and whether similar developmental mechanisms underlie variation in the expressivity of a single allele across conditions and among alleles is poorly understood. Besides their fundamental biological significance, these questions have important implications for the interpretation of functional genetic analyses, for example, if these factors alter the ordering of allelic series or patterns of complementation. We examined the impact of genetic background and rearing environment for a series of mutations spanning the range of phenotypic effects for both the scalloped and vestigial genes, which influence wing development in Drosophila melanogaster. Genetic background and rearing environment influenced the phenotypic outcome of mutations, including intra-genic interactions, particularly for mutations of moderate expressivity. We examined whether cellular correlates (such as cell proliferation during development) of these phenotypic effects matched the observed phenotypic outcome. While cell proliferation decreased with mutations of increasingly severe effects, surprisingly it did not co-vary strongly with the degree of background dependence. We discuss these findings and propose a phenomenological model to aid in understanding the biology of genes, and how this influences our interpretation of allelic effects in genetic analysis.

Severe obesity and bariatric surgery alter the platelet mRNA profile.

Mechanisms explaining the relationship between obesity and cardiovascular disease (CVD) are needed. Despite growing recognition of the importance of the anucleate platelet transcriptome, low levels of RNA in platelets make assessment difficult. We sought to perform unbiased platelet RNA profiling in obesity by performing a prospective study of severe obesity and weight loss via bariatric surgery on platelet characteristics and mRNA profile in 26 pre-menopausal, non-diabetic women (31.6 ± 8.4 years; BMI 43.0 ± 6.5 kg/m2) who underwent sleeve gastrectomy. Totally, 10 women of similar age with normal BMI served as controls. Platelet activation via flow cytometry was assessed before and after surgery. RNA-sequencing (RNAseq) was performed on platelet isolates from a subset of 13 subjects (eight obese women and five normal-BMI subjects). Platelet count, size, and age did not differ between control and obese women. However, platelet surface P-selectin and CD40 were higher in obesity. RNAseq demonstrated 629 differentially abundant transcripts in obesity. Notably, S100A9 and AGER, established markers of cardiovascular risk, were two of the most highly upregulated transcripts (each > 2.5 fold). At 6 months post-operatively, subjects lost 26.1 ± 5.8% body weight and inducible platelet P-selectin expression was reduced. Expression of 170 transcripts was affected by surgery, but only a small fraction (46/629) were genes found altered in obesity. We demonstrate that obesity is associated with an altered platelet transcriptome and increased platelet activation, which is partly attenuated by bariatric surgery. These observations suggest that platelets may contribute to increased cardiovascular risk in obesity through a variety of mechanisms.

Emergence and dissemination of SARS-CoV-2 XBB.1.5 in New York.

The recombinant SARS-CoV-2 Omicron XBB.1.5 variant was first detected in New York City (NYC) and rapidly became the predominant variant in the area by early 2023. The increased occurrence of circulating variants within the SARS-CoV-2 XBB-sublineage prompted the modification of COVID-19 mRNA vaccines by Moderna and Pfizer-BioNTech. This update, implemented in mid-September 2023, involved the incorporation of a monovalent XBB.1.5 component. Considering that NYC probably played a central role in the emergence of the XBB.1.5 variant, we conducted phylogeographic analysis to investigate the emergence and spread of this variant in the metropolitan area. Our analysis confirms that XBB.1.5 emerged within or near the NYC area and indicates that XBB.1.5 had a diffusion velocity similar to that of the variant Alpha in the same study area. Additionally, the analysis of 2,392 genomes collected in the context of the genomic surveillance program at NYU Langone Health system showed that there was no increased proportion of XBB.1.5, relative to all cocirculating variants, in the boosted compared to unvaccinated individuals. This study provides a comprehensive description of the emergence and dissemination of XBB.1.5.

Generation of quality-controlled SARS-CoV-2 variant stocks.

One of the main challenges in the fight against coronavirus disease 2019 (COVID-19) stems from the ongoing evolution of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) into multiple variants. To address this hurdle, research groups around the world have independently developed protocols to isolate these variants from clinical samples. These isolates are then used in translational and basic research-for example, in vaccine development, drug screening or characterizing SARS-CoV-2 biology and pathogenesis. However, over the course of the COVID-19 pandemic, we have learned that the introduction of artefacts during both in vitro isolation and subsequent propagation to virus stocks can lessen the validity and reproducibility of data. We propose a rigorous pipeline for the generation of high-quality SARS-CoV-2 variant clonal isolates that minimizes the acquisition of mutations and introduces stringent controls to detect them. Overall, the process includes eight stages: (i) cell maintenance, (ii) isolation of SARS-CoV-2 from clinical specimens, (iii) determination of infectious virus titers by plaque assay, (iv) clonal isolation by plaque purification, (v) whole-virus-genome deep-sequencing, (vi and vii) amplification of selected virus clones to master and working stocks and (viii) sucrose purification. This comprehensive protocol will enable researchers to generate reliable SARS-CoV-2 variant inoculates for in vitro and in vivo experimentation and will facilitate comparisons and collaborative work. Quality-controlled working stocks for most applications can be generated from acquired biorepository virus within 1 month. An additional 5-8 d are required when virus is isolated from clinical swab material, and another 6-7 d is needed for sucrose-purifying the stocks.

Selective adaptation of SARS-CoV-2 Omicron under booster vaccine pressure: a multicentre observational study.

BACKGROUND: High rates of vaccination and natural infection drive immunity and redirect selective viral adaptation. Updated boosters are installed to cope with drifted viruses, yet data on adaptive evolution under increasing immune pressure in a real-world situation are lacking. METHODS: Cross-sectional study to characterise SARS-CoV-2 mutational dynamics and selective adaptation over >1 year in relation to vaccine status, viral phylogenetics, and associated clinical and demographic variables. FINDINGS: The study of >5400 SARS-CoV-2 infections between July 2021 and August 2022 in metropolitan New York portrayed the evolutionary transition from Delta to Omicron BA.1-BA.5 variants. Booster vaccinations were implemented during the Delta wave, yet booster breakthrough infections and SARS-CoV-2 re-infections were almost exclusive to Omicron. In adjusted logistic regression analyses, BA.1, BA.2, and BA.5 had a significant growth advantage over co-occurring lineages in the boosted population, unlike BA.2.12.1 or BA.4. Selection pressure by booster shots translated into diffuse adaptive evolution in Delta spike, contrasting with strong, receptor-binding motif-focused adaptive evolution in BA.2-BA.5 spike (Fisher Exact tests; non-synonymous/synonymous mutation rates per site). Convergent evolution has become common in Omicron, engaging spike positions crucial for immune escape, receptor binding, or cleavage. INTERPRETATION: Booster shots are required to cope with gaps in immunity. Their discriminative immune pressure contributes to their effectiveness but also requires monitoring of selective viral adaptation processes. Omicron BA.2 and BA.5 had a selective advantage under booster vaccination pressure, contributing to the evolution of BA.2 and BA.5 sublineages and recombinant forms that predominate in 2023. FUNDING: The study was supported by NYU institutional funds and partly by the Cancer Center Support Grant P30CA016087 at the Laura and Isaac Perlmutter Cancer Center.

Delta-Omicron recombinant escapes therapeutic antibody neutralization.

The emergence of recombinant viruses is a threat to public health, as recombination may integrate variant-specific features that together result in escape from treatment or immunity. The selective advantages of recombinant SARS-CoV-2 isolates over their parental lineages remain unknown. We identified a Delta-Omicron (AY.45-BA.1) recombinant in an immunosuppressed transplant recipient treated with monoclonal antibody Sotrovimab. The single recombination breakpoint is located in the spike N-terminal domain adjacent to the Sotrovimab binding site. While Delta and BA.1 are sensitive to Sotrovimab neutralization, the Delta-Omicron recombinant is highly resistant. To our knowledge, this is the first described instance of recombination between circulating SARS-CoV-2 variants as a functional mechanism of resistance to treatment and immune escape.

Remdesivir resistance in transplant recipients with persistent COVID-19.

The medical community currently lacks robust data regarding the incidence, prevalence, and clinical significance of mutations associated with resistance to severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) therapeutics. This report describes two renal transplant recipients who, after remdesivir exposure, developed a de novo V792I RNA-dependent RNA polymerase (RdRp) mutation that has recently been found to confer resistance to remdesivir in vitro . To the best of our knowledge, this publication is the first to document the emergence of V792I in patients treated with remdesivir. Our work underscores the critical need for augmented efforts to identify concerning mutations and address their clinical implications.

Clinical and genomic signatures of SARS-CoV-2 Delta breakthrough infections in New York.

BACKGROUND: In 2021, Delta became the predominant SARS-CoV-2 variant worldwide. While vaccines have effectively prevented COVID-19 hospitalization and death, vaccine breakthrough infections increasingly occurred. The precise role of clinical and genomic determinants in Delta infections is not known, and whether they contributed to increased rates of breakthrough infections compared to unvaccinated controls. METHODS: We studied SARS-CoV-2 variant distribution, dynamics, and adaptive selection over time in relation to vaccine status, phylogenetic relatedness of viruses, full genome mutation profiles, and associated clinical and demographic parameters. FINDINGS: We show a steep and near-complete replacement of circulating variants with Delta between May and August 2021 in metropolitan New York. We observed an increase of the Delta sublineage AY.25 (14% in vaccinated, 7% in unvaccinated), its spike mutation S112L, and AY.44 (8% in vaccinated, 2% in unvaccinated) with its nsp12 mutation F192V in breakthroughs. Delta infections were associated with younger age and lower hospitalization rates than Alpha. Delta breakthrough infections increased significantly with time since vaccination, and, after adjusting for confounders, they rose at similar rates as in unvaccinated individuals. INTERPRETATION: We observed a modest adaptation of Delta genomes in breakthrough infections in New York, suggesting an improved genomic framework to support Delta's epidemic growth in times of waning vaccine protection despite limited impact on vaccine escape. FUNDING: The study was supported by NYU institutional funds. The NYULH Genome Technology Center is partially supported by the Cancer Center Support Grant P30CA016087 at the Laura and Isaac Perlmutter Cancer Center.

Delta-Omicron recombinant escapes therapeutic antibody neutralization.

Background: The emergence of recombinant viruses is a threat to public health. Recombination of viral variants may combine variant-specific features that together catalyze viral escape from treatment or immunity. The selective advantages of recombinant SARS-CoV-2 isolates over their parental lineages remain unknown. Methods: Multi-method amplicon and metagenomic sequencing of a clinical swab and the in vitro grown virus allowed for high-confidence detection of a novel recombinant variant. Mutational, phylogeographic, and structural analyses determined features of the recombinant genome and spike protein. Neutralization assays using infectious as well as pseudotyped viruses and point mutants thereof defined the recombinant's sensitivity to a panel of monoclonal antibodies and sera from vaccinated and/or convalescent individuals. Results: A novel Delta-Omicron SARS-CoV-2 recombinant was identified in an unvaccinated, immunosuppressed kidney transplant recipient treated with monoclonal antibody Sotrovimab. The recombination breakpoint is located in the spike N-terminal domain, adjacent to the Sotrovimab quaternary binding site, and results in a 5'-Delta AY.45 and a 3'-Omicron BA.1 mosaic spike protein. Delta and BA.1 are sensitive to Sotrovimab neutralization, whereas the Delta-Omicron recombinant is highly resistant to Sotrovimab, both with and without the RBD resistance mutation E340D. Conclusions: Recombination between circulating SARS-CoV-2 variants can functionally contribute to immune escape. It is critical to validate phenotypes of mosaic viruses and monitor immunosuppressed COVID-19 patients treated with monoclonal antibodies for the selection of recombinant and immune escape variants. (Funded by NYU, the National Institutes of Health, and others).

Dominance of Alpha and Iota variants in SARS-CoV-2 vaccine breakthrough infections in New York City.

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. The median age of patients in the study was 48 years; 7 required hospitalization and 1 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 χ2 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.

Clinical and genomic signatures of rising SARS-CoV-2 Delta breakthrough infections in New York.

In 2021, Delta has become the predominant SARS-CoV-2 variant worldwide. While vaccines effectively prevent COVID-19 hospitalization and death, vaccine breakthrough infections increasingly occur. The precise role of clinical and genomic determinants in Delta infections is not known, and whether they contribute to increased rates of breakthrough infections compared to unvaccinated controls. Here, we show a steep and near complete replacement of circulating variants with Delta between May and August 2021 in metropolitan New York. We observed an increase of the Delta sublineage AY.25, its spike mutation S112L, and nsp12 mutation F192V in breakthroughs. Delta infections were associated with younger age and lower hospitalization rates than Alpha. Delta breakthroughs increased significantly with time since vaccination, and, after adjusting for confounders, they rose at similar rates as in unvaccinated individuals. Our data indicate a limited impact of vaccine escape in favor of Delta's increased epidemic growth in times of waning vaccine protection.