Next Article in Journal
Alterations of Plasma Biochemical and Immunological Parameters and Spatiotemporal Expression of TLR2 and TLR9 in Gibel Carp (Carassius auratus gibelio) after CyHV-2 Infection
Previous Article in Journal
Correction: Chen et al. Identification of the Causal Agent of Brown Leaf Spot on Kiwifruit and Its Sensitivity to Different Active Ingredients of Biological Fungicides. Pathogens 2022, 11, 673
Previous Article in Special Issue
Novel Molecular Consortia of Cannabidiol with Nonsteroidal Anti-Inflammatory Drugs Inhibit Emerging Coronaviruses’ Entry
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Genomic Mutations in SARS-CoV-2 Genome following Infection in Syrian Golden Hamster and Associated Lung Pathologies

by
Gudepalya Renukaiah Rudramurthy
*,
Chakenahalli N. Naveenkumar
,
Kumaraswamy Bharathkumar
,
Radha K. Shandil
and
Shridhar Narayanan
Foundation for Neglected Disease Research (FNDR), Plot No. 20A, KIADB Industrial Area, Bengaluru 561203, Karnataka, India
*
Author to whom correspondence should be addressed.
Pathogens 2023, 12(11), 1328; https://doi.org/10.3390/pathogens12111328
Submission received: 22 October 2023 / Revised: 6 November 2023 / Accepted: 6 November 2023 / Published: 8 November 2023
(This article belongs to the Special Issue SARS-CoV-2: From Virus Replication Cycle to Antiviral Strategies)

Abstract

:
The continuous evolution of the SARS-CoV-2 virus led to constant developments and efforts in understanding the significance and impacts of SARS-CoV-2 variants on human health. Our study aimed to determine the accumulation of genetic mutations and associated lung pathologies in male and female hamsters infected with the ancestral Wuhan strain of SARS-CoV-2. The present study showed no significant difference in the viral load between male and female hamsters and peak infection was found to be on day four post infection in both sexes of the animals. Live virus particles were detected up to 5 days post infection (dpi) through the TCID-50 assay, while qRT-PCR could detect viral RNA up to 14 dpi from all the infected animals. Further, the determination of the neutralizing antibody titer showed the onset of the humoral immune response as early as 4 dpi in both sexes against SARS-CoV-2, and a significant cross-protection against the delta variant of SARS-CoV-2 was observed. Histopathology showed edema, inflammation, inflammatory cell infiltration, necrosis, and degeneration of alveolar and bronchial epithelium cells from 3 dpi to 14 dpi in both sexes. Furthermore, next-generation sequencing (NGS) showed up to 10 single-nucleotide polymorphisms (SNPs) in the SARS-CoV-2 (ancestral Wuhan strain) genome isolated from both male and female hamsters. The mutation observed at the 23014 position (Glu484Asp) in the SARS-CoV-2 genome isolated from both sexes of the hamsters plays a significant role in the antiviral efficacy of small molecules, vaccines, and the Mabs-targeting S protein. The present study shows that either of the genders can be used in the pre-clinical efficacy of antiviral agents against SARS-CoV-2 in hamsters. However, considering the major mutation in the S protein, the understanding of the genetic mutation in SARS-CoV-2 after passing through hamsters is crucial in deciding the efficacy of the antiviral agents targeting the S protein. Importance: Our study findings indicate the accumulation of genomic mutations in SARS-CoV-2 after passing through the Syrian golden hamsters. Understanding the genomic mutations showed that either of the hamster genders can be used in the pre-clinical efficacy of antiviral agents and vaccines.

1. Introduction

Coronaviruses (CoVs) belonging to the family Coronaviridae and the order Nidovirales are genetically diverse and enveloped viruses, consisting of four genera, alpha, beta, gamma, and delta. Coronaviruses are capable of infecting multiple animal species and several cross-species transmission incidences have been identified [1]. Humans are known to be infected by only alpha and beta CoVs, and the majority of these infections are thought to be due to zoonotic transmission [2,3,4,5]. The ongoing pandemic caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), known as coronavirus disease 2019 (COVID-19), was reported in early December 2019, in Wuhan, China [6,7]. The SARS-CoV-2 belonging to beta coronaviruses is more contagious [8] and continuous mutations in the spike protein led to the emergence of multiple variants such as B.1.1.7 (alpha), B.1.351 (beta), B.1.1.28.1/P.1 (gamma), B.1.617.2 (delta), B.1.617.1 (or kappa), and B.1.1.529 (omicron), which are circulating globally in the ongoing COVID-19 pandemic [9,10], and many different lethal variants will keep emerging in the future. The World Health Organization designated the B.1.617.2 (delta) variant as a variant of concern (VOC) and mutations in the spike protein have increased the binding efficiency to the ACE2 receptor, the transmissibility, secondary attack rates, and the escape from immunity and neutralization [11].
The clinical symptoms and severity of SARS-CoV-2 infection are highly variable with a significant mortality, and populations with underlying medical conditions are highly vulnerable [6]. Millions of human lives have been claimed since the emergence of COVID-19 with a catastrophic impact, and this has pushed the economy of the world into a significantly bad shape. Several different vaccines have been developed at a rapid pace and have been given approval for emergency use by the respective government agencies. The swift vaccination programs and management strategies across the globe helped in managing the severity of SARS-CoV-2 in subsequent infection waves. However, the development of specific antivirals is another area of focus and urgent need, and further, the specific antivirals or therapies against any human CoVs are not yet available [12,13].
A key strategy to minimize the effect of the COVID-19 pandemic on human health and to protect from the ongoing pandemic is the development of effective antivirals and vaccines. Preclinical studies involving animals play a critical role in the determination of efficacy of the biologics and antivirals and help further take up the clinical studies based on the preclinical data. Hence, developing and validating suitable animal infection models play a crucial role in the preclinical study of antiviral inhibitors and vaccine efficacy. However, the absence of appropriate receptors for many human viruses in animals impedes the development of suitable animal infection models. Small mammals such as Syrian hamsters (Mesocricetus auratus) are known to express ACE-2 receptors and interact with the receptor binding domain (RBD) of SARS-CoV-2 [14] and mediate entry [15]. Further, hamsters are susceptible to several respiratory virus infections such as SARS-CoV, influenza, and SARS-CoV-2, and have been suitably used as infection models for antiviral and vaccine efficacy testing [16,17,18]. The pathogenicity, severity of lung disease, and tissue tropism during SARS-CoV-2 infection in hamsters have been evaluated and hamsters serve as suitable mammalian infection models for COVID-19 [19].
In the present study, we evaluated the SARS-CoV-2 replicative ability, course of infection, and lung pathology in male and female Syrian golden hamsters with an ancestral Wuhan Strain of SARS-CoV-2. Further, the mutation capability of SARS-CoV-2 after passing through male and female hamsters was evaluated to determine the accumulation of genomic mutations in SARS-CoV-2 post-animal infection. The viral load in the lungs was determined through a standard plaque assay and quantitative real-time PCR. This study included the identification of differences in histopathological changes in the lungs of both males and females infected with the Wuhan strain of SARS-CoV-2. Furthermore, the RNA isolated from the virus stock (SARS-CoV-2) used in animal infection and the viral RNAs extracted from the infected animal tissues were sequenced to determine any mutations in the SARS-CoV-2 genome. The SARS-CoV-2 infection study was carried out in a high-containment (BSL-3) facility.

2. Materials and Methods

2.1. Cell Lines and Virus Propagation

The Vero-E6 cell (C1008) (African green monkey kidney epithelial cell) was obtained from Elabscience Biotechnology Inc. (Cat no. EP-CL-0491) and maintained in Dulbecco’s modified Eagle’s medium (DMEM) supplemented with 10% fetal bovine serum (FBS), an antibiotic antimycotic solution containing penicillin (100 units/ml), streptomycin (100 µg/ml), and amphotericin B (0.25 µg/ml).
SARS-CoV-2 isolates, USA-WA1/2020 (ancestral Wuhan) (Cat No; NR-52281) and B.1.617.2 (delta) (Cat No; NR- NR-55611), used in the present study were obtained from BEI Resources. The virus stocks were prepared by propagating the virus strains in Vero E6 cells by following the standard protocol [20,21]. The virus stocks were titrated using a gold-standard plaque assay [20] to determine the plaque forming units (PFU/ml) and 50% tissue culture infective dose (TCID50/ml) and stored at −80 °C in small aliquots [22]. The TCID50/ml of the stocks was determined using the Reed and Muench method [23]. Furthermore, the viral RNA was extracted from both stocks and subjected to next-generation sequencing (NGS) to confirm the genome sequence of SARS-CoV-2.

2.2. SARS-CoV-2 Infection in Syrian Golden Hamsters

The male and female hamsters were used in the present study to understand the replicative capability, pathology, and mutation rates of SARS-CoV-2 (ancestral Wuhan strain). The Syrian golden hamsters (males and females) used in the present study were obtained from ACTREC (Advanced Centre for Treatment, Research and Education in Cancer, Mumbai). The health status of the animals and body weights were regularly monitored during the quarantine period for 7 days before being used in this study. The studies were conducted in our BSL-3 containment facility. The males and females were separately grouped (n = 5) and acclimatized at the BSL-3 facility for 7 days before the study initiation. The animals were anesthetized with ketamine at 150 mg/kg body weight and Xylazine at 10 mg/Kg body weight through the intraperitoneal route and intranasally infected (100 µL) with SARS-CoV-2 virus, USA-WA1/2020 (ancestral Wuhan) at 105 TCID50 per animal. The uninfected group (Mock control) in each study was intranasally given 100 µL of DMEM to mimic the infection group.
The animals were sacrificed on different days post-infection (dpi) of 2, 3, 4, 5, 7, and 14 dpi to understand the pathogenesis and immune response during SARS-CoV-2 infection. The body weight of the animals was recorded before the sacrifice at each time point. The whole lungs were collected from the animals and the right lobe of the lungs was used to determine the viral load, while the left lobes were fixed in 10% buffered formalin for histopathology study. The mock control groups were terminated at the last time point (14 dpi). The animals were euthanized and sacrificed by following the standard approved protocol and the institutional CPCSEA number is FNDR-FB-83-PD-2021.

2.3. Determination of Viral Load from the Lung Samples

The infectious virus load in the lung tissues of the control and infected groups was determined using standard TCID-50 assay and real-time quantitative PCR. The tissue samples were added with sterile PBS at a 1:2 ratio (w/v) and homogenized using a Tissue Homogenizer (BIOSPEC PRODUCTS INC., USA Cat No; 985370EUR-07) for 10 seconds and then the homogenate was clarified through centrifugation. The viral load from the supernatant was determined through standard TCID-50 assay and is expressed as TCID-50/lung tissue.
Further, the viral RNA was extracted from all the tissue homogenates using QIAamp Viral RNA Mini Kit (Qiagen) by following the manufacturer’s instructions. Quantitative real-time PCR (qRT-PCR) was performed to determine the genomic viral RNA using US-2019-nCoV Assays for the detection of SARS-CoV-2 RNA (Sigma Darmstadt, Germany, Cat no.: CDA00011). The viral RNA per lung tissue was calculated and compared to determine the course of infection in male and female hamsters with two different variants of SARS-CoV-2.

2.4. Next-Generation Sequencing (NGS) to Study Mutation in the Viral RNA

The NGS study was conducted to determine any mutations in the viral RNA after passing through hamsters in comparison to virus stocks used in animal inoculation. For this study, the viral RNA was extracted from the virus stock used in the infection of hamsters and from the lung homogenates of infected male (n = 1) and female (n = 1) groups at 4 dpi (Table 1). The extracted samples were confirmed for the presence of SARS-CoV-2 viral RNA using qRT-PCR. The whole-genome sequencing of the isolated RNA samples was conducted using the NextSeq500 platform (https://www.illumina.com/systems/sequencing-platforms/nextseq.html, accessed on 5 November 2023). The consensus sequence was extracted from each sample, and a single-nucleotide polymorphism (SNP) was identified using the mpileup utility of Samtools (version; v 0.1.18). The most likely lineages to the sample consensus sequences were assigned using the Pangolin tool.

2.5. Immune Response against SARS-CoV-2 Infection in Hamsters

The neutralizing antibodies in the serum samples of infected hamsters were determined through the Microneutralization Test (MNT) to understand the immune response elicited by hamsters against SARS-CoV-2 infection. The blood samples were collected at time points of 2, 3, 4, 5, 7, and 14 dpi from all the animals from different groups, and serum was separated, heat-inactivated at 56 °C for 1 h, and preserved at −80 °C until further use. The neutralizing antibody titer in the sera sample was determined using MNT against SARS-CoV-2, USA-WA1/2020 (ancestral Wuhan) to determine the immunological protection. Further, the neutralizing antibodies against SARS-CoV-2, delta variant (B.1.617.2), were estimated to determine the cross-protection elicited by SARS-CoV-2 (ancestral Wuhan) infection in hamsters. In brief, the serum samples were serially diluted and incubated with the respective SARS-CoV-2 variants at 100 TCID-50 for 1 hour at 37 °C. Upon completion of the incubation, the virus load was determined to calculate the 50% virus neutralization. To summarize, the ancestral strain (Wuhan strain)-infected animal sera samples were tested against the SARS-CoV-2 (ancestral Wuhan and delta variant) to understand the earliest protection and cross-protection elicited during SARS-CoV-2 infection in animals.

2.6. Lung Histopathology

The lungs were collected from infected groups and mock control group animals on the day of sacrifice at different time points. The images of whole lungs were captured to observe any gross pathological changes in comparison to the control group. The left lobes of the lungs were perfused with 10% buffered formalin and subjected to histopathological analysis. In brief, the samples were processed for paraffin embedding after fixing in 10% buffered formalin. Subsequently, thin sections (3 µm) of paraffin blocks were prepared and mounted on a glass slide. The sections were then stained with hematoxylin and eosin and the images were captured at 10× and 60× magnification.

2.7. Statistical Analysis

The data collected from this study such as body weight and viral load were analyzed using Graph pad Prism software (Version; 9.3.1, Boston, MA, USA). The data were statistically analyzed using a two-tailed unpaired t-test to determine the significance between the genders and variants of SARS-CoV-2 infection.

3. Results

3.1. Change in Body Weight Remained Similar in Both Sexes

The body weight of the study animals was measured before the sacrifice at every time point to determine the change in body weight. The schematic diagram depicting the study design and time points is shown in Figure 1. The mock control group animals showed a decrease in body weight on day one post-infection followed by recovery in the body weight on subsequent days. The infected group showed a reduced body weight (up to 5%) until 5 dpi and subsequently recovered the body weight. The change in body weight showed a significant difference between the mock and infected groups and a similar trend was observed in both sexes of the animals (Figure 2A,B).

3.2. Whole-Genome Sequencing of Viral RNA Showed Mutations in the Spike Protein

The viral RNA isolated from the virus stock (used in the infection of hamsters) and from the infected lung samples of male (n = 1) and female (n = 1) hamsters were subjected to whole-genome sequencing using next-generation sequencing (NGS), to identify mutations in the SARS-CoV-2 genome. The whole genome of the stock virus (ancestral Wuhan) used in the infection of hamsters was sequenced and compared with the reference sequence SARS-CoV- 2, isolate Wuhan-Hu-1, complete genome (NC_045512.2). The comparison showed 8 SNPs, 5 of which were in the spike (S) protein gene, 2 of which were in the ORF1ab, and 1 of which was in the ORF8 (Table 1). The comparison details of SNPs and respective changes in the amino acid sequences are given in Table 2. The viral RNA genomes isolated from the lungs of hamsters (male and female) were compared with the stock virus (used in infection) to identify the SNPs. The comparison showed 10 SNPs in both male and female hamsters (Table 1). The mutation sites remain similar between the viral genomes isolated from male and female hamsters; however, at position 22296, viral RNA from the male hamster showed a mutation from G to A with a protein change from Arg245His, while the viral RNA from the female hamster did not show any mutation (Table 3).

3.3. Lung Viral Load Revealed No Significant Difference between Male and Female Hamsters

The infectious SARS-CoV-2 virus particles and SARS-CoV-2 RNA copies in the lung samples of infected animals were determined using the TCID-50 assay and qRT-PCR, respectively. The TCID-50 assay showed a peak average viral load of log10 6.0 and qRT-PCR showed average viral RNA copies of log10 8.0 at 4 dpi in both male and female hamsters. The viral load/RNA copies’ data showed no significant difference among male and female hamsters (Figure 3A,B). The TCID-50 assay could not detect SRAS-CoV-2 virus particles at 7 dpi and 14 dpi from both males and females. However, the qRT-PCR assay could detect SARS-CoV-2 RNA both at 7 dpi and 14 dpi with an approximate viral RNA copy of log10 5.0 (Figure 3C,D).

3.4. Immune Response against SARS-CoV-2 Has Been Detected in Infected Hamsters as Early as 4 dpi

The antibody titers in the serum samples of infected animals at different time points were estimated against SARS-CoV-2 variants to determine the early time point at which the immune response was elicited in the hamsters. The neutralizing antibodies in serum samples were determined using the microneutralization assay. The neutralizing antibodies against SARS-CoV-2 were observed at 4 dpi against both variants in both males and females. The hamsters (male and female) infected with the SARS-CoV-2 (ancestral Wuhan strain) showed a similar neutralizing antibody titer against the delta variant from 4 dpi and onward. Further, no significant difference was observed between male and female hamsters in the antibody neutralization capability against both the ancestral Wuhan strain and the delta variant (Figure 4A,B).

3.5. Significant Lung Pathology Has Been Observed as Early as 2 dpi in SARS-CoV-2-Infected Hamsters

The SARS-CoV-2-induced progression of respiratory disease was detected through the evaluation of lung histopathological changes. The gross pathological examination of the lungs showed no abnormal gross lesions until 2 dpi in both male and female hamsters infected with the SARS-CoV-2 (ancestral Wuhan strain). At 3 dpi, slight generalized edema with areas of focal congestion was observed in all the groups. The lung pathology was found to be severe including marked inflammation, multifocal areas of hemorrhage, congestion, and edema at 4 dpi and onward to 7 dpi in all the groups. However, a reduced lung pathology was observed at 14 dpi in all the groups (Figure 5).
The histopathological examination revealed minimal to mild inflammation with focal inflammatory cell infiltration, mild necrosis, degeneration of alveolar and bronchial epithelium cells, and peribronchial and perivascular inflammatory cell infiltration with vascular congestion at 2 dpi and 3 dpi in both male and female hamsters (Figure 6).
At 4 dpi and 5 dpi, the lung samples showed multifocal areas of inflammation, moderate to marked broncho-interstitial pneumonia, and an accumulation of mononuclear inflammatory cells like neutrophils, lymphocytes, and macrophages. Furthermore, moderate alveolar and bronchial epithelium degeneration, necrosis with cell debris, and inflammatory cells in the alveolar and bronchial spaces were prominent. Vascular changes like edema, mild endotheliosis, congestion, and areas of hemorrhage were observed in alveoli, and in peribronchial and perivascular areas, mild to moderate bronchial and alveolar epithelial hyperplasia with hyperplasia of type II alveolar epithelia were observed (Figure 6).
The lung pathology at 7 dpi in both males and females showed a similar lung pathology such as marked to severe multifocal broncho-interstitial pneumonia with a marked increase in lung cellularity. Further, the inflammatory cell infiltrates and exudate identified at 3 to 5 dpi were largely replaced by proliferative bronchiolar and alveolar epithelia. Moreover, marked alveolar and bronchial epithelium degeneration and marked hyperplasia of Type II alveolar cells were observed (Figure 6). At 14 dpi, the lungs showed a minimal focal infiltration of inflammatory cells and congestion. However, hyperplasia of the bronchial epithelium cells was observed in all the animals at 14 dpi (Figure 6).

4. Discussion

Syrian hamsters (Mesocricetus auratus) have been used as animal infection models for many respiratory and non-respiratory viruses such as influenza virus, SARS-CoV, and Japanese encephalitis virus [24]. Golden Syrian hamsters that are susceptible to almost every strain of SARS-CoV-2 are being explored in the study of SARS-CoV-2-related pathogenesis, vaccine, and antiviral efficacy testing [11,24,25]. The present study is designed to determine the viral replication and pathogenesis associated with SARS-CoV-2 infection in male and female hamsters. Further, this study was aimed at the determination of the accumulation of genetic mutations in the ancestral Wuhan strain of SARS-CoV-2 in male and female hamsters.
The body weight measurement at different time points showed a body weight loss of up to 5% in both sexes of the animals up to 4 dpi followed by a gain in body weight. Previous studies have shown that the SARS-CoV-2 infection in hamsters is associated with clinical manifestations such as reduced body weight up to 7 dpi followed by recovery [19,26,27]. Further, it has been observed that the change in body weight in hamsters is associated with virus inoculum dose rather than the age of the animal [28,29]. The estimation of SARS-CoV-2 viral RNA through qRT-PCR and viral particles through the TCID-50 assay showed peak infection at 4 dpi in both male and female hamsters, and the present results agree with the previous studies [19,26]. None of the groups showed live virus particles at 7 dpi and 14 dpi through the TCID-50 assay; however, qRT-PCR could detect the viral RNA at 7 dpi and 14 dpi, and this is due to the presence of viral RNA in the lungs even after the virus clearance [24]. The TCID-50 assay detects live virus particles, whereas qRT-PCR detects viral RNA from both live and dead virus particles.
The immune response against SARS-CoV-2 infection has been evaluated through the determination of the neutralizing antibody titer from the serum samples collected at different time points. The neutralizing antibody titers were detected as early as 4 dpi in both males and females from infected hamsters’ sera (Figure 4). The neutralizing antibody titer remains slightly higher against homologous infected sera compared to cross-reactive neutralizing antibodies. Further, the difference in antibody titer was not observed between the sexes of the animals. In the earlier study, no significant difference in antibody titers was observed amongst variants from B.1.1.7-infected serum samples at 14 and 28 dpi [30]. The present study showed the early onset of humoral immune response and the cross-protection elicited against different variants of SARS-CoV-2.
The gross lung pathology showed inflammation, multiple hemorrhages, and edema as early as 3 dpi up to 14 dpi in all the animals except mock control animals with severe gross lung pathology. Previous studies have shown that the SARS-CoV-2-infection-related lung pathologies in hamsters resemble those of human lung pathologies including inflammation, alveolar epithelial injury, and pneumonitis [30,31,32]. The histopathology study showed inflammation as early as 2 dpi in both male and female hamsters. The onset of inflammatory lesions as early as 2 dpi indicates acute lung damage and increased severity by 4 dpi up to 7 dpi (Figure 6). Previous studies have shown widespread necrosis and inflammation at 2 dpi and 4 dpi [33,34]; further, alveolar and interstitial macrophage infiltration and type II hyperplasia were observed at 7 dpi [35]. The evaluation of lung histopathology revealed no significant differences among the sexes of the animals. A study by Plunkard et al. showed largely consistent lung pathologies across different variants of SARS-CoV-2 [27].
The NGS of the SARS-CoV-2 genome was conducted to determine the SNPs between the virus stock used in infection versus the virus passed through both male and female hamsters. The NGS data showed 8 SNPs in the SARS-CoV-2 (ancestral Wuhan strain) virus stock propagated in the Vero E6 cells and used in the infection of hamsters. Out of 8 SNPs, 5 were in the spike protein gene, 2 were in the ORF1ab, and 1 was in ORF8 from female and male hamsters (Table 2). Earlier studies showed the accumulation of adaptive mutations in the spike protein of the SARS-CoV-2 genome after different passages in Vero E6 cells [36,37]. The NGS data showed 10 SNPs in the SRAS-CoV-2 (ancestral Wuhan strain) genome isolated from both female and male hamsters compared to the virus stock used in the infection of hamsters. The number of mutations in the spike protein was found to be higher compared to other proteins such as ORF1ab. Mutations in the S protein of SARS-CoV-2 at key residues play a significant role in the interaction with the ACE2 receptor. The E484 of the Receptor binding domain (RBD) has been identified as an important position, and a change in the amino acid from E to K, Q, or P is known to reduce the antibody neutralization and change receptor affinity [38,39,40,41]. Further, the exposure of SARS-CoV-2 to mAbs has identified E484K as an escape mutation and only this mutation reduced the neutralizing ability of the combination of monoclonal antibodies (Mabs) [42,43]. Several earlier studies demonstrated a partial or complete loss of antiviral activity due to mutations in the SARS-CoV-2 genome [44]. An earlier study by Sonnleitner et al. reported the immune escape mutations in the SARS-CoV-2 genome from immunocompromised patients using NGS and the results from this study suggest that specific mutations may serve as targets for the development of vaccines and treatments for COVID-19 in the future [45]. Furthermore, a majority (five SNPs) of the mutations that were found in the virus stock (propagated in Vero E6 cells and used in the infection of hamsters) were replaced by the wildtype variant (NC_045512.2) in the viral genome isolated from hamsters in both sexes (Table S1). A significant difference in mutation and associated protein changes was not observed in the SARS-CoV-2 genome isolated from male and female hamsters, except for a mutation at position 22296 (Arg245His) in males, while in females, no mutation was observed (Table 1 and Table 2). The lack of significant difference in the genetic mutations in the SARS-CoV-2 (ancestral Wuhan) genome in male and female hamsters along with similar lung pathologies provides the basis for using either of the sexes in pre-clinical efficacy testing involving small molecules, vaccines, and MAbs. However, the identification of the mutation sites in the S protein after passing through hamsters is crucial in the preclinical efficacy study of antiviral agents targeting the S protein including MAbs and small-molecule inhibitors in hamsters.
Major features such as sequence generation, assembly, alignment, variant identification, annotation, prioritization, and visualization are all shared by NGS [46,47,48] and require stringent quality control to warrant the validity of results. Earlier studies have shown variations in the frequency of SNPs ranging from 2.72% to 96.74% and variable mapping rates in the viral genome and this depends on the quality of the sample/storage conditions [49,50]. Further, the low-frequency variants reflect contaminations between samples, intra-host heterogeneity, coinfection with different viral strains, or PCR polymerase errors [51,52]. Although the present study provides information on the mutation sites in SARS-CoV-2 after passing through the hamsters, understanding the genome mutations in the virus used in the preclinical efficacy study is crucial in determining the antiviral effect of vaccines, small molecules, and Mabs including antiviral drug resistance.

5. Conclusions

The continued emergence of new variants of SARS-CoV-2 emphasizes the concerted efforts and need to understand the SARS-CoV-2 pathology. Our study findings indicate the accumulation of genomic mutations in SARS-CoV-2 after passing through the Syrian golden hamsters with no differences among the sexes of the hamsters. Hence, understanding the genomic mutations in SARS-CoV-2 following infection in animals plays a significant role in testing the efficacy of antiviral compounds and vaccines in hamster animal models.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/pathogens12111328/s1, Table S1: Comparison of the SNPs identified in the virus stock used in infection and virus genome from infected female and male hamsters with the reference sequence (NC_045512.2).

Author Contributions

Conceptualization, G.R.R., R.K.S. and S.N.; Investigation, G.R.R., C.N.N. and K.B.; Methodology, G.R.R., C.N.N. and K.B.; Writing—Original Draft, G.R.R.; Writing—Review and Editing, G.R.R., R.K.S. and S.N. All authors have read and agreed to the published version of the manuscript.

Funding

This research received funding from the Biotechnology Industry Research Assistance Council (BIRAC), India; Grant number: BT/CTH0004/CTH/01/20.

Institutional Review Board Statement

The animal study protocol was approved by the Institutional Review Board (or Ethics Committee) of the Foundation for Neglected Disease Research; FNDR-FB-071-PD-2021 (22-06-2021) and FNDR-FB-083-PD-2021 (23-12-2021).

Informed Consent Statement

Not applicable.

Data Availability Statement

The datasets generated in the current study are available from the corresponding author upon reasonable request.

Acknowledgments

The following reagent was deposited by the Centers for Disease Control and Prevention and obtained through BEI Resources, NIAID, NIH: SARS-Related Coronavirus 2, Isolate USA-WA1/2020, NR-52281. The following reagent was obtained through BEI Resources, NIAID, NIH: SARS-Related Coronavirus 2, Isolate hCoV-19/USA/PHC658/2021 (Lineage B.1.617.2; Delta Variant), NR-55611, contributed by Richard Webby and Anami Patel.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Graham, R.L.; Baric, R.S. Recombination, Reservoirs, and the Modular Spike: Mechanisms of Coronavirus Cross-Species Transmission. J. Virol. 2010, 84, 3134–3146. [Google Scholar] [CrossRef]
  2. Anthony, S.J.; Gilardi, K.; Menachery, V.D.; Goldstein, T.; Ssebide, B.; Mbabazi, R.; Navarrete-Macias, I.; Liang, E.; Wells, H.; Hicks, A.; et al. Further Evidence for Bats as the Evolutionary Source of Middle East Respiratory Syndrome Coronavirus. mBio 2017, 8, e00373-17. [Google Scholar] [CrossRef]
  3. Hu, B.; Ge, X.; Wang, L.-F.; Shi, Z. Bat Origin of Human Coronaviruses. Virol. J. 2015, 12, 221. [Google Scholar] [CrossRef] [PubMed]
  4. Huynh, J.; Li, S.; Yount, B.; Smith, A.; Sturges, L.; Olsen, J.C.; Nagel, J.; Johnson, J.B.; Agnihothram, S.; Gates, J.E.; et al. Evidence Supporting a Zoonotic Origin of Human Coronavirus Strain NL63. J. Virol. 2012, 86, 12816–12825. [Google Scholar] [CrossRef]
  5. Murgolo, N.; Therien, A.G.; Howell, B.; Klein, D.; Koeplinger, K.; Lieberman, L.A.; Adam, G.C.; Flynn, J.; McKenna, P.; Swaminathan, G.; et al. SARS-CoV-2 Tropism, Entry, Replication, and Propagation: Considerations for Drug Discovery and Development. PLoS Pathog. 2021, 17, e1009225. [Google Scholar] [CrossRef] [PubMed]
  6. Huang, C.; Wang, Y.; Li, X.; Ren, L.; Zhao, J.; Hu, Y.; Zhang, L.; Fan, G.; Xu, J.; Gu, X.; et al. Clinical Features of Patients Infected with 2019 Novel Coronavirus in Wuhan, China. Lancet 2020, 395, 497–506. [Google Scholar] [CrossRef] [PubMed]
  7. Zhou, Z.; Ren, L.; Zhang, L.; Zhong, J.; Xiao, Y.X.; Jia, Z.; Guo, L.; Yang, J.; Wang, C.; Jiang, S.; et al. Heightened Innate Immune Responses in the Respiratory Tract of COVID-19 Patients. Cell Host Microbe 2020, 27, 883–890.e2. [Google Scholar] [CrossRef] [PubMed]
  8. Sun, G.; Cui, Q.; Garcia, G.; Wang, C.; Zhang, M.; Arumugaswami, V.; Riggs, A.D.; Shi, Y. Comparative Transcriptomic Analysis of SARS-CoV-2 Infected Cell Model Systems Reveals Differential Innate Immune Responses. Sci. Rep. 2021, 11, 17146. [Google Scholar] [CrossRef] [PubMed]
  9. Krause, P.R.; Fleming, T.R.; Longini, I.M.; Peto, R.; Briand, S.; Heymann, D.L.; Beral, V.; Snape, M.D.; Rees, H.; Ropero, A.-M.; et al. SARS-CoV-2 Variants and Vaccines. N. Engl. J. Med. 2021, 385, 179–186. [Google Scholar] [CrossRef] [PubMed]
  10. WHO. Update on Omicron. 2021. Available online: https://www.who.int/fr/news/item/28-11-2021-update-on-omicron (accessed on 5 November 2023).
  11. Mohandas, S.; Dhruv Yadav, P.; Shete, A.; Nyayanit, D.; Sapkal, G.; Lole, K.; Gupta, N. SARS-CoV-2 Delta Variant Pathogenesis and Host Response in Syrian Hamsters. Viruses 2021, 13, 1773. [Google Scholar] [CrossRef] [PubMed]
  12. Sheahan, T.; Sims, A.; Zhou, S.; Graham, R.; Hill, C.; Leist, S.; Schäfer, A.; Dinnon, K.; Montgomery, S.; Agostini, M.; et al. An Orally Bioavailable Broad-Spectrum Antiviral Inhibits SARS-CoV-2 and Multiple Endemic, Epidemic and Bat Coronavirus. bioRxiv 2020. [Google Scholar] [CrossRef]
  13. Kaptein, S.J.; Jacobs, S.; Langendries, L.; Seldeslachts, L.; ter Horst, S.; Liesenborghs, L.; Hens, B.; Vergote, V.; Heylen, E.; Maas, E.; et al. Antiviral Treatment of SARS-CoV-2-Infected Hamsters Reveals a Weak Effect of Favipiravir and a Complete Lack of Effect for Hydroxychloroquine. bioRxiv 2020. [Google Scholar] [CrossRef]
  14. Muñoz-Fontela, C.; Dowling, W.E.; Funnell, S.G.P.; Gsell, P.S.; Riveros-Balta, A.X.; Albrecht, R.A.; Andersen, H.; Baric, R.S.; Carroll, M.W.; Cavaleri, M.; et al. Animal Models for COVID-19. Nature 2020, 586, 509–515. [Google Scholar] [CrossRef] [PubMed]
  15. Liu, Y.; Hu, G.; Wang, Y.; Ren, W.; Zhao, X.; Ji, F.; Zhu, Y.; Feng, F.; Gong, M.; Ju, X.; et al. Functional and Genetic Analysis of Viral Receptor ACE2 Orthologs Reveals a Broad Potential Host Range of SARS-CoV-2. Proc. Natl. Acad. Sci. USA 2021, 118, e2025373118. [Google Scholar] [CrossRef]
  16. Roberts, A.; Vogel, L.; Guarner, J.; Hayes, N.; Murphy, B.; Zaki, S.; Subbarao, K. Severe Acute Respiratory Syndrome Coronavirus Infection of Golden Syrian Hamsters. J. Virol. 2005, 79, 503–511. [Google Scholar] [CrossRef]
  17. Iwatsuki-Horimoto, K.; Nakajima, N.; Ichiko, Y.; Sakai-Tagawa, Y.; Noda, T.; Hasegawa, H.; Kawaoka, Y. Syrian Hamster as an Animal Model for the Study of Human Influenza Virus Infection. J. Virol. 2018, 92, 1–14. [Google Scholar] [CrossRef] [PubMed]
  18. Mehla, R.; Kokate, P.; Bhosale, S.R.; Vaidya, V.; Narayanan, S.; Shandil, R.K.; Singh, M.; Rudramurthy, G.R.; Naveenkumar, C.N.; Bharathkumar, K.; et al. A Live Attenuated COVID-19 Candidate Vaccine for Children: Protection against SARS-CoV-2 Challenge in Hamsters. Vaccines 2023, 11, 255. [Google Scholar] [CrossRef] [PubMed]
  19. Chan, J.F.W.; Zhang, A.J.; Yuan, S.; Poon, V.K.M.; Chan, C.C.S.; Lee, A.C.Y.; Chan, W.M.; Fan, Z.; Tsoi, H.W.; Wen, L.; et al. Simulation of the Clinical and Pathological Manifestations of Coronavirus Disease 2019 (COVID-19) in a Golden Syrian Hamster Model: Implications for Disease Pathogenesis and Transmissibility. Clin. Infect. Dis. 2020, 71, 2428–2446. [Google Scholar] [CrossRef]
  20. Case, J.B.; Bailey, A.L.; Kim, A.S.; Chen, R.E.; Diamond, M.S. Growth, Detection, Quantification, and Inactivation of SARS-CoV-2. Virology 2020, 548, 39–48. [Google Scholar] [CrossRef]
  21. Jureka, A.S.; Silvas, J.A.; Basler, C.F. Propagation, Inactivation, and Safety Testing of SARS-CoV-2. Viruses 2020, 12, 622. [Google Scholar] [CrossRef]
  22. Rudramurthy, G.R.; Shandil, R.K.; Narayanana, S. In-Vitro Screening of Repurposed Drug Library against Severe Acute Respiratory Syndrome Coronavirus-2. Med. Res. Arch. 2023, 11, 1–13. [Google Scholar] [CrossRef]
  23. Reed, L.J.; Muench, H. A Simple Method of Estimating Fifty per Cent Endpoints. Am. J. Epidemiol. 1938, 27, 493–497. [Google Scholar] [CrossRef]
  24. Dhakal, S.; Ruiz-Bedoya, C.A.; Zhou, R.; Creisher, P.S.; Villano, J.S.; Littlefield, K.; Castillo, J.R.; Marinho, P.; Jedlicka, A.E.; Ordonez, A.A.; et al. Sex Differences in Lung Imaging and SARS-CoV-2 Antibody Responses in a COVID-19 Golden Syrian Hamster Model. mBio 2021, 12, e0097421. [Google Scholar] [CrossRef]
  25. Yuan, S.; Yin, X.; Meng, X.; Chan, J.F.W.; Ye, Z.W.; Riva, L.; Pache, L.; Chan, C.C.Y.; Lai, P.M.; Chan, C.C.S.; et al. Clofazimine Broadly Inhibits Coronaviruses Including SARS-CoV-2. Nature 2021, 593, 418–423. [Google Scholar] [CrossRef] [PubMed]
  26. Osterrieder, N.; Bertzbach, L.D.; Dietert, K.; Abdelgawad, A.; Vladimirova, D.; Kunec, D.; Hoffmann, D.; Beer, M.; Gruber, A.D.; Trimpert, J. Age-Dependent Progression of SARS-CoV-2 Infection in Syrian Hamsters. Viruses 2020, 12, 779. [Google Scholar] [CrossRef]
  27. Plunkard, J.; Mulka, K.; Zhou, R.; Tarwater, P.; Zhong, W.; Lowman, M.; Wong, A.; Pekosz, A.; Villano, J. SARS-CoV-2 Variant Pathogenesis Following Primary Infection and Reinfection in Syrian Hamsters. mBio 2023, 14, e0007823. [Google Scholar] [CrossRef] [PubMed]
  28. Lee, A.C.Y.; Zhang, A.J.; Chan, J.F.W.; Li, C.; Fan, Z.; Liu, F.; Chen, Y.; Liang, R.; Sridhar, S.; Cai, J.P.; et al. Oral SARS-CoV-2 Inoculation Establishes Subclinical Respiratory Infection with Virus Shedding in Golden Syrian Hamsters. Cell Rep. Med. 2020, 1, 100121. [Google Scholar] [CrossRef] [PubMed]
  29. Imai, M.; Iwatsuki-Horimoto, K.; Hatta, M.; Loeber, S.; Halfmann, P.J.; Nakajima, N.; Watanabe, T.; Ujie, M.; Takahashi, K.; Ito, M.; et al. Syrian Hamsters as a Small Animal Model for SARS-CoV-2 Infection and Countermeasure Development. Proc. Natl. Acad. Sci. USA 2020, 117, 16587–16595. [Google Scholar] [CrossRef] [PubMed]
  30. O’Donnell, K.L.; Pinski, A.N.; Clancy, C.S.; Gourdine, T.; Shifflett, K.; Fletcher, P.; Messaoudi, I.; Marzi, A. Pathogenic and Transcriptomic Differences of Emerging SARS-CoV-2 Variants in the Syrian Golden Hamster Model. EBioMedicine 2021, 73, 103675. [Google Scholar] [CrossRef]
  31. Boudewijns, R.; Thibaut, H.J.; Kaptein, S.J.F.; Li, R.; Vergote, V.; Seldeslachts, L.; Van Weyenbergh, J.; De Keyzer, C.; Bervoets, L.; Sharma, S.; et al. STAT2 Signaling Restricts Viral Dissemination but Drives Severe Pneumonia in SARS-CoV-2 Infected Hamsters. Nat. Commun. 2020, 11, 5838. [Google Scholar] [CrossRef]
  32. Sia, S.F.; Yan, L.M.; Chin, A.W.H.; Fung, K.; Choy, K.T.; Wong, A.Y.L.; Kaewpreedee, P.; Perera, R.A.P.M.; Poon, L.L.M.; Nicholls, J.M.; et al. Pathogenesis and Transmission of SARS-CoV-2 in Golden Hamsters. Nature 2020, 583, 834–838. [Google Scholar] [CrossRef] [PubMed]
  33. Gruber, A.D.; Firsching, T.C.; Trimpert, J.; Dietert, K. Hamster Models of COVID-19 Pneumonia Reviewed: How Human Can They Be? Vet. Pathol. 2022, 59, 528–545. [Google Scholar] [CrossRef]
  34. Tostanoski, L.H.; Wegmann, F.; Martinot, A.J.; Loos, C.; McMahan, K.; Mercado, N.B.; Yu, J.; Chan, C.N.; Bondoc, S.; Starke, C.E.; et al. Ad26 Vaccine Protects against SARS-CoV-2 Severe Clinical Disease in Hamsters. Nat. Med. 2020, 26, 1694–1700. [Google Scholar] [CrossRef] [PubMed]
  35. Mulka, K.R.; Beck, S.E.; Solis, C.V.; Johanson, A.L.; Queen, S.E.; McCarron, M.E.; Richardson, M.R.; Zhou, R.; Marinho, P.; Jedlicka, A.; et al. Progression and Resolution of Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) Infection in Golden Syrian Hamsters. Am. J. Pathol. 2022, 192, 195–207. [Google Scholar] [CrossRef]
  36. Ogando, N.S.; Dalebout, T.J.; Zevenhoven-Dobbe, J.C.; Limpens, R.W.A.L.; van der Meer, Y.; Caly, L.; Druce, J.; de Vries, J.J.C.; Kikkert, M.; Barcena, M.; et al. SARS-Coronavirus-2 Replication in Vero E6 Cells: Replication Kinetics, Rapid Adaptation and Cytopathology. J. Gen. Virol. 2020, 101, 925–940. [Google Scholar] [CrossRef] [PubMed]
  37. Chen, Y.; Liu, M.-Q.; Luo, Y.; Jiang, R.-D.; Si, H.-R.; Zhu, Y.; Li, B.; Shen, X.-R.; Lin, H.-F.; Zhao, K.; et al. Genetic Mutation of SARS-CoV-2 during Consecutive Passages in Permissive Cells. Virol. Sin. 2021, 36, 1073–1076. [Google Scholar] [CrossRef] [PubMed]
  38. Huang, Y.; Yang, C.; Xu, X.F.; Xu, W.; Liu, S.W. Structural and Functional Properties of SARS-CoV-2 Spike Protein: Potential Antivirus Drug Development for COVID-19. Acta Pharmacol. Sin. 2020, 41, 1141–1149. [Google Scholar] [CrossRef]
  39. Wrapp, D.; Wang, N.; Corbett, K.S.; Goldsmith, J.A.; Hsieh, C.L.; Abiona, O.; Graham, B.S.; McLellan, J.S. Cryo-EM Structure of the 2019-NCoV Spike in the Prefusion Conformation. Science 2020, 367, 1264–1269. [Google Scholar] [CrossRef]
  40. Harvey, W.T.; Carabelli, A.M.; Jackson, B.; Gupta, R.K.; Thomson, E.C.; Harrison, E.M.; Ludden, C.; Reeve, R.; Rambaut, A.; Peacock, S.J.; et al. SARS-CoV-2 Variants, Spike Mutations and Immune Escape. Nat. Rev. Microbiol. 2021, 19, 409–424. [Google Scholar] [CrossRef]
  41. Greaney, A.J.; Loes, A.N.; Crawford, K.H.D.; Starr, T.N.; Malone, K.D.; Chu, H.Y.; Bloom, J.D. Comprehensive Mapping of Mutations in the SARS-CoV-2 Receptor-Binding Domain That Affect Recognition by Polyclonal Human Plasma Antibodies. Cell Host Microbe 2021, 29, 463–476.e6. [Google Scholar] [CrossRef]
  42. Weisblum, Y.; Schmidt, F.; Zhang, F.; DaSilva, J.; Poston, D.; Lorenzi, J.C.C.; Muecksch, F.; Rutkowska, M.; Hoffmann, H.H.; Michailidis, E.; et al. Escape from Neutralizing Antibodies 1 by SARS-CoV-2 Spike Protein Variants. eLife 2020, 9, e61312. [Google Scholar] [CrossRef]
  43. Baum, A.; Fulton, B.O.; Wloga, E.; Copin, R.; Pascal, K.E.; Russo, V.; Giordano, S.; Lanza, K.; Negron, N.; Ni, M.; et al. Antibody Cocktail to SARS-CoV-2 Spike Protein Prevents Rapid Mutational Escape Seen with Individual Antibodies. Science 2020, 369, 1014–1018. [Google Scholar] [CrossRef]
  44. European Centre for Disease Prevention and Control. SARS-CoV-2 Variant Mutations Conferring Reduced Susceptibility to Antiviral Drugs and Monoclonal Antibodies: A Non-Systematic Literature Review for Surveillance Purposes; European Centre for Disease Prevention and Control: Stockholm, Sweden, 2023. [Google Scholar]
  45. Sonnleitner, S.T.; Prelog, M.; Sonnleitner, S.; Hinterbichler, E.; Halbfurter, H.; Kopecky, D.B.C.; Almanzar, G.; Koblmüller, S.; Sturmbauer, C.; Feist, L.; et al. Cumulative SARS-CoV-2 Mutations and Corresponding Changes in Immunity in an Immunocompromised Patient Indicate Viral Evolution within the Host. Nat. Commun. 2022, 13, 2560. [Google Scholar] [CrossRef] [PubMed]
  46. John, G.; Sahajpal, N.S.; Mondal, A.K.; Ananth, S.; Williams, C.; Chaubey, A.; Rojiani, A.M.; Kolhe, R. Next-Generation Sequencing (Ngs) in Covid-19: A Tool for Sars-Cov-2 Diagnosis, Monitoring New Strains and Phylodynamic Modeling in Molecular Epidemiology. Curr. Issues Mol. Biol. 2021, 43, 845–867. [Google Scholar] [CrossRef] [PubMed]
  47. Roy, S.; Coldren, C.; Karunamurthy, A.; Kip, N.S.; Klee, E.W.; Lincoln, S.E.; Leon, A.; Pullambhatla, M.; Temple-Smolkin, R.L.; Voelkerding, K.V.; et al. Standards and Guidelines for Validating Next-Generation Sequencing Bioinformatics Pipelines: A Joint Recommendation of the Association for Molecular Pathology and the College of American Pathologists. J. Mol. Diagn. 2018, 20, 4–27. [Google Scholar] [CrossRef] [PubMed]
  48. Wadapurkar, R.M.; Vyas, R. Computational Analysis of next Generation Sequencing Data and Its Applications in Clinical Oncology. Inform. Med. Unlocked 2018, 11, 75–82. [Google Scholar] [CrossRef]
  49. Bin Mahmood, T.; Saha, A.; Hossan, M.I.; Mizan, S.; Arman, S.M.A.S.; Chowdhury, A.S. A next Generation Sequencing (NGS) Analysis to Reveal Genomic and Proteomic Mutation Landscapes of SARS-CoV-2 in South Asia. Curr. Res. Microb. Sci. 2021, 2, 100065. [Google Scholar] [CrossRef]
  50. Liu, T.; Chen, Z.; Chen, W.; Chen, X.; Hosseini, M.; Yang, Z.; Li, J.; Ho, D.; Turay, D.; Gheorghe, C.P.; et al. A Benchmarking Study of SARS-CoV-2 Whole-Genome Sequencing Protocols Using COVID-19 Patient Samples. iScience 2021, 24, 102892. [Google Scholar] [CrossRef]
  51. Jacot, D.; Pillonel, T.; Greub, G.; Bertelli, C. Assessment of SARS-CoV-2 Genome Sequencing: Quality Criteria and Low-Frequency Variants. J. Clin. Microbiol. 2021, 59, e0094421. [Google Scholar] [CrossRef] [PubMed]
  52. Kebschull, J.M.; Zador, A.M. Sources of PCR-Induced Distortions in High-Throughput Sequencing Data Sets. Nucleic Acids Res. 2015, 43, e143. [Google Scholar] [CrossRef]
Figure 1. Schematic representation depicting the study design; the Syrian golden hamsters (male and female) were infected with SARS-CoV-2 (ancestral Wuhan). The animals were sacrificed on different days post-infection (dpi) to collect the lung and blood samples for viral load estimation (lung) and neutralizing antibody estimation (serum).
Figure 1. Schematic representation depicting the study design; the Syrian golden hamsters (male and female) were infected with SARS-CoV-2 (ancestral Wuhan). The animals were sacrificed on different days post-infection (dpi) to collect the lung and blood samples for viral load estimation (lung) and neutralizing antibody estimation (serum).
Pathogens 12 01328 g001
Figure 2. Percentage body weight change. (A,B): male and female hamsters infected with the SARS-CoV-2 (ancestral Wuhan strain), respectively. The body weight of the animals was recorded from day 0 (before infection) until day 14 post-infection. The infected animals showed up to a 5% decrease in body weight up to 5 dpi; however, mock animals showed an increase in body weight. Later, the animals gained body weight with a significant difference in body weight between mock and infected groups in both sexes of the animals. *** p < 0.001 and ** p < 0.01 using a two-tailed unpaired t-test.
Figure 2. Percentage body weight change. (A,B): male and female hamsters infected with the SARS-CoV-2 (ancestral Wuhan strain), respectively. The body weight of the animals was recorded from day 0 (before infection) until day 14 post-infection. The infected animals showed up to a 5% decrease in body weight up to 5 dpi; however, mock animals showed an increase in body weight. Later, the animals gained body weight with a significant difference in body weight between mock and infected groups in both sexes of the animals. *** p < 0.001 and ** p < 0.01 using a two-tailed unpaired t-test.
Pathogens 12 01328 g002
Figure 3. Determination of infectious viral particles and viral RNA from lung homogenates: male and female hamsters infected with SARS-CoV-2 (ancestral Wuhan strain). Both TCID-50 and qRTPCR assay showed peak viral load at 4 dpi regardless of the sex of the animal with no significant difference (A,B). Infectious virus particles were not detected using TCID50 assay at 7 dpi and onward, while qRTPCR detected viral RNA up to 14 dpi with a significant difference between the assays (C,D). * p < 0.1, and ns, non-significant using a two-tailed unpaired t-test.
Figure 3. Determination of infectious viral particles and viral RNA from lung homogenates: male and female hamsters infected with SARS-CoV-2 (ancestral Wuhan strain). Both TCID-50 and qRTPCR assay showed peak viral load at 4 dpi regardless of the sex of the animal with no significant difference (A,B). Infectious virus particles were not detected using TCID50 assay at 7 dpi and onward, while qRTPCR detected viral RNA up to 14 dpi with a significant difference between the assays (C,D). * p < 0.1, and ns, non-significant using a two-tailed unpaired t-test.
Pathogens 12 01328 g003
Figure 4. Determination of neutralizing antibodies from infectious hamster serum samples: (A,B): male and female hamsters infected with SARS-CoV-2 (ancestral Wuhan strain), respectively. Neutralizing antibodies were detected as early as 4 dpi regardless of the sex of the animal with no significant difference. Further, the cross-protecting antibodies were observed regardless of the sex of the animals with no significant difference. ns, non-significant using a two-tailed unpaired t-test.
Figure 4. Determination of neutralizing antibodies from infectious hamster serum samples: (A,B): male and female hamsters infected with SARS-CoV-2 (ancestral Wuhan strain), respectively. Neutralizing antibodies were detected as early as 4 dpi regardless of the sex of the animal with no significant difference. Further, the cross-protecting antibodies were observed regardless of the sex of the animals with no significant difference. ns, non-significant using a two-tailed unpaired t-test.
Pathogens 12 01328 g004
Figure 5. Gross pathology of the lungs. Gross pathological changes were observed at 3 dpi and onward in all the groups. The gross pathology was found to be severe at 4 dpi and onward to 7 dpi and reduced at 14 dpi in all the groups. Pathogens 12 01328 i001 Represent the areas of congestion and hemorrhage and Pathogens 12 01328 i002 represent lung consolidation.
Figure 5. Gross pathology of the lungs. Gross pathological changes were observed at 3 dpi and onward in all the groups. The gross pathology was found to be severe at 4 dpi and onward to 7 dpi and reduced at 14 dpi in all the groups. Pathogens 12 01328 i001 Represent the areas of congestion and hemorrhage and Pathogens 12 01328 i002 represent lung consolidation.
Pathogens 12 01328 g005
Figure 6. Histopathology of the lungs. Thin sections of the tissue samples were stained by hematoxylin and eosin and then the images were captured at 10× and 60× magnification, respectively, using an eyepiece and objective. The histopathological features such as inflammation and necrosis started as early as 2 dpi, showing an acute inflammatory response with diffuse alveolar damage, hemorrhages, congestion (Pathogens 12 01328 i003), and an accumulation of eosinophilic edematous exudate (Pathogens 12 01328 i004). Moderate to marked broncho-interstitial pneumonia with alveolar damage and inflammatory cell infiltration (Pathogens 12 01328 i005), hemorrhage (Pathogens 12 01328 i006), and hyperplasia of the bronchial epithelium (Pathogens 12 01328 i007) at 3 and 4 dpi in both male and female hamsters. Severe broncho-interstitial pneumonia with a marked increase in lung cellularity with hyperplasia of Type II alveolar cells (Pathogens 12 01328 i008) was observed at 7 dpi. Mild congestion (Pathogens 12 01328 i009), infiltration of inflammatory cells (Pathogens 12 01328 i010), and mild hyperplasia of the bronchial epithelium (Pathogens 12 01328 i011) were observed at 14 dpi.
Figure 6. Histopathology of the lungs. Thin sections of the tissue samples were stained by hematoxylin and eosin and then the images were captured at 10× and 60× magnification, respectively, using an eyepiece and objective. The histopathological features such as inflammation and necrosis started as early as 2 dpi, showing an acute inflammatory response with diffuse alveolar damage, hemorrhages, congestion (Pathogens 12 01328 i003), and an accumulation of eosinophilic edematous exudate (Pathogens 12 01328 i004). Moderate to marked broncho-interstitial pneumonia with alveolar damage and inflammatory cell infiltration (Pathogens 12 01328 i005), hemorrhage (Pathogens 12 01328 i006), and hyperplasia of the bronchial epithelium (Pathogens 12 01328 i007) at 3 and 4 dpi in both male and female hamsters. Severe broncho-interstitial pneumonia with a marked increase in lung cellularity with hyperplasia of Type II alveolar cells (Pathogens 12 01328 i008) was observed at 7 dpi. Mild congestion (Pathogens 12 01328 i009), infiltration of inflammatory cells (Pathogens 12 01328 i010), and mild hyperplasia of the bronchial epithelium (Pathogens 12 01328 i011) were observed at 14 dpi.
Pathogens 12 01328 g006
Table 1. Mapping and consensus statistics of NGS data: The SARS-CoV2 isolate of the present study used in the infection of hamsters was mapped with reference sequence NC_045512.2 and the virus isolated from the hamsters was mapped with the stock virus (used in infection) to identify the single-nucleotide polymorphisms (SNPs) and lineage. The analysis identified that the ancestral Wuhan strain (used in infection) and virus isolated from hamsters belong to lineage A.
Table 1. Mapping and consensus statistics of NGS data: The SARS-CoV2 isolate of the present study used in the infection of hamsters was mapped with reference sequence NC_045512.2 and the virus isolated from the hamsters was mapped with the stock virus (used in infection) to identify the single-nucleotide polymorphisms (SNPs) and lineage. The analysis identified that the ancestral Wuhan strain (used in infection) and virus isolated from hamsters belong to lineage A.
SampleReference SequenceTotal Length of
Consensus (bp)
No. of SNPsLineage Identified
Virus stock used in the infection of HamstersNC_045512.2 299038A
Virus isolated from the infected hamster (Male) Virus stock (used in the infection of hamsters)2990310A
Virus isolated from the infected hamster (Female) 2990310A
Table 2. Comparison of SNPs identified in the stock virus (ancestral Wuhan strain) with the reference sequence (NC_045512.2): The SNPs identified in the viral RNA isolated from the virus stock (used in infection) were compared to the reference sequence (NC_045512.2) and the corresponding protein change is mentioned.
Table 2. Comparison of SNPs identified in the stock virus (ancestral Wuhan strain) with the reference sequence (NC_045512.2): The SNPs identified in the viral RNA isolated from the virus stock (used in infection) were compared to the reference sequence (NC_045512.2) and the corresponding protein change is mentioned.
Position in
Reference
Sequence (NC_045512.2)
Reported
Base in
Reference Sequence (NC_045512.2)
Identified
Alternate Base
in Virus Stock (Used in the Infection of Hamsters)
Gene NameProtein
Change
8782CTORF1abSer2839Ser
18060CTORF1abLeu5932Leu
21759AGSHis66Arg
22296AGSHis245Arg
22482CTSThr307Ile
23606CTSArg682Trp
23607GTSArg682Leu
28144TCORF8Leu84Ser
Table 3. Comparison of SNPs identified in the virus isolated from hamsters (male and female) with the virus stock (ancestral Wuhan strain) used in the infection of hamsters. Protein change remains the same in both males and females unless otherwise mentioned.
Table 3. Comparison of SNPs identified in the virus isolated from hamsters (male and female) with the virus stock (ancestral Wuhan strain) used in the infection of hamsters. Protein change remains the same in both males and females unless otherwise mentioned.
Position in
Stock Virus (Used in Infection)
Reported
Base in Stock Virus (Used in Infection)
Identified
Alternate Base
from an Infected Animal (Male)
Identified
Alternate Base
from Infected Animal (Female)
Gene NameProtein
Change
17827CAAORF1abGln5855Lys
21801AGGSAsp80Gly
21849A-CSGlu96Ala
22206AGGSAsp215Gly
23014ACCSGlu484Asp
23525CTTSHis655Tyr
24734CTTSHis1058Tyr
21759GAASArg66His
22296GAGSArg245His (Male)
23606TCCSTrp682Arg
23607TGGSLeu682Arg
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Rudramurthy, G.R.; Naveenkumar, C.N.; Bharathkumar, K.; Shandil, R.K.; Narayanan, S. Genomic Mutations in SARS-CoV-2 Genome following Infection in Syrian Golden Hamster and Associated Lung Pathologies. Pathogens 2023, 12, 1328. https://doi.org/10.3390/pathogens12111328

AMA Style

Rudramurthy GR, Naveenkumar CN, Bharathkumar K, Shandil RK, Narayanan S. Genomic Mutations in SARS-CoV-2 Genome following Infection in Syrian Golden Hamster and Associated Lung Pathologies. Pathogens. 2023; 12(11):1328. https://doi.org/10.3390/pathogens12111328

Chicago/Turabian Style

Rudramurthy, Gudepalya Renukaiah, Chakenahalli N. Naveenkumar, Kumaraswamy Bharathkumar, Radha K. Shandil, and Shridhar Narayanan. 2023. "Genomic Mutations in SARS-CoV-2 Genome following Infection in Syrian Golden Hamster and Associated Lung Pathologies" Pathogens 12, no. 11: 1328. https://doi.org/10.3390/pathogens12111328

APA Style

Rudramurthy, G. R., Naveenkumar, C. N., Bharathkumar, K., Shandil, R. K., & Narayanan, S. (2023). Genomic Mutations in SARS-CoV-2 Genome following Infection in Syrian Golden Hamster and Associated Lung Pathologies. Pathogens, 12(11), 1328. https://doi.org/10.3390/pathogens12111328

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop