key: cord-1018453-nv6s49zq
authors: Smith, Darci R.; Singh, Christopher; Green, Jennetta; Lueder, Matthew R.; Arnold, Catherine E.; Voegtly, Logan J.; Long, Kyle A.; Rice, Gregory K.; Luquette, Andrea E.; Miner, Haven L.; Glang, Lindsay; Bennett, Andrew J.; Miller, Robin H.; Malagon, Francisco; Cer, Regina Z.; Bishop-Lilly, Kimberly A.
title: Genomic and Virological Characterization of SARS-CoV-2 Variants in a Subset of Unvaccinated and Vaccinated U.S. Military Personnel
date: 2022-01-27
journal: Front Med (Lausanne)
DOI: 10.3389/fmed.2021.836658
sha: 82930bf60e5a16884ce6803c6c3a2fe51a3bf379
doc_id: 1018453
cord_uid: nv6s49zq

The emergence of SARS-CoV-2 variants complicates efforts to control the COVID-19 pandemic. Increasing genomic surveillance of SARS-CoV-2 is imperative for early detection of emerging variants, to trace the movement of variants, and to monitor effectiveness of countermeasures. Additionally, determining the amount of viable virus present in clinical samples is helpful to better understand the impact these variants have on viral shedding. In this study, we analyzed nasal swab samples collected between March 2020 and early November 2021 from a cohort of United States (U.S.) military personnel and healthcare system beneficiaries stationed worldwide as a part of the Defense Health Agency's (DHA) Global Emerging Infections Surveillance (GEIS) program. SARS-CoV-2 quantitative real time reverse-transcription PCR (qRT-PCR) positive samples were characterized by next-generation sequencing and a subset was analyzed for isolation and quantification of viable virus. Not surprisingly, we found that the Delta variant is the predominant strain circulating among U.S. military personnel beginning in July 2021 and primarily represents cases of vaccine breakthrough infections (VBIs). Among VBIs, we found a 50-fold increase in viable virus in nasal swab samples from Delta variant cases when compared to cases involving other variants. Notably, we found a 40-fold increase in viable virus in nasal swab samples from VBIs involving Delta as compared to unvaccinated personnel infected with other variants prior to the availability of approved vaccines. This study provides important insight about the genomic and virological characterization of SARS-CoV-2 isolates from a unique study population with a global presence.

The emergence of SARS-CoV-2 variants complicates efforts to control the COVID-19 pandemic. Increasing genomic surveillance of SARS-CoV-2 is imperative for early detection of emerging variants, to trace the movement of variants, and to monitor effectiveness of countermeasures. Additionally, determining the amount of viable virus present in clinical samples is helpful to better understand the impact these variants have on viral shedding. In this study, we analyzed nasal swab samples collected between March 2020 and early November 2021 from a cohort of United States (U.S.) military personnel and healthcare system beneficiaries stationed worldwide as a part of the Defense Health Agency's (DHA) Global Emerging Infections Surveillance (GEIS) program. SARS-CoV-2 quantitative real time reverse-transcription PCR (qRT-PCR) positive samples were characterized by next-generation sequencing and a subset was analyzed for isolation and quantification of viable virus. Not surprisingly, we found that the Delta variant is the predominant strain circulating among U.S. military personnel beginning in July 2021 and primarily represents cases of vaccine breakthrough infections (VBIs). Among VBIs, we found a 50-fold increase in viable virus in nasal swab samples from Delta variant cases when compared to cases involving other variants. Notably, we found a 40-fold increase in viable virus in nasal swab samples from VBIs involving Delta as compared to unvaccinated personnel infected with other variants prior to the availability of approved vaccines. This study provides important insight about the genomic and virological characterization of SARS-CoV-2 isolates from a unique study population with a global presence.

The COVID-19 pandemic is currently a leading cause of death globally and new SARS-CoV-2 variants continue to emerge. Genomic surveillance of variants is necessary to characterize mutations that could affect transmissibility and spread, antigenicity or virulence. Furthermore, wet lab studies are necessary to evaluate how genetic differences may affect viral fitness and transmissibility. The Delta variant is currently the predominant strain globally and determining factors that drive its ability to spread rapidly is important. In this study, we characterized SARS-CoV-2 positive samples from U.S. military personnel and their beneficiaries by next-generation sequencing and isolation of viable virus. Consistent with other studies, we found higher levels of infectious virus from Delta samples when compared to non-Delta infections. Strikingly, we found the difference in titer between Delta and other strains to be so profound as to be unaffected by vaccination status, suggesting that increased transmissibility of the Delta variant is in part due to higher amounts of virus shedding. This helps explain the rapid spread of the Delta variant and provides the impetus to increase control measures such as vaccination, boosters, masking and distancing requirements. It will be necessary to continue genomic and virological characterization of new variants, such as Omicron.

The emergence of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), the causative agent of coronavirus disease 2019 , has led to a historic global pandemic. Multiple SARS-CoV-2 genetic variants are circulating globally, some of which have mutations of concern or mutations of interest that could impact transmissibility and spread, antigenicity or virulence (1). The World Health Organization (WHO) classifies many of these as variants of concern (VOC) or variants of interest (VOI) based on mutations that may impact SARS-CoV-2 countermeasures, including vaccines, therapeutics, and diagnostics. Current VOC as defined by the WHO include B. (3) . The Delta variant is estimated to be more transmissible than the Alpha variant, is less sensitive to neutralizing antibodies, and causes increased hospitalization rates mainly in unvaccinated individuals (4) (5) (6) . There are currently multiple vaccines in use either with FDA approval or emergency use authorization (EUA) that have shown efficacy, however vaccine breakthrough infections (VBIs) do occur and may increase with the emergence of novel variants (7) .

Quantitative real time reverse-transcription PCR (qRT-PCR) is the benchmark method for COVID-19 diagnosis, but an important limitation of this approach is the assay only detects RNA and does not detect infectious virus. This can be problematic since several studies have described the persistence of SARS-CoV-2 RNA within different body sites (8) (9) (10) and viral RNA levels as a measurement of viral load does not always correlate to the amount of viable virus in a given sample. Quantitative assays that measure replication competent virus are important for determining infectious viral titers in clinical samples, which is critical to informing infection prevention and control guidelines for current and future SARS-CoV-2 variants.

Detection of infectious virus is determined by in vitro viral culture on susceptible cell lines. However, the ability of viral culture methods to provide insight on infectious shedding is hampered by the fact that this method is labor-intensive and requires biosafety level 3 (BSL-3) facilities for SARS-CoV-2. There are also assay variables such as the cell line and methods used to quantify the amount of virus. Vero cells and their subclones or derivatives are commonly used to culture SARS-CoV-2 and Vero E6 cells expressing human transmembrane serine protease 2 (TMPRSS2) have demonstrated enhanced susceptibility and isolation of SARS-CoV-2 (11) . Plaque or median tissue culture infectious dose (TCID 50 ) assays are commonly used to quantify viable virus, but new methods have been described like the S-Fuse assay that uses a cell reporter system and is based on the detection of GFP positive syncytia formation between infected and neighboring cells (12) .

While some studies have described the amount of viable virus shed by individuals infected with the various SARS-CoV-2 variants (12) (13) (14) , more studies are needed to better characterize infectious virus levels in vaccinated and unvaccinated individuals. This is especially important because each study uses different methods and focuses on different populations of individuals from different geographic locations. In this study, we analyzed nasal swab samples collected between March 2020 and early November 2021 from a cohort of U.S. military personnel and beneficiaries stationed worldwide as a part of the DHA's GEIS program. SARS-CoV-2 qRT-PCR positive samples were characterized by next-generation sequencing and a subset was analyzed for isolation and quantification of viable virus. Following the EUA approval of the COVID-19 vaccines, samples were collected from vaccinated individuals to characterize VBIs. This study provides unique insight about the genomic and virological characteristics of SARS-CoV-2 isolates from a generally healthy population with a global presence that spans from the beginning of the pandemic to early November 2021.

Viral genome sequencing was conducted on the Illumina MiSeq platform using the ARTIC nCoV-2019 Sequencing protocol (15) and the YouSeq SARS-CoV-2 Coronavirus NGS Library prep kit (YouSeq). Approximately 100 ng of RNA was reverse-transcribed as in the protocol, however the YouSeq reverse transcriptase was replaced with SuperScript IV (ThermoFisher Scientific). cDNA was amplified using multiplex PCR and either the associated ARTIC primer pools or YouSeq primer pools. Libraries prepared via the ARTIC protocol were cleaned using 1× AMPure XP beads (Beckman Coulter) and re-suspended in nuclease free molecular grade water. The samples were then processed following the QiaSeq FX protocol (Qiagen) and libraries were completed. Libraries were quality-checked using an Agilent Bioanalyzer High sensitivity kit (Agilent) and quantified using the Qubit DNA High Sensitivity assay (ThermoFisher Scientific) prior to sequencing using Illumina MiSeq v3 2 × 300 chemistry (Illumina).

Viral Amplicon Illumina Workflow (VAIW) was used to collate SARS-CoV-2 consensus genomes (16) . Briefly, Illumina reads were quality trimmed and filtered to >Q20 and minimum length of 50 bp using bbmap (17) . Paired reads were merged using bbmerge with default settings (18) . Trimmed, filtered, and merged reads were aligned to the Wuhan reference genome (NCBI GenBank accession NC_045512.2/MN908947.3) using bbmap with local alignment and maximum insertion/deletion of 500 bp (17) . Primers were trimmed from sequences using align_trim from ARTIC workflow/pipeline (19) . Once a high quality consensus genome was obtained, Single Nucleotide Variants (SNVs) were determined using SAMtools mpileup (20) and iVar (Intrahost variant analysis of replicates) (21) 

Patient samples were cultured for SARS-CoV-2 using a standard plaque assay on Vero E6 cells expressing TMPRSS2 [4] in six-well plates and a cytopathic effect (CPE) assay on Vero E6/TMPRSS2 cells in T25 cm 2 flasks. For the plaque assay, duplicate wells were infected with 0.2 mL aliquots from serial 10fold dilutions in Minimum Essential Medium (MEM) followed by a 1 h incubation at 37 • C, 5% CO 2 to allow virus adsorption to occur. After incubation, cells were overlaid with MEM containing 0.5% agar supplemented with 5% heat-inactivated FBS, 1% Penicillin/Streptomycin and incubated for 48 h at 37 • C, 5% CO 2 . Cells were fixed in 10% formalin prior to staining with crystal violet. For the CPE assay, cells were seeded in T25 cm 2 flasks and each flask was infected with 0.5 mL aliquots from a 1:10 dilution in Dulbecco's Modified Eagle Medium (DMEM) followed by a 1 h incubation at 37 • C, 5% CO 2 . After incubation, 5 mL of DMEM supplemented with 5% heatinactivated FBS, 1% Penicillin/Streptomycin was incubated for 4 days at 37 • C, 5% CO 2 . CPE was monitored daily. After 4 days, supernatant was passed onto fresh cells to allow additional time to amplify. The presence of SARS-CoV-2 was confirmed by lateral flow immunoassay using the Quidel Sofia 2 analyzer and/or qRT-PCR.

Statistical analysis was performed using GraphPad Prism. Tukey's multiple comparisons test was used to compare the amount of viable virus between the unvaccinated and vaccinated groups.

We have sequenced over 2,300 nasal swab samples collected from March 2020 until early November 2021 and of those, we produced 1,304 coding complete genomes for which we then determined the Pango lineages (Figure 1 and Supplementary Table 1 ). Smaller numbers of samples were received early in the pandemic compared to after the variants began to emerge when there was more impetus to sequence a greater number of samples. The first samples sequenced in March and April 2020 were similar to strains at the root of the pandemic (lineage A and sublineages) and are classified as A. (Figures 1B,D) . The "other" category for Figures 1B,D was used for any sublineage with 5 or less representatives. The Omicron variant was not detected during the timeframe of our current study.

A subset of samples were selected to determine if viable virus was present where we compared viral titers in samples collected from unvaccinated individuals and cases of VBI (Figure 2) . A summary of the results broken down by Pango lineage is shown in Tables 1, 2 for unvaccinated and vaccinated individuals, respectively. We detected viable virus in 25 samples received from unvaccinated personnel collected between March 2020 and January 2021, which is prior to the emergence of VOCs. The mean amount of virus detected in unvaccinated personnel was 3.2 log 10 PFU/mL, which is similar to the amount of virus detected from VBIs associated with all variants except for Delta (3.1 log 10 PFU/mL). Significantly more virus was detected in VBIs associated with the Delta variant when compared to unvaccinated personnel. The mean amount of virus detected in VBIs associated with the Delta variant was 4.6 log 10 PFU/mL, which was statistically significant by Tukey's multiple comparisons test (p = 0.0001) when compared to the amount of virus detected in unvaccinated personnel prior to the emergence of VOCs. In contrast to the unvaccinated samples, almost all the VBI samples associated with the Delta variant had detectable virus that was above the level of detection for the plaque assay and one sample had a titer as high as 7 log 10 PFU/mL. Samples from unvaccinated personnel infected with the Delta variant were not available for this study since the majority of military personnel were fully vaccinated when the Delta variant became the predominant strain. The amount of virus detected in VBIs associated with the Delta variant compared to VBIs with other variants was also statistically significant by Tukey's multiple comparisons test (p = 0.0132) although with less power because of the low number of VBI samples associated with other variants. The results are further broken down by vaccine manufacturer (Figure 2B) where we continued to observe a statistically significant difference in viral titers for unvaccinated individuals with non-Delta variants compared to samples from VBI associated with the Delta variant from individuals that received the Pfizer (p = 0.0106) and Johnson & Johnson (p = 0.0427) vaccines. Interestingly, VBI associated with the Delta variant that received the Moderna vaccine was not statistically significant when compared to the unvaccinated group. This may be a reflection of the limited number of samples in our analyses. 

The Delta variant displaced the Alpha variant and was identified as the predominant strain in military personnel starting in July 2021. Currently, multiple sublineages of B.1.617.2 continue to be identified. A greater proportion of the VBI cases in our dataset have been associated with the Delta variant, commensurate with its current preponderance in the population at this point in time.

We quantified the amount of viable virus in a subset of samples, and our results suggest that the increased transmissibility of the Delta variant is in part due to higher amounts of virus shedding. We detected 50-fold more infectious virus in VBIs involving Delta compared to VBIs with other variants. We also detected 40-fold more infectious virus in VBIs involving Delta compared to unvaccinated personnel infected with non-Delta variants.

Our results are similar to what has been reported for the rest of the population in the U.S. where the Delta variant is currently the predominant strain and is more than twice as infectious as previous strains. The improved viral fitness of the Delta variant has been associated with mutations in the furin cleavage site that allow the virus to enter cells more efficiently (27, 28) . Most research has focused on analyzing mutations in the spike (S) protein to better understand viral fitness of emerging SARS-CoV-2 variants. However, a recent study found that four mutations in the nucleocapsid (N) protein, which is an important protein for viral replication, improve the viruses' ability to form infectious particles (29) . These four mutations in N, which are universally found in more transmissible SARS-CoV-2 variants, increased mRNA delivery and expression and two of the mutations (S202R and R203M) produced more infectious virus (29) . Interestingly, the Delta variant contains the R203M mutation and the Omicron variant contains the R203K mutation, which provide a possible explanation for the increased fitness of these variants. This is consistent with other studies that found that samples from people infected with the Delta variant have higher viral RNA levels and an increased duration of viral shedding compared to samples collected from people infected with other variants (14, (30) (31) (32) . Furthermore, the levels of viral RNA from infected individuals have been shown to be similar for vaccinated and unvaccinated individuals (14, 33) . More studies are needed to determine the level and duration of viral shedding with the newly emerging Omicron variant.

Numerous studies measure viral load by the amount of viral RNA in a clinical sample, but do not determine the amount of viable virus. This is likely due to the cumbersome nature of the assays used to quantify levels of infectious virus coupled with the requirement to perform these assays in a BSL-3 laboratory. However, a recent study found a significant amount of variability in the relationship between viral titers for infectious virus and viral RNA levels (13) . Their results indicate that determining viral load by RNA levels in clinical samples may have limitations and assays for viable virus should be included for emerging variants. Additionally, they detected six times as much infectious virus for the same amount of RNA for Delta variant samples compared to Alpha variant samples. Our results are similar where we found that the Delta variant cases shed viable virus at higher levels. Similar levels of virus were detected from unvaccinated personnel with non-Delta infections prior to the availability of vaccines granted EUA and FDA approval compared to VBI due to other variants besides Delta. In contrast, 50-fold more infectious virus was detected in samples from VBI due to Delta compared to VBIs with other variants. When we analyzed our results by vaccine manufacturer, we continued to detect significantly higher levels of infectious virus in samples from individuals who received the Pfizer and Johnson & Johnson vaccines and experienced a VBI due to Delta when compared to the unvaccinated personnel with non-Delta infections. Interestingly, VBI associated with the Delta variant that received the Moderna vaccine was not statistically significant when compared to the unvaccinated group. This may be a reflection of the limited number of samples in our analyses. Another limitation of our study is the lack of access to metadata for the clinical specimens, which affects our interpretation of the results by differences in age, pre-existing conditions, and days from exposure or symptom onset. However, most of our samples did come from active duty military personnel stationed worldwide who are considered young (age , are required to maintain a physical fitness regimen and pass routine physicals, and are therefore considered relatively healthy compared to the general population.

Other studies have also found that viable virus is more likely to be isolated from clinical samples from Delta compared to Alpha (14) . The authors do point out that there was a significantly longer period between being considered fully vaccinated and when VBI with the Delta variant occurred, which is similar to our results. However, our data demonstrating that we detected 40-fold more infectious virus in VBIs involving Delta compared to unvaccinated personnel with non-Delta infections prior to the availability of vaccines granted EUA and FDA approval suggests that the increase in the amount of viable virus is due to the variant and not the vaccination status. This is consistent with the study by Luo et al. (14) that found the recovery of infectious virus was significantly higher with the Delta variant compared to Alpha in both vaccinated and unvaccinated groups.

Our results highlight the importance of characterizing the SARS-CoV-2 genome and infectious viral titers of emerging variants isolated from unvaccinated and vaccinated personnel. This work contributes to the growing body of evidence that the increased fitness of Delta is primarily due to viral genetic determinants that increase infectivity or replication. This helps to explain the rapid spread of the Delta variant and provides the impetus to increase control measures such as vaccination, boosters, masking and distancing requirements to prevent ongoing transmission. Furthermore, this study demonstrates the value of pairing genomic and virological characterization when assessing novel SARS-CoV-2 variants. These parallel approaches should be applied to newly emerging variants, such as Omicron, and other variants that may arise in the future.

The datasets presented in this study can be found in online repositories: GenBank (https://www.ncbi.nlm.nih.gov/genbank/) and GISAID (https://www.gisaid.org/). The names of the repository/repositories and accession number(s) can be found in the Supplementary Material. 

We thank Justice Akuoku Frimpong, Kevin Schully, and Christina Farris for technical assistance with the viral cultures.

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmed. 2021.836658/full#supplementary-material

Covid-19 vaccines and variants of concern: a review

Estimated transmissibility and impact of SARS-CoV-2 lineage B.1.1.7 in England

GISAID. hCov19 Variants. Available online at

Reduced sensitivity of SARS-CoV-2 variant Delta to antibody neutralization

Effect of Delta variant on viral burden and vaccine effectiveness against new SARS-CoV-2 infections in the UK

Hospital admission and emergency care attendance risk for SARS-CoV-2 delta (B.1.617.2) compared with alpha (B.1.1.7) variants of concern: a cohort study

SARS-CoV-2 breakthrough infections in vaccinated individuals: measurement, causes and impact

Temporal dynamics in viral shedding and transmissibility of COVID-19

Virological assessment of hospitalized patients with COVID-2019

Persistent detection and infectious potential of sars-cov-2 virus in clinical specimens from COVID-19 patients

Enhanced isolation of SARS-CoV-2 by TMPRSS2-expressing cells

Release of infectious virus and cytokines in nasopharyngeal swabs from individuals infected with non-alpha or alpha SARS-CoV-2 variants: an observational retrospective study

Quantitative measurement of infectious virus in SARS-CoV-2 Alpha, Delta and Epsilon variants reveals higher infectivity (viral titer:RNA ratio) in clinical samples containing the Delta and Epsilon variants. medRxiv

Infection with the SARS-CoV-2 Delta variant is associated with higher recovery of infectious virus compared to the Alpha variant in both unvaccinated and vaccinated individuals

nCoV-2019 sequencing protocol

A Custom Pipeline to Analyze the SARS-CoV-2 Genomes Prepared With an Amplicon (ARTIC (v3) and YouSeq (v2)) Based Library Protocols

BBMAP: A Fast, Accurate, Splice-Aware Aligner

BBMerge -accurate paired shotgun read merging via overlap

Available online at

The sequence alignment/map format and SAMtools

An amplicon-based sequencing framework for accurately measuring intrahost virus diversity using PrimalSeq and iVar

Available online at

MAFFT: a novel method for rapid multiple sequence alignment based on fast Fourier transform

IQ-TREE 2: new models and efficient methods for phylogenetic inference in the genomic era

Available online at

A dynamic nomenclature proposal for SARS-CoV-2 lineages to assist genomic epidemiology

Spike protein cleavage-activation mediated by the SARS-CoV-2 P681R mutation: a case-study from its first appearance in variant of interest (VOI) A.23.1 identified in Uganda. bioRxiv

The SARS-CoV-2 variants associated with infections in INdia, B.1.617, show enhanced spike cleavage by furin

Rapid assessment of SARS-CoV-2 evolved variants using virus-like particles

Viral Infection and Transmission in a large well-traced outbreak caused by the Delta SARS-CoV-2 variant. medRxiv

Clinical and virological features of SARS-CoV-2 variants of concern: a retrospective cohort study comparing B.1.1.7 (Alpha), B.1.315 (Beta), and B.1.617.2 (Delta)

The Delta SARS-CoV-2 variant has a higher viral load than the Beta and the historical variants in nasopharyngeal samples from newly diagnosed COVID-19 patients

Viral dynamics of acute SARS-CoV-2 infection and applications to diagnostic and public health strategies

The remaining authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.