Viruses, Vol. 15, Pages 73: Evolutionary and Genetic Recombination Analyses of Coxsackievirus A6 Variants Associated with Hand, Foot, and Mouth Disease Outbreaks in Thailand between 2019 and 2022

Conceptualization, J.P., N.W. and Y.P.; methodology, N.S., N.K. and S.K.; software, J.P. and J.C.; validation, N.S., N.K. and S.K.; formal analysis, J.P. and J.C.; investigation, J.P. and N.W.; resources, R.Y., N.K. and S.K.; data curation, J.P. and Y.P.; writing—original draft preparation, J.P.; writing—review and editing, J.P., J.C., N.W. and Y.P.; visualization, N.W. and Y.P.; supervision, N.W. and Y.P.; project administration, N.W. and Y.P.; funding acquisition, N.W. and Y.P. All authors have read and agreed to the published version of the manuscript.

Figure 1. Temporal distribution of enteroviral infections, by epidemiological month, Thailand, 2019–2022 (n = 755).

Figure 1. Temporal distribution of enteroviral infections, by epidemiological month, Thailand, 2019–2022 (n = 755).

Viruses 15 00073 g001

Figure 2. Phylogenetic comparisons of sequences from partial VP1 and 3Dpol. Phylogenetic comparisons of sequences from (A) partial VP1 (nt positions 2630–3286, 657 bp) and (B) partial 3Dpol (nt positions 6188–7388, 1200 bp) using samples and sequences of other coxsackievirus (CV)-A6 variants obtained from GenBank. Maximum-likelihood trees were reconstructed using 1000 bootstrap re-samples to demonstrate the robustness of groupings; values ≥ 70% are shown.

Figure 2. Phylogenetic comparisons of sequences from partial VP1 and 3Dpol. Phylogenetic comparisons of sequences from (A) partial VP1 (nt positions 2630–3286, 657 bp) and (B) partial 3Dpol (nt positions 6188–7388, 1200 bp) using samples and sequences of other coxsackievirus (CV)-A6 variants obtained from GenBank. Maximum-likelihood trees were reconstructed using 1000 bootstrap re-samples to demonstrate the robustness of groupings; values ≥ 70% are shown.

Viruses 15 00073 g002

Figure 3. Unrooted phylogenetic trees of coxsackievirus (CV)-A6 strains based on different genomic regions. Maximum-likelihood trees were constructed based on the 5′-UTR and the P1, P2, and 3Dpol regions. Sequences are numbered according to the numbering of the Gdula prototype strain (GenBank accession number AY421764) and previously referenced sequences. Bootstrap values at key nodes are shown as percentages of 1000 replicates. The nearly complete genome sequences of the 20 CV-A6/TH strains, obtained in this study, are highlighted using colored dots. The scale bars represent the number of substitutions per site.

Figure 3. Unrooted phylogenetic trees of coxsackievirus (CV)-A6 strains based on different genomic regions. Maximum-likelihood trees were constructed based on the 5′-UTR and the P1, P2, and 3Dpol regions. Sequences are numbered according to the numbering of the Gdula prototype strain (GenBank accession number AY421764) and previously referenced sequences. Bootstrap values at key nodes are shown as percentages of 1000 replicates. The nearly complete genome sequences of the 20 CV-A6/TH strains, obtained in this study, are highlighted using colored dots. The scale bars represent the number of substitutions per site.

Viruses 15 00073 g003

Figure 4. Time-scaled phylogenetic tree of VP1 sequences of coxsackievirus (CV)-A6 variants. The phylogenetic tree was generated using CV-A6 VP1 sequences characterized in this study (red) and previously published sequences. A maximum clade credibility tree was constructed from 10,000 trees, sampled from the posterior distribution with mean node ages. Clades described in previous studies (A, B, C, D1, D2, and D3) are shown. Several subclades were associated with severe hand, foot, and mouth outbreaks in Thailand; posterior probability support is provided. Branch colors denote different countries.

Figure 4. Time-scaled phylogenetic tree of VP1 sequences of coxsackievirus (CV)-A6 variants. The phylogenetic tree was generated using CV-A6 VP1 sequences characterized in this study (red) and previously published sequences. A maximum clade credibility tree was constructed from 10,000 trees, sampled from the posterior distribution with mean node ages. Clades described in previous studies (A, B, C, D1, D2, and D3) are shown. Several subclades were associated with severe hand, foot, and mouth outbreaks in Thailand; posterior probability support is provided. Branch colors denote different countries.

Viruses 15 00073 g004

Figure 5. (A) Schematic diagram of the enteroviral (EV) genome structure. (B) Recombination analysis of putative recombinant coxsackievirus (CV)-A6 strains. SimPlot analysis of the recombinant strain CVA6-CU3494. The vertical axis represents the nucleotide sequence similarity (%) between the putative recombinant strain and other EV-A lineages, constructed using SimPlot v3.5.1 with a window size of 200 nt and a step size of 20 nt. (C) Bootscan plot of recombination events between the daughter strain EVA71-CU3494 and major (CV-A6) or minor (CV-A10) parental strains. The bootstrapping support value was computed using the RDP4 program with a window size of 200 nt, a step size of 10 nt, and 1000 bootstrap replicates. (D) Amino acid changes between the CV-A6 subclade and the D3/A and D3/Y strains.

Figure 5. (A) Schematic diagram of the enteroviral (EV) genome structure. (B) Recombination analysis of putative recombinant coxsackievirus (CV)-A6 strains. SimPlot analysis of the recombinant strain CVA6-CU3494. The vertical axis represents the nucleotide sequence similarity (%) between the putative recombinant strain and other EV-A lineages, constructed using SimPlot v3.5.1 with a window size of 200 nt and a step size of 20 nt. (C) Bootscan plot of recombination events between the daughter strain EVA71-CU3494 and major (CV-A6) or minor (CV-A10) parental strains. The bootstrapping support value was computed using the RDP4 program with a window size of 200 nt, a step size of 10 nt, and 1000 bootstrap replicates. (D) Amino acid changes between the CV-A6 subclade and the D3/A and D3/Y strains.

Viruses 15 00073 g005

Table 1. Demographic characteristics and characterization of 20 CV-A6/TH strains based on structural sequences.

Table 1. Demographic characteristics and characterization of 20 CV-A6/TH strains based on structural sequences.

CV-A6/TH StrainsAge (years)SexSampleCollection DateSubcladeNearest
CV-A6 Strain in Structural Region (Accession No.)Similarity with the
Nearest CV-A6 Strain (% nt)GenBank Accession No.CU34821MStool9-Sep-19D3/AChina/2018 (MN845848)99.5OP896719CU34946FThroat swab18-Sep-19D3/YChina/2017 (MN845834)98.5OP896715CU34953FThroat swab21-Sep-19D3/AChina/2018 (OL839942)99.2OP896717CU34963MThroat swab23-Sep-19D3/YChina/2017 (MN845834)98.4OP896716CU35108MThroat swab8-Oct-19D3/AChina/2018 (OL839942)99.2OP896718CU35141MThroat swab18-Oct-19D3/YChina/2017 (MN845834)98.6OP896714CU35473MStool24-Jan-20D3/YChina/2017 (MN845834)98.3OP896713CU37741MStool22-Jun-22D3/AChina/2018 (OL839942)98.0OP896722CU37803MThroat swab1-Jul-22D3/AChina/2018 (OL839942)98.0OP896725CU37864FThroat swab1-Jul-22D3/AChina/2018 (OL839942)98.0OP896729CU37921MNasal swab5-Jul-22D3/AChina/2018 (OL839942)98.0OP896728CU38122MThroat swab1-Aug-22D3/AChina/2018 (OL839942)98.0OP896726CU38353M Throat swab31-Aug-22D3/AChina/2018 (OL839942)97.9OP896727CU38462MThroat swab8-Sep-22D3/AChina/2018 (OL839942)97.9OP896724CU38473FThroat swab8-Sep-22D3/AChina/2018 (OL839942)97.7OP896723CU39133FThroat swab20-Sep-22D3/AChina/2018 (OL839942)97.9OP896730CU39196MThroat swab20-Sep-22D3/AChina/2018 (OL839942)97.7OP896731CU39244MThroat swab21-Sep-22D3/AChina/2018 (OL839942)97.7OP896721CU39251MThroat swab21-Sep-22D3/AChina/2018 (MN845848)98.5OP896720CU39334MThroat swab21-Sep-22D3/AChina/2018 (OL839942)97.7OP896732

留言 (0)

沒有登入
gif