SARS-CoV-2 lineage assignments using phylogenetic placement/UShER are superior to pangoLEARN machine-learning method.

With the rapid spread and evolution of SARS-CoV-2, the ability to monitor its transmission and distinguish among viral lineages is critical for pandemic response efforts. The most commonly used software for the lineage assignment of newly isolated SARS-CoV-2 genomes is pangolin, which offers two methods of assignment, pangoLEARN and pUShER. PangoLEARN rapidly assigns lineages using a machine-learning algorithm, while pUShER performs a phylogenetic placement to identify the lineage corresponding to a newly sequenced genome. In a preliminary study, we observed that pangoLEARN (decision tree model), while substantially faster than pUShER, offered less consistency across different versions of pangolin v3. Here, we expand upon this analysis to include v3 and v4 of pangolin, which moved the default algorithm for lineage assignment from pangoLEARN in v3 to pUShER in v4, and perform a thorough analysis confirming that pUShER is not only more stable across versions but also more accurate. Our findings suggest that future lineage assignment algorithms for various pathogens should consider the value of phylogenetic placement.

Implementation of California COVIDNet - a multi-sector collaboration for statewide SARS-CoV-2 genomic surveillance.

Introduction: The SARS-CoV-2 pandemic represented a formidable scientific and technological challenge to public health due to its rapid spread and evolution. To meet these challenges and to characterize the virus over time, the State of California established the California SARS-CoV-2 Whole Genome Sequencing (WGS) Initiative, or "California COVIDNet". This initiative constituted an unprecedented multi-sector collaborative effort to achieve large-scale genomic surveillance of SARS-CoV-2 across California to monitor the spread of variants within the state, to detect new and emerging variants, and to characterize outbreaks in congregate, workplace, and other settings. Methods: California COVIDNet consists of 50 laboratory partners that include public health laboratories, private clinical diagnostic laboratories, and academic sequencing facilities as well as expert advisors, scientists, consultants, and contractors. Data management, sample sourcing and processing, and computational infrastructure were major challenges that had to be resolved in the midst of the pandemic chaos in order to conduct SARS-CoV-2 genomic surveillance. Data management, storage, and analytics needs were addressed with both conventional database applications and newer cloud-based data solutions, which also fulfilled computational requirements. Results: Representative and randomly selected samples were sourced from state-sponsored community testing sites. Since March of 2021, California COVIDNet partners have contributed more than 450,000 SARS-CoV-2 genomes sequenced from remnant samples from both molecular and antigen tests. Combined with genomes from CDC-contracted WGS labs, there are currently nearly 800,000 genomes from all 61 local health jurisdictions (LHJs) in California in the COVIDNet sequence database. More than 5% of all reported positive tests in the state have been sequenced, with similar rates of sequencing across 5 major geographic regions in the state. Discussion: Implementation of California COVIDNet revealed challenges and limitations in the public health system. These were overcome by engaging in novel partnerships that established a successful genomic surveillance program which provided valuable data to inform the COVID-19 public health response in California. Significantly, California COVIDNet has provided a foundational data framework and computational infrastructure needed to respond to future public health crises.

Trends of Enterovirus D68 Concentrations in Wastewater, California, USA, February 2021-April 2023.

In this retrospective study, we measured enterovirus D68 (EV-D68) genomic RNA in wastewater solids longitudinally at 2 California, USA, wastewater treatment plants twice per week for 26 months. EV-D68 RNA was undetectable except when concentrations increased from mid-July to mid-December 2022, which coincided with a peak in confirmed EV-D68 cases.

Pathogen genomics in public health laboratories: successes, challenges, and lessons learned from California's SARS-CoV-2 Whole-Genome Sequencing Initiative, California COVIDNet.

The capacity for pathogen genomics in public health expanded rapidly during the coronavirus disease 2019 (COVID-19) pandemic, but many public health laboratories did not have the infrastructure in place to handle the vast amount of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) sequence data generated. The California Department of Public Health, in partnership with Theiagen Genomics, was an early adopter of cloud-based resources for bioinformatics and genomic epidemiology, resulting in the creation of a SARS-CoV-2 genomic surveillance system that combined the efforts of more than 40 sequencing laboratories across government, academia and industry to form California COVIDNet, California's SARS-CoV-2 Whole-Genome Sequencing Initiative. Open-source bioinformatics workflows, ongoing training sessions for the public health workforce, and automated data transfer to visualization tools all contributed to the success of California COVIDNet. While challenges remain for public health genomic surveillance worldwide, California COVIDNet serves as a framework for a scaled and successful bioinformatics infrastructure that has expanded beyond SARS-CoV-2 to other pathogens of public health importance.

Severe Acute Respiratory Syndrome Coronavirus 2 Delta Variant Genomic Variation Associated With Breakthrough Infection in Northern California: A Retrospective Cohort Study.

BACKGROUND: The association between severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) genomic variation and breakthrough infection is not well defined among persons with Delta variant SARS-CoV-2 infection. METHODS: In a retrospective cohort, we assessed whether individual nonlineage defining mutations and overall genomic variation (including low-frequency alleles) were associated with breakthrough infection, defined as SARS-CoV-2 infection after coronavirus disease 2019 primary vaccine series. We identified all nonsynonymous single-nucleotide polymorphisms, insertions, and deletions in SARS-CoV-2 genomes with ≥5% allelic frequency and population frequency of ≥5% and ≤95%. Using Poisson regression, we assessed the association with breakthrough infection for each individual mutation and a viral genomic risk score. RESULTS: Thirty-six mutations met our inclusion criteria. Among 12 744 persons infected with Delta variant SARS-CoV-2, 5949 (47%) were vaccinated and 6795 (53%) were unvaccinated. Viruses with a viral genomic risk score in the highest quintile were 9% more likely to be associated with breakthrough infection than viruses in the lowest quintile, but including the risk score improved overall predictive model performance (measured by C statistic) by only +0.0006. CONCLUSIONS: Genomic variation within SARS-CoV-2 Delta variant was weakly associated with breakthrough infection, but several potential nonlineage defining mutations were identified that might contribute to immune evasion by SARS-CoV-2.

Adeno-associated virus type 2 in US children with acute severe hepatitis.

As of August 2022, clusters of acute severe hepatitis of unknown aetiology in children have been reported from 35 countries, including the USA1,2. Previous studies have found human adenoviruses (HAdVs) in the blood from patients in Europe and the USA3-7, although it is unclear whether this virus is causative. Here we used PCR testing, viral enrichment-based sequencing and agnostic metagenomic sequencing to analyse samples from 16 HAdV-positive cases from 1 October 2021 to 22 May 2022, in parallel with 113 controls. In blood from 14 cases, adeno-associated virus type 2 (AAV2) sequences were detected in 93% (13 of 14), compared to 4 (3.5%) of 113 controls (P < 0.001) and to 0 of 30 patients with hepatitis of defined aetiology (P < 0.001). In controls, HAdV type 41 was detected in blood from 9 (39.1%) of the 23 patients with acute gastroenteritis (without hepatitis), including 8 of 9 patients with positive stool HAdV testing, but co-infection with AAV2 was observed in only 3 (13.0%) of these 23 patients versus 93% of cases (P < 0.001). Co-infections by Epstein-Barr virus, human herpesvirus 6 and/or enterovirus A71 were also detected in 12 (85.7%) of 14 cases, with higher herpesvirus detection in cases versus controls (P < 0.001). Our findings suggest that the severity of the disease is related to co-infections involving AAV2 and one or more helper viruses.

Travel history among persons infected with SARS-CoV-2 variants of concern in the United States, December 2020-February 2021.

The first three SARS-CoV-2 phylogenetic lineages classified as variants of concern (VOCs) in the United States (U.S.) from December 15, 2020 to February 28, 2021, Alpha (B.1.1.7), Beta (B.1.351), and Gamma (P.1) lineages, were initially detected internationally. This investigation examined available travel history of coronavirus disease 2019 (COVID-19) cases reported in the U.S. in whom laboratory testing showed one of these initial VOCs. Travel history, demographics, and health outcomes for a convenience sample of persons infected with a SARS-CoV-2 VOC from December 15, 2020 through February 28, 2021 were provided by 35 state and city health departments, and proportion reporting travel was calculated. Of 1,761 confirmed VOC cases analyzed, 1,368 had available data on travel history. Of those with data on travel history, 1,168 (85%) reported no travel preceding laboratory confirmation of SARS-CoV-2 and only 105 (8%) reported international travel during the 30 days preceding a positive SARS-CoV-2 test or symptom onset. International travel was reported by 92/1,304 (7%) of persons infected with the Alpha variant, 7/55 (22%) with Beta, and 5/9 (56%) with Gamma. Of the first three SARS-CoV-2 lineages designated as VOCs in the U.S., international travel was common only among the few Gamma cases. Most persons infected with Alpha and Beta variant reported no travel history, therefore, community transmission of these VOCs was likely common in the U.S. by March 2021. These findings underscore the importance of global surveillance using whole genome sequencing to detect and inform mitigation strategies for emerging SARS-CoV-2 VOCs.

Multiple Mutations Associated with Emergent Variants Can Be Detected as Low-Frequency Mutations in Early SARS-CoV-2 Pandemic Clinical Samples.

Genetic analysis of intra-host viral populations provides unique insight into pre-emergent mutations that may contribute to the genotype of future variants. Clinical samples positive for SARS-CoV-2 collected in California during the first months of the pandemic were sequenced to define the dynamics of mutation emergence as the virus became established in the state. Deep sequencing of 90 nasopharyngeal samples showed that many mutations associated with the establishment of SARS-CoV-2 globally were present at varying frequencies in a majority of the samples, even those collected as the virus was first detected in the US. A subset of mutations that emerged months later in consensus sequences were detected as subconsensus members of intra-host populations. Spike mutations P681H, H655Y, and V1104L were detected prior to emergence in variant genotypes, mutations were detected at multiple positions within the furin cleavage site, and pre-emergent mutations were identified in the nucleocapsid and the envelope genes. Because many of the samples had a very high depth of coverage, a bioinformatics pipeline, "Mappgene", was established that uses both iVar and LoFreq variant calling to enable identification of very low-frequency variants. This enabled detection of a spike protein deletion present in many samples at low frequency and associated with a variant of concern.

Evaluation of co-circulating pathogens and microbiome from COVID-19 infections.

Co-infections or secondary infections with SARS-CoV-2 have the potential to affect disease severity and morbidity. Additionally, the potential influence of the nasal microbiome on COVID-19 illness is not well understood. In this study, we analyzed 203 residual samples, originally submitted for SARS-CoV-2 testing, for the presence of viral, bacterial, and fungal pathogens and non-pathogens using a comprehensive microarray technology, the Lawrence Livermore Microbial Detection Array (LLMDA). Eighty-seven percent of the samples were nasopharyngeal samples, and 23% of the samples were oral, nasal and oral pharyngeal swabs. We conducted bioinformatics analyses to examine differences in microbial populations of these samples, as a proxy for the nasal and oral microbiome, from SARS-CoV-2 positive and negative specimens. We found 91% concordance with the LLMDA relative to a diagnostic RT-qPCR assay for detection of SARS-CoV-2. Sixteen percent of all the samples (32/203) revealed the presence of an opportunistic bacterial or frank viral pathogen with the potential to cause co-infections. The two most detected bacteria, Streptococcus pyogenes and Streptococcus pneumoniae, were present in both SARS-CoV-2 positive and negative samples. Human metapneumovirus was the most prevalent viral pathogen in the SARS-CoV-2 negative samples. Sequence analysis of 16S rRNA was also conducted to evaluate bacterial diversity and confirm LLMDA results.

Neutralizing Immunity Induced Against the Omicron BA.1 and BA.2 Variants in Vaccine Breakthrough Infections.

BACKGROUND: As of early 2022, the Omicron variants are the predominant circulating lineages globally. Understanding neutralizing antibody responses against Omicron BA.1 and BA.2 after vaccine breakthrough infections will provide insights into BA.2 infectivity and susceptibility to subsequent reinfection. METHODS: Live virus neutralization assays were used to study immunity against Delta and Omicron BA.1 and BA.2 variants in samples from 86 individuals, 24 unvaccinated (27.9%) and 62 vaccinated (72.1%), who were infected with Delta (n = 42, 48.8%) or BA.1 (n = 44, 51.2%). Among the 62 vaccinated individuals, 39 were unboosted (62.9%), whereas 23 were boosted (37.1%). RESULTS: In unvaccinated infections, neutralizing antibodies (nAbs) against the three variants were weak or undetectable, except against Delta for Delta-infected individuals. Both Delta and BA.1 breakthrough infections resulted in strong nAb responses against ancestral wild-type and Delta lineages, but moderate nAb responses against BA.1 and BA.2, with similar titers between unboosted and boosted individuals. Antibody titers against BA.2 were generally higher than those against BA.1 in breakthrough infections. CONCLUSIONS: These results underscore the decreased immunogenicity of BA.1 compared to BA.2, insufficient neutralizing immunity against BA.2 in unvaccinated individuals, and moderate to strong neutralizing immunity induced against BA.2 in Delta and BA.1 breakthrough infections.

Limited cross-variant immunity from SARS-CoV-2 Omicron without vaccination.

SARS-CoV-2 Delta and Omicron are globally relevant variants of concern. Although individuals infected with Delta are at risk of developing severe lung disease, infection with Omicron often causes milder symptoms, especially in vaccinated individuals1,2. The question arises of whether widespread Omicron infections could lead to future cross-variant protection, accelerating the end of the pandemic. Here we show that without vaccination, infection with Omicron induces a limited humoral immune response in mice and humans. Sera from mice overexpressing the human ACE2 receptor and infected with Omicron neutralize only Omicron, but not other variants of concern, whereas broader cross-variant neutralization was observed after WA1 and Delta infections. Unlike WA1 and Delta, Omicron replicates to low levels in the lungs and brains of infected animals, leading to mild disease with reduced expression of pro-inflammatory cytokines and diminished activation of lung-resident T cells. Sera from individuals who were unvaccinated and infected with Omicron show the same limited neutralization of only Omicron itself. By contrast, Omicron breakthrough infections induce overall higher neutralization titres against all variants of concern. Our results demonstrate that Omicron infection enhances pre-existing immunity elicited by vaccines but, on its own, may not confer broad protection against non-Omicron variants in unvaccinated individuals.

Neutralizing immunity in vaccine breakthrough infections from the SARS-CoV-2 Omicron and Delta variants.

Virus-like particle (VLP) and live virus assays were used to investigate neutralizing immunity against Delta and Omicron SARS-CoV-2 variants in 259 samples from 128 vaccinated individuals. Following Delta breakthrough infection, titers against WT rose 57-fold and 3.1-fold compared with uninfected boosted and unboosted individuals, respectively, versus only a 5.8-fold increase and 3.1-fold decrease for Omicron breakthrough infection. Among immunocompetent, unboosted patients, Delta breakthrough infections induced 10.8-fold higher titers against WT compared with Omicron (p = 0.037). Decreased antibody responses in Omicron breakthrough infections relative to Delta were potentially related to a higher proportion of asymptomatic or mild breakthrough infections (55.0% versus 28.6%, respectively), which exhibited 12.3-fold lower titers against WT compared with moderate to severe infections (p = 0.020). Following either Delta or Omicron breakthrough infection, limited variant-specific cross-neutralizing immunity was observed. These results suggest that Omicron breakthrough infections are less immunogenic than Delta, thus providing reduced protection against reinfection or infection from future variants.

Predominance of antibody-resistant SARS-CoV-2 variants in vaccine breakthrough cases from the San Francisco Bay Area, California.

Associations between vaccine breakthrough cases and infection by different SARS coronavirus 2 (SARS-CoV-2) variants have remained largely unexplored. Here we analysed SARS-CoV-2 whole-genome sequences and viral loads from 1,373 persons with COVID-19 from the San Francisco Bay Area from 1 February to 30 June 2021, of which 125 (9.1%) were vaccine breakthrough infections. Vaccine breakthrough infections were more commonly associated with circulating antibody-resistant variants carrying ≥1 mutation associated with decreased antibody neutralization (L452R/Q, E484K/Q and/or F490S) than infections in unvaccinated individuals (78% versus 48%, P = 1.96 × 10-8). Differences in viral loads were non-significant between unvaccinated and fully vaccinated cases overall (P = 0.99) and according to lineage (P = 0.09-0.78). Symptomatic vaccine breakthrough infections had comparable viral loads (P = 0.64), whereas asymptomatic breakthrough infections had decreased viral loads (P = 0.023) compared with infections in unvaccinated individuals. In 5 cases with serial samples available for serologic analyses, vaccine breakthrough infections were found to be associated with low or undetectable neutralizing antibody levels attributable to an immunocompromised state or infection by an antibody-resistant lineage. Taken together, our results show that vaccine breakthrough infections are overrepresented by antibody-resistant SARS-CoV-2 variants, and that symptomatic breakthrough infections may be as efficient in spreading COVID-19 as unvaccinated infections, regardless of the infecting lineage.

Severe Acute Respiratory Syndrome Coronavirus 2 and Respiratory Virus Sentinel Surveillance, California, USA, May 10, 2020-June 12, 2021.

State and local health departments established the California Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) and Respiratory Virus Sentinel Surveillance System to conduct enhanced surveillance for SARS-CoV-2 and other respiratory pathogens at sentinel outpatient testing sites in 10 counties throughout California, USA. We describe results obtained during May 10, 2020‒June 12, 2021, and compare persons with positive and negative SARS-CoV-2 PCR results by using Poisson regression. We detected SARS-CoV-2 in 1,696 (19.6%) of 8,662 specimens. Among 7,851 specimens tested by respiratory panel, rhinovirus/enterovirus was detected in 906 (11.5%) specimens and other respiratory pathogens in 136 (1.7%) specimens. We also detected 23 co-infections with SARS-CoV-2 and another pathogen. SARS-CoV-2 positivity was associated with male participants, an age of 35-49 years, Latino race/ethnicity, obesity, and work in transportation occupations. Sentinel surveillance can provide useful virologic and epidemiologic data to supplement other disease monitoring activities and might become increasingly useful as routine testing decreases.

Effectiveness of Abbott BinaxNOW Rapid Antigen Test for Detection of SARS-CoV-2 Infections in Outbreak among Horse Racetrack Workers, California, USA.

The Abbott BinaxNOW rapid antigen test is cheaper and faster than real-time reverse transcription PCR (rRT-PCR) for detecting severe acute respiratory syndrome coronavirus 2. We compared BinaxNOW with rRT-PCR in 769 paired specimens from 342 persons during a coronavirus disease outbreak among horse racetrack workers in California, USA. We found positive percent agreement was 43.3% (95% CI 34.6%-52.4%), negative percent agreement 100% (95% CI 99.4%-100%), positive predictive value 100% (95% CI 93.5%-100%), and negative predictive value 89.9% (95% CI 87.5%-92.0%). Among 127 rRT-PCR-positive specimens, the 55 with paired BinaxNOW-positive results had a lower mean cycle threshold than the 72 with paired BinaxNOW-negative results (17.8 vs. 28.5; p<0.001). Of 100 specimens with cycle threshold <30, a total of 51 resulted in positive virus isolation; 45 (88.2%) of those were BinaxNOW-positive. Our comparison supports immediate isolation for BinaxNOW-positive persons and confirmatory testing for negative persons.

Transmission, infectivity, and neutralization of a spike L452R SARS-CoV-2 variant.

We identified an emerging severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) variant by viral whole-genome sequencing of 2,172 nasal/nasopharyngeal swab samples from 44 counties in California, a state in the western United States. Named B.1.427/B.1.429 to denote its two lineages, the variant emerged in May 2020 and increased from 0% to >50% of sequenced cases from September 2020 to January 2021, showing 18.6%-24% increased transmissibility relative to wild-type circulating strains. The variant carries three mutations in the spike protein, including an L452R substitution. We found 2-fold increased B.1.427/B.1.429 viral shedding in vivo and increased L452R pseudovirus infection of cell cultures and lung organoids, albeit decreased relative to pseudoviruses carrying the N501Y mutation common to variants B.1.1.7, B.1.351, and P.1. Antibody neutralization assays revealed 4.0- to 6.7-fold and 2.0-fold decreases in neutralizing titers from convalescent patients and vaccine recipients, respectively. The increased prevalence of a more transmissible variant in California exhibiting decreased antibody neutralization warrants further investigation.

Introduction, Transmission Dynamics, and Fate of Early Severe Acute Respiratory Syndrome Coronavirus 2 Lineages in Santa Clara County, California.

We combined viral genome sequencing with contact tracing to investigate introduction and evolution of severe acute respiratory syndrome coronavirus 2 lineages in Santa Clara County, California, from 27 January to 21 March 2020. From 558 persons with coronavirus disease 2019, 101 genomes from 143 available clinical samples comprised 17 lineages, including SCC1 (n = 41), WA1 (n = 9; including the first 2 reported deaths in the United States, with postmortem diagnosis), D614G (n = 4), ancestral Wuhan Hu-1 (n = 21), and 13 others (n = 26). Public health intervention may have curtailed the persistence of lineages that appeared transiently during February and March. By August, only D614G lineages introduced after 21 March were circulating in Santa Clara County.

Transmission, infectivity, and antibody neutralization of an emerging SARS-CoV-2 variant in California carrying a L452R spike protein mutation.

We identified a novel SARS-CoV-2 variant by viral whole-genome sequencing of 2,172 nasal/nasopharyngeal swab samples from 44 counties in California. Named B.1.427/B.1.429 to denote its 2 lineages, the variant emerged around May 2020 and increased from 0% to >50% of sequenced cases from September 1, 2020 to January 29, 2021, exhibiting an 18.6-24% increase in transmissibility relative to wild-type circulating strains. The variant carries 3 mutations in the spike protein, including an L452R substitution. Our analyses revealed 2-fold increased B.1.427/B.1.429 viral shedding in vivo and increased L452R pseudovirus infection of cell cultures and lung organoids, albeit decreased relative to pseudoviruses carrying the N501Y mutation found in the B.1.1.7, B.1.351, and P.1 variants. Antibody neutralization assays showed 4.0 to 6.7-fold and 2.0-fold decreases in neutralizing titers from convalescent patients and vaccine recipients, respectively. The increased prevalence of a more transmissible variant in California associated with decreased antibody neutralization warrants further investigation.

A diagnostic host response biosignature for COVID-19 from RNA profiling of nasal swabs and blood.

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), which causes coronavirus disease-19 (COVID-19), has emerged as the cause of a global pandemic. We used RNA sequencing to analyze 286 nasopharyngeal (NP) swab and 53 whole-blood (WB) samples from 333 patients with COVID-19 and controls. Overall, a muted immune response was observed in COVID-19 relative to other infections (influenza, other seasonal coronaviruses, and bacterial sepsis), with paradoxical down-regulation of several key differentially expressed genes. Hospitalized patients and outpatients exhibited up-regulation of interferon-associated pathways, although heightened and more robust inflammatory responses were observed in hospitalized patients with more clinically severe illness. Two-layer machine learning-based host classifiers consisting of complete (>1000 genes), medium (<100), and small (<20) gene biomarker panels identified COVID-19 disease with 85.1-86.5% accuracy when benchmarked using an independent test set. SARS-CoV-2 infection has a distinct biosignature that differs between NP swabs and WB and can be leveraged for COVID-19 diagnosis.