Method versatility in RNA extraction-free PCR detection of SARS-CoV-2 in saliva samples.

Early in the pandemic, a simple, open-source, RNA extraction-free RT-qPCR protocol for SARS-CoV-2 detection in saliva was developed and made widely available. This simplified approach (SalivaDirect) requires only sample treatment with proteinase K prior to PCR testing. However, feedback from clinical laboratories highlighted a need for a flexible workflow that can be seamlessly integrated into their current health and safety requirements for the receiving and handling of potentially infectious samples. To address these varying needs, we explored additional pre-PCR workflows. We built upon the original SalivaDirect workflow to include an initial incubation step (95 °C for 30 min, 95 °C for 5 min or 65 °C for 15 min) with or without addition of proteinase K. The limit of detection for the workflows tested did not significantly differ from that of the original SalivaDirect workflow. When tested on de-identified saliva samples from confirmed COVID-19 individuals, these workflows also produced comparable virus detection and assay sensitivities, as determined by RT-qPCR analysis. Exclusion of proteinase K did not negatively affect the sensitivity of the assay. The addition of multiple heat pretreatment options to the SalivaDirect protocol increases the accessibility of this cost-effective SARS-CoV-2 test as it gives diagnostic laboratories the flexibility to implement the workflow which best suits their safety protocols.

Nasal host response-based screening for undiagnosed respiratory viruses: a pathogen surveillance and detection study.

BACKGROUND: Symptomatic patients who test negative for common viruses are an important possible source of unrecognised or emerging pathogens, but metagenomic sequencing of all samples is inefficient because of the low likelihood of finding a pathogen in any given sample. We aimed to determine whether nasopharyngeal CXCL10 screening could be used as a strategy to enrich for samples containing undiagnosed viruses. METHODS: In this pathogen surveillance and detection study, we measured CXCL10 concentrations from nasopharyngeal swabs from patients in the Yale New Haven health-care system, which had been tested at the Yale New Haven Hospital Clinical Virology Laboratory (New Haven, CT, USA). Patients who tested negative for a panel of respiratory viruses using multiplex PCR during Jan 23-29, 2017, or March 3-14, 2020, were included. We performed host and pathogen RNA sequencing (RNA-Seq) and analysis for viral reads on samples with CXCL10 higher than 1 ng/mL or CXCL10 testing and quantitative RT-PCR (RT-qPCR) for SARS-CoV-2. We used RNA-Seq and cytokine profiling to compare the host response to infection in samples that were virus positive (rhinovirus, seasonal coronavirus CoV-NL63, or SARS-CoV-2) and virus negative (controls). FINDINGS: During Jan 23-29, 2017, 359 samples were tested for ten viruses on the multiplex PCR respiratory virus panel (RVP). 251 (70%) were RVP negative. 60 (24%) of 251 samples had CXCL10 higher than 150 pg/mL and were identified for further analysis. 28 (47%) of 60 CXCL10-high samples were positive for seasonal coronaviruses. 223 (89%) of 251 samples were PCR negative for 15 viruses and, of these, CXCL10-based screening identified 32 (13%) samples for further analysis. Of these 32 samples, eight (25%) with CXCL10 concentrations higher than 1 ng/mL and sufficient RNA were selected for RNA-Seq. Microbial RNA analysis showed the presence of influenza C virus in one sample and revealed RNA reads from bacterial pathobionts in four (50%) of eight samples. Between March 3 and March 14, 2020, 375 (59%) of 641 samples tested negative for 15 viruses on the RVP. 32 (9%) of 375 samples had CXCL10 concentrations ranging from 100 pg/mL to 1000 pg/mL and four of those were positive for SARS-CoV-2. CXCL10 elevation was statistically significant, and a distinguishing feature was found in 28 (8%) of 375 SARS-CoV-2-negative samples versus all four SARS-CoV-2-positive samples (p=4·4 × 10-5). Transcriptomic signatures showed an interferon response in virus-positive samples and an additional neutrophil-high hyperinflammatory signature in samples with high amounts of bacterial pathobionts. The CXCL10 cutoff for detecting a virus was 166·5 pg/mL for optimal sensitivity and 1091·0 pg/mL for specificity using a clinic-ready automated microfluidics-based immunoassay. INTERPRETATION: These results confirm CXCL10 as a robust nasopharyngeal biomarker of viral respiratory infection and support host response-based screening followed by metagenomic sequencing of CXCL10-high samples as a practical approach to incorporate clinical samples into pathogen discovery and surveillance efforts. FUNDING: National Institutes of Health, the Hartwell Foundation, the Gruber Foundation, Fast Grants for COVID-19 research from the Mercatus Center, and the Huffman Family Donor Advised Fund.

Lineage abundance estimation for SARS-CoV-2 in wastewater using transcriptome quantification techniques.

Effectively monitoring the spread of SARS-CoV-2 mutants is essential to efforts to counter the ongoing pandemic. Predicting lineage abundance from wastewater, however, is technically challenging. We show that by sequencing SARS-CoV-2 RNA in wastewater and applying algorithms initially used for transcriptome quantification, we can estimate lineage abundance in wastewater samples. We find high variability in signal among individual samples, but the overall trends match those observed from sequencing clinical samples. Thus, while clinical sequencing remains a more sensitive technique for population surveillance, wastewater sequencing can be used to monitor trends in mutant prevalence in situations where clinical sequencing is unavailable.

Combining genomic and epidemiological data to compare the transmissibility of SARS-CoV-2 variants Alpha and Iota.

SARS-CoV-2 variants shaped the second year of the COVID-19 pandemic and the discourse around effective control measures. Evaluating the threat posed by a new variant is essential for adapting response efforts when community transmission is detected. In this study, we compare the dynamics of two variants, Alpha and Iota, by integrating genomic surveillance data to estimate the effective reproduction number (Rt) of the variants. We use Connecticut, United States, in which Alpha and Iota co-circulated in 2021. We find that the Rt of these variants were up to 50% larger than that of other variants. We then use phylogeography to show that while both variants were introduced into Connecticut at comparable frequencies, clades that resulted from introductions of Alpha were larger than those resulting from Iota introductions. By monitoring the dynamics of individual variants throughout our study period, we demonstrate the importance of routine surveillance in the response to COVID-19.

Sequencing SARS-CoV-2 genomes from saliva.

Genomic sequencing is crucial to understanding the epidemiology and evolution of Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). Often, genomic studies rely on remnant diagnostic material, typically nasopharyngeal (NP) swabs, as input into whole-genome SARS-CoV-2 next-generation sequencing pipelines. Saliva has proven to be a safe and stable specimen for the detection of SARS-CoV-2 RNA via traditional diagnostic assays; however, saliva is not commonly used for SARS-CoV-2 sequencing. Using the ARTIC Network amplicon-generation approach with sequencing on the Oxford Nanopore MinION, we demonstrate that sequencing SARS-CoV-2 from saliva produces genomes comparable to those from NP swabs, and that RNA extraction is necessary to generate complete genomes from saliva. In this study, we show that saliva is a useful specimen type for genomic studies of SARS-CoV-2.

An outbreak of SARS-CoV-2 on a transplant unit in the early vaccination era.

BACKGROUND: Solid organ transplant recipients are at increased risk of COVID-19-associated morbidity and mortality. AIMS: We describe a nosocomial outbreak investigation on an immunocompromised inpatient unit. METHODS: Patients positive for SARS-CoV-2 were identified. An epidemiologic investigation was assisted with whole genome sequencing of positive samples. RESULTS: Two patients were identified as potential index cases; one presented with diarrhea and was initially not isolated, and the other developed hypoxemia on hospital day 18 before testing positive. Following identification of a SARS-CoV-2 cluster, the unit was closed and all patients and staff received surveillance testing revealing eight additional positive patients and staff members. Whole genome sequencing confirmed an outbreak. Enhanced infection prevention practices mitigated further spread. Asymptomatic patients with COVID-19 were successfully treated with bamlanivimab. DISCUSSION: Preventing SARS-CoV-2 outbreaks in transplant units poses unique challenges as patients may have atypical presentations of COVID-19. Immunocompromised patients who test positive for SARS-CoV-2 while asymptomatic may benefit from monoclonal antibody therapy to prevent disease progression. All hospital staff members working with immunocompromised patients should be promptly encouraged to follow infection prevention behaviors and receive SARS-CoV-2 vaccination. CONCLUSION: SARS-CoV-2 outbreaks on immunocompromised units can be mitigated through prompt identification of cases and robust infection prevention practices.

Longitudinal Immune Profiling of a Severe Acute Respiratory Syndrome Coronavirus 2 Reinfection in a Solid Organ Transplant Recipient.

BACKGROUND: The underlying immunologic deficiencies enabling severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) reinfection are currently unknown. We describe deep longitudinal immune profiling of a transplant recipient hospitalized twice for coronavirus disease 2019 (COVID-19). METHODS: A 66-year-old male renal transplant recipient was hospitalized with COVID-19 March 2020 then readmitted to the hospital with COVID-19 233 days after initial diagnosis. Virologic and immunologic investigations were performed on samples from the primary and secondary infections. RESULTS: Whole viral genome sequencing and phylogenetic analysis revealed that viruses causing both infections were caused by distinct genetic lineages without evidence of immune escape mutations. Longitudinal comparison of cellular and humoral responses during primary SARS-CoV-2 infection revealed that this patient responded to the primary infection with low neutralization titer anti-SARS-CoV-2 antibodies that were likely present at the time of reinfection. CONCLUSIONS: The development of neutralizing antibodies and humoral memory responses in this patient failed to confer protection against reinfection, suggesting that they were below a neutralizing titer threshold or that additional factors may be required for efficient prevention of SARS-CoV-2 reinfection. Development of poorly neutralizing antibodies may have been due to profound and relatively specific reduction in naive CD4 T-cell pools. Seropositivity alone may not be a perfect correlate of protection in immunocompromised patients.

Multiple Transmission Chains within COVID-19 Cluster, Connecticut, USA, 20201.

In fall 2020, a coronavirus disease cluster comprising 16 cases occurred in Connecticut, USA. Epidemiologic and genomic evidence supported transmission among persons at a school and fitness center but not a workplace. The multiple transmission chains identified within this cluster highlight the necessity of a combined investigatory approach.

Variant abundance estimation for SARS-CoV-2 in wastewater using RNA-Seq quantification.

Effectively monitoring the spread of SARS-CoV-2 variants is essential to efforts to counter the ongoing pandemic. Wastewater monitoring of SARS-CoV-2 RNA has proven an effective and efficient technique to approximate COVID-19 case rates in the population. Predicting variant abundances from wastewater, however, is technically challenging. Here we show that by sequencing SARS-CoV-2 RNA in wastewater and applying computational techniques initially used for RNA-Seq quantification, we can estimate the abundance of variants in wastewater samples. We show by sequencing samples from wastewater and clinical isolates in Connecticut U.S.A. between January and April 2021 that the temporal dynamics of variant strains broadly correspond. We further show that this technique can be used with other wastewater sequencing techniques by expanding to samples taken across the United States in a similar timeframe. We find high variability in signal among individual samples, and limited ability to detect the presence of variants with clinical frequencies <10%; nevertheless, the overall trends match what we observed from sequencing clinical samples. Thus, while clinical sequencing remains a more sensitive technique for population surveillance, wastewater sequencing can be used to monitor trends in variant prevalence in situations where clinical sequencing is unavailable or impractical.

Sequencing SARS-CoV-2 Genomes from Saliva.

Genomic sequencing is crucial to understanding the epidemiology and evolution of SARS-CoV-2. Often, genomic studies rely on remnant diagnostic material, typically nasopharyngeal swabs, as input into whole genome SARS-CoV-2 next-generation sequencing pipelines. Saliva has proven to be a safe and stable specimen for the detection of SARS-CoV-2 RNA via traditional diagnostic assays, however saliva is not commonly used for SARS-CoV-2 sequencing. Using the ARTIC Network amplicon-generation approach with sequencing on the Oxford Nanopore MinION, we demonstrate that sequencing SARS-CoV-2 from saliva produces genomes comparable to those from nasopharyngeal swabs, and that RNA extraction is necessary to generate complete genomes from saliva. In this study, we show that saliva is a useful specimen type for genomic studies of SARS-CoV-2.

Combining genomic and epidemiological data to compare the transmissibility of SARS-CoV-2 lineages.

Emerging SARS-CoV-2 variants have shaped the second year of the COVID-19 pandemic and the public health discourse around effective control measures. Evaluating the public health threat posed by a new variant is essential for appropriately adapting response efforts when community transmission is detected. However, this assessment requires that a true comparison can be made between the new variant and its predecessors because factors other than the virus genotype may influence spread and transmission. In this study, we develop a framework that integrates genomic surveillance data to estimate the relative effective reproduction number (R t ) of co-circulating lineages. We use Connecticut, a state in the northeastern United States in which the SARS-CoV-2 variants B.1.1.7 and B.1.526 co-circulated in early 2021, as a case study for implementing this framework. We find that the R t of B.1.1.7 was 6-10% larger than that of B.1.526 in Connecticut in the midst of a COVID-19 vaccination campaign. To assess the generalizability of this framework, we apply it to genomic surveillance data from New York City and observe the same trend. Finally, we use discrete phylogeography to demonstrate that while both variants were introduced into Connecticut at comparable frequencies, clades that resulted from introductions of B.1.1.7 were larger than those resulting from B.1.526 introductions. Our framework, which uses open-source methods requiring minimal computational resources, may be used to monitor near real-time variant dynamics in a myriad of settings.

Multiplex qPCR discriminates variants of concern to enhance global surveillance of SARS-CoV-2.

With the emergence of Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) variants that may increase transmissibility and/or cause escape from immune responses, there is an urgent need for the targeted surveillance of circulating lineages. It was found that the B.1.1.7 (also 501Y.V1) variant, first detected in the United Kingdom, could be serendipitously detected by the Thermo Fisher TaqPath COVID-19 PCR assay because a key deletion in these viruses, spike Δ69-70, would cause a "spike gene target failure" (SGTF) result. However, a SGTF result is not definitive for B.1.1.7, and this assay cannot detect other variants of concern (VOC) that lack spike Δ69-70, such as B.1.351 (also 501Y.V2), detected in South Africa, and P.1 (also 501Y.V3), recently detected in Brazil. We identified a deletion in the ORF1a gene (ORF1a Δ3675-3677) in all 3 variants, which has not yet been widely detected in other SARS-CoV-2 lineages. Using ORF1a Δ3675-3677 as the primary target and spike Δ69-70 to differentiate, we designed and validated an open-source PCR assay to detect SARS-CoV-2 VOC. Our assay can be rapidly deployed in laboratories around the world to enhance surveillance for the local emergence and spread of B.1.1.7, B.1.351, and P.1.

Longitudinal immune profiling of a SARS-CoV-2 reinfection in a solid organ transplant recipient.

The underlying immunologic deficiencies enabling SARS-CoV-2 reinfections are currently unknown. Here we describe a renal-transplant recipient who developed recurrent, symptomatic SARS-CoV-2 infection 7 months after primary infection. To elucidate the immunological mechanisms responsible for reinfection, we performed longitudinal profiling of cellular and humoral responses during both primary and recurrent SARS-CoV-2 infection. We found that the patient responded to the primary infection with transient, poor-quality adaptive immune responses that was further compromised by intervening treatment for acute rejection of the renal allograft prior to reinfection. Importantly, we identified the development of neutralizing antibodies and humoral memory responses prior to SARS-CoV-2 reinfection. However, these neutralizing antibodies failed to confer protection against reinfection, suggesting that additional factors are required for efficient prevention of SARS-CoV-2 reinfection. Further, we found no evidence supporting viral evasion of primary adaptive immune responses, suggesting that susceptibility to reinfection may be determined by host factors rather than pathogen adaptation.

Case Study: Longitudinal immune profiling of a SARS-CoV-2 reinfection in a solid organ transplant recipient.

Prior to the emergence of antigenically distinct SARS-CoV-2 variants, reinfections were reported infrequently - presumably due to the generation of durable and protective immune responses. However, case reports also suggested that rare, repeated infections may occur as soon as 48 days following initial disease onset. The underlying immunologic deficiencies enabling SARS-CoV-2 reinfections are currently unknown. Here we describe a renal transplant recipient who developed recurrent, symptomatic SARS-CoV-2 infection - confirmed by whole virus genome sequencing - 7 months after primary infection. To elucidate the immunological mechanisms responsible for SARS-CoV-2 reinfection, we performed longitudinal profiling of cellular and humoral responses during both primary and recurrent SARS-CoV-2 infection. We found that the patient responded to the primary infection with transient, poor-quality adaptive immune responses. The patient's immune system was further compromised by intervening treatment for acute rejection of the renal allograft prior to reinfection. Importantly, we also identified the development of neutralizing antibodies and the formation of humoral memory responses prior to SARS-CoV-2 reinfection. However, these neutralizing antibodies failed to confer protection against reinfection, suggesting that additional factors are required for efficient prevention of SARS-CoV-2 reinfection. Further, we found no evidence supporting viral evasion of primary adaptive immune responses, suggesting that susceptibility to reinfection may be determined by host factors rather than pathogen adaptation in this patient. In summary, our study suggests that a low neutralizing antibody presence alone is not sufficient to confer resistance against reinfection. Thus, patients with solid organ transplantation, or patients who are otherwise immunosuppressed, who recover from infection with SARS-CoV-2 may not develop sufficient protective immunity and are at risk of reinfection.

Early introductions and transmission of SARS-CoV-2 variant B.1.1.7 in the United States.

The emergence and spread of SARS-CoV-2 lineage B.1.1.7, first detected in the United Kingdom, has become a global public health concern because of its increased transmissibility. Over 2,500 COVID-19 cases associated with this variant have been detected in the United States (US) since December 2020, but the extent of establishment is relatively unknown. Using travel, genomic, and diagnostic data, we highlight that the primary ports of entry for B.1.1.7 in the US were in New York, California, and Florida. Furthermore, we found evidence for many independent B.1.1.7 establishments starting in early December 2020, followed by interstate spread by the end of the month. Finally, we project that B.1.1.7 will be the dominant lineage in many states by mid- to late March. Thus, genomic surveillance for B.1.1.7 and other variants urgently needs to be enhanced to better inform the public health response.

PCR assay to enhance global surveillance for SARS-CoV-2 variants of concern.

With the emergence of SARS-CoV-2 variants that may increase transmissibility and/or cause escape from immune responses 1-3 , there is an urgent need for the targeted surveillance of circulating lineages. It was found that the B.1.1.7 (also 501Y.V1) variant first detected in the UK 4,5 could be serendipitously detected by the ThermoFisher TaqPath COVID-19 PCR assay because a key deletion in these viruses, spike Δ69-70, would cause a "spike gene target failure" (SGTF) result. However, a SGTF result is not definitive for B.1.1.7, and this assay cannot detect other variants of concern that lack spike Δ69-70, such as B.1.351 (also 501Y.V2) detected in South Africa 6 and P.1 (also 501Y.V3) recently detected in Brazil 7 . We identified a deletion in the ORF1a gene (ORF1a Δ3675-3677) in all three variants, which has not yet been widely detected in other SARS-CoV-2 lineages. Using ORF1a Δ3675-3677 as the primary target and spike Δ69-70 to differentiate, we designed and validated an open source PCR assay to detect SARS-CoV-2 variants of concern 8 . Our assay can be rapidly deployed in laboratories around the world to enhance surveillance for the local emergence spread of B.1.1.7, B.1.351, and P.1.

Early introductions and community transmission of SARS-CoV-2 variant B.1.1.7 in the United States.

The emergence and spread of SARS-CoV-2 lineage B.1.1.7, first detected in the United Kingdom, has become a global public health concern because of its increased transmissibility. Over 2500 COVID-19 cases associated with this variant have been detected in the US since December 2020, but the extent of establishment is relatively unknown. Using travel, genomic, and diagnostic data, we highlight the primary ports of entry for B.1.1.7 in the US and locations of possible underreporting of B.1.1.7 cases. Furthermore, we found evidence for many independent B.1.1.7 establishments starting in early December 2020, followed by interstate spread by the end of the month. Finally, we project that B.1.1.7 will be the dominant lineage in many states by mid to late March. Thus, genomic surveillance for B.1.1.7 and other variants urgently needs to be enhanced to better inform the public health response.

COVID-19 Outcomes and Genomic Characterization of SARS-CoV-2 Isolated From Veterans in New England States: Retrospective Analysis.

BACKGROUND: Clinical and virologic characteristics of COVID-19 infections in veterans in New England have not been described. The average US veteran is a male older than the general US population. SARS-CoV-2 infection is known to cause poorer outcomes among men and older adults, making the veteran population an especially vulnerable group for COVID-19. OBJECTIVE: This study aims to evaluate clinical and virologic factors impacting COVID-19 outcomes. METHODS: This retrospective chart review included 476 veterans in six New England states with confirmed SARS-CoV-2 infection between April and September 2020. Whole genome sequencing was performed on SARS-CoV-2 RNA isolated from these veterans, and the correlation of genomic data to clinical outcomes was evaluated. Clinical and demographic variables were collected by manual chart review and were correlated to the end points of peak disease severity (based on oxygenation requirements), hospitalization, and mortality using multivariate regression analyses. RESULTS: Of 476 veterans, 274 had complete and accessible charts. Of the 274 veterans, 92.7% (n=254) were men and 83.2% (n=228) were White, and the mean age was 63 years. In the multivariate regression, significant predictors of hospitalization (C statistic 0.75) were age (odds ratio [OR] 1.05, 95% CI 1.03-1.08) and non-White race (OR 2.39, 95% CI 1.13-5.01). Peak severity (C statistic 0.70) also varied by age (OR 1.07, 95% CI 1.03-1.11) and O2 requirement on admission (OR 45.7, 95% CI 18.79-111). Mortality (C statistic 0.87) was predicted by age (OR 1.06, 95% CI 1.01-1.11), dementia (OR 3.44, 95% CI 1.07-11.1), and O2 requirement on admission (OR 6.74, 95% CI 1.74-26.1). Most (291/299, 97.3%) of our samples were dominated by the spike protein D614G substitution and were from SARS-CoV-2 B.1 lineage or one of 37 different B.1 sublineages, with none representing more than 8.7% (26/299) of the cases. CONCLUSIONS: In a cohort of veterans from the six New England states with a mean age of 63 years and a high comorbidity burden, age was the largest predictor of hospitalization, peak disease severity, and mortality. Non-White veterans were more likely to be hospitalized, and patients who required oxygen on admission were more likely to have severe disease and higher rates of mortality. Multiple SARS-CoV-2 lineages were distributed in patients in New England early in the COVID-19 era, mostly related to viruses from New York State with D614G mutation.

Widespread Transcriptional Readthrough Caused by Nab2 Depletion Leads to Chimeric Transcripts with Retained Introns.

Nascent RNA sequencing has revealed that pre-mRNA splicing can occur shortly after introns emerge from RNA polymerase II (RNA Pol II). Differences in co-transcriptional splicing profiles suggest regulation by cis- and/or trans-acting factors. Here, we use single-molecule intron tracking (SMIT) to identify a cohort of regulators by machine learning in budding yeast. Of these, Nab2 displays reduced co-transcriptional splicing when depleted. Unexpectedly, these splicing defects are attributable to aberrant "intrusive" transcriptional readthrough from upstream genes, as revealed by long-read sequencing. Transcripts that originate from the intron-containing gene's own transcription start site (TSS) are efficiently spliced, indicating no direct role of Nab2 in splicing per se. This work highlights the coupling between transcription, splicing, and 3' end formation in the context of gene organization along chromosomes. We conclude that Nab2 is required for proper 3' end processing, which ensures gene-specific control of co-transcriptional RNA processing.