Optimisation and evaluation of viral genomic sequencing of SARS-CoV-2 rapid diagnostic tests: a laboratory and cohort-based study.

BACKGROUND: Sequencing of SARS-CoV-2 from rapid diagnostic tests (RDTs) can bolster viral genomic surveillance efforts; however, approaches to maximise and standardise pathogen genome recovery from RDTs remain underdeveloped. We aimed to systematically optimise the elution of genetic material from RDT components and to evaluate the efficacy of RDT sequencing for outbreak investigation. METHODS: In this laboratory and cohort-based study we seeded RDTs with inactivated SARS-CoV-2 to optimise the elution of genomic material from RDT lateral flow strips. We measured the effect of changes in buffer type, time in buffer, and rotation on PCR cycle threshold (Ct) value. We recruited individuals older than 18 years residing in the greater Boston area, MA, USA, from July 18 to Nov 5, 2022, via email advertising to students and staff at Harvard University, MA, USA, and via broad social media advertising. All individuals recruited were within 5 days of a positive diagnostic test for SARS-CoV-2; no other relevant exclusion criteria were applied. Each individual completed two RDTs and one PCR swab. On Dec 29, 2022, we also collected RDTs from a convenience sample of individuals who were positive for SARS-CoV-2 and associated with an outbreak at a senior housing facility in MA, USA. We extracted all returned PCR swabs and RDT components (ie, swab, strip, or buffer); samples with a Ct of less than 40 were subject to amplicon sequencing. We compared the efficacy of elution and sequencing across RDT brands and components and used RDT-derived sequences to infer transmission links within the outbreak at the senior housing facility. We conducted metagenomic sequencing of negative RDTs from symptomatic individuals living in the senior housing facility. FINDINGS: Neither elution duration of greater than 10 min nor rotation during elution impacted viral titres. Elution in Buffer AVL (Ct=31·4) and Tris-EDTA Buffer (Ct=30·8) were equivalent (p=0·34); AVL outperformed elution in lysis buffer and 50% lysis buffer (Ct=40·0, p=0·0029 for both) as well as Universal Viral Transport Medium (Ct=36·7, p=0·079). Performance of RDT strips was poorer than that of matched PCR swabs (mean Ct difference 10·2 [SD 4·3], p<0·0001); however, RDT swabs performed similarly to PCR swabs (mean Ct difference 4·1 [5·2], p=0·055). No RDT brand significantly outperformed another. Across sample types, viral load predicted the viral genome assembly length. We assembled greater than 80% complete genomes from 12 of 17 RDT-derived swabs, three of 18 strips, and four of 11 residual buffers. We generated outbreak-associated SARS-CoV-2 genomes using both amplicon and metagenomic sequencing and identified multiple introductions of the virus that resulted in downstream transmission. INTERPRETATION: RDT-derived swabs are a reasonable alternative to PCR swabs for viral genomic surveillance and outbreak investigation. RDT-derived lateral flow strips yield accurate, but significantly fewer, viral reads than matched PCR swabs. Metagenomic sequencing of negative RDTs can identify viruses that might underlie patient symptoms. FUNDING: The National Science Foundation, the Hertz Foundation, the National Institute of General Medical Sciences, Harvard Medical School, the Howard Hughes Medical Institute, the US Centers for Disease Control and Prevention, the Broad Institute and the National Institute of Allergy and Infectious Diseases.

Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) Household Transmission during the Omicron Era in Massachusetts: A Prospective, Case-Ascertained Study using Genomic Epidemiology.

Background: Households are a major setting for SARS-CoV-2 infections, but there remains a lack of knowledge regarding the dynamics of viral transmission, particularly in the setting of widespread pre-existing SARS-CoV-2 immunity and evolving variants. Methods: We conducted a prospective, case-ascertained household transmission study in the greater Boston area in March-July 2022. Anterior nasal swabs, along with clinical and demographic data, were collected for 14 days. Nasal swabs were tested for SARS-CoV-2 by PCR. Whole genome sequencing was performed on high-titer samples. Results: We enrolled 33 households in a primary analysis set, with a median age of participants of 25 years old (range 2-66); 98% of whom had received at least 2 doses of a COVID-19 vaccine. 58% of households had a secondary case during follow up and the secondary attack rate (SAR) for contacts infected was 39%. We further examined a strict analysis set of 21 households that had only 1 PCR+ case at baseline, finding an SAR of 22.5%. Genomic epidemiology further determined that there were multiple sources of infection for household contacts, including the index case and outside introductions. When limiting estimates to only highly probable transmissions given epidemiologic and genomic data, the SAR was 18.4%. Conclusions: Household contacts of a person newly diagnosed with COVID-19 are at high risk for SARS-CoV-2 infection in the following 2 weeks. This is, however, not only due to infection from the household index case, but also because the presence of an infected household member implies increased SARS-CoV-2 community transmission. Further studies to understand and mitigate household transmission are needed.

Antigen test swabs are comparable to nasopharyngeal swabs for sequencing of SARS-CoV-2.

Viral genomic surveillance has been integral in the global response to the SARS-CoV-2 pandemic. Surveillance efforts rely on the availability of representative clinical specimens from ongoing testing activities. However, testing practices have recently shifted due to the widespread availability and use of rapid antigen tests, which could lead to gaps in future monitoring efforts. As such, genomic surveillance strategies must adapt to include laboratory workflows that are robust to sample type. To that end, we compare the results of RT-qPCR and viral genome sequencing using samples from positive BinaxNOW COVID-19 Antigen Card swabs (N = 555) to those obtained from nasopharyngeal (NP) swabs used for nucleic acid amplification testing (N = 135). We show that swabs obtained from antigen cards are comparable in performance to samples from NP swabs, providing a viable alternative and allowing for the potential expansion of viral genomic surveillance to outpatient clinic as well as other settings where rapid antigen tests are often used.

The 2022 RSV surge was driven by multiple viral lineages.

The US experienced an early and severe respiratory syncytial virus (RSV) surge in autumn 2022. Despite the pressure this has put on hospitals and care centers, the factors promoting the surge in cases are unknown. To investigate whether viral characteristics contributed to the extent or severity of the surge, we sequenced 105 RSV-positive specimens from symptomatic patients diagnosed with RSV who presented to the Massachusetts General Hospital (MGH) and its outpatient practices in the Greater Boston Area. Genomic analysis of the resulting 77 genomes (54 with >80% coverage, and 23 with >5% coverage) demonstrated that the surge was driven by multiple lineages of RSV-A (91%; 70/77) and RSV-B (9%; 7/77). Phylogenetic analysis of all US RSV-A revealed 12 clades, 4 of which contained Massachusetts and Washington genomes. These clades individually had times to most recent common ancestor (tMRCA) between 2014 and 2017, and together had a tMRCA of 2009, suggesting that they emerged well before the COVID-19 pandemic. Similarly, the RSV-B genomes had a tMRCA between 2016 and 2019. We found that the RSV-A and RSV-B genomes in our sample did not differ statistically from the estimated clock rate of the larger phylogenetic tree (10.6 and 12.4 substitutions per year, respectively). In summary, the polyphyletic nature of viral genomes sequenced in the US during the autumn 2022 surge is inconsistent with the emergence of a single, highly transmissible causal RSV lineage.

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.

Multimodal surveillance of SARS-CoV-2 at a university enables development of a robust outbreak response framework.

BACKGROUND: Universities are vulnerable to infectious disease outbreaks, making them ideal environments to study transmission dynamics and evaluate mitigation and surveillance measures. Here, we analyze multimodal COVID-19-associated data collected during the 2020-2021 academic year at Colorado Mesa University and introduce a SARS-CoV-2 surveillance and response framework. METHODS: We analyzed epidemiological and sociobehavioral data (demographics, contact tracing, and WiFi-based co-location data) alongside pathogen surveillance data (wastewater and diagnostic testing, and viral genomic sequencing of wastewater and clinical specimens) to characterize outbreak dynamics and inform policy. We applied relative risk, multiple linear regression, and social network assortativity to identify attributes or behaviors associated with contracting SARS-CoV-2. To characterize SARS-CoV-2 transmission, we used viral sequencing, phylogenomic tools, and functional assays. FINDINGS: Athletes, particularly those on high-contact teams, had the highest risk of testing positive. On average, individuals who tested positive had more contacts and longer interaction durations than individuals who never tested positive. The distribution of contacts per individual was overdispersed, although not as overdispersed as the distribution of phylogenomic descendants. Corroboration via technical replicates was essential for identification of wastewater mutations. CONCLUSIONS: Based on our findings, we formulate a framework that combines tools into an integrated disease surveillance program that can be implemented in other congregate settings with limited resources. FUNDING: This work was supported by the National Science Foundation, the Hertz Foundation, the National Institutes of Health, the Centers for Disease Control and Prevention, the Massachusetts Consortium on Pathogen Readiness, the Howard Hughes Medical Institute, the Flu Lab, and the Audacious Project.

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.

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.

Proteomic profiling dataset of chemical perturbations in multiple biological backgrounds.

While gene expression profiling has traditionally been the method of choice for large-scale perturbational profiling studies, proteomics has emerged as an effective tool in this context for directly monitoring cellular responses to perturbations. We previously reported a pilot library containing 3400 profiles of multiple perturbations across diverse cellular backgrounds in the reduced-representation phosphoproteome (P100) and chromatin space (Global Chromatin Profiling, GCP). Here, we expand our original dataset to include profiles from a new set of cardiotoxic compounds and from astrocytes, an additional neural cell model, totaling 5300 proteomic signatures. We describe filtering criteria and quality control metrics used to assess and validate the technical quality and reproducibility of our data. To demonstrate the power of the library, we present two case studies where data is queried using the concept of "connectivity" to obtain biological insight. All data presented in this study have been deposited to the ProteomeXchange Consortium with identifiers PXD017458 (P100) and PXD017459 (GCP) and can be queried at https://clue.io/proteomics .

The Origins and Future of Sentinel: An Early-Warning System for Pandemic Preemption and Response.

While investigating a signal of adaptive evolution in humans at the gene LARGE, we encountered an intriguing finding by Dr. Stefan Kunz that the gene plays a critical role in Lassa virus binding and entry. This led us to pursue field work to test our hypothesis that natural selection acting on LARGE-detected in the Yoruba population of Nigeria-conferred resistance to Lassa Fever in some West African populations. As we delved further, we conjectured that the "emerging" nature of recently discovered diseases like Lassa fever is related to a newfound capacity for detection, rather than a novel viral presence, and that humans have in fact been exposed to the viruses that cause such diseases for much longer than previously suspected. Dr. Stefan Kunz's critical efforts not only laid the groundwork for this discovery, but also inspired and catalyzed a series of events that birthed Sentinel, an ambitious and large-scale pandemic prevention effort in West Africa. Sentinel aims to detect and characterize deadly pathogens before they spread across the globe, through implementation of its three fundamental pillars: Detect, Connect, and Empower. More specifically, Sentinel is designed to detect known and novel infections rapidly, connect and share information in real time to identify emerging threats, and empower the public health community to improve pandemic preparedness and response anywhere in the world. We are proud to dedicate this work to Stefan Kunz, and eagerly invite new collaborators, experts, and others to join us in our efforts.

Phylogenetic analysis of SARS-CoV-2 in Boston highlights the impact of superspreading events.

Analysis of 772 complete severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) genomes from early in the Boston-area epidemic revealed numerous introductions of the virus, a small number of which led to most cases. The data revealed two superspreading events. One, in a skilled nursing facility, led to rapid transmission and significant mortality in this vulnerable population but little broader spread, whereas other introductions into the facility had little effect. The second, at an international business conference, produced sustained community transmission and was exported, resulting in extensive regional, national, and international spread. The two events also differed substantially in the genetic variation they generated, suggesting varying transmission dynamics in superspreading events. Our results show how genomic epidemiology can help to understand the link between individual clusters and wider community spread.

Avant-garde: an automated data-driven DIA data curation tool.

Several challenges remain in data-independent acquisition (DIA) data analysis, such as to confidently identify peptides, define integration boundaries, remove interferences, and control false discovery rates. In practice, a visual inspection of the signals is still required, which is impractical with large datasets. We present Avant-garde as a tool to refine DIA (and parallel reaction monitoring) data. Avant-garde uses a novel data-driven scoring strategy: signals are refined by learning from the dataset itself, using all measurements in all samples to achieve the best optimization. We evaluate the performance of Avant-garde using benchmark DIA datasets and show that it can determine the quantitative suitability of a peptide peak, and reach the same levels of selectivity, accuracy, and reproducibility as manual validation. Avant-garde is complementary to existing DIA analysis engines and aims to establish a strong foundation for subsequent analysis of quantitative mass spectrometry data.

Phylogenetic analysis of SARS-CoV-2 in the Boston area highlights the role of recurrent importation and superspreading events.

SARS-CoV-2 has caused a severe, ongoing outbreak of COVID-19 in Massachusetts with 111,070 confirmed cases and 8,433 deaths as of August 1, 2020. To investigate the introduction, spread, and epidemiology of COVID-19 in the Boston area, we sequenced and analyzed 772 complete SARS-CoV-2 genomes from the region, including nearly all confirmed cases within the first week of the epidemic and hundreds of cases from major outbreaks at a conference, a nursing facility, and among homeless shelter guests and staff. The data reveal over 80 introductions into the Boston area, predominantly from elsewhere in the United States and Europe. We studied two superspreading events covered by the data, events that led to very different outcomes because of the timing and populations involved. One produced rapid spread in a vulnerable population but little onward transmission, while the other was a major contributor to sustained community transmission, including outbreaks in homeless populations, and was exported to several other domestic and international sites. The same two events differed significantly in the number of new mutations seen, raising the possibility that SARS-CoV-2 superspreading might encompass disparate transmission dynamics. Our results highlight the failure of measures to prevent importation into MA early in the outbreak, underscore the role of superspreading in amplifying an outbreak in a major urban area, and lay a foundation for contact tracing informed by genetic data.

Chasing Tails: Cathepsin-L Improves Structural Analysis of Histones by HX-MS.

The N-terminal regions (tails) of histone proteins are dynamic elements that protrude from the nucleosome and are involved in many aspects of chromatin organization. Their epigenetic role is well-established, and post-translational modifications present on these regions contribute to transcriptional regulation. Considering their biological significance, relatively few structural details have been established for histone tails, mainly because of their inherently disordered nature. Although hydrogen/deuterium exchange mass spectrometry (HX-MS) is well-suited for the analysis of dynamic structures, it has seldom been employed in this context, presumably because of the poor N-terminal coverage provided by pepsin. Inspired from histone-clipping events, we profiled the activity of cathepsin-L under HX-MS quench conditions and characterized its specificity employing the four core histones (H2A, H2B, H3 and H4). Cathepsin-L demonstrated cleavage patterns that were substrate- and pH-dependent. Cathepsin-L generated overlapping N-terminal peptides about 20 amino acids long for H2A, H3, and H4 proving its suitability for the analysis of histone tails dynamics. We developed a comprehensive HX-MS method in combination with pepsin and obtained full sequence coverage for all histones. We employed our method to analyze histones H3 and H4. We observe rapid deuterium exchange of the N-terminal tails and cooperative unfolding (EX1 kinetics) in the histone-fold domains of histone monomers in-solution. Overall, this novel strategy opens new avenues for investigating the dynamic properties of histones that are not apparent from the crystal structures, providing insights into the structural basis of the histone code.

Tales of diversity: Genomic and morphological characteristics of forty-six Arthrobacter phages.

The vast bacteriophage population harbors an immense reservoir of genetic information. Almost 2000 phage genomes have been sequenced from phages infecting hosts in the phylum Actinobacteria, and analysis of these genomes reveals substantial diversity, pervasive mosaicism, and novel mechanisms for phage replication and lysogeny. Here, we describe the isolation and genomic characterization of 46 phages from environmental samples at various geographic locations in the U.S. infecting a single Arthrobacter sp. strain. These phages include representatives of all three virion morphologies, and Jasmine is the first sequenced podovirus of an actinobacterial host. The phages also span considerable sequence diversity, and can be grouped into 10 clusters according to their nucleotide diversity, and two singletons each with no close relatives. However, the clusters/singletons appear to be genomically well separated from each other, and relatively few genes are shared between clusters. Genome size varies from among the smallest of siphoviral phages (15,319 bp) to over 70 kbp, and G+C contents range from 45-68%, compared to 63.4% for the host genome. Although temperate phages are common among other actinobacterial hosts, these Arthrobacter phages are primarily lytic, and only the singleton Galaxy is likely temperate.