Advertisement
Letter to the Editor| Volume 86, ISSUE 4, e84-e86, April 2023

Download started.

Ok

Genomic surveillance of SARS-CoV-2 in mainland China after ending the zero-COVID policy, December 2022–January 2023

  • Pengcheng Liu
    Affiliations
    Department of Clinical Laboratory, Children’s Hospital of Fudan University, National Children’s Medical Center, Shanghai 201102, China
    Search for articles by this author
  • Jin Xu
    Correspondence
    Correspondence to: Department of Clinical Laboratory, Children’s Hospital of Fudan University, 399 Wanyuan Rd, Minhang District, Shanghai 201102, China.
    Affiliations
    Department of Clinical Laboratory, Children’s Hospital of Fudan University, National Children’s Medical Center, Shanghai 201102, China
    Shanghai Institute of Infectious Disease and Biosecurity, Fudan University, Shanghai, China
    Search for articles by this author
Published:March 01, 2023DOI:https://doi.org/10.1016/j.jinf.2023.02.040
      Dear Editor,
      After implementing the dynamic zero-COVID policy for nearly three years, China has recently begun to adjust its COVID-19 prevention and control strategies with the attenuated pathogenicity of omicron subvariants and increasing vaccination coverage, notably by announcing the 10 measures on 7 December 2022.

      Xinhua. China Focus: COVID-19 response further optimized with 10 new measures; 2023. Available from: https://english.news.cn/20221207/ca014c043bf24728b8dcbc0198565fdf/c.html [Accessed February 15, 2023].

      A recent publication in this journal by Jingyi Liang et al. showed that the cumulative SARS-CoV-2 infection rate rose rapidly to 70% within three weeks after the ending of zero-COVID policy in Macao, China.
      • Liang J.
      • Liu R.
      • He W.
      • Zeng Z.
      • Wang Y.
      • Wang B.
      • et al.
      Infection rates of 70% of the population observed within 3 weeks after release of COVID-19 restrictions in Macao, China.
      Herein, we report the genomic characteristics and phylodynamic of SARS-CoV-2 complete genome sequences in mainland China after ending zero-COVID policy based on the viral genomic surveillance data deposited in the Global Initiative on Sharing All Influenza Data (GISAID) database.
      A total of 8254 complete SARS-CoV-2 genomic sequences from mainland China were extracted from the GISAID database during the period of December 2022 and January 2023 (date of access, February 6, 2023). Sequences were typed using the Phylogenetic Assignment of Named Global Outbreak Lineages (PANGOLIN; version 4.2) web tool.
      • O’Toole A.
      • Scher E.
      • Underwood A.
      • Jackson B.
      • Hill V.
      • McCrone J.T.
      • et al.
      Assignment of epidemiological lineages in an emerging pandemic using the pangolin tool.
      These genomes corresponded to 88 Pango-nomenclature-system-named subvariants, with the dominant lineages of BA.5.2.48 (4881/8254, 59.1%) and BF.7.14 (2223/8254, 26.9%), followed by BA.5.2.49 (525/8254, 6.4%), BA.5.2 (205/8254, 2.5%) and BF.7 (93/8254, 1.1%) (Fig. 1A). The proportion of these dominant lineages were not significantly changed over the 2 months of outbreak (Fig. 1B). While the other 83 lineages accounted for only 4.0%, such as BQ.1.1, BN.1.3 and BN.1.2 (Fig. 1A). What is more, the proportion of these lineages had fallen consistently over the time (Fig. 1B). All these findings demonstrated that the emerging BA.5.2.48 and BF.7.14 were the absolutely dominant drivers of the current COVID-19 outbreak, and the lineages found to be dominant internationally during the same period (i.e., XBB and BQ.1) were quickly cleared and did not prevail in China.

      CDC. COVID data tracker: variant proportions; 2023. Available from: https://covid.cdc.gov/covid-data-tracker/#variant-proportions [Accessed February 15, 2023].

      Fig. 1
      Fig. 1The composition of SARS-CoV-2 lineages circulating in mainland China during December 2022 to January 2023. (A) The percentage distribution of SARS-CoV-2 lineages during the whole period. (B) Temporal distributions of the SARS-CoV-2 lineages.
      In order to investigate the phylodynamics of BA.5.2.48 and BF.7.14 subvariants, we randomly chose complete genomic sequences of the subvariants from mainland China collected on or before January 31, 2023 in the GISAID database, respectively (date of access, February 16, 2023). To test for temporal structure in the datasets, we performed linear regressions of root-to-tip genetic distances against sampling dates using TempEst (version 1.5.3). The results indicated that both the BA.5.2.48 (n = 350) and BF.7.14 (n = 357) datasets were considered to have sufficient temporal signals for molecular dating after discarding several outliers (Supplemental Fig. S1). Then, Bayesian time-scaled phylogenetic analyses were performed using Markov Chain Monte Carlo (MCMC) methods implemented in BEAST software (version 1.10.4) to estimate the time of the most recent common ancestor (tMRCA), rate of evolution, as well as to infer the effective population size for both BA.5.2.48 and BF.7.14 subvariants. A GTR +Γ substitution model was used along with a strict clock model. A Bayesian skyline plot was used as the tree prior to estimate the median effective population size through time with a 95% highest posterior density (95% HPD). The MCMC chain length was set at 350 million and 300 million states for the BA.5.2.48 and BF.7.14, respectively. Convergence was assessed using Tracer (version 1.7.2), and 10% of states were removed to account for burn-in. Finally, maximum clade credibility (MCC) trees were generated using TreeAnnotator (version 1.10.4).
      Our Bayesian phylodynamic analysis revealed that the tMRCA for the BA.5.2.48 and BF.7.14 subvariants were estimated to be 16 July 2022 (95% HPD, 14 June–11 August 2022) and 9 August 2022 (95% HPD, 28 July–15 August 2022) with posterior probability (PP) = 1, respectively (Fig. 2AB). The estimated emerging dates of the subvariants were slight earlier than the first complete genomic sequences collected in GISAID database (15 August 2022 for BA.5.2.48 and BF.7.14). The mean substitution rates of BA.5.2.48 and BF.7.14 were estimated to be 5.353 × 10−4 (95% HPD, 4.141 × 10−4 to 6.893 × 10−4) substitution/site/year and 5.332 × 10−4 (95% HPD, 4.562 × 10−4 to 6.129 × 10−4) substitution/site/year, respectively. The results were in line with a previous study that estimated the average substitution rate of global SARS-CoV-2 (all lineages) was 5.83 × 10−4 substitution/site/year.
      • Chaguza C.
      • Hahn A.M.
      • Petrone M.E.
      • Zhou S.
      • Ferguson D.
      • Breban M.I.
      • et al.
      Accelerated SARS-CoV-2 intrahost evolution leading to distinct genotypes during chronic infection.
      Indeed, at the beginning of the pandemic, the evolutionary rate of SARS-CoV-2 was about 1.1 × 10−3 substitution /site/year, which means that the new variants present about 2 times slower than the original variants.
      • Candido D.S.
      • Claro I.M.
      • de Jesus J.G.
      • Souza W.M.
      • Moreira F.R.R.
      • Dellicour S.
      • et al.
      Evolution and epidemic spread of SARS-CoV-2 in Brazil.
      • Duchene S.
      • Featherstone L.
      • Haritopoulou-Sinanidou M.
      • Rambaut A.
      • Lemey P.
      • Baele G.
      Temporal signal and the phylodynamic threshold of SARS-CoV-2.
      The Bayesian skyline plot (BSP) showed that the effective population size for both BA.5.2.48 and BF.7.14 underwent expansions during the early phase of the outbreak in November 2022, matching the change of prevention and control policies by announcing the 20 measures during the same period in China.

      Xinhua. China Focus: China releases measures to optimize COVID-19 response; 2023. Available from: https://english.news.cn/20221111/d4399114a082438eaac32d08a02bf58d/c.html [Accessed February 15, 2023].

      For BA.5.2.48, the BSP analysis revealed one timeframe with a rapid exponential growth from around 7 November to 25 November 2022, continuing at a high average until the end of the sampling time (Fig. 2C). However, the effective population size for BF.7.14 increased earlier than BA.5.2.48 and showed two growth intervals over the time (Fig. 2D). The difference in the growth patterns of the effective population size of these two subvariants might be due to the different viral fitness and the adjustment of COVID-19 control strategies. Besides, which was unexpected, the difference possibly arose from bias due to the limited sample size.
      Fig. 2
      Fig. 2Phylodynamic analyses of BA.5.2.48 and BF.7.14 in mainland China by 31 January 2023. (A) Time-scaled maximum clade credibility (MCC) tree of BA.5.2.48. (B) Time-scaled maximum clade credibility (MCC) tree of BF.7.14. (C) Bayesian skyline plot (BSP) of BA.5.2.48. The colored area corresponds to the 95% credibility intervals of the highest probability density (95% HPD) region. (D) Bayesian skyline plot (BSP) of BF.7.14. The colored area corresponds to the 95% credibility intervals of the highest probability density (95% HPD) region. tMRCA, time of the most recent common ancestor; PP, posterior probability; Ne, effective population size.
      This study provides insights into the genomic epidemiology of SARS-CoV-2 in mainland China after ending the zero-COVID policy. We report that BA.5.2.48 and BF.7.14 were the dominant drivers of the current outbreak. Continuous genomic surveillance of SARS-CoV-2 is crucial for monitoring emerging variants with the potential to cause new COVID-19 outbreaks in China.

      Funding

      This study was funded by Key Development Program of Children's Hospital of Fudan University (No: EK2022ZX05).

      Declaration of Competing Interest

      All authors claim no conflict of interest.

      Acknowledgment

      We gratefully acknowledge all data contributors, i.e., the authors and their originating laboratories responsible for obtaining the specimens, and their submitting laboratories for generating the genetic sequence and metadata and sharing via the GISAID Initiative, on which this research is based.

      Appendix A. Supplementary material

      References

      1. Xinhua. China Focus: COVID-19 response further optimized with 10 new measures; 2023. Available from: https://english.news.cn/20221207/ca014c043bf24728b8dcbc0198565fdf/c.html [Accessed February 15, 2023].

        • Liang J.
        • Liu R.
        • He W.
        • Zeng Z.
        • Wang Y.
        • Wang B.
        • et al.
        Infection rates of 70% of the population observed within 3 weeks after release of COVID-19 restrictions in Macao, China.
        J Infect. 2023; 4453 (00035-X): S0163
        • O’Toole A.
        • Scher E.
        • Underwood A.
        • Jackson B.
        • Hill V.
        • McCrone J.T.
        • et al.
        Assignment of epidemiological lineages in an emerging pandemic using the pangolin tool.
        Virus Evol. 2021; 7: veab064
      2. CDC. COVID data tracker: variant proportions; 2023. Available from: https://covid.cdc.gov/covid-data-tracker/#variant-proportions [Accessed February 15, 2023].

        • Chaguza C.
        • Hahn A.M.
        • Petrone M.E.
        • Zhou S.
        • Ferguson D.
        • Breban M.I.
        • et al.
        Accelerated SARS-CoV-2 intrahost evolution leading to distinct genotypes during chronic infection.
        Cell Rep Med. 2023; 4: 100943
        • Candido D.S.
        • Claro I.M.
        • de Jesus J.G.
        • Souza W.M.
        • Moreira F.R.R.
        • Dellicour S.
        • et al.
        Evolution and epidemic spread of SARS-CoV-2 in Brazil.
        Science. 2020; 369: 1255-1260
        • Duchene S.
        • Featherstone L.
        • Haritopoulou-Sinanidou M.
        • Rambaut A.
        • Lemey P.
        • Baele G.
        Temporal signal and the phylodynamic threshold of SARS-CoV-2.
        Virus Evol. 2020; 6: veaa061
      3. Xinhua. China Focus: China releases measures to optimize COVID-19 response; 2023. Available from: https://english.news.cn/20221111/d4399114a082438eaac32d08a02bf58d/c.html [Accessed February 15, 2023].