Characteristics of SARS-CoV-2 Omicron BA.5 variants in Shanghai after ending the zero-COVID policy in December 2022: a clinical and genomic analysis

Introduction: An unprecedented surge of Omicron infections appeared nationwide in China in December 2022 after the adjustment of the COVID-19 response policy. Here, we report the clinical and genomic characteristics of SARS-CoV-2 infections among children in Shanghai during this outbreak.

Methods: A total of 64 children with symptomatic COVID-19 were enrolled. SARS-CoV-2 whole genome sequences were obtained using next-generation sequencing (NGS) technology. Patient demographics and clinical characteristics were compared between variants. Phylogenetic tree, mutation spectrum, and the impact of unique mutations on SARS-CoV-2 proteins were analysed in silico.

Results: The genomic monitoring revealed that the emerging BA.5.2.48 and BF.7.14 were the dominant variants. The BA.5.2.48 infections were more frequently observed to experience vomiting/diarrhea and less frequently present cough compared to the BF.7.14 infections among patients without comorbidities in the study. The high-frequency unique non-synonymous mutations were present in BA.5.2.48 (N:Q241K) and BF.7.14 (nsp2:V94L, nsp12:L247F, S:C1243F, ORF7a:H47Y) with respect to their parental lineages. Of these mutations, S:C1243F, nsp12:L247F, and ORF7a:H47Y protein were predicted to have a deleterious effect on the protein function. Besides, nsp2:V94L and nsp12:L247F were predicted to destabilize the proteins.

Discussion: Further in vitro to in vivo studies are needed to verify the role of these specific mutations in viral fitness. In addition, continuous genomic monitoring and clinical manifestation assessments of the emerging variants will still be crucial for the effective responses to the ongoing COVID-19 pandemic.

1 Introduction

The coronavirus disease 2019 (COVID-19) pandemic caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) resulted in a global emergence during the past three years since cases were first reported. The emergence of new variants will impose a risk of future surges (Wang et al., 2020; Zhu et al., 2020). Given the widespread and continuous evolution of SARS-CoV-2, numerous variants of concern (VOCs) have emerged and successively dominated multiple waves of the COVID-19 pandemic globally (Boehm et al., 2021). The new VOCs are often associated with increased transmissibility and/or immune evasion properties, which led to their rapid spread globally (Karim and Karim, 2022). Currently, the Omicron variant (B.1.1.529) is the predominant VOC around the world, since its emergence in South Africa in November 2021 (Petersen et al., 2022). From a genomic perspective, it shares several mutations with the previously identified VOCs, such as Alpha (B.1.1.7), Beta (B.1.351), Gamma (P.1), and Delta (B.1.617.2), but it also harbors a large number of specific mutations (Mohapatra et al., 2023). Up to now, a series of Omicron sub-lineages including BA.1 (original Omicron), BA.2, BA.3, BA.4, BA.5, and XBB have emerged and then caused the waves of COVID-19 globally due to further neutralization escape (Ai et al., 2022; Khan et al., 2022; Planas et al., 2022; Qu et al., 2022; Mohapatra et al., 2023). This highlights the importance of continuous genomic monitoring of SARS-CoV-2 variants.

Since the outbreak of COVID-19 in late 2019, China has adhered to policies of zero-COVID for almost three years with strictly enforced lockdowns and other restrictive measures, including social distancing, school closure, mask use, and case isolation (Lai et al., 2020; Liu et al., 2022). Given the attenuated pathogenicity of omicron subvariants and increasing vaccination coverage, China lifted the zero-COVID strategies, notably by announcing the ‘10 measures’ about the optimization of COVID-19 rules on 7 December 2022 (Xinhua, 2022). After that, China experiences a nationwide outbreak of COVID-19. Leung et al. (2023) estimated that the cumulative infection attack rate in Beijing was 75.7% (95% credible interval (CrI): 60.7–84.4) on 22 December 2022 and 92.3% (95% CrI: 91.4–93.1) on 31 January 2023. A recent study by Liang et al. (2023) showed that the cumulative SARS-CoV-2 infection rate rose rapidly to 70% within three weeks after the ending of the zero-COVID policy in Macao. A study conducted in Guangzhou also revealed that the infection attack ratio reached to 80.7% (95% CrI: 72.2–86.8) at 30 days after easing the zero-COVID policy (Huang et al., 2023). Such an unprecedented epidemic raised concerns about specific and real-time data on the viral genetic sequencing, monitoring of variants, and disease impact (World Health Organization, 2022).

Shanghai, with a population of 25 million, is a leading economic center in China. The prevalence of SARS-CoV-2 variants in Shanghai can be considered a snapshot of China. Herein, we report the clinical and genetic characteristics of SARS-CoV-2 infections among children with COVID-19 in Shanghai after ending the zero-COVID policy in December 2022, based on viral genetic sequencing and clinical data.

2 Materials and methods

2.1 Study population and data collection

This study randomly selected and enrolled 64 pediatric cases with symptomatic COVID-19, who were admitted to the Children’s Hospital of Fudan University in late December 2022. Clinical data were collected via electronic medical charts, including demographic information, clinical symptoms, laboratory findings, and outcomes.

2.2 Sample selection and sequencing

Nasopharyngeal swabs obtained from the enrolled cases and confirmed as SARS-COV-2 positive by real-time PCR with cycle threshold (Ct) < 30 were selected for genome sequencing. Viral RNA from the swabs was extracted using an automatic magnetic extraction device and accompanying kit (Daan Gene Co., Ltd) following the manufacturer’s instructions. The SARS-CoV-2 amplicon libraries were generated with a 15 μL viral RNA template, by using the VAHTS RNA Multi-PCR Library Prep Kit according to the manufacturer’s protocol. Libraries were then sequenced on a Nova Seq instrument (Illumina, San Diego, CA, United States) with 2 × 150-bp paired-endreads. Raw reads were trimmed for adapters and filtered for quality (average q20 threshold and read length > 50 nt) using Trimmomatic (version 0.39). The last 8 nucleotides were also removed from all reads. Reference-based assembly was performed with Bowtie2 (version 2.3.5), aligning against the GenBank reference genome MN908947.3. SNPs variants were called through a pipeline based on GATK (version 4.0.6.0), and all SNPs having a minimum supporting read frequency of 50% were retained.

2.3 Mutation identification and phylogenetic analysis

Lineages of the SARS-CoV-2 consensus sequences obtained were assigned using Pangolin (version 4.2) (O'Toole et al., 2021). Mutations were identified using Nextclade (version 2.11.0) (Hadfield et al., 2018). All sequences were aligned using the MAFFT (version 7.511) (Katoh et al., 2002) and a Maximum likelihood (ML) phylogenetic tree was reconstructed by IQ-TREE2 (version 2.1.2, COVID-edition) (Minh et al., 2020) using the best-fit model of nucleotide substitution TIM + F + I + I + R2 inferred by ModelFinder (Kalyaanamoorthy et al., 2017) and rooted using Wuhan Hu-1 (MN908947.3) as an outgroup. The phylogenetic tree was visualized and modified with the FigTree (version 1.4.4). The illustrated figures of protein structures with the high-frequency (>50%) unique mutations were prepared using UCSF Chimera (version 1.16) (Pettersen et al., 2004).

2.4 Analysis of mutation impact on viral protein function

Preliminary functional analysis of the high-frequency unique non-synonymous mutations in proteins was performed using PredictSNP¹ (Bendl et al., 2014). Wuhan Hu-1 (MN908947.3) was selected as canonical protein sequence for the analysis. PredictSNP comprises scores from different predictors (MAPP, PhD-SNP, PolyPhen-1, PolyPhen-2, SIFT, SNAP) and uses the information of them to create its own score. PredictSNP then classifies mutations as “neutral” or “deleterious” and transforms the individual confidence scores of each predictor into one comparable scale ranging from 0 to 100%, which represents the percentage of expected accuracy.

2.5 Analysis of mutation impact on viral protein stability

The I-Mutant3.0² and DynaMut2³ web servers were used for predicting SARS-COV-2 protein conformational stability changes upon the identified high-frequency unique mutations by default pH and temperature. I-Mutant3.0 is a suite of Support Vector Machine (SVM) based predictors and offers the opportunity to predict automatically protein stability changes upon single-site mutations starting from protein sequence alone (Capriotti et al., 2005). DynaMut2 is a structure-based approach for assessing mutation effects on protein stability by using normal mode analysis (NMA) approaches with graph-based distance matrix (Rodrigues et al., 2018). In this study, the wide-type 3D structures of the nucleocapsid (N) protein (PDB ID: 8FD5), spike (S) glycoprotein (PDB ID: 6VXX), nsp2 (PDB ID: 7MSW), nsp12 (PDB ID: 7C2K) and ORF7a protein (PDB ID: 7CI3) were retrieved from Protein Data Bank (PDB). Both tools classify each mutation as stabilizing or destabilizing by providing the predicted Gibbs free energy change (ΔΔG). A positive ΔΔG value corresponds to the mutation predicted to be stabilizing, and a negative value suggests that the mutation can destabilize the protein.

2.6 Statistical analysis

Categorical variables were expressed as numbers (%) and compared using Pearson chi-squared or Fisher’s exact tests. Continuous variables were expressed as medians (interquartile range) and compared using Mann–Whitney U tests. All of the tests were two-tailed, and a p value <0.05 represented statistical significance. The statistical analyses were conducted in SPSS version 26.0 software (IBM, New York, United States).

3 Results

3.1 Demographic and clinical characteristics

A total of 64 children with COVID-19 were included in this study. The median age of the patients was 17 months (IQR: 3–43 months), and 64.1% were male. Fever was the most common symptom, observed in 98.4% of the patients, followed by cough (64.1%), shortness of breath (23.4%), and vomiting/diarrhea (21.9%). Abnormal chest imaging was observed in 54.7% of the patients. The median duration of hospitalization was 3 days (IQR: 1.8–4.4 months). Among those, 15 patients (23.4%) had at least one comorbidity. The most common comorbidities were chronic neurological disorders (n = 9), followed by chronic haematological disorders (n = 4), chronic hepatic disorders (n = 2), and solid tumor (n = 1). Children with comorbidities were more likely to present shortness of breath, abnormal chest imaging, PICU admission, and longer duration of hospitalization than those without comorbidities (Table 1).

Table 1

Table 1. Clinical characteristics and laboratory findings of children with SARS-CoV-2 infection.

3.2 Genotype and phylogenetic analysis

SARS-CoV-2 genome assemblies were obtained from the enrolled patients. Our genomic monitoring indicated that all cases clustered into the Omicron BA.5.2* lineage, with the dominant lineages of BA.5.2.48 (33/64, 51.6%) and BF.7.14 (26/64, 40.6%). Besides, sporadic detection of other sub-lineages such as BA.5.2.49 (3/64, 4.7%) was also observed (Figure 1A). As expected, the ML tree based on the complete genomes showed that all sequences could be classified into two main clades and clustered tightly with BA.5 lineage (Figure 1B). Thirty-seven sequences belonged to clade 1, including 33 BA.5.2.48, 3 BA.5.2.49, and 1 BA.5.2. The other 27 sequences belonged to clade 2, including 26 BF.7.14 and 1 BF.7 (Figure 1B).

Figure 1

Figure 1. Genetic and clinical characterization of the Omicron variant driving the wave of SARS-CoV-2 outbreak in Shanghai after ending the zero-COVID policy in December 2022. (A) The composition of SARS-CoV-2 sub-lineages in this study. BA.5.2.48 and BF.7.14 were the dominant sub-lineages. (B) Maximum likelihood (ML) tree of 64 genomes sequenced in this study. The tree was rooted with Wuhan Hu-1 (MN908947.3). (C) Comparison of clinical characteristics of infections with BF.7.14 and BA.5.2.48 among children without comorbidities. Statistical evaluations were made with Pearson chi-squared or Fisher’s exact tests. ns, not significant; * significant level at p value <0.05.

Further, we compared the demographic and clinical characteristics of patients with BA.5.2.48 and BF.7.14 infections. To avoid biases caused by comorbidities, children with and without comorbidities were analyzed individually. The results showed that the BA.5.2.48 infections were more frequently observed to experience vomiting/diarrhea and less frequently present cough compared to the BF.7.14 infections among patients without comorbidities (Table 2) (Figure 1C). Further, to rule out the possibility of other confounding factors, common viruses including norovirus, adenovirus, and rotavirus were detected in the stool samples from the patients with vomiting/diarrhea. The results were all negative for these common viruses. However, there were no significant differences in clinical characteristics and laboratory findings between BA.5.2.48 and BF.7.14 infections among children with comorbidities (Table 3).

Table 2

Table 2. Clinical characteristics and laboratory findings of children without comorbidities according to the infection with different Omicron lineages.

Table 3

Table 3. Clinical characteristics and laboratory findings of children with comorbidities according to the infection with different Omicron lineages.

3.3 Mutational analysis

Compared with the reference genome of Wuhan-Hu-1 strain, a total of 194 mutations and deletions were spotted across different genome regions; however, only 111 variant sites with a prevalence of ≥2 sequences were presented. Among the 111 variant sites, there were 70 (63.1%) non-synonymous mutations, 30 (27.0%) synonymous mutations, 5 (4.5%) deletions, and 6 (5.4%) mutations in the untranslated region (UTR) (Table 4). The synonymous mutations were mainly located in non-structural protein (nsp) regions (Table 4).

Table 4

Table 4. Mutations and deletions in 64 sequences of SARS-CoV-2 isolated in this study.

Most of these variant sites were in common with the known mutations in their parental lineages (BA.5.2 for BA.5.2.48, and BF.7 for BF.7.14). However, we found 5 unique non-synonymous mutations in BA.5.2.48 and BF.7.14, respectively, with high frequency (>50%) (Table 5). Among these unique mutations, N:Q241K is a characteristic mutation of BA.5.2.48 linage and nsp2:V94L, nsp12:L247F, S:C1243F, and ORF7a:H47Y are the characteristic mutations of BF.7.14 linage. Figure 2 is a graphical representation that shows the location of these unique mutations in each region of the complete SARS-CoV-2 genome and the high frequency unique mutations in the protein crystal structures.

Table 5

Table 5. Unique non-synonymous mutations in BA.5.2.48 and BF.7.14 compared with their parental lineages.

Figure 2

Figure 2. Schematic representation of the unique amino acid mutations of the variants in this study. The schematic diagram of BA.5.2.48 (A) and BF.7.14 (B) unique amino acid mutation sites on the SARS-CoV-2 genome. The high frequency unique mutations are shown in red. (C) The location of the high-frequency unique mutations in the protein crystal structures. V94L mutation is located in the nsp2 N-terminal. L247F mutation is located in the NiRAN domain, which lies at the N terminal end of the RdRp (nsp12) domain. H47Y mutation is located in the ectodomain of the ORF7a protein. Q241K mutation is located in an intrinsically disordered region of the N protein, which connects the N-terminal domain and the C-terminal domain. Mutation positions are framed in red circles, and mutant residues are represented in stick form. C1243F mutation is located at the cytoplasmic region of the spike protein without a resolution crystal structure.

Further, we assessed the associations of the mutations with the symptoms of vomiting/diarrhea and cough. Interestingly, we found that the frequency of nsp12:L247F, S:C1243F, and ORF7a:H47Y (characteristic mutations of BF.7.14) was significantly higher among those with cough than without (Figure 3A). Besides, the frequency of these three mutations was lower among those with vomiting/diarrhea than without (Figure 3B).

Figure 3

Figure 3. The associations of the mutations with the symptoms of cough (A) and vomiting/diarrhea (B). The mutation frequencies between groups were compared using Pearson chi-squared or Fisher’s exact tests. Only mutations with statistically significant difference in frequencies between groups are presented. Only children without comorbidities were included. Color gradient indicates mutation frequencies. Synonymous mutations are colored in green.

3.4 Effect of mutations on protein function

The predicted effects of pathogenicity for the high-frequency mutations are show in Figure 4A. A mutation was classified as deleterious only when it was predicted as deleterious by more than three tools. Our results revealed that the Q241K mutation in N protein and the V94L mutation in nsp2 were predicted to have a neutral effect on the protein function. However, L247F mutation in nsp12, C1243F mutation in S protein, and H47Y mutation in ORF7a protein were predicted to be deleterious by the consensus classifier.

Figure 4

Figure 4. In silico analysis of protein function and structural stability changes upon the high-frequency unique mutations using different tools. (A) Effect of the mutations on protein function. The red color stands for deleterious mutation, whereas the green color represents neutral mutation. (B) Effect of the mutations on protein structural stability. Site 1,243 without resolution in the spike (S) protein crystallographic structure, located at the cytoplasmic region was not analyzed for energy estimation by the structure-based tool Dynamut2. NA, not available.

3.5 Effect of mutations on protein stability

Structural stability of proteins due to the high-frequency mutations were analyzed using I-Mutant and DynaMut stability predictors. The tools provided almost consensus predicted results. The V94L mutation in nsp2, as well as the L247F mutation in nsp12, destabilized the proteins. While the H47Y mutation increased the stability of structure of the ORF7a protein. The Q241K mutation in the N protein and C1243F mutation in the S protein had little effect on the protein stability due to the low absolute ΔΔG values (Figure 4B).

4 Discussion

In this study, we have reported the first comprehensive clinical and viral genomic analysis of SARS-CoV-2 infections among children hospitalized with COVID-19 in Shanghai after ending the zero-COVID policy. Whole genome sequencing of samples obtained from the enrolled 64 pediatric patients revealed that all cases clustered into the Omicron BA.5.2* lineage, with the dominant omicron sub-lineages of BA.5.2.48 and BF.7.14. This result was in line with our previous study based on 8,254 SARS-CoV-2 complete genomes available on the GISAID database from the Chinese mainland during December 2022 and January 2023, which indicated that the genomes corresponded to 88 Pango-nomenclature-system-named subvariants, with the dominant lineages of BA.5.2.48 (4,881/8254, 59.1%) and BF.7.14 (2,223/8254, 26.9%), and the proportion of these dominant lineages were not significant changed over the two months of outbreak (Liu and Xu, 2023). In addition, BA.5.2.48 and BF.7.14 were also found to be the dominant lineages among SARS-CoV-2 positive passengers on flights from China to Italy in late December 2022 (Novazzi et al., 2023). However, the lineages found to be dominant internationally during the same period, such as BQ.1, BQ.1.1, and XBB.1.5 were quickly cleared and did not prevail in China (Liu and Xu, 2023). Taken together, all these results demonstrated that the emerging BA.5.2.48 and BF.7.14 were the absolutely dominant drivers of the current COVID-19 outbreak after ending the zero-COVID policy, which could be attributed to the high fitness of lineages or a random founder effect in China.

Fever and cough were the most common symptoms among children with COVID-19 in this study. This result was consistent with the earlier community outbreak in Shanghai driven by BA.2.2.1 sub-lineage in spring 2022 (Ao et al., 2022; Ling et al., 2022; Shen et al., 2023). However, the demographic characteristics, clinical symptoms, laboratory findings, and outcomes might vary between different SARS-CoV-2 variants in children hospitalized with COVID-19 (Boncuoglu et al., 2022; Quintero et al., 2022; Tagarro et al., 2022; Sahin et al., 2023). Consequently, we further compared the clinical features of children with BA.5.2.48 and BF.7.14 infections. We found that the BA5.2.48 infections were more frequently observed to experience vomiting/diarrhea and less frequently present cough compared to the BF.7.14 infections among patients without comorbidities. To figure out susceptible mutations related to the variation of symptoms, we assessed the associations of the mutations with the symptoms of vomiting/diarrhea and cough. We found that the frequency of the characteristic mutation nsp12:L247F, S:C1243F, and ORF7a:H47Y was significantly varied among those with cough or vomiting/diarrhea than without. These observations suggest that these three characteristic mutations might contribute to the variation of symptoms between children infected with BF.7.14 and BA.5.2.48 by affecting the tissue tropism of the variants or other mechanisms. However, possibly arose bias due to the limited sample size in this study. Besides, the results might also be affected by the demographic characteristics, especially age. Thus, the interpretation of the results must be cautious and might not precisely reflect the broader population. Given the rapid evolution of SARS-CoV-2, further studies are still needed to elucidate the clinical features and severity of different variants with specific mutations.

Adaptive mutations in the SARS-CoV-2 genome could alter its pathogenic potential, and at the same time would increase the infectivity and immune escape capacity. Single amino acid changes are worth monitoring because they can be phenotypically relevant. Perhaps one of the best exemplars of the impacts of amino acid changes in the SARS-CoV-2 is the D614G mutation in S protein. D614G substitution was first identified in early 2020 and rapidly spread throughout the global population by increasing the infectivity and stability of virion (Korber et al., 2020; Zhang et al., 2020; Plante et al., 2021). With this background, we further investigated the unique amino acid changes in the BA.5.2.48 and BF.7.14 sub-lineages, and predicted the effects of these mutations on the stability and function of viral proteins. Stability is a parameter which is crucial to judge the functional and structural activity of a protein. Protein stability dictates the conformational structure of the protein, thereby determining its function. Any change in protein stability may cause misfolding, degradation or aberrant conglomeration of proteins. Understanding the stability changes in SARS-CoV-2 proteins is essential for predicting virus infectivity. Changes in Gibbs free energy of unfolding (ΔΔG) between the wild-type and mutant proteins could predict the effects of mutations on the stability of protein structure (Pan et al., 2022). In order to ascertain the significance of mutations on protein function, we analyzed the pathogenicity of the mutations as deleterious or neutral. It is important to note that in case of protein, damaging mostly defines instability. Generally, this is used for human proteins. As a consequence, if the human protein is damaging in nature because of mutations, then the human protein–protein interactions may occur with high or low binding affinity. Now in case of virus, similar consequences may happen, which means if the virus protein is damaged because of mutations, it may interact with human proteins with similar binding affinity. As a result, the virus may acquire characteristics like transmissibility, escaping antibodies. For example, the D614G was predicted to be deleterious and instable using I-mutant and PredictSNP servers. Thus, the basic premise for the study was that mutations will be contributing to the viral evolution only if they are deleterious and neutral mutations would not be affecting the protein function (Laskar and Ali, 2021).

We found 5 high-frequency amino acid changes (N: Q241K, nsp2: V94L, nsp12: L247F, S: C1243F, ORF7a: H47Y) in the BA.5.2.48 and BF.7.14 sub-lineages. Among these mutations, the C1243F mutation in S protein, L247F mutation in nsp12, and H47Y mutation in ORF7a protein were predicted to have a deleterious effect on the protein function. S protein decorates the surface of coronavirus and plays a critical role in viral entry (Gallagher and Buchmeier, 2001). It comprises two functional subunits responsible for binding to the host cell receptor (S1 subunit) and membrane fusion (S2 subunit) (Walls et al., 2020). In the S1 subunit, there is an N-terminal domain (14–305 residues) and a receptor-binding domain (RBD, 319–541 residues); the fusion peptide (FP) (788–806 residues), heptapeptide repeat sequence 1 (HR1) (912–984 residues), HR2 (1,163–1,213 residues), transmembrane domain (1,213–1,237 residues), and intracellular domain (1,237–1,273 residues) comprise the S2 subunit (Huang et al., 2020). The C1243F mutation is located in the intracellular domain of S2 subunit. The mutations in the intracellular domain are unlikely to drive immune evasion. However, the mutations in this domain may affect the S protein expression at the cell surface and syncytia formation by mediating intracellular trafficking and membrane location of S protein (Cattin-Ortola et al., 2021; Li et al., 2022).

nsp12, also named RNA-dependent RNA polymerase (RdRp), catalyzes the synthesis of viral RNA and thus plays a central role in the replication and transcription cycle, with the assistance of nsp7 and nsp8 as cofactors. The structure of the nsp12 contains a right-hand RdRp domain (367–920 residues) and a nidovirus RdRp-associated nucleotidyltransferase domain (NiRAN, 60–249 residues) (Gao et al., 2020). The L247F mutation is located in the NiRAN domain, which lies at the N terminal end of the RdRp. Although the NiRAN domain is essential for viral propagation, its functions during the viral life cycle remain unclear. A recent study revealed that the NiRAN domain catalyzes the covalent link of RNA 5′ end to the first residue of nsp9, thus being an intermediate to form cap core (GpppA) with GTP catalyzed again by NiRAN (Yan et al., 2022). Therefore, the L247F mutation in nsp12 may affect the SARS-CoV-2 replication in the host cells.

ORF7a protein is a type-I transmembrane protein, consisting of an N-terminal signaling region (1–15 residues), an immunoglobulin-like ectodomain (16–96 residues), a hydrophobic transmembrane domain (97–116 residues), and a typical endoplasmic reticulum retention motif (117–121 residues) (Zhou et al., 2021). The H47Y mutation is located in the ectodomain. A recent study suggested that the Immunoglobulin-like fold ectodomain of the ORF7a interacts with high efficiency to the CD14+ monocytes in human peripheral blood, and ORF7a may also suppress the antigen-presenting ability of these monocytes and trigger the significant upregulation of multiple proinflammatory cytokines (Zhou et al., 2021). Further in vitro and in vivo studies are needed to verify the role of ORF7a: H47Y mutation in viral fitness.

5 Conclusion

Our results revealed that the current large-scale COVID-19 outbreak in Shanghai after ending the zero-COVID policy was driven by the emerging BA.5.2.48 and BF.7.14 variants with unique deleterious mutations. In addition, this study described the clinical characteristics of pediatric cases infected with BA.5.2.48 and BF.7.14. Continuous genomic monitoring and clinical manifestation assessments of the emerging variants will be crucial for countering the ongoing COVID-19 pandemic.

Data availability statement

The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found at: https://gisaid.org, EPI_ISL_17371203 to EPI_ISL_17371266.

Ethics statement

The studies involving humans were approved by Children’s Hospital of Fudan University. The studies were conducted in accordance with the local legislation and institutional requirements. Written informed consent for participation was not required from the participants or the participants’ legal guardians/next of kin in accordance with the national legislation and institutional requirements.

Author contributions

PL: Conceptualization, Formal analysis, Methodology, Writing – original draft, Writing – review & editing. JC: Conceptualization, Investigation, Resources, Writing – review & editing. HT: Data curation, Investigation, Writing – review & editing. JL: Investigation, Validation, Writing – review & editing. LL: Resources, Validation, Writing – review & editing. MX: Investigation, Resources, Writing – review & editing. XZ: Investigation, Resources, Writing – review & editing. XF: Investigation, Resources, Writing – review & editing. XW: Methodology, Resources, Writing – review & editing. HZ: Investigation, Resources, Writing – review & editing. RJ: Investigation, Resources, Writing – review & editing. YG: Investigation, Resources, Writing – review & editing. YZ: Investigation, Resources, Writing – review & editing. MZ: Conceptualization, Project administration, Supervision, Writing – review & editing. JX: Conceptualization, Project administration, Supervision, Writing – review & editing.

Funding

The author(s) declare financial support was received for the research, authorship, and/or publication of this article. This study was supported by a 3-Year Action Plan for Strengthening Public Health System in Shanghai (No: GWVI-11.2-YQ43), a grant from Shanghai Municipal Health Commission (No: 202040099), and a Key Development Program of Children’s Hospital of Fudan University (No: EK2022ZX05).

Conflict of interest

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

The reviewer XZ declared a shared parent affiliation with the authors and a shared affiliation with the author JX to the handling editor at the time of review.

Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Footnotes

1. ^https://loschmidt.chemi.muni.cz/predictsnp

2. ^http://gpcr2.biocomp.unibo.it/cgi/predictors/I-Mutant3.0/I-Mutant3.0.cgi

3. ^http://biosig.unimelb.edu.au/dynamut

References

Ai, J., Wang, X., He, X., Zhao, X., Zhang, Y., Jiang, Y., et al. (2022). Antibody evasion of SARS-CoV-2 omicron BA.1, BA.1.1, BA.2, and BA.3 sub-lineages. Cell Host Microbe 30, 1077–1083.e4. doi: 10.1016/j.chom.2022.05.001

PubMed Abstract | Crossref Full Text | Google Scholar

Ao, Y., Li, J., Wei, Z., Wang, Z., Tian, H., Qiu, Y., et al. (2022). Clinical and virological characteristics of SARS-CoV-2 omicron BA.2.2 variant outbreaks during April to May, 2022, Shanghai, China. J. Infect. 85, 573–607. doi: 10.1016/j.jinf.2022.07.027

PubMed Abstract | Crossref Full Text | Google Scholar

Bendl, J., Stourac, J., Salanda, O., Pavelka, A., Wieben, E. D., Zendulka, J., et al. (2014). PredictSNP: robust and accurate consensus classifier for prediction of disease-related mutations. PLoS Comput. Biol. 10:e1003440. doi: 10.1371/journal.pcbi.1003440

PubMed Abstract | Crossref Full Text | Google Scholar

Boehm, E., Kronig, I., Neher, R. A., Eckerle, I., Vetter, P., and Kaiser, L. (2021). Novel SARS-CoV-2 variants: the pandemics within the pandemic. Clin. Microbiol. Infect. 27, 1109–1117. doi: 10.1016/j.cmi.2021.05.022

PubMed Abstract | Crossref Full Text | Google Scholar

Boncuoglu, E., Kiymet, E., Sahinkaya, S., Cem, E., Yilmaz Celebi, M., Gulderen, M., et al. (2022). Did hospitalization age decrease in children in the omicron (B.1.1.529) era? Pediatr. Infect. Dis. J. 41:e403. doi: 10.1097/INF.0000000000003600

PubMed Abstract | Crossref Full Text | Google Scholar

Capriotti, E., Fariselli, P., and Casadio, R. (2005). I-Mutant2.0: predicting stability changes upon mutation from the protein sequence or structure. Nucleic Acids Res. 33, W306–W310. doi: 10.1093/nar/gki375

PubMed Abstract | Crossref Full Text | Google Scholar

Cattin-Ortola, J., Welch, L. G., Maslen, S. L., Papa, G., James, L. C., and Munro, S. (2021). Sequences in the cytoplasmic tail of SARS-CoV-2 spike facilitate expression at the cell surface and syncytia formation. Nat. Commun. 12:5333. doi: 10.1038/s41467-021-25589-1

PubMed Abstract | Crossref Full Text | Google Scholar

Gallagher, T. M., and Buchmeier, M. J. (2001). Coronavirus spike proteins in viral entry and pathogenesis. Virology 279, 371–374. doi: 10.1006/viro.2000.0757

PubMed Abstract | Crossref Full Text | Google Scholar

Gao, Y., Yan, L., Huang, Y., Liu, F., Zhao, Y., Cao, L., et al. (2020). Structure of the RNA-dependent RNA polymerase from COVID-19 virus. Science 368, 779–782. doi: 10.1126/science.abb7498

PubMed Abstract | Crossref Full Text | Google Scholar

Hadfield, J., Megill, C., Bell, S. M., Huddleston, J., Potter, B., Callender, C., et al. (2018). Nextstrain: real-time tracking of pathogen evolution. Bioinformatics 34, 4121–4123. doi: 10.1093/bioinformatics/bty407

Crossref Full Text | Google Scholar

Huang, Y., Yang, C., Xu, X. F., Xu, W., and Liu, S. W. (2020). Structural and functional properties of SARS-CoV-2 spike protein: potential antivirus drug development for COVID-19. Acta Pharmacol. Sin. 41, 1141–1149. doi: 10.1038/s41401-020-0485-4

PubMed Abstract | Crossref Full Text | Google Scholar

Huang, J., Zhao, S., Chong, K. C., Zhou, Y., Lu, W., Fang, F., et al. (2023). Infection rate in Guangzhou after easing the zero-COVID policy: seroprevalence results to ORF8 antigen. Lancet Infect. Dis. 23, 403–404. doi: 10.1016/S1473-3099(23)00112-3

PubMed Abstract | Crossref Full Text | Google Scholar

Kalyaanamoorthy, S., Minh, B. Q., Wong, T. K. F., von Haeseler, A., and Jermiin, L. S. (2017). ModelFinder: fast model selection for accurate phylogenetic estimates. Nat. Methods 14, 587–589. doi: 10.1038/nmeth.4285

PubMed Abstract | Crossref Full Text | Google Scholar

Karim, S. S. A., and Karim, Q. A. (2022). Omicron SARS-CoV-2 variant: a new chapter in the COVID-19 pandemic. Lancet 399:142. doi: 10.1016/S0140-6736(21)02758-6

Crossref Full Text | Google Scholar

Katoh, K., Misawa, K., Kuma, K., and Miyata, T. (2002). MAFFT: a novel method for rapid multiple sequence alignment based on fast Fourier transform. Nucleic Acids Res. 30, 3059–3066. doi: 10.1093/nar/gkf436

PubMed Abstract | Crossref Full Text | Google Scholar

Khan, K., Karim, F., Ganga, Y., Bernstein, M., Jule, Z., Reedoy, K., et al. (2022). Omicron BA.4/BA.5 escape neutralizing immunity elicited by BA.1 infection. Nat. Commun. 13:4686. doi: 10.1038/s41467-022-32396-9

Crossref Full Text | Google Scholar

Korber, B., Fischer, W. M., Gnanakaran, S., Yoon, H., Theiler, J., Abfalterer, W., et al. (2020). Tracking changes in SARS-CoV-2 spike: evidence that D614G increases infectivity of the COVID-19 virus. Cell 182, 812–827.e19. doi: 10.1016/j.cell.2020.06.043

PubMed Abstract | Crossref Full Text | Google Scholar

Lai, S., Ruktanonchai, N. W., Zhou, L., Prosper, O., Luo, W., Floyd, J. R., et al. (2020). Effect of non-pharmaceutical interventions to contain COVID-19 in China. Nature 585, 410–413. doi: 10.1038/s41586-020-2293-x

PubMed Abstract | Crossref Full Text | Google Scholar

Laskar, R., and Ali, S. (2021). Differential mutation profile of SARS-CoV-2 proteins across deceased and asymptomatic patients. Chem. Biol. Interact. 347:109598. doi: 10.1016/j.cbi.2021.109598

PubMed Abstract | Crossref Full Text | Google Scholar

Leung, K., Lau, E. H. Y., Wong, C. K. H., Leung, G. M., and Wu, J. T. (2023). Estimating the transmission dynamics of SARS-CoV-2 omicron BF.7 in Beijing after adjustment of the zero-COVID policy in November-December 2022. Nat. Med. 29, 579–582. doi: 10.1038/s41591-023-02212-y

PubMed Abstract | Crossref Full Text | Google Scholar

Li, Q., Liu, Y., and Zhang, L. (2022). Cytoplasmic tail determines the membrane trafficking and localization of SARS-CoV-2 spike protein. Front. Mol. Biosci. 9:1004036. doi: 10.3389/fmolb.2022.1004036

PubMed Abstract | Crossref Full Text | Google Scholar

Liang, J., Liu, R., He, W., Zeng, Z., Wang, Y., Wang, B., et al. (2023). Infection rates of 70% of the population observed within 3 weeks after release of COVID-19 restrictions in Macao, China. J. Infect. 86, 402–404. doi: 10.1016/j.jinf.2023.01.029

PubMed Abstract | Crossref Full Text | Google Scholar

Ling, Y., Lu, G., Liu, F., Tan, Y., Xu, X., Wei, D., et al. (2022). The omicron BA.2.2.1 subvariant drove the wave of SARS-CoV-2 outbreak in Shanghai during spring 2022. Cell Discov. 8:97. doi: 10.1038/s41421-022-00468-1

PubMed Abstract | Crossref Full Text | Google Scholar

Liu, J., Liu, M., and Liang, W. (2022). The dynamic COVID-zero strategy in China. China CDC Wkly. 4, 74–75. doi: 10.46234/ccdcw2022.015

PubMed Abstract | Crossref Full Text | Google Scholar

Liu, P., and Xu, J. (2023). Genomic surveillance of SARS-CoV-2 in mainland China after ending the zero-COVID policy, December 2022-January 2023. J. Infect. 86, e84–e86. doi: 10.1016/j.jinf.2023.02.040

PubMed Abstract | Crossref Full Text | Google Scholar

Minh, B. Q., Schmidt, H. A., Chernomor, O., Schrempf, D., Woodhams, M. D., Haeseler, V., et al. (2020). IQ-TREE 2: new models and efficient methods for phylogenetic inference in the genomic era. Mol. Biol. Evol. 37, 1530–1534. doi: 10.1093/molbev/msaa015

PubMed Abstract | Crossref Full Text | Google Scholar

Mohapatra, R. K., Verma, S., Kandi, V., Sarangi, A. K., Seidel, V., Das, S. N., et al. (2023). The SARS-CoV-2 omicron variant and its multiple sub-lineages: transmissibility, vaccine development, antiviral drugs, monoclonal antibodies, and strategies for infection control – a review. ChemistrySelect 8:e20220138. doi: 10.1002/slct.202201380

Crossref Full Text | Google Scholar

Novazzi, F., Giombini, E., Rueca, M., Baj, A., Fabeni, L., Genoni, A., et al. (2023). Genomic surveillance of SARS-CoV-2 positive passengers on flights from China to Italy, December 2022. Euro. Surveill. 28:2300008. doi: 10.2807/1560-7917.ES.2023.28.2.2300008

Crossref Full Text | Google Scholar

O'Toole, A., Scher, E., Underwood, A., Jackson, B., Hill, V., McCrone, J. T., et al. (2021). Assignment of epidemiological lineages in an emerging pandemic using the pangolin tool. Virus Evol. 7:veab064. doi: 10.1093/ve/veab064

PubMed Abstract | Crossref Full Text | Google Scholar

Pan, Q., Nguyen, T. B., Ascher, D. B., and Pires, D. E. V. (2022). Systematic evaluation of computational tools to predict the effects of mutations on protein stability in the absence of experimental structures. Brief. Bioinform. 23:bbac025. doi: 10.1093/bib/bbac025

PubMed Abstract | Crossref Full Text | Google Scholar

Petersen, E., Ntoumi, F., Hui, D. S., Abubakar, A., Kramer, L. D., Obiero, C., et al. (2022). Emergence of new SARS-CoV-2 variant of concern omicron (B.1.1.529) – highlights Africa's research capabilities, but exposes major knowledge gaps, inequities of vaccine distribution, inadequacies in global COVID-19 response and control efforts. Int. J. Infect. Dis. 114, 268–272. doi: 10.1016/j.ijid.2021.11.040

PubMed Abstract | Crossref Full Text | Google Scholar

Pettersen, E. F., Goddard, T. D., Huang, C. C., Couch, G. S., Greenblatt, D. M., Meng, E. C., et al. (2004). UCSF chimera – a visualization system for exploratory research and analysis. J. Comput. Chem. 25, 1605–1612. doi: 10.1002/jcc.20084

PubMed Abstract | Crossref Full Text | Google Scholar

Planas, D., Saunders, N., Maes, P., Guivel-Benhassine, F., Planchais, C., Buchrieser, J., et al. (2022). Considerable escape of SARS-CoV-2 omicron to antibody neutralization. Nature 602, 671–675. doi: 10.1038/s41586-021-04389-z

PubMed Abstract | Crossref Full Text | Google Scholar

Plante, J. A., Liu, Y., Liu, J., Xia, H., Johnson, B. A., Lokugamage, K. G., et al. (2021). Spike mutation D614G alters SARS-CoV-2 fitness. Nature 592, 116–121. doi: 10.1038/s41586-020-2895-3

PubMed Abstract | Crossref Full Text | Google Scholar

Qu, P., Faraone, J., Evans, J. P., Zou, X., Zheng, Y. M., Carlin, C., et al. (2022). Neutralization of the SARS-CoV-2 omicron BA.4/5 and BA.2.12.1 subvariants. New Engl. J. Med. 386:2526–2528. doi: 10.1056/NEJMc2206725

Crossref Full Text | Google Scholar

Quintero, A. M., Eisner, M., Sayegh, R., Wright, T., Ramilo, O., Leber, A. L., et al. (2022). Differences in SARS-CoV-2 clinical manifestations and disease severity in children and adolescents by infecting variant. Emerg. Infect. Dis. 28, 2270–2280. doi: 10.3201/eid2811.220577

PubMed Abstract | Crossref Full Text | Google Scholar

Rodrigues, C. H. M., Pires, D. E. V., and Ascher, D. B. (2018). DynaMut: predicting the impact of mutations on protein conformation, flexibility and stability. Nucleic Acids Res. 46, W350–W355. doi: 10.1093/nar/gky300

Crossref Full Text | Google Scholar

Sahin, A., Karadag-Oncel, E., Buyuksen, O., Ekemen-Keles, Y., Ustundag, G., Elvan-Tuz, A., et al. (2023). The diversity in the clinical features of children hospitalized with COVID-19 during the nonvariant, alpha (B.1.1.7), Delta (B.1.617.2), and omicron (B.1.1.529) variant periods of SARS CoV-2: caution for neurological symptoms in omicron variant. J. Med. Virol. 95:e28628. doi: 10.1002/jmv.28628

Crossref Full Text | Google Scholar

Shen, N., Wu, Y. F., Chen, Y. W., Fang, X. Y., Zhou, M., Wang, W. Y., et al. (2023). Clinical characteristics of pediatric cases infected with the SARS-CoV-2 omicron variant in a tertiary children's medical center in Shanghai, China. World J. Pediatr. 19, 87–95. doi: 10.1007/s12519-022-00621-6

PubMed Abstract | Crossref Full Text | Google Scholar

Tagarro, A., Coya, O. N., Perez-Villena, A., Iglesias, B., Navas, A., Aguilera-Alonso, D., et al. (2022). Features of COVID-19 in children during the omicron wave compared with previous waves in Madrid, Spain. Pediatr. Infect. Dis. J. 41, e249–e251. doi: 10.1097/INF.0000000000003482

PubMed Abstract | Crossref Full Text | Google Scholar

Walls, A. C., Park, Y. J., Tortorici, M. A., Wall, A., McGuire, A. T., and Veesler, D. (2020). Structure, function, and antigenicity of the SARS-CoV-2 spike glycoprotein. Cell 183:1735. doi: 10.1016/j.cell.2020.11.032

PubMed Abstract | Crossref Full Text | Google Scholar

Wang, C., Horby, P. W., Hayden, F. G., and Gao, G. F. (2020). A novel coronavirus outbreak of global health concern. Lancet 395, 470–473. doi: 10.1016/S0140-6736(20)30185-9

PubMed Abstract | Crossref Full Text | Google Scholar

World Health Organization. (2022). WHO meets with Chinese officials on current COVID-19 situation. Available at: https://www.who.int/news/item/30-12-2022-who-meets-with-chinese-officials-on-current-covid-19-situation (Accessed May 10, 2023).

Google Scholar

Xinhua. (2022). China Focus: COVID-19 response further optimized with 10 new measures. Available at: https://english.news.cn/20221207/ca014c043bf24728b8dcbc0198565fdf/c.html (Accessed May 10, 2023).

Google Scholar

Yan, L., Huang, Y., Ge, J., Liu, Z., Lu, P., Huang, B., et al. (2022). A mechanism for SARS-CoV-2 RNA capping and its inhibition by nucleotide analog inhibitors. Cell 185, 4347–4360.e4317. doi: 10.1016/j.cell.2022.09.037

Crossref Full Text | Google Scholar

Zhang, L., Jackson, C. B., Mou, H., Ojha, A., Peng, H., Quinlan, B. D., et al. (2020). SARS-CoV-2 spike-protein D614G mutation increases virion spike density and infectivity. Nat. Commun. 11:6013. doi: 10.1038/s41467-020-19808-4

PubMed Abstract | Crossref Full Text | Google Scholar

Zhou, Z. L., Huang, C. L., Zhou, Z. C., Huang, Z. X., Su, L. L., Kang, S. S., et al. (2021). Structural insight reveals SARS-CoV-2 ORF7a as an immunomodulating factor for human CD14(+) monocytes. Iscience. 24:102187. doi: 10.1016/j.isci.2021.102187

PubMed Abstract | Crossref Full Text | Google Scholar

Zhu, N., Zhang, D., Wang, W., Li, X., Yang, B., Song, J., et al. (2020). A novel coronavirus from patients with pneumonia in China, 2019. N. Engl. J. Med. 382, 727–733. doi: 10.1056/NEJMoa2001017

Crossref Full Text | Google Scholar