Five-year (2017-2022) evolutionary dynamics of human coronavirus OC43 in southern France based on whole genome next-generation sequencing

HCoV-OC43 genomes and their evolution are scarcely studied worldwide and in France with only 361 genomes available as of October 2023. Here we implemented an in-house PCR amplification system to obtain retrospectively by next-generation sequencing then analyze HCoV-OC43 genomes for infections diagnosed with this virus in southern France between 02/2017 and 10/2022. Multiplex PCR amplification using a set of in-house primers designed using the Gemi software was carried out on residues of HCoV-OC43 RNA-positive nasopharyngeal samples, before next-generation sequencing (NGS) using Illumina technology on a NovaSeq 6000 instrument. HCoV-OC43 genome assembly, bioinformatic analyses, and phylogeny reconstruction were then carried out using CLC Genomics, Mafft, BioEdit, Nextstrain, Nextclade, MEGA, iTOL and RDP4 softwares. A total of 34 PCR primer pairs were designed for amplification then NGS of HCoV-OC43 genomes. A total of 185 genomes were obtained, 17, 79 and 89 belonging to genotypes G, J and K, respectively. These three genotypes circulated exclusively or co-circulated according to the year. A total of 303, 940, and 1,300 amino acid substitutions were detected in genotype G, J and K, respectively, compared with reference genomes of the same genotype dating back to 2017-2018. Possible recombinations were detected for 10 HCoV-OC43 genomes classified in genotypes K or J. Overall, the present study more than doubled the set of HCoV-OC43 genomes available worldwide for the 2017-2022 period, and contributed to the monitoring of the HCoV-OC43 evolutionary dynamics.