A total of 172 archived samples from clinically suspected HFMD cases were included in this study. Of which, 148 samples were collected from secondary care hospitals in Kozhikode, followed by 18 from Wayanad, 4 from Malappuram districts of Kerala state, and 2 from Shimoga district of Karnataka state in South India. Enteroviruses were identified in 107 (62%) cases out of 172 clinically suspected cases of HFMD. Among these 172 samples, 107 were pan enterovirus RNA positive by real-time RT-PCR. Of the 107 cases, 89 (83%) were below five years of age, and 18 (17%) were between 5 and 65 years (Table 1). Eleven (10.3%) adults showed HFMD infection in the study. The number of affected males were 68 (64%), while 39 (36%) were females. The mean age of the patients with HFMD was 4.4 years, ranging from 9 months to 65 years. Meanwhile, the mean age of the male and female population was 5.4 and 2.8, respectively. The most common clinical features observed among the enterovirus associated HFMD cases were fever (90, 84%) and vesicular rash (107, 100%). Respiratory symptoms such as cough (18, 16%), coryza (6, 6%), breathlessness (1, 1%) and sore throat (9, 8%) were also observed among patients. Some of the HFMD cases had oral ulcers (3, 3%), headache (4, 4%), vomiting (2, 2%), myalgia (1, 1%), abdominal pain (1, 1%) and lymphadenopathy (1, 1%). Rash in patients with hand, foot and mouth disease was illustrated in Fig. 2. We also noticed other rash types such as papule (25, 23%), maculopapular rash (19, 18%), crusted lesion (11, 10%), macule (7, 5%), erythema (8, 8%), pustule (3, 3%), erosion (3, 3%) and bullae (1, 1%) in patients. Rash was seen in hand (87, 81%), foot (91, 85%), mouth (51, 48%), soles (39, 36%), palms (41, 38%), buttocks (32, 30%), and other parts of the body (< 20%).
Table 1
Demographic and clinical characteristics of enterovirus associated HFMD
Variable
|
CVA6
(n = 69)
|
CVA16
(n = 16)
|
CVA10
(n = 1)
|
EV types not identified
(n = 21)
|
Total EVs
(n = 107)
|
|
n (%)
|
n (%)
|
n (%)
|
n (%)
|
n (%)
|
Gender
|
|
|
|
|
|
Male
|
44 (64)
|
10 (63)
|
1 (100)
|
12 (57)
|
67 (63)
|
Female
|
25 (36)
|
6 (37)
|
0
|
9 (43)
|
40 (37)
|
Age
|
|
|
|
|
|
Mean (± SD)
|
2.56 ± 3.9
|
12.68 ± 17.3
|
NA
|
4.6 ± 5.9
|
4.4 ± 8.4
|
Range (in years)
|
0.8 to 30
|
1 to 65
|
NA
|
0.8 to 22
|
0.8 to 65
|
Age group
(in years)
|
|
|
|
|
|
≤ 5
|
63 (91)
|
9 (56)
|
1 (100)
|
16 (76)
|
89 (83)
|
> 5
|
6 (9)
|
7 (44)
|
0
|
5 (24)
|
18 (17)
|
Clinical Parameters
|
|
|
|
|
|
Cutaneous rash
|
69 (100)
|
16 (100)
|
1 (100)
|
21 (100)
|
107 (100)
|
Fever
|
58 (84)
|
11 (69)
|
1 (100)
|
20 (95)
|
90 (84)
|
Cough
|
10 (15)
|
3 (19)
|
0
|
5 (19)
|
18 (16)
|
Sore throat
|
3 (4)
|
4 (25)
|
0
|
2 (10)
|
9 (8)
|
Coryza
|
3 (4)
|
0
|
0
|
3 (14)
|
6 (6)
|
Headache
|
1 (1)
|
0
|
0
|
3 (14)
|
4 (4)
|
Oral ulcers
|
1 (1)
|
2 (13)
|
0
|
0
|
3 (3)
|
Vomiting
|
0
|
0
|
0
|
2 (10)
|
2 (2)
|
Breathlessness
|
1 (1)
|
0
|
0
|
0
|
1 (1)
|
Abdominal pain
|
0
|
0
|
0
|
1 (5)
|
1 (1)
|
Lymphadenopathy
|
1 (1)
|
0
|
0
|
0
|
1 (1)
|
Myalgia
|
0
|
0
|
0
|
1 (1)
|
1 (1)
|
Types of Rash
|
|
|
|
|
|
Vesicle
|
69 (100)
|
16 (100)
|
1 (100)
|
21 (100)
|
107 (100)
|
Papule
|
17 (25)
|
3 (19)
|
0
|
5 (25)
|
25 (23)
|
Maculopapular rash
|
12 (17)
|
2 (13)
|
0
|
5 (26)
|
19 (18)
|
Crusted lesion
|
11 (16)
|
0
|
0
|
0
|
11 (10)
|
Erythema
|
5 (7)
|
0
|
1 (100)
|
2 (5)
|
8 (8)
|
Macule
|
5 (7)
|
1 (6)
|
0
|
1 (10)
|
7 (5)
|
Pustule
|
1 (1)
|
0
|
0
|
2 (10)
|
3 (3)
|
Erosion
|
1 (1)
|
2 (13)
|
0
|
0
|
3 (3)
|
Bullae
|
0
|
1 (6)
|
0
|
0
|
1 (1)
|
Rash Distribution
|
|
|
|
|
|
Hand
|
55 (80)
|
13 (81)
|
0
|
19 (91)
|
87 (81)
|
Foot
|
61 (88)
|
10 (63)
|
1 (100)
|
19 (91)
|
91 (85)
|
Mouth
|
33 (48)
|
9 (56)
|
0
|
9 (43)
|
51 (48)
|
Palms
|
30 (44)
|
5 (31)
|
1 (100)
|
5 (24)
|
41 (38)
|
Soles
|
31 (45)
|
5 (31)
|
1 (100)
|
2 (10)
|
39 (36)
|
Buttocks
|
21 (30)
|
6 (38)
|
1 (100)
|
4 (19)
|
32 (30)
|
Trunk
|
13 (19)
|
1 (6)
|
0
|
5 (24)
|
19 (17)
|
Knees
|
11 (16)
|
3 (19)
|
0
|
2 (10)
|
16 (15)
|
Thigh
|
8 (12)
|
1 (6)
|
0
|
3 (14)
|
12 (11)
|
Elbow
|
8 (12)
|
1 (6)
|
0
|
2 (19)
|
13 (12)
|
Face
|
6 (9)
|
0
|
0
|
2 (10)
|
8 (8)
|
Enterovirus RNA was detected in vesicular (97, 91%) and throat swabs (10, 9%) of 107 cases of HFMD. Of the typed enterovirus cases (86, 80%), 69 (64%) were CVA6, 16 (15%) were CVA16 and 1 (1%) was CVA10. However, 21 (20%) cases of enteroviruses were not typable by any of the molecular typing strategies used in this study. In 2015, 24 HFMD cases were caused by CVA6 while one was CVA16. There were 20 HFMD cases due to CVA6 infection, while nine were CVA16 in 2016. In the following year 2017, HFMD was identified in 25 cases by CVA6, six by CVA16 and one by CVA10. HFMD cases were predominantly found in the rainy season from June to November.
The PCR and sequencing enabled CVA6 typing in 27 samples based on 5’VP1, followed by VP2 based typing in five samples and 5’NCR based typing in nine samples. Meanwhile, CVA6 specific nested RT-PCR based on 3’VP1 region enabled the typing of CVA6 strains in the other 28 samples, of which sequencing was possible only in 15 samples. Out of 16 (15%) samples positive for CVA16, PCR and sequencing could identify CVA16 in six samples based on 5’VP1, followed by VP2 targeted typing in two samples and 5’NCR based typing in one sample. Meanwhile, CVA16 specific real-time RT-PCR could type CVA16 in the other seven samples. CVA10 strain was identified only in one sample (1%) by 5’VP1 based PCR and sequencing. However, 21 samples were not typed by any of the enterovirus typing methods described in the study. A total of 66 sequences obtained from the study were deposited to GenBank database under the following accession numbers: MG885750, MG885751, MG840455 - MG840478, MG840480 - MG840487, MG869706 - MG869710, MG885752 - MG885768, MH160045 - MH160049, and MH160051 - MH160055.
The phylogenetic dendrogram of CVA6 based on 5’VP1 (398 bp), 3’VP1 (475 bp), VP2 (540 bp) and 5’NCR (301 bp) partial sequences are shown in Fig. 3. Based on partial 5’VP1 (n = 74), 22 CVA6 strains clustered together with other Indian strain (MH539787, 2016) with 96–100% nucleotide identity, followed by five strains showed 97–97.5% similarity with Australian strain (MH111055, 2017) in subclade D3 (Fig. 3A). Phylogeny of 3’VP1 sequences (n = 67) showed 14 CVA6 strains clustered together with Indian strain (MH539787, 2016) in 96.2–98.7% identity and one with 98.1% similarity with Australian strain (MH111055, 2017) in subclade D3 (Fig. 3B). VP2 sequence (n = 39) based phylogeny showed two CVA6 subclade D3 strains in high nucleotide similarity (99.1–99.3%) with Indian strain (MH539787, 2016), while the other three strains formed a cluster with Australian strain (MH111055, 2017) in 96.9–97.78% identity (Fig. 3C). Of the 5’NCR (301 bp) targeted phylogenetic tree of 41 sequences, four CVA6 strains showed high similarity (99.0%) with Indian strain (MH539787, 2016), while the other five strains showed 98–99% identity to Australian strain (MH111055, 2017) in subclade D3 (Fig. 3D). Phylogenetic analysis of CVA16 based on 5’VP1 (413 bp), VP2 (533 bp) and 5’NCR (301 bp) regions depicted in Fig. 4. The phylogenic tree based on 5’VP1 (n = 34) showed five CVA16 strains in high similarity (96.4–97.6%) with Indian strain (KY792583, 2013), followed by one with 98.3% identity to Indian strain (KY792578, 2013) in B1c subclade (Fig. 4A). Meanwhile, two CVA16 strains showed high nucleotide identity (97.2%) with Indian strain (KY792584, 2015) based on VP2 targeted phylogeny (n = 33) (Fig. 4B). One CVA16 strain was 99% identical to Indian strain (KY792577, 2012) based on 5’NCR sequences (n = 42) (Fig. 4C). CVA10 strain from the study showed 95.6% similarity with Indian strain (MH118041, 2017) in subclade D based on 5’VP1 region (410 bp, n = 33), depicted in Fig. 5. Nucleotide divergence among the inter-clades of CVA6 showed 13.6–30.9% for 5’VP1, 13.1–24.0% for 3’VP1, 15.6–19.2% for VP2 and 4.7–11.3% for 5’NCR. Inter-clade nucleotide divergence of CVA16 was 24.9–27.9% for 5’VP1, 21.4–24.8% for VP2 and 5.7–10.0% for 5’NCR, whereas nucleotide divergence of CVA10 inter-clades was 14.4–27.8% based on 5’VP1.