Objectives HCoV-HKU1 diversity and evolution was scarcely studied. We performed next-generation sequencing (NGS) and analysis of HCoV-HKU1 genomes over five years. Methods NGS used Illumina technology on NovaSeq 6000 following whole genome PCR amplification by a in-house set of primers designed using Gemi and PrimalScheme. Genome assembly and analyses used CLC Genomics, Mafft, BioEdit, Nextstrain, Nextclade, MEGA and iTol bioinformatic tools. Spike molecular modelling and dynamics simulations used Molegro Molecular Viewer/Hyperchem. Results Twenty-eight PCR systems allowed obtaining 158 HCoV-HKU1 genomes including 69 and 89 of genotype A and B, respectively. Both genotypes co-circulated during the study period but one predominated each year. 1,683 amino acid substitutions including 80 in ≥10 genomes were detected in genotype A relatively to a 2004 reference. H512R in spike, first detected in 2009 and reported as involved in antibody neutralization, was found in all genotype A, almost always with V387I and K478N, and was predicted here to significantly improve cellular TMPRSS2 protein binding. 1,802 amino acid substitutions including 64 in ≥10 genomes were detected in genotype B relatively to a 2005 reference. Conclusion This study substantially expands the global set of HCoV-HKU1 genomes. Genomics with protein structural analyses contributed to our understanding of HCoV-HKU1 evolution.

The diversity and evolution of the genomes of Human Bocavirus (HBoV), which causes respiratory diseases, have been scarcely studied. Here, we aimed to obtain and characterize HBoV genomes from patients’s nasopharyngeal samples collected between 2017-2022 period (five years). Next-generation sequencing (NGS) used Illumina technology after having implemented using GEMI an in-house multiplex PCR amplification strategy. Genomes were assembled and analyzed with CLC Genomics, Mafft, BioEdit, MeV, Nextclade, MEGA, and iTol. A total of 213 genomes were obtained. Phylogeny classified them all as of Bocavirus 1 (HBoV1) species. Five HBoV1 genotypes determined by hierarchical clustering analysis of 27 variable genome positions were scattered over the study period although with differences in yearly prevalence. A total of 167 amino acid substitutions were detected. Besides, coinfection was observed for 52% of the samples, rhinoviruses then adenoviruses (HAdVs) being the most common viruses. Principal component analysis showed that HBoV1 genotype α tended to be significantly correlated with HAdV co-infection. Subsequent HAdV typing for HBoV1-positive samples and negative controls demonstrated that HAdVC species predominated but HAdVB was that significantly HBoV1-associated. Overall, we described here the first HBoV1 genomes sequenced for France. HBoV1 and HAdVB association deserves further investigation.