Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) Mutational Pattern in the Fourth Pandemic Phase in Greece

The aim of this study is to investigate the circulating variants of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) from Athens and from rural areas in Greece during July and August 2021. We also present a rapid review of literature regarding significant SARS-CoV-2 mutations and their impact on public health. A total of 2500 nasopharyngeal swab specimens were collected from suspected COVID-19 cases (definition by WHO 2021b). Viral nucleic acid extraction was implemented using an automatic extractor and the RNA recovered underwent qRT-PCR in order to characterize the specimens as positive or negative for SARS-CoV-2. The positive specimens were then used to identify specific Spike gene mutations and characterize the emerging SARS-CoV-2 variants. For this step, various kits were utilized. From the 2500 clinical specimens, 220 were tested positive for SARS-CoV-2 indicating a prevalence of 8.8% among suspected cases. The RT-PCR Ct (Cycle threshold) Value ranged from 19 to 25 which corresponds to medium to high copy numbers of the virus in the positive samples. From the 220 positive specimens 148 (67.3%) were from Athens and 72 (32.7%) from Greek rural areas. As far as the Spike mutations investigated: N501Y appeared in all the samples, D614G mutation appeared in 212 (96.4%) samples with a prevalence of 87.2% in Athens and 98.6% in the countryside, E484K had a prevalence of 10.8% and 12.5% in Athens and the rural areas, respectively. K417N was found in 18 (12.2%) samples from Athens and four (5.6%) from the countryside, P681H was present in 51 (34.5%) Athenian specimens and 14 (19.4%) specimens from rural areas, HV69-70 was carried in 32.4% and 19.4% of the samples from Athens and the countryside, respectively. P681R had a prevalence of 87.2% in Athens and 98.6% in rural areas, and none of the specimens carried the L452R mutation. 62 (28.2%) samples carried the N501Y, P681H, D614G and HV69-70 mutations simultaneously and the corresponding variant was characterized as the Alpha (UK) variant (B 1.1.7). Only six (2.7%) samples from the center of Athens had the N501Y, E484K, K417N and D614G mutations simultaneously and the virus responsible was characterized as the Beta (South African) variant (B 1.351). Our study explored the SARS-CoV-2 variants using RT-PCR in a representative cohort of samples collected from Greece in July and August 2021. The prevalent mutations identified were N501Y (100%), D614G (96.4%), P681R (90.1%) and the variants identified were the Delta (90.1%), Alpha (28.2%) and Beta (2.7%).     

D614G substitution at the hinge region enhances the stability of trimeric SARS-CoV-2 spike protein

Mutations in the spike protein of SARS-CoV-2 are the major causes for the modulation of ongoing COVID-19 infection. Currently, the D614G substitution in the spike protein has become dominant worldwide. It is associated with higher infectivity than the ancestral (D614)variant. We demonstrate using Gaussian network model-based normal mode analysis that the D614G substitution occurs at the hinge region that facilitates domain-domain motions between receptor binding domain and S2 region of the spike protein. Computer-aided mutagenesis and inter-residue energy calculations reveal that contacts involving D614 are energetically frustrated. However, contacts involving G614 are energetically favourable, implying the substitution strengthens residue contacts that are formed within as well as between protomers. We also find that the free energy difference (ΔΔG) between two variants is -2.6 kcal/mol for closed and -2.0 kcal/mol for 1-RBD up conformation. Thus, the thermodynamic stability has increased upon D614G substitution. Whereas the reverse mutation in spike protein structures having G614 substitution has resulted in the free energy differences of 6.6 kcal/mol and 6.3 kcal/mol for closed and 1-RBD up conformations, respectively, indicating that the overall thermodynamic stability has decreased. These results suggest that the D614G substitution modulates the flexibility of spike protein and confers enhanced thermodynamic stability irrespective of conformational states. This data concurs with the known information demonstrating increased availability of the functional form of spikeprotein trimer upon D614G substitution.     

B.1.617.2 (Delta) Variant of SARS-CoV-2: features,transmission and potential strategies

Coronavirus disease 2019 (COVID-19) caused by the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has become a pandemic. With the continuous evolution of the viral genome, SARS-CoV-2 has evolved many variants. B.1.617.2, also called Delta, is one of the most concerned variants. The Delta variant was first reported in India at the end of 2020 but has spread globally, by now, to 135 countries and is not stand still. Delta shared some mutations with other variants, and owned its special mutations on spike proteins, which may be responsible for its strong transmission and increasing virulence. Under these circumstances, a systematic summary of Delta is necessary. This review will focus on the Delta variant. We will describe all the characteristics of Delta (including biological features and clinical characteristics), analyze potential reasons for its strong transmission, and provide potential protective ways for combating Delta.     

Cross-variant protection against SARS-CoV-2 infection in hamsters immunized with monovalent and bivalent inactivated vaccines

Rapid development and successful use of vaccines against SARS-CoV-2 might hold the key to curb the ongoing pandemic of COVID-19. Emergence of vaccine-evasive SARS-CoV-2 variants of concern (VOCs) has posed a new challenge to vaccine design and development. One urgent need is to determine what types of variant-specific and bivalent vaccines should be developed. Here, we compared homotypic and heterotypic protection against SARS-CoV-2 infection of hamsters with monovalent and bivalent whole-virion inactivated vaccines derived from representative VOCs. In addition to the ancestral SARS-CoV-2 Wuhan strain, Delta (B.1.617.2; δ) and Theta (P.3; θ) variants were used in vaccine preparation. Additional VOCs including Omicron (B.1.1.529) and Alpha (B.1.1.7) variants were employed in the challenge experiment. Consistent with previous findings, Omicron variant exhibited the highest degree of immune evasion, rendering all different forms of inactivated vaccines substantially less efficacious. Notably, monovalent and bivalent Delta variant-specific inactivated vaccines provided optimal protection against challenge with Delta variant. Yet, some cross-variant protection against Omicron and Alpha variants was seen with all monovalent and bivalent inactivated vaccines tested. Taken together, our findings support the notion that an optimal next-generation inactivated vaccine against SARS-CoV-2 should contain the predominant VOC in circulation. Further investigations are underway to test whether a bivalent vaccine for Delta and Omicron variants can serve this purpose.     

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

Severe COVID‐19 is characterized by lung abnormalities, including the presence of syncytial pneumocytes. Syncytia form when SARS‐CoV‐2 spike protein expressed on the surface of infected cells interacts with the ACE2 receptor on neighboring cells. The syncytia forming potential of spike variant proteins remain poorly characterized. Here, we first assessed Alpha (B.1.1.7) and Beta (B.1.351) spread and fusion in cell cultures, compared with the ancestral D614G strain. Alpha and Beta replicated similarly to D614G strain in Vero, Caco‐2, Calu‐3, and primary airway cells. However, Alpha and Beta formed larger and more numerous syncytia. Variant spike proteins displayed higher ACE2 affinity compared with D614G. Alpha, Beta, and D614G fusion was similarly inhibited by interferon‐induced transmembrane proteins (IFITMs). Individual mutations present in Alpha and Beta spikes modified fusogenicity, binding to ACE2 or recognition by monoclonal antibodies. We further show that Delta spike also triggers faster fusion relative to D614G. Thus, SARS‐CoV‐2 emerging variants display enhanced syncytia formation.     

Mutational sensitivity of D614G in spike protein of SARS-CoV-2 in Jordan

BackgroundSpike protein is the surface glycoprotein of the severe acute respiratory syndrome-coronavirus-2 (SARS-CoV-2) necessary for the entry of the virus via the transmembrane receptors of the human respiratory cells causing COVID-19 disease.AimHere, we aimed to predict the three-dimensional monomer structure of spike protein of SARS-CoV-2 from 20 Jordanian nasopharyngeal samples and to determine the percentage of single amino acid variants (SAV) in the spike protein of SARS-CoV-2.MethodsThe output of the Protein Homology/analogY Recognition Engine V 2.0 (Phyre2) found four single amino acid variants in the spike gene.ResultsThe first variant represented by 5% of samples that showed tyrosine deletion at Y144 located in the N terminal domain. The second and the dominant variant, represented by 62%, showed aspartate a coil amino acid substitution to glycine an extracellular amino acid at D614G located in the spike recognition binding site. The third variant, represented by 5%, showed aspartate substitution to tyrosine at D1139Y, and the fourth variant, represented by 5% glycine substitution to serine at G1167S.ConclusionOur results have shown low mutational sensitivity in all variants except to D614G the one with the most likely neutral mutational sensitivity that all variants might not explicitly affect the function of spike glycoprotein. However, D614G might change the viral conformational plasticity and hence a potential viral fitness gain but one must be cautious about drawing any concrete conclusions about the severity of symptoms and viral transmission from genomic data only.General significanceStudying mutations such as D614G in deep is essential to control the pandemic in terms of immune systems, antibodies, or even vaccines.     

Detection of the SARS-CoV-2 D614G mutation using engineered Cas12a guide RNA

Detection of pathogens with single-nucleotide variations is indispensable for the disease tracing, but remains technically challenging. The D614G mutation in the SARS-CoV-2 spike protein is known to markedly enhance viral infectivity but is difficult to detect. Here, we report an effective approach called “synthetic mismatch integrated crRNA guided Cas12a detection” (symRNA-Cas12a) to detect the D614 and G614 variants effectively. Using this method, we systemically screened a pool of crRNAs that contain all the possible nucleotide substitutions covering the -2 to +2 positions around the mutation and identify one crRNA that can efficiently increase the detection specificity by 13-fold over the ancestral crRNA. With this selected crRNA, the symRNA-Cas12a assay can detect as low as 10 copies of synthetic mutant RNA and the results are confirmed to be accurate by Sanger sequencing. Overall, we have developed the symRNA-Cas12a method to specifically, sensitively and rapidly detect the SARS-CoV-2 D614G mutation.     

High-resolution melting analysis to discriminate between the SARS-CoV-2 Omicron variants BA.1 and BA.2

High-resolution melting (HRM) analysis was conducted to discriminate between SARS-CoV-2 Omicron variant BA.1 (B.1.1.529.1) and subvariant BA.2 (B.1.1.529.2). We performed two-step PCR consisting of the first PCR and the second nested PCR to prepare the amplicon for HRM analysis, which detected G339D, N440K, G446S and D796Y variations in the SARS-CoV-2 spike protein. The melting temperatures (Tms) of the amplicons from the cDNA of the Omicron variant BA.1 and subvariant BA.2 receptor binding domain (RBD) in spike protein were the same: 75.2 °C (G339D variation) and 73.4 °C (D796Y variation). These Tms were distinct from those of SARS-CoV-2 isolate Wuhan-Hu-1, and were specific to the Omicron variant. In HRM analyses that detected the N440K and G446S variations, the Tms of amplicons from the cDNA of the Omicron variant BA.1 and subvariant BA.2 RBDs were 73.0 °C (N440K and G446S variations) and 73.5 °C (G446S variation). This difference indicates that the SARS-CoV-2 Omicron variants BA.1 and BA.2 can be clearly discriminated. Our study demonstrates the usefulness of HRM analysis after two-step PCR for the discrimination of SARS-CoV-2 variants.     

Non-synonymous mutations of SARS-CoV-2 leads epitope loss and segregates its variants

The non-synonymous mutations of SARS-CoV-2 isolated from across the world have been identified during the last few months. The surface glycoprotein spike of SARS-CoV-2 forms the most important hotspot for amino acid alterations followed by the ORF1a/ORF1ab poly-proteins. It is evident that the D614G mutation in spike glycoprotein and P4715L in RdRp is the important determinant of SARS-CoV-2 evolution since its emergence. P4715L in RdRp, G251V in ORF3a and S1498F of Nsp3 is associated with the epitope loss that may influence pathogenesis caused by antibody escape variants. The phylogenomics distinguished the ancestral viral samples from China and most part of Asia, isolated since the initial outbreak and the later evolved variants isolated from Europe and Americas. The evolved variants have been found to predominant globally with the loss of epitopes from its proteins. These have implications for SARS-CoV-2 transmission, pathogenesis and immune interventions.     

Population Dynamics and Structural Effects at Short and Long Range Support the Hypothesis of the Selective Advantage of the G614 SARS-CoV-2 Spike Variant

SARS-CoV-2 epidemics quickly propagated worldwide, sorting virus genomic variants in newly established propagules of infections. Stochasticity in transmission within and between countries or an actual selective advantage could explain the global high frequency reached by some genomic variants. Using statistical analyses, demographic reconstructions, and molecular dynamics simulations, we show that the globally invasive G614 spike variant 1) underwent a significant demographic expansion in most countries explained neither by stochastic effects nor by overrepresentation in clinical samples, 2) increases the spike S1/S2 furin-like site conformational plasticity (short-range effect), and 3) modifies the internal motion of the receptor-binding domain affecting its cross-connection with other functional domains (long-range effect). Our results support the hypothesis of a selective advantage at the basis of the spread of the G614 variant, which we suggest may be due to structural modification of the spike protein at the S1/S2 proteolytic site, and provide structural information to guide the design of variant-specific drugs.     

Increased pathogenicity and aerosol transmission for one SARS-CoV-2 B.1.617.2 Delta variant over the wild-type strain in hamsters

During the two-year pandemic of coronavirus disease 2019 (COVID-19), its causative agent, severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), has been evolving. SARS-CoV-2 Delta, a variant of concern, has become the dominant circulating strain worldwide within just a few months. Here, we performed a comprehensive analysis of a new B.1.617.2 Delta strain (Delta630) compared with the early WIV04 strain (WIV04) in vitro and in vivo, in terms of replication, infectivity, pathogenicity, and transmission in hamsters. When inoculated intranasally, Delta630 led to more pronounced weight loss and more severe disease in hamsters. Moreover, 40% mortality occurred about one week after infection with 10⁴ PFU of Delta630, whereas no deaths occurred even after infection with 10⁵ PFU of WIV04 or other strains belonging to the Delta variant. Moreover, Delta630 outgrew over WIV04 in the competitive aerosol transmission experiment. Taken together, the Delta630 strain showed increased replication ability, pathogenicity, and transmissibility over WIV04 in hamsters. To our knowledge, this is the first SARS-CoV-2 strain that causes death in a hamster model, which could be an asset for the efficacy evaluation of vaccines and antivirals against infections of SARS-CoV-2 Delta strains. The underlying molecular mechanisms of increased virulence and transmission await further analysis.     

Reduced neutralisation of the Delta (B.1.617.2) SARS-CoV-2 variant of concern following vaccination

Vaccines are proving to be highly effective in controlling hospitalisation and deaths associated with SARS-CoV-2 infection but the emergence of viral variants with novel antigenic profiles threatens to diminish their efficacy. Assessment of the ability of sera from vaccine recipients to neutralise SARS-CoV-2 variants will inform the success of strategies for minimising COVID19 cases and the design of effective antigenic formulations. Here, we examine the sensitivity of variants of concern (VOCs) representative of the B.1.617.1 and B.1.617.2 (first associated with infections in India) and B.1.351 (first associated with infection in South Africa) lineages of SARS-CoV-2 to neutralisation by sera from individuals vaccinated with the BNT162b2 (Pfizer/BioNTech) and ChAdOx1 (Oxford/AstraZeneca) vaccines. Across all vaccinated individuals, the spike glycoproteins from B.1.617.1 and B.1.617.2 conferred reductions in neutralisation of 4.31 and 5.11-fold respectively. The reduction seen with the B.1.617.2 lineage approached that conferred by the glycoprotein from B.1.351 (South African) variant (6.29-fold reduction) that is known to be associated with reduced vaccine efficacy. Neutralising antibody titres elicited by vaccination with two doses of BNT162b2 were significantly higher than those elicited by vaccination with two doses of ChAdOx1. Fold decreases in the magnitude of neutralisation titre following two doses of BNT162b2, conferred reductions in titre of 7.77, 11.30 and 9.56-fold respectively to B.1.617.1, B.1.617.2 and B.1.351 pseudoviruses, the reduction in neutralisation of the delta variant B.1.617.2 surpassing that of B.1.351. Fold changes in those vaccinated with two doses of ChAdOx1 were 0.69, 4.01 and 1.48 respectively. The accumulation of mutations in these VOCs, and others, demonstrate the quantifiable risk of antigenic drift and subsequent reduction in vaccine efficacy. Accordingly, booster vaccines based on updated variants are likely to be required over time to prevent productive infection. This study also suggests that two dose regimes of vaccine are required for maximal BNT162b2 and ChAdOx1-induced immunity.     

A variant-proof SARS-CoV-2 vaccine targeting HR1 domain in S2 subunit of spike protein

The emerging SARS-CoV-2 variants, commonly with many mutations in S1 subunit of spike (S) protein are weakening the efficacy of the current vaccines and antibody therapeutics. This calls for the variant-proof SARS-CoV-2 vaccines targeting the more conserved regions in S protein. Here, we designed a recombinant subunit vaccine, HR121, targeting the conserved HR1 domain in S2 subunit of S protein. HR121 consisting of HR1–linker1–HR2–linker2–HR1, is conformationally and functionally analogous to the HR1 domain present in the fusion intermediate conformation of S2 subunit. Immunization with HR121 in rabbits and rhesus macaques elicited highly potent cross-neutralizing antibodies against SARS-CoV-2 and its variants, particularly Omicron sublineages. Vaccination with HR121 achieved near-full protections against prototype SARS-CoV-2 infection in hACE2 transgenic mice, Syrian golden hamsters and rhesus macaques, and effective protection against Omicron BA.2 infection in Syrian golden hamsters. This study demonstrates that HR121 is a promising candidate of variant-proof SARS-CoV-2 vaccine with a novel conserved target in the S2 subunit for application against current and future SARS-CoV-2 variants.Subject terms: Mechanisms of disease, Molecular modelling     

Discovery of potential anti-SARS-CoV-2 drugs based on large-scale screening in vitro and effect evaluation in vivo

The coronavirus disease 2019 (COVID-19) pandemic caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is a global crisis. Clinical candidates with high efficacy, ready availability, and that do not develop resistance are in urgent need. Despite that screening to repurpose clinically approved drugs has provided a variety of hits shown to be effective against SARS-CoV-2 infection in cell culture, there are few confirmed antiviral candidates in vivo. In this study, 94 compounds showing high antiviral activity against SARS-CoV-2 in Vero E6 cells were identified from 2,580 FDA-approved small-molecule drugs. Among them, 24 compounds with low cytotoxicity were selected, and of these, 17 compounds also effectively suppressed SARS-CoV-2 infection in HeLa cells transduced with human ACE2. Six compounds disturb multiple processes of the SARS-CoV-2 life cycle. Their prophylactic efficacies were determined in vivo using Syrian hamsters challenged with SARS-CoV-2 infection. Seven compounds reduced weight loss and promoted weight regain of hamsters infected not only with the original strain but also the D614G variant. Except for cisatracurium, six compounds reduced hamster pulmonary viral load, and IL-6 and TNF-α mRNA when assayed at 4 d postinfection. In particular, sertraline, salinomycin, and gilteritinib showed similar protective effects as remdesivir in vivo and did not induce antiviral drug resistance after 10 serial passages of SARS-CoV-2 in vitro, suggesting promising application for COVID-19 treatment.

     

In-silico genomic landscape characterization and evolution of SARS-CoV-2 variants isolated in India shows significant drift with high frequency of mutations

In-silico studies on SARS-CoV-2 genome are considered important to identify the significant pattern of variations and its possible effects on the structural and functional characteristics of the virus. The current study determined such genetic variations and their possible impact among SARS-CoV-2 variants isolated in India. A total of 546 SARS-CoV-2 genomic sequences (India) were retrieved from the gene bank (NCBI) and subjected to alignment against the Wuhan variant (NC_045512.2), the corresponding amino acid changes were analyzed using NCBI Protein-BLAST. These 546 variants revealed 841 mutations; most of these were non-synonymous 464/841 (55.1%), there was no identical variant compared to the original strain. All genes; coding and non-coding showed nucleotide changes, most of the structural genes showed frequent nonsynonymous mutations. The most affected genes were ORF1a/b followed by the S gene which showed 515/841 (61.2%) and 120/841 (14.3%) mutations, respectively. The most frequent non-synonymous mutation 486/546 (89.01%) occurred in the S gene (structural gene) at position 23,403 where A changed to G leading to the replacement of aspartic acid by glycine in position (D614G). Interestingly, four variants also showed deletion. The variants MT800923 and MT800925 showed 12 consecutive nucleotide deletion in position 21982–21993 resulting in 4 consecutive amino acid deletions that were leucine, glycine, valine, and tyrosine in positions 141, 142, 143, and 144 respectively. The present study exhibited a higher mutations rate per variant compared to other studies carried out in India.     

Vaccine-escape and fast-growing mutations in the United Kingdom,the United States,Singapore, Spain,India, and other COVID-19-devastated countries

Recently, the SARS-CoV-2 variants from the United Kingdom (UK), South Africa, and Brazil have received much attention for their increased infectivity, potentially high virulence, and possible threats to existing vaccines and antibody therapies. The question remains if there are other more infectious variants transmitted around the world. We carry out a large-scale study of 506,768 SARS-CoV-2 genome isolates from patients to identify many other rapidly growing mutations on the spike (S) protein receptor-binding domain (RBD). We reveal that essentially all 100 most observed mutations strengthen the binding between the RBD and the host angiotensin-converting enzyme 2 (ACE2), indicating the virus evolves toward more infectious variants. In particular, we discover new fast-growing RBD mutations N439K, S477N, S477R, and N501T that also enhance the RBD and ACE2 binding. We further unveil that mutation N501Y involved in United Kingdom (UK), South Africa, and Brazil variants may moderately weaken the binding between the RBD and many known antibodies, while mutations E484K and K417N found in South Africa and Brazilian variants, L452R and E484Q found in India variants, can potentially disrupt the binding between the RBD and many known antibodies. Among these RBD mutations, L452R is also now known as part of the California variant B.1.427. Finally, we hypothesize that RBD mutations that can simultaneously make SARS-CoV-2 more infectious and disrupt the existing antibodies, called vaccine escape mutations, will pose an imminent threat to the current crop of vaccines. A list of most likely vaccine escape mutations is given, including S494P, Q493L, K417N, F490S, F486L, R403K, E484K, L452R, K417T, F490L, E484Q, and A475S. Mutation T478K appears to make the Mexico variant B.1.1.222 the most infectious one. Our comprehensive genetic analysis and protein-protein binding study show that the genetic evolution of SARS-CoV-2 on the RBD, which may be regulated by host gene editing, viral proofreading, random genetic drift, and natural selection, gives rise to more infectious variants that will potentially compromise existing vaccines and antibody therapies.     

Nucleocapsid mutations in SARS-CoV-2 augment replication and pathogenesis

While SARS-CoV-2 continues to adapt for human infection and transmission, genetic variation outside of the spike gene remains largely unexplored. This study investigates a highly variable region at residues 203–205 in the SARS-CoV-2 nucleocapsid protein. Recreating a mutation found in the alpha and omicron variants in an early pandemic (WA-1) background, we find that the R203K+G204R mutation is sufficient to enhance replication, fitness, and pathogenesis of SARS-CoV-2. The R203K+G204R mutant corresponds with increased viral RNA and protein both in vitro and in vivo. Importantly, the R203K+G204R mutation increases nucleocapsid phosphorylation and confers resistance to inhibition of the GSK-3 kinase, providing a molecular basis for increased virus replication. Notably, analogous alanine substitutions at positions 203+204 also increase SARS-CoV-2 replication and augment phosphorylation, suggesting that infection is enhanced through ablation of the ancestral ‘RG’ motif. Overall, these results demonstrate that variant mutations outside spike are key components in SARS-CoV-2’s continued adaptation to human infection.     

Characteristics of replication and pathogenicity of SARS-CoV-2 Alpha and Delta isolates

The continuously arising of SARS-CoV-2 variants has been posting a great threat to public health safety globally, from B.1.17 (Alpha), B.1.351 (Beta), P.1 (Gamma), B.1.617.2 (Delta) to B.1.1.529 (Omicron). The emerging or reemerging of the SARS-CoV-2 variants of concern is calling for the constant monitoring of their epidemics, pathogenicity and immune escape. In this study, we aimed to characterize replication and pathogenicity of the Alpha and Delta variant strains isolated from patients infected in Laos. The amino acid mutations within the spike fragment of the isolates were determined via sequencing. The more efficient replication of the Alpha and Delta isolates was documented than the prototyped SARS-CoV-2 in Calu-3 and Caco-2 cells, while such features were not observed in Huh-7, Vero E6 and HPA-3 cells. We utilized both animal models of human ACE2 (hACE2) transgenic mice and hamsters to evaluate the pathogenesis of the isolates. The Alpha and Delta can replicate well in multiple organs and cause moderate to severe lung pathology in these animals. In conclusion, the spike protein of the isolated Alpha and Delta variant strains was characterized, and the replication and pathogenicity of the strains in the cells and animal models were also evaluated.