Background: Coxsackievirus A16 (CVA16) is one of the major etiological agents of hand, foot and mouth disease (HFMD). This study aimed to investigate the molecular epidemiology and evolutionary characteristics of CVA16. Methods: Throat swabs were collected from children with HFMD and suspected HFMD during 2010–2019. Enteroviruses (EVs) were detected and typed by real-time reverse transcription-polymerase chain reaction (RT-PCR) and RT-PCR. The genotype, evolutionary rate, the most recent common ancestor, population dynamics and selection pressure of CVA16 were analyzed based on viral protein gene (VP1) by bioinformatics software. Results: A total of 4709 throat swabs were screened. EVs were detected in 3180 samples and 814 were CVA16 positive. More than 81% of CVA16-positive children were under 5 years old. The prevalence of CVA16 showed obvious periodic fluctuations with a high level during 2010–2012 followed by an apparent decline during 2013–2017. However, the activities of CVA16 increased gradually during 2018–2019. All the Beijing CVA16 strains belonged to sub-genotype B1, and B1b was the dominant strain. One B1c strain was detected in Beijing for the first time in 2016. The estimated mean evolutionary rate of VP1 gene was 4.49 × 10^–3 substitution/site/year. Methionine gradually fixed at site-23 of VP1 since 2012. Two sites were detected under episodic positive selection, one of which (site-223) located in neutralizing linear epitope PEP71. Conclusions: The dominant strains of CVA16 belonged to clade B1b and evolved in a fast evolutionary rate during 2010–2019 in Beijing. To provide more favorable data for HFMD prevention and control, it is necessary to keep attention on molecular epidemiological and evolutionary characteristics of CVA16.

Recombination events in human adenovirus (HAdV) have led to some new highly pathogenic or infectious types. It is vital to monitor recombinant HAdVs, especially in children with acute respiratory tract infections (ARIs). In the retrospective study, HAdV positive specimens were collected from pediatric patients with ARIs during 2015 to 2021, then typed by sequence analysis of the penton base, hexon and fiber gene sequence. For those with inconsistent typing results, a modified method with species-specific primer sets of a fiber gene sequence was developed to distinguish co-infections of different types from recombinant HAdV infections. Then, plaque assays combined with meta-genomic next-generation sequencing (mNGS) were used to reveal the HAdV genomic characteristics. There were 466 cases positive for HAdV DNA (2.89%, 466/16,097) and 350 (75.11%, 350/466) successfully typed with the most prevalent types HAdV-B3 (56.57%, 198/350) and HAdV-B7 (32.00%, 112/350), followed by HAdV-C1 (6.00%, 21/350). Among 35 cases (7.51%, 35/466) with inconsistent typing results, nine cases were confirmed as co-infections by different types of HAdVs, and 26 cases as recombinant HAdVs in six genetic patterns primarily clustered to species C (25 cases) in pattern 1–5, or species D (1 case) in pattern 6. The novel recombinant HAdV of species D was identified with multiple recombinant events among HAdV-D53, HAdV-D64, and HAdV-D8, and officially named as HAdV-D115. High-frequency recombination of HAdVs in six genetic recombination patterns were identified among children with ARIs in Beijing. Specifically, there is a novel Adenovirus D human/CHN/S8130/2023/115[P22H8F8] designed as HAdV D115.