key: cord-0852174-ekj9kar6
authors: Rowe, Lori A.; Beddingfield, Brandon J.; Goff, Kelly; Killeen, Stephanie Z.; Chirichella, Nicole R.; Melton, Alexandra; Roy, Chad J.; Maness, Nicholas J.
title: Intra-Host SARS-CoV-2 Evolution in the Gut of Mucosally-Infected Chlorocebus aethiops (African Green Monkeys)
date: 2022-01-01
journal: Viruses
DOI: 10.3390/v14010077
sha: eec1bfaccf4262bede41d940940f6c4e994e876c
doc_id: 852174
cord_uid: ekj9kar6

In recent months, several SARS-CoV-2 variants have emerged that enhance transmissibility and escape host humoral immunity. Hence, the tracking of viral evolutionary trajectories is clearly of great importance. Little is known about SARS-CoV-2 evolution in nonhuman primate models used to test vaccines and therapies and to model human disease. Viral RNA was sequenced from rectal swabs from Chlorocebus aethiops (African green monkeys) after experimental respiratory SARS-CoV-2 infection. Two distinct patterns of viral evolution were identified that were shared between all collected samples. First, mutations in the furin cleavage site that were initially present in the virus as a consequence of VeroE6 cell culture adaptation were not detected in viral RNA recovered in rectal swabs, confirming the necessity of this motif for viral infection in vivo. Three amino acid changes were also identified; ORF 1a S2103F, and spike D215G and H655Y, which were detected in rectal swabs from all sampled animals. These findings are demonstrative of intra-host SARS-CoV-2 evolution and may identify a host-adapted variant of SARS-CoV-2 that would be useful in future primate models involving SARS-CoV-2 infection.

The COVID-19 pandemic, caused by the coronavirus SARS-CoV-2, has killed more than 4 million people to date. Despite the development and rollout of safe and highly effective vaccines on an unparalleled time scale, nearly unchecked viral spread has continued due to vaccine refusal and hesitancy in nations with adequate vaccine supply coupled with inadequate supply in many regions of the world. As a result, several new variants have emerged with enhanced replicative or infectious capacity and immune escape [1] . Since the discovery of the D614G mutation [2, 3] , noted early in the pandemic and found to enable enhanced infection in cells and is now present in all sequenced isolates, several variants of interest and concern (as defined by the CDC) have arisen [1, [4] [5] [6] [7] [8] [9] [10] [11] . These include B.1.1.7 (alpha), which was originally detected in the United Kingdom and rapidly spread globally, B.1.351 (beta), originally detected in South Africa, and P.1 (gamma), originally detected in Japan in a traveler from Brazil. More recently, the B.1.617.2 (delta) variant rapidly became globally dominant, which is of great concern due to its greatly enhanced transmissibility relative to other variants and its ability to infect vaccinated individuals [12] . At the time of this submission, the B.1.1.529 (omicron) is rapidly increasing in frequency and may overtake delta as the globally dominant variant [13, 14] . Collectively, these data suggest a concerning scenario wherein continued SARS-CoV-2 evolution may facilitate persistence and spread in populations, including those with preexisting immunity to the virus.

Several animal species have been tested in efforts to develop a model that faithfully recapitulates human disease and associated pathological consequences, ultimately to be used in the evaluation of vaccines and therapies for COVID-19. Multiple nonhuman primate (NHP) species have been explored to this end, including Macaca mulatta (rhesus macaque) (RhM), Papio anubis (baboons) [15] , Macaca nemestrina (pigtail macaques) [16] , Macaca fascicularis (cynomolgus macaque), and Chlorocebus aethiops (African green monkeys; AGM) [17] [18] [19] . Although some data suggest enhanced disease in AGM relative to other species [17] , SARS-CoV-2 replicates to high titer in all of these species. However, in most cases, viral replication in upper respiratory sites is restricted to the first several days after infection. Thus, detection of intrahost evolution in these sites may be limited. Importantly, it is now becoming clear that the virus can persist and continue to replicate in gastrointestinal sites [20, 21] , which provides an opportunity to examine viral evolution beyond the first several days of infection and may illuminate tissue-specific compartmentalized evolution.

We recently conducted a study to compare viral and host dynamics in RhM and AGM. Four of each species were infected via the combined intratracheal/intranasal route, and four of each were infected via aerosolized virus. All animals were euthanized four weeks after infection. Viral dynamics were largely the same regardless of route, and AGM tended, on average, to have detectable viral RNA for a longer period, including multiple animals that maintained replicating virus (as evidenced by qRT PCR detection of subgenomic viral RNA) in rectal swab samples for the entirety of the study [22] . Unfortunately, most samples from that study were exhausted and could not be used for sequencing.

In this study, a focused set of gastrointestinal samples from a previous study [22] was used to examine viral dynamics in SARS-CoV-2 infected African green monkeys (AGM). Consistent patterns of evolution were identified, suggestive of adaptive evolution. The available samples were exclusively rectal swabs, which suggest viral replication in the gastrointestinal system may be an important source of viral variants. Notably, these results have important implications for the detection of SARS-CoV-2 variants in wastewater [23, 24] .

Eight RhM and eight AGM were challenged with SARS-CoV-2, WA1/2020 isolate, as described [22] . Briefly, four of each species were challenged via the combined intratracheal/intranasal route (2 × 10 6 TCID50) and four via aerosolized virus from the same stock. Aerosolization inherently leads to a lower effective dose, so animals challenged via this route were estimated to receive a dose between 100 and 1000 times lower than those that received direct mucosal instillation. Viral titers were quantified in several respiratory sites using primers and probes specific to the N gene (for total viral RNA) and specific for both the N and E subgenomic RNAs to approximate replicating virus.

Rectal Swabs were collected and stored in RNA/DNA Shield (Zymo Research, Irvine, CA, USA). RNA was isolated using the Zymo Quick-RNA Viral kit and converted to cDNA using Protoscript II (New England Biolabs, Ipswich, MA, USA) as follows: 10 µL Template RNA, 1 µL 10 µM random hexamers, and 1 µL 10 mM dNTPs were incubated at 65 • C for 5 min and then placed directly on ice for 1 min. The following was then added: 4 µL PSII buffer, 2 µL 100 mM DTT, 1 µL RNase inhibitor, and 1 µL PSII reverse transcriptase and incubated at 42 • C for 50 min, then 70 • C for 10 min, followed by a hold at 4 • C.

DNA libraries were made using the standard SWIFT Normalase Amplicon Panels protocol (SWIFT Biosciences, Ann Arbor, MI, USA) utilizing the SNAP UD indexing primers.

The libraries were normalized to 4 nM and pooled. Paired-end sequencing (2 × 150) was performed on the Illumina (San Diego, CA, USA) MiSeq platform.

Primer sequences were trimmed, and sequence reads were aligned to the SARS-CoV-2 genome (WA1/2020 isolate, accession MN985325) using the built-in mapping function in Geneious Prime software. Variants were called that were present at greater than 10% of reads at that site. Variants detected in the SARS-CoV-2 spike structure were visualized with UCSF Chimera, developed by the Resources for Biocomputing, Visualization, and Informatics at the University of California, San Francisco, CA, USA [25] .

A high-density overlapping amplicon approach was used to amplify and sequence the entire SARS-CoV-2 genome from the viral stock used to experimentally infect the AGM in our previous study [23] . Rectal swabs isolated from three AGM infections at 2 or 4 weeks post-infection were analyzed in the same fashion. The challenge stock harbored several mutations of various frequencies relative to the WA1/2020 patient isolate ( Figure 1A ). These included two mutations in spike, F79L, and a high-frequency mutation that changed the arginine (R) at position 682 to a leucine (R682L) ( Figure 1A ). This residue is the second arginine and one of the key residues in the furin cleavage site, most likely the result of adaptation to VeroE6 cells. This stock was used to infect 16 NHP (8 RhM and 8 AGM). Samples were acquired from nasal, oropharyngeal, and rectal swabs from time points spanning four weeks post-infection. RNA extracted from swabs at early time points was exhausted from extensive qRT PCR characterization, and focus was shifted to rectal swabs collected from late time points, including 14 and 28 days post-infection in three of the AGMs. Sequencing of these samples revealed a consistent pattern of evolution in vivo with three mutations arising in all three animals; ORF1a/b S2103F, spike D215G, and spike H655Y ( Figure 1B ). In contrast, F79L and R682L mutations, present at approximately 60 and 80% in the virus stock, respectively, were completely absent in all rectal samples. Importantly, we note that since available samples used in this report were solely from rectal swabs and solely from AGM, we cannot draw conclusions regarding whether the mutations we identified are specific to this compartment or to this species. Our data simply provide evidence that ongoing viral replication in this compartment can lead to the accumulation of viral mutations, which has not been noted in animal models of COVID-19. More work will be necessary to assess whether these or other mutations occur in respiratory sites, which would be critical for viral transmission.

The ORF1a/b residue 2103 lies in nonstructural protein 3 (nsp3), which is the largest protein in the SARS-CoV-2 proteome and localizes with other viral proteins to the cell membrane [27, 28] . The S2103F mutation in this protein has been detected at low frequency in human samples, including in India [29] , but is not restricted to a particular clade and, through the course of the pandemic, has never been detected in greater than 1% of all sequenced isolates (nextstrain.org accessed on 29 November 2021 [30, 31] ).

The loss of the F79L mutation, present at a relatively low frequency in the inoculum and lost in vivo, may simply reflect the outgrowth of viral variants containing the wild type residue at this position rather than selection. However, the R682L mutation has been noted before in stocks of this virus expanded in Vero E6 cells [32, 33] and likely eliminates the furin cleavage site in spike. Reversion to the wild-type amino acid, therefore, likely confirms the necessity of the furin cleavage site for infection in vivo, which has been previously suggested [34] . This finding, although not unexpected, confirms that the loss of the furin cleavage site in vitro can be ascribed to replication in VeroE6 cells [35] . Since VeroE6 cells are sourced from AGM kidneys, this finding confirms the loss of the furin cleavage site to be an in vitro phenomenon and not specific to replication in AGM cells per se. The D215G mutation rose from approximately 15% prevalence in the stock to near 100% in all three animals. This mutation lies in the N-terminal domain of the spike protein ( Figure 1C) and is detected at low frequency in multiple viral lineages and is a defining variant of the B.1.351 (beta) lineage (nextstrain.org). One report showed this mutation has a weak but detectable impact on augmenting cell to cell fusion [36] . The H655Y mutation lies immediately adjacent to the furin cleavage site ( Figure 1C ) and has been detected in multiple viral lineages, with the P.1 (gamma) and B.1.1.529 (omicron) lineages being most prominent (nextstrain.org). This variant has also been detected in other animal species disease models, including cats and mink [37, 38] . A recent report also showed that virus harboring a 655Y amino acid showed enhanced replicative capacity in vitro and spike protein cleavage [38] , likely explaining the H655Ymutation in the AGM species.

In this study, we define a consistent pattern of viral evolution in SARS-CoV-2-infected AGM. Despite a limited sample size, the ubiquity with which this pattern was detected strongly suggests a selective advantage to the identified mutations. All samples used in this study were derived from rectal rather than respiratory samples, and the data do not clearly suggest a possible mechanism for this selection. It is plausible that replication in the GI tract presents a unique set of selective forces on the virus driving evolution to the mutations identified. Long-term persistence of the virus in the gastrointestinal tract, as has been suggested for SARS-CoV-2 in humans and is a well-described phenomenon in other coronaviruses [39] , may also allow selected mutations to accumulate. Even in the latter case, these data suggest the virus continues to replicate and evolve throughout the gastrointestinal system, which may further reveal novel aspects of the virus-host relationship and may have important implications for the detection of variants in wastewater.

These data and study support two important conclusions. First, previous reports that an intact furin cleavage site is important for infection in vivo were confirmed in this study. Second, the combination of the three mutations identified (ORF1a/b S2103F, spike D215G, and spike H655Y) may represent an adaptation to AGM or to NHP generally. Future experiments may include infection of AGM (or other nonhuman primate species) with this virus to assess the possible enhancement of pathogenesis and to optimize infectious challenges in product evaluation studies. Funding: This work was supported by NIH/NIAID contract HHSN27220170033I (CJR) and grant P51OD01110459, the TNPRC base grant.

The use of nonhuman primates in this project was reviewed and approved by the Tulane Institu-tional Animal Care and Use Committee (IACUC), protocol P0447.

Sequence data obtained in this study are available on the Global Initiative on Sharing All Influenza Data (GISAID) database under accession numbers EPI_ISL_8207137, EPI_ISL_8207138, EPI_ISL_8207139, and EPI_ISL_8207140.

The authors declare no conflict of interest.

The Spike of Concern-The Novel Variants of SARS-CoV-2. Viruses

SARS-CoV-2 spike-protein D614G mutation increases virion spike density and infectivity

Tracking Changes in SARS-CoV-2 Spike: Evidence that D614G Increases Infectivity of the COVID-19 Virus

Resistance of SARS-CoV-2 variants to neutralization by monoclonal and serum-derived polyclonal antibodies

Correlates of Neutralization Against SARS-CoV-2 Variants of Concern by Early Pandemic Sera

The Emerging Concern and Interest SARS-CoV-2 Variants

Implications of the Novel Mutations in the SARS-CoV-2 Genome for Transmission, Disease Severity, and the Vaccine Development. Front. Med. 2021, 8, 636532

Neutralising antibody escape of SARS-CoV-2 spike protein: Risk assessment for antibody-based Covid-19 therapeutics and vaccines

SARS-CoV-2 variants of concern partially escape humoral but not T-cell responses in COVID-19 convalescent donors and vaccinees

N-terminal domain antigenic mapping reveals a site of vulnerability for SARS-CoV-2

The emergence and ongoing convergent evolution of the SARS-CoV-2 N501Y lineages

Comparison of two highly-effective mRNA vaccines for COVID-19 during periods of Alpha and Delta variant prevalence

SARS-CoV-2 B.1.1.529 (Omicron) Variant-United States

SARS-CoV-2 Omicron has extensive but incomplete escape of Pfizer BNT162b2 elicited neutralization and requires ACE2 for infection

Responses to acute infection with SARS-CoV-2 in the lungs of rhesus macaques, baboons and marmosets

The pigtail macaque (Macaca nemestrina) model of COVID-19 reproduces diverse clinical outcomes and reveals new and complex signatures of disease

Acute Respiratory Distress in Aged, SARS-CoV-2-Infected African Green Monkeys but Not Rhesus Macaques

Respiratory disease in rhesus macaques inoculated with SARS-CoV-2

Characteristic and quantifiable COVID-19-like abnormalities in CT-and PET/CT-imaged lungs of SARS-CoV-2-infected crab-eating macaques (Macaca fascicularis)

Evolution of antibody immunity to SARS-CoV-2

SARS-CoV-2 infection of African green monkeys results in mild respiratory disease discernible by PET/CT imaging and shedding of infectious virus from both respiratory and gastrointestinal tracts

Myeloid cell-driven nonregenerative pulmonary scarring is conserved in multiple 1 nonhuman primate species regardless of SARS-CoV-2 infection modality

Tracking COVID-19 with wastewater to understand asymptomatic transmission

Changing composition of SARS-CoV-2 lineages and rise of Delta variant

UCSF Chimera-A visualization system for exploratory research and analysis

SARS-CoV-2 neutralizing antibody structures inform therapeutic strategies

Proteomics Analysis Unravels the Functional Repertoire of Coronavirus Nonstructural Protein 3

Nsp3 of coronaviruses: Structures and functions of a large multi-domain protein

SARS-CoV-2 has observably higher propensity to accept uracil as nucleotide 1 substitution: Prevalence of amino acid substitutions and their predicted functional 2 implications in circulating SARS-CoV-2 in India up to

Maximum-likelihood phylodynamic analysis

Nextstrain: Real-time tracking of pathogen evolution

Human airway cells prevent SARS-CoV-2 multibasic cleavage site cell culture adaptation

High-throughput, single-copy sequencing reveals SARS-CoV-2 spike variants coincident with mounting humoral immunity during acute COVID-19

Loss of furin cleavage site attenuates SARS-CoV-2 pathogenesis

SARS-CoV-2 growth, furin-cleavage-site adaptation and neutralization using serum from acutely infected, hospitalized COVID-19 patients

SARS-CoV-2 Alpha, Beta, and Delta variants display enhanced Spike-mediated syncytia formation

Transmission of SARS-CoV-2 in domestic cats imposes a narrow bottleneck

SARS-CoV-2 variants of concern have acquired mutations associated with an 1 increased spike cleavage

Naturally Occurring Animal Coronaviruses as Models for Studying Highly Pathogenic Human Coronaviral Disease