Keyi Zhang1
Mei Hong2
Yong Zhang1,3*
Zhenzhi Han1
Jinbo Xiao1Huanhuan Lu1Yang Song1Dongmei Yan1Dongyan Wang1Shuangli Zhu1
Wenbo Xu1,3Guizhen Wu1*- 1WHO WPRO Regional Polio Reference Laboratory, NHC Key Laboratory for Biosafety, NHC Key Laboratory for Medical Virology, National Institute for Viral Disease Control and Prevention, Chinese Center for Disease Control and Prevention, Beijing, China
- 2Tibet Center for Disease Control and Prevention, Lhasa, China
- 3Center for Biosafety Mega-Science, Chinese Academy of Sciences, Wuhan, China
Enterovirus B75 (EV-B75) is a newly identified serotype of the enterovirus B species. To date, only 112 cases related to EV-B75 have been reported worldwide, and research on EV-B75 is still limited with only two full-length genome sequences available in GenBank. The present study reported seven EV-B75 sequences from a child with acute flaccid paralysis and six asymptomatic close contacts in Shigatse, Tibet. Phylogenetic analysis revealed that the Tibetan strain was possibly imported from neighboring India. Seroepidemiological analyses indicated that EV-B75 has not yet caused a large-scale epidemic in Tibet. Similarity plots and boot scanning analyses revealed frequent intertypic recombination in the non-structural region of all seven Tibet EV-B75 strains. All seven Tibetan strains were temperature-sensitive, suggesting their poor transmissibility in the environment. Overall, though the seven Tibetan strains did not cause large-scale infection, prevention and control of the novel enterovirus cannot be underestimated.
Introduction
Enteroviruses (genus Enterovirus, family Picornaviridae) are divided into 15 species assigned to enteric enterovirus (EV) A–L and respiratory rhinovirus A–C (Walker et al., 2019). Human pathogenic enteroviruses are classified under EV-A–D, including poliovirus (PV) (I–III), coxsackievirus (CV), echovirus, and novel enteroviruses. These diverse serotypes cause a diverse array of clinical features such as aseptic meningitis, hand foot and mouth disease (HFMD), neonatal sepsis-like disease, pancreatitis, encephalitis, myocarditis, pericarditis, acute flaccid paralysis (AFP), and severe respiratory disease. As per the International Committee on Taxonomy of Viruses (ICTV) classification, enteroviruses newly discovered after 1976 are named as novel enteroviruses, namely EV-D68 to the current EV-A121, with a total of 53 genotypes.
Most novel enteroviruses are not the main pathogens of intestinal-transmitted diseases and rarely cause concentrated outbreaks worldwide. However, increasing novel enteroviruses have been found and defined in recent years using innovations in molecular typing methods. Occasionally, a previously rare enterovirus genotype may cause large-scale outbreaks or continuous epidemics. EV-A71, a member of the novel enterovirus family, was defined in 1974. Following the outbreak of HFMD caused by EV-A71 infection in Linyi, Shandong Province, in 2007 (Zhang et al., 2009) and in Fuyang, Anhui Province, in 2008 (Zhang et al., 2010), laboratory surveillance of HFMD was conducted nationwide in China. Since then, the reported cases of EV-A71 have increased annually, resulting in EV-A71 being the predominant HFMD pathogen and one of the leading pathogens for severe and fatal cases of HFMD. EV-D68, another rarely reported enterovirus, caused an endemic outbreak of severe respiratory illness in the United States, with more than 1,150 laboratory-confirmed cases (Greninger et al., 2015). Thus, comprehensive surveillance and precise investigation of viral prevalence and transmission is critical for the control and intervention of enterovirus diseases.
Currently, based on the novel enterovirus isolates available in GenBank database, most novel enterovirus genotypes except for EV-D68 and EV-A71 have been isolated in the mainland of China, and the clinical manifestation of the vast majority of cases related to novel enterovirus is AFP (Oberste et al., 2004; Bingjun et al., 2008; Tao et al., 2010, 2013, 2014a; Xu et al., 2011a,b; Sun et al., 2013; Hu et al., 2014; Liu et al., 2014; Lu et al., 2014; Tang et al., 2014; Zhang et al., 2014, 2016; Fan et al., 2015; Song et al., 2017; Han et al., 2018; Huang et al., 2018; Wang et al., 2019). A considerable number of strains have been isolated from healthy people (Wu et al., 2013; Hu et al., 2014; Tao et al., 2014b; Fan et al., 2015; Tang et al., 2016; Zhang et al., 2016; Song et al., 2017; Han et al., 2018; Zhang et al., 2019), directly indicating that the population carries the novel enteroviruses and is at risk of developing the related diseases. Of these, a large number of cases arising from some novel enterovirus serotypes have already been reported in China, such as EV-B85 in Xinjiang (Sun et al., 2013) (33 cases, 2011, AFP and his asymptomatic contacts), EV-A90 in Xinjiang (Huang et al., 2018) (5 cases, 2011, AFP) and Shandong (Tao et al., 2013) (3 cases, 2001, AFP), EV-C96 in Yunnan (Bingjun et al., 2008) (18 cases, 1999-2010, AFP), and Jiangxi Province (13 cases, 2016, unknown), and EV-C104 in Shandong (Xiang et al., 2013) (five cases, 2011–2012, lower respiratory tract infection). However, due to the lack of full-length genome sequences and detailed epidemiological information, the study of novel enteroviruses has not attracted sufficient attention.
Enterovirus B75 (EV-B75) is a novel enterovirus serotype within EV-B, identified in 2004 (Oberste et al., 2004). In this study, we report the isolation of EV-B75 from a case of AFP and six asymptomatic close contacts from Shigatse, Tibet, China, in 2007. We determined the full-length genomic sequence of these seven Tibetan isolates, and used phylogenetic analysis to track the infection source and to determine recombination events incorporating sequences from the current outbreaks, from public databases. Serum neutralization tests were used to evaluate the population prevalence of EV-B75 in Tibet, and temperature sensitivity tests were used to determine the environmental transmission capacity of the seven strains.
Materials and Methods
Sample Collection
The collection, use, and analysis of all samples involved in this study were approved by the Ethics Review Committee of the National Institute for Viral Disease Control and Prevention, Chinese Center for Disease Control and Prevention. Written informed consent for the use of their clinical samples was acquired from the parents of all the children included in this study. All experimental protocols were approved by the National Institute for Viral Disease Control and Prevention, and the methods were carried out in accordance with the approved guidelines.
Seven strains of EV-B75 were isolated from the feces of a patient with AFP and his close contacts during an epidemiological investigation in Shigatse, Tibet, by the Center for Disease Control and Prevention of Tibet Autonomous Region in 2007. All children from which the EV-B75 virus had been isolated, were under 5 years of age. We named the seven EV-B75 Tibetan strains as Y16/XZ/CHN/2007, Y17/XZ/CHN/2007, Y18/XZ/CHN/2007, Y20/XZ/CHN/2007, Y24/XZ/CHN/2007, Y25/XZ/CHN/2007, and Y26/XZ/CHN/2007, hereafter abbreviated as Y16, Y17, Y18, Y20, Y24, Y25, and Y26. In addition, 112 serum samples from Tibetan children used in this study were also collected by the Tibet Autonomous Region Center for Disease Control and Prevention from December 2010 to January 2011, and sent to the National Polio Reference Laboratory for monitoring antibody levels after polio vaccination. All samples were collected from healthy children (aged 1–14 years) in Tibet, including 36 from Shigatse and 76 from Lhasa.
Virus Isolation and Plaque Purification
All collected fecal specimens were pre-treated according to the standard operating method of the National Polio Reference Laboratory to obtain fecal suspensions, which were then, respectively, inoculated onto human rhabdomyosarcoma (RD) cells for continuous culture for three generations to isolate the virus. The cell line was provided by the WHO Global Poliovirus Specialized Laboratory in the United States and was originally purchased from the American Type Culture Collection (Manassas, VA, United States). When a complete EV-like cytopathic effect (CPE) was observed, the cell cultures were harvested.
All seven viral strains (Y16, Y17, Y18, Y20, Y24, Y25, and Y26) were further purified using plaque tests and then amplified on RD cells. Briefly, RD cells were cultured in a 6-well plate. After the cells were confluent, 1 ml of virus solution diluted 105 times was added and the cells were incubated for 1.5 h (36°C, 5% CO2). The 6-well plate was slightly shaken every 30 min to avoid accumulation of the virus solution. Next, a 2.8% low melting point agarose (melting point 65°C) solution, sterilized at high temperature and mixed with the same amount of 2x MEM to form a maintenance medium, was slowly poured into the culture to fix the virus, and the plate was then cultured in an inverted position. On the next day, neutral red dye solution was added, followed by culturing for 2–8 days. When spots were observed, they were picked and further amplified on RD cells.
Molecular Typing and Full-Length Genome Sequencing
We used Primscript One-Step RT-PCR Kit Ver.2 (Takara, Dali, China) for RNA fragment amplification. The VP1 region sequence was amplified by using primers E490/E492 (E490: TGIGTIYTITGYRTICATTGGAT, E492: GGRTTIGTIGWYTGCCA) (Oberste et al., 2006), and the full-length genomic sequence was amplified by the primers designed by the only two full-length genome sequences in GenBank.
The PCR reaction process was divided into three stages: the first stage was a reverse transcription stage, which was maintained at 50°C for 30 min to reverse transcribe RNA into cDNA, and at 94°C for 2 min to inactivate the reverse transcriptase; the second stage was “denaturation-annealing-extension” amplification stage. The DNA was helically released to form single-stranded DNA at 94°C for 30 s, and primers were bound to the template at 50°C for 30 s and dNTP at 72°C for 80 s, 40 cycles in total. The third stage was the extension stage, which only required 72°C for 7 min and the product was finally preserved at 4°C.
All amplification-enriched products were purified for sequencing using the QIAQUICK gel extract kit (Qiagen, Hillen, Germany), and the amplicons were sequenced on an ABI 3130 Genetic Analyzer (Applied Biosystems, Foster City, CA, United States) as described above.
The results indicated that the VP1 sequences of the seven Tibet strains showed the highest similarity with those of the EV-B75 strain in the GenBank databases, which reached a percentage of as high as 90.95%. We identified all seven strains as EV-B75 based on the enterovirus molecular typing criteria.
Phylogenetic and Recombination Analysis
For VP1 phylogenetic analysis, the seven Tibetan EV-B75 sequences generated in this study were combined with all the publicly available VP1 gene sequences in GenBank database with known sampling dates. Maximum likelihood trees were constructed using the GTR + G model as suggested by ModelGenerator0.85 and were implemented in MEGA7 with 1,000 bootstrap replicates. Potential inter-typic and intra-typic recombination was analyzed using SimPlot (version 3.5.1) with a 200-nt window moving with 20-nt steps. To investigate the recombination relationship among the seven Tibetan EV-B75 strains and other enterovirus genotypes, the EV-B75 genome was spliced at the suggested recombination break points. BLAST searches were separately performed using the P2 and P3 coding regions. Closely related enterovirus genome sequences with similarity higher than 85% were selected.
Test of Neutralization
The Y24 strain was selected as the attack virus in the neutralizing test because it showed the highest titer among the seven Tibetan EV-B75 strains. In total, 112 serum samples were inactivated at 56°C for 30 min, and then diluted at five gradients from 1:4 to 1:1024 (1:4, 1:16, 1:64, 1:256, and 1:1024) for detection. The mixtures comprising a virus culture infection dose (CCID50) of 100 (50 μl), human RD cell suspension (100 μl), and serum dilution (50 μl) were then incubated at 36°C in a 5% CO2 incubator. With observation for 7 days, the highest dilution of serum that protected 50% of the cultures was recorded based on an EV-like CPE. If the neutralization antibody titer was observed at a dilution higher than 1:8, the serum sample was considered positive and the geometric mean titer was subsequently calculated.
Assay of Temperature Sensitivity
The temperature sensitivity of the seven plaque-purified EV-B75 strains and two selected control strains (HTYT-ARL-AFP02F/XJ/CHN/2011, showing no temperature sensitivity and HTYPS-QDH11F/XJ/CHN/2011, showing temperature sensitivity) were assayed on monolayer RD cells in 24-well plates. The 24-well plates were inoculated with 50 μL of undiluted virus stocks. Two different incubation plates were used: one was placed at an optimum temperature of 36°C, and the other was placed at the supra-optimal temperature of 39.5°C for virus propagation. After adsorption at 36°C or 39.5°C for 1 h, non-adsorbed virus cultures were discarded and 100 μL of maintenance fluid was added to each well. The 24-well plates were then incubated at 36 and 39.5°C, respectively, and sampled at six time points after infection (4, 8, 16, 24, 48, and 72 h). The CCID50 was calculated by the end-point dilution of monolayer RD cells in a 96-well plate at 36°C. Virus isolates that showed a drop of more than 2 log values at different temperatures were considered temperature-sensitive.
Results
Global Epidemiology of Novel EV-B75
To date, there are 112 gene sequences in the GenBank database, associated with various clinical manifestations including AFP (Oberste et al., 2004; Bingjun et al., 2008; Rao et al., 2012; Laxmivandana et al., 2013; Sadeuh-Mba et al., 2013; Tao et al., 2014a; Fernandez-Garcia et al., 2017; Maan et al., 2019; Sadeuh-Mba et al., 2019) (47 isolates worldwide), aseptic meningitis (Avellon et al., 2006) (8 isolates in Spain, 2005), encephalitis (Lewthwaite et al., 2010; Kumar et al., 2012) (13 isolates in India, 2006–2010), diarrhea (Rao et al., 2013; Patil et al., 2015) (7 isolates in India, 2005–2011), respiratory diseases (Opanda et al., 2016; Hellferscee et al., 2017) (5 isolates in total), and others (Oberste et al., 2004) (Figure 1). With the addition of the seven Tibetan EV-B75 isolates from this study, no significant or typical temporal and spatial aggregation of EV-B75 was observed worldwide. However, after EV-B75 was molecularly finalized in 2004, the number of isolates was significantly increased in the following years compared with those before finalization. However, the number of related reports was significantly reduced in recent years, possibly due to a lag in reporting. Among these, the case reports of EV-B75 in Spain were relatively concentrated, with one report of six aseptic meningitis-related cases isolated from children’s hospitals in 2005, which was the only report of EV-B75 related to aseptic meningitis worldwide, and another report of the eight encephalitis-associated strains isolated the next year (Cabrerizo et al., 2008). In contrast, there were relatively more isolates in India, with a total of 25. Cases were reported every year from 2005 to 2011, and the types of cases involved were relatively diverse, with AFP (Rao et al., 2012) and encephalitis (Kumar et al., 2012) as the main presentations. Notably, only seven EV-B75 isolates from patients with acute diarrhea worldwide are available in India (Rao et al., 2013; Patil et al., 2015).
Figure 1. Time, region, and clinical distribution of EV-B75 worldwide.
Notably, there are only two EV-B75 full-length genome sequences in the GenBank database, one of which is the prototype strain (AY556070, USA-OK85-10362, hereafter referred as strain OK85-10362) and the other is isolated from a child with AFP in Shandong Province, China in 1997 (GQ329840, 102-SD-CHN-97).
Full-Length Genome Analysis of the Seven Tibetan EV-B75 Strains
Human enteroviruses are small, non-enveloped icosahedral viruses, with a positive-sense single stranded RNA genome, composed of approximately 7,500 nucleotides (nt), including a well-known long open reading frame (ORF) and a newly discovered short second ORF (Guo et al., 2019). The full-length genome sequences of the seven Tibet EV-B75 strains were determined and analyzed. The results showed that the genomes of all strains were 7,425 nucleotides in length. The overall nucleotide composition of the seven strains was as follows: 27.87–28.05% A, 24.08–24.34% T, 23.06–23.35% C, and 24.57–24.69% G. Aligned with the prototype strain of EV-B75 from the United States (OK-10362), all seven strains contained some deletions and insertions, including one deletion at position 16, one deletion at position 24, three deletions at position 107, one insertion at position 572, one deletion at position 7,355, and one insertion at position 7,421 in that order.
The long ORF, franked by a 5′ untranslated region (UTR) and a 3′UTR, was initiated by AUG744 and ended at position 7323, totaling 6579 nucleotides in length, encoding a single polypeptide of 2,193 amino acid, which was cleaved into 3 polyprotein precursors P1, P2, and P3, encoding the structural proteins (VP1 to VP4) and non-structural proteins. The second ORF, located in the VI domain of the viral internal ribosomal entry site (IRES) was initiated by AUG593 with the WIGHPV (TGG ATT GGC CAT CCG GTG) domain, which is the specific domain of ORF2p and is highly conserved among different species of HEVs; it is important to trigger viral replication and particle release from intestinal epithelial cells (Guo et al., 2019). Non-coding 3′UTRs were 103 nucleotides in length, preceding a long poly(A) tail.
Pairwise comparisons of nucleotide sequences and the deduced amino acid sequences were conducted among the seven Tibetan EV-B75 strains. The full-length genome nucleotide and amino acid similarities were 94.4–99.9 and 99.7–99.9%, respectively, suggesting that the seven Tibetan isolates might have a common origin or ancestor and that their isolation time was relatively concentrated. Furthermore, the nucleotide sequence and amino acid sequence similarities of the VP1 region and P1 region with the EV-B75 prototype strain were 85.3–85.5 and 98.3–98.6%, and 82.3–85.7 and 98.3–100%, respectively. Considering the cutoff value of 75% VP1 similarity for EV typing, these data verified the molecular typing results. As expected, the nucleotide sequence and amino acid sequence similarities of the P2 and P3 regions with the EV-B75 prototype strain were 79.3–81.7 and 97.3–99.0%, and 77.9–82.5 and 95.3–100%, respectively, but they shared a distinctly high identity with other EV-B prototype strains other than EV-B75, suggesting that recombination might occur in these coding regions.
Phylogenetic Analysis of the Seven Tibet EV-B75 Strains
A maximum likelihood phylogenetic tree was constructed for the complete VP1 sequences of the seven Tibet EV-B75 strains and all sequences available in GenBank to study their evolutionary relationships and epidemic patterns (Figure 2).
Figure 2. Maximum likelihood phylogenetic tree based on the entire VP1 coding region sequences of EV-B75 available from the GenBank database. The seven Tibetan EV-B75 strains isolated in this study are indicated by red circles, and the prototype of EV-B75 is indicated by a green circle. The scale bars indicate the substitution per site per year. The numbers at the nodes indicate the bootstrap support for the node (percentage of 1,000 bootstrap replicates).
All the seven Tibet EV-B75 strains were grouped into one branch (with a bootstrap support value of 100%), which verified the high similarity of their VP1 region. They formed a single lineage and were clustered with the Indian strain (JN204061, India-N-694), which was isolated from a patient with AFP in 2007–2009 in India, indicating that they had a closer genetic relationship with the Indian strain. In addition, Shigatse is located at the southwest border of China, adjacent to India. Due to the geographical advantage and the coincidence of onset time and case characteristics, we suspect that the seven Tibetan strains may have originated from India.
Shandong strain (GQ329840, 102-SD-CHN-97) isolated from a patient with AFP in 1997 formed a distinctive lineage with another Chinese strain (AY570801, CHN-FJ96-71) isolated from Fujian Province in 1996, occurred on a separate branch and were grouped separately. Considering that the geographic distance between the two provinces is more than 1000 km, the close genetic relationship possibly reflects a long-distance transmission. This result suggests that more than two different EV-B75 genotypes are existing in China.
Phylogenetic trees were also constructed based on the entire VP1, P1, P2, and P3 coding region nucleotide sequences of the seven Tibet EV-B75 strains along with the prototype of EV-B from GenBank database. Consistent with the phylogenetic tree of the P1 coding region (Figure 3B), the VP1 phylogenetic tree (Figure 3A) indicated that the seven Tibet EV-B75 strains formed a single lineage and were clustered with the EV-B75 prototype (OK-10362) as expected, verifying the primary molecular typing results. Unlike the VP1 and P1 phylogenetic trees, the P2 and P3 coding regions phylogenetic trees showed that the seven Tibet EV-B75 strains were clustered together and shared the highest similarity with the prototype strains of EV-B88, EV-B77, and EV-B86, respectively, rather than with the prototype of EV-B75, suggesting the occurrence of recombination between EV-B75 and other EV-B serotypes (Figures 3C,D).
Figure 3. Maximum likelihood phylogenetic tree based on the VP1, P1, P2, and P3 coding regions of the prototype sequence of all EV-B in the GenBank database and the seven Tibet EV-B75 strains in this study. The seven Tibet EV-B75 strains in this study are indicated by red circles, and the prototype of EV-B75 is indicated by a green circle. The scale bars indicate the substitution per site per year. The numbers at the nodes indicate the bootstrap support for the node (percentage of 1,000 bootstrap replicates). Coding sequences of (A) VP1, (B) P1, (C) P2, and (D) P3 are shown.
Recombination Events Between the Seven Tibet EV-B75 Strains and Other EV-B Strains
Potential evidence of recombination in the genome of the seven Tibet EV-B75 strains was investigated via similarity plot and boot scanning analysis. Screening of closely related sequences available in GenBank was performed using BLAST from NCBI. The P2 and P3 regions of the seven strains were tested separately. All sequences downloaded from NCBI showing over 85% similarity with the query sequence were used to construct a phylogenetic tree for the P3 region and to set a series in SimPlot (Version3.5.1) if clustered in the same branch and the same serotype. In the P1 region, all seven Tibet EV-B75 strains showed the highest similarity to the EV-B75 prototype strain. In other regions, multiple recombination events were revealed with relatively higher similarities to other EV-B viruses. The similarity plots and boot scanning analysis results indicated that the seven Tibetan EV-B75 strains recombined with the EV-107 strain (AB426609-TN94-0349) at position 4620–5240 and with the CVB1 strain (JN797615-1167438) at position 5340–7425 as shown in Figure 4.
Figure 4. Similarity and boot scanning analysis of the seven Tibetan EV-B75 strains with the prototype of EV-B75, a field EV-107 strain (AB426609-TN-94-0349), and a field CVB1 strain (JN797615-1167438). All seven Tibetan EV-B75 strains were set in a series as a query sequence.
Seroprevalence of EV-B75 in Tibet
Considering that the first case (Y16) of EV-B75 infection had caused six asymptomatic infections in the same region, we examined the prevalence of EV-B75 among Tibetan population. A total of 112 serum samples from children aged 1–14 years old in Tibet were tested. These were collected from December 2010 to January 2011, including 36 serum samples from Shigatse (27 males and 9 females) and 76 serum samples from Lhasa (44 males and 32 females). Of these, only 3 serum samples were seropositive for EV-B75 (≥1:8) with a total positive rate of 2.61%, of which one was from Lhasa and two were from Shigatse. The geometric mean titer of the positive serum was only 1:4.52. Compared with the seroprevalence rate of EV-71 or even some rarely isolated novel enteroviruses in China, this rate was rather low (Zhu et al., 2010). The low serum positive rate suggested that EV-B75 was limited in the population of Shigatse and Lhasa regions, and that the proportion of people who had been infected or in contact with this serotype was low, due to which there was a limited spread and no epidemic in the region.
Environmental Tolerance of the Seven Tibetan EV-B75 Strains
The tolerance of these seven strains to the environment determines their ability to spread. The viral replication ability of these seven EV-B75 strains was compared at different temperatures to determine whether they were sensitive to high temperature (39.5°C). The virus titer of all seven Tibet EV-B75 strains cultured at 36°C for 72 h was higher than 2 logs compared to that at the elevated temperature (39.5°C) (Figure 5). These results indicated that all seven Tibetan EV-B75 strains were temperature-sensitive compared to the EV-B85 strain (Sun et al., 2013) which showed non-temperature sensitivity and the EV-B106 strain (Song et al., 2017) which showed temperature sensitivity.
Figure 5. Temperature sensitivity test curves of the seven Tibetan EV-B75 strains. Blue and orange lines represent the growth trends of the viruses on RD cells at 36 and 39.5°C, respectively. The Xinjiang EV-B85 strain (HTYT-ARL-AFP02F/XJ/CHN/2011, showing no temperature sensitivity) and the EV-B106 strain (KS-MGTH90F/XJ/CHN/2011, showing temperature sensitivity) were used as experimental controls. (A) Strain Y16/XZ/CHN/2007; (B) strain Y17/XZ/CHN/2007; (C) strain Y18/XZ/CHN/2007; (D) strain Y20/XZ/CHN/2007; (E) strain Y24/XZ/CHN/2007; (F) strain Y25/XZ/CHN/2007; (G) strain Y26/XZ/CHN/2007; (H) strain HTYT-ARL-AFP02F/XJ/CHN/2011 (EV-B85); (I) strain KS-MGTH90F/XJ/CHN/2011 (EV-B106).
Discussion
Based on the sequence similarity among the seven Tibetan EV-B75 strains from Shigatse, Tibet, we speculated that they had the same origin or ancestor. Based on the clustering analysis of phylogenetic trees, combined with the geographical location and clinical features of EV-B75 cases isolated in India in 2007, we speculated that the Tibetan strains were imported from India. Combined with the phylogenetic analysis of Indian strains, the AFP and encephalitis isolates in India are at least clustered into three different clusters, suggesting that EV-B75 has more than one undergoing coevolution, and that its epidemic may be far beyond the currently reported epidemic situation. This further confirmed our speculation that the Tibetan case was imported from India. The Shandong strain and Fujian strain from China were clustered in the phylogenetic tree, separate from other branches, suggesting that there might be two existing branches of EV-B75 in China.
We also investigated the seroepidemiology of EV-B75 in Tibet, China and found a very low seropositive rate, which might indicate the seven strains isolated in Tibet, China in 2007 to be the first batch of EV-B75 serotype isolates found, and that these have not widely circulated in the population. This strengthens the hypothesis that the seven Tibet EV-B75 strains were from India. The results of temperature sensitivity also suggested that the transmissibility of EV-B75 was limited and that it was not resistant to high temperature, which also contributed to the relatively limited spread of the Tibetan strain.
Although the transmission range of EV-B75 in Tibet was speculated to be limited in this study, EV-B75 has two existing circulations in China, and the transmission potential of the novel enterovirus in China cannot be underestimated. Further, due to the seroprevalence of other novel enteroviruses reported in China, even though there are only a few isolates, the seropositive rate of novel enteroviruses in China is not low, indicating the potential for an outbreak (Hu et al., 2014; Zhang et al., 2014; Fan et al., 2015; Song et al., 2017; Han et al., 2018; Huang et al., 2018; Lu et al., 2020; Xiao et al., 2020). The proportion of people who have been exposed to the virus is relatively high, and it thus has great transmission potential in China.
Recombinant analysis revealed that EV-B75 strains were extensively recombined with other EV-B enteroviruses in the P2 and P3 regions, and that recombination was indeed common in enteroviruses, especially in these two regions. This serotype might thus undergo gene sequence exchange and evolution with other enteroviruses in a human or environmental epidemic; further, this genotype has existed in the human population for some time and has been epidemic for some time, rather than being a new-born serotype. Thus, it needs to be monitored.
At present, the clinical symptoms of EV-B75 serotype-related strains in GenBank are relatively diverse, involving AFP, aseptic meningitis, encephalitis, acute diarrhea, and some flu-like symptoms, and these have been collected from feces, cerebrospinal fluid, and pharyngeal swabs. Although there are currently 112 EV-B75-related cases, limited studies have targeted EV-B75 as the causative agent. The relationship between pathogens and diseases is far from being conclusive and needs to be proved using more detailed epidemiological information, more full-length genome sequences, and observations of pathological changes in animals. In this study, we explored the possible sources of this outbreak of EV-B75 in Tibet, investigated the population prevalence of EV-B75 and the transmission ability of these strains in Tibet, and answered some epidemiological questions about this outbreak. However, due to the limited case information in GenBank, the gene sequences with different clinical manifestations were not compared. Afterward, based on the existing experiments, we will conduct animal experiments to study the dynamic distribution characteristics and pathological changes of the virus in different tissues and organs of animals after the attack, in order to try to explain the potential relationship between the pathogen and the disease. Our study on EV-B75 is just a stepping stone, and additional studies are needed to supplement.
China is a country with a long border shared with many countries, which are facing great pressure on the health system. In addition, with the progress of globalization and increasingly frequent communication between China and other countries, the pressure on the health system at China’s border areas will also increase. For instance, in 2011, China experienced an outbreak caused by wild poliovirus type I imported from Pakistan, affecting 10 young children and 11 adults, resulting in 2 deaths (Luo et al., 2013). As a lesson learned from these cases, more attention should be paid for the prevention of epidemics along the border.
Data Availability Statement
The full-length genomic sequences detected in this study have been deposited to China Virus Identification Net under accession numbers CVIN_AA002552–CVIN_AA002558 and to GenBank under accession numbers MW183133–MW183139.
Author Contributions
KZ, YZ, GW, and WX: conceptualization. KZ, YZ, ZH, JX, HL, YS, DY, DW, and SZ: formal analysis. MH: resources. KZ: writing—original draft preparation. YZ: writing—review and editing and funding acquisition. YZ and GW: project administration. All authors have read and agreed to the published version of the manuscript.
Funding
This study was supported by the National Key Technology R&D Program of China (Project Nos. 2017ZX10104001 and 2018ZX10711001), and Beijing Natural Science Foundation (Project No. L192014). We also acknowledge the funding received from the Key Technologies R&D Program of the National Ministry of Science (Project Nos. 2018ZX10713002 and 2018ZX10713001-003). The funding body was not involved in the study design, clinical sample collection, data analysis, and interpretation or writing of the manuscript.
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
References
Avellon, A., Rubio, G., Palacios, G., Casas, I., Rabella, N., Reina, G., et al. (2006). Enterovirus 75 and aseptic meningitis, Spain, 2005. Emerg. Infect. Dis. 12, 1609–1611.
Bingjun, T., Yoshida, H., Yan, W., Lin, L., Tsuji, T., Shimizu, H., et al. (2008). Molecular typing and epidemiology of non-polio enteroviruses isolated from Yunnan Province, the People’s Republic of China. J. Med. Virol. 80, 670–679. doi: 10.1002/jmv.21122
Cabrerizo, M., Echevarria, J. E., Gonzalez, I., de Miguel, T., and Trallero, G. (2008). Molecular epidemiological study of HEV-B enteroviruses involved in the increase in meningitis cases occurred in Spain during 2006. J. Med. Virol. 80, 1018–1024. doi: 10.1002/jmv.21197
Fan, Q., Zhang, Y., Hu, L., Sun, Q., Cui, H., Yan, D., et al. (2015). A novel recombinant enterovirus type EV-A89 with low epidemic strength in Xinjiang. China. Sci. Rep. 5:18558. doi: 10.1038/srep18558
Fernandez-Garcia, M. D., Kebe, O., Fall, A. D., and Ndiaye, K. (2017). Identification and molecular characterization of non-polio enteroviruses from children with acute flaccid paralysis in West Africa, 2013-2014. Sci. Rep. 7:3808. doi: 10.1038/s41598-017-03835-3831
Greninger, A. L., Naccache, S. N., Messacar, K., Clayton, A., Yu, G., Somasekar, S., et al. (2015). A novel outbreak enterovirus D68 strain associated with acute flaccid myelitis cases in the USA (2012–14): a retrospective cohort study. Lancet Infect. Dis. 15, 671–682. doi: 10.1016/s1473-3099(15)70093-70099
Guo, H., Li, Y., Liu, G., Jiang, Y., Shen, S., Bi, R., et al. (2019). A second open reading frame in human enterovirus determines viral replication in intestinal epithelial cells. Nat. Commun. 10:4066. doi: 10.1038/s41467-019-12040-12049
Han, Z., Zhang, Y., Huang, K., Cui, H., Hong, M., Tang, H., et al. (2018). Genetic characterization and molecular epidemiological analysis of novel enterovirus EV-B80 in China. Emerg. Microbes Infect. 7:193. doi: 10.1038/s41426-018-0196-199
Hellferscee, O., Tempia, S., Walaza, S., Variava, E., Dawood, H., Wolter, N., et al. (2017). Enterovirus genotypes among patients with severe acute respiratory illness, influenza-like illness, and asymptomatic individuals in South Africa, 2012-2014. J. Med. Virol. 89, 1759–1767. doi: 10.1002/jmv.24869
Hu, L., Zhang, Y., Hong, M., Zhu, S., Yan, D., Wang, D., et al. (2014). Phylogenetic evidence for multiple intertypic recombinations in enterovirus B81 strains isolated in Tibet. China. Sci. Rep. 4:6035. doi: 10.1038/srep06035
Huang, K., Zhang, Y., Song, Y., Cui, H., Yan, D., Zhu, S., et al. (2018). Antigenic characteristics and genomic analysis of novel EV-A90 enteroviruses isolated in Xinjiang. China. Sci. Rep. 8:10247. doi: 10.1038/s41598-018-28469-28469
Kumar, A., Shukla, D., Kumar, R., Idris, M. Z., Misra, U. K., and Dhole, T. N. (2012). Molecular epidemiological study of enteroviruses associated with encephalitis in children from India. J. Clin. Microbiol. 50, 3509–3512. doi: 10.1128/JCM.01483-1412
Laxmivandana, R., Yergolkar, P., Gopalkrishna, V., and Chitambar, S. D. (2013). Characterization of the non-polio enterovirus infections associated with acute flaccid paralysis in South-Western India. PLoS One 8:e61650. doi: 10.1371/journal.pone.0061650
Lewthwaite, P., Perera, D., Ooi, M. H., Last, A., Kumar, R., Desai, A., et al. (2010). Enterovirus 75 encephalitis in children, southern India. Emerg. Infect. Dis. 16, 1780–1782. doi: 10.3201/eid1611.100672
Liu, X., Tao, Z., Wang, H., Lin, X., Song, L., Li, Y., et al. (2014). Complete genome analysis of human enterovirus B73 isolated from an acute flaccid paralysis patient in Shandong. China. Virus Genes 49, 38–44. doi: 10.1007/s11262-014-1077-1075
Lu, H., Hong, M., Zhang, Y., Xiao, J., Zhang, M., Zhang, K., et al. (2020). A novel interspecies recombinant enterovirus (Enterovirus A120) isolated from a case of acute flaccid paralysis in China. Emerg. Microbes Infect. 9, 1733–1743. doi: 10.1080/22221751.2020.1796527
Lu, J., Zheng, H., Zhang, Y., Guo, X., Wu, D., Li, H., et al. (2014). Whole genomic sequence and replication kinetics of a new enterovirus C96 isolated from Guangdong. China with a different cell tropism. PLoS One 9:e86877. doi: 10.1371/journal.pone.0086877
Luo, H. M., Zhang, Y., Wang, X. Q., Yu, W. Z., Wen, N., Yan, D. M., et al. (2013). Identification and control of a poliomyelitis outbreak in Xinjiang. China. N. Engl. J. Med. 369, 1981–1990. doi: 10.1056/NEJMoa1303368
Maan, H. S., Dhole, T. N., and Chowdhary, R. (2019). Identification and characterization of nonpolio enterovirus associated with nonpolio-acute flaccid paralysis in polio endemic state of Uttar Pradesh. Northern India. PLoS One 14:e0208902. doi: 10.1371/journal.pone.0208902
Oberste, M. S., Maher, K., Williams, A. J., Dybdahl-Sissoko, N., Brown, B. A., Gookin, M. S., et al. (2006). Species-specific RT-PCR amplification of human enteroviruses: a tool for rapid species identification of uncharacterized enteroviruses. J. Gen. Virol. 87(Pt 1), 119–128.
Oberste, M. S., Michele, S. M., Maher, K., Schnurr, D., Cisterna, D., Junttila, N., et al. (2004). Molecular identification and characterization of two proposed new enterovirus serotypes. EV74 and EV75. J. Gen. Virol. 85(Pt 11), 3205–3212. doi: 10.1099/vir.0.80148-80140
Opanda, S. M., Wamunyokoli, F., Khamadi, S., Coldren, R., and Bulimo, W. D. (2016). Genotyping of enteroviruses isolated in Kenya from pediatric patients using partial VP1 region. Springerplus 5:158. doi: 10.1186/s40064-016-1834-1830
Patil, P. R., Chitambar, S. D., and Gopalkrishna, V. (2015). Molecular surveillance of non-polio enterovirus infections in patients with acute gastroenteritis in Western India: 2004-2009. J. Med. Virol. 87, 154–161. doi: 10.1002/jmv.23992
Rao, C. D., Yergolkar, P., and Shankarappa, K. S. (2012). Antigenic diversity of enteroviruses associated with nonpolio acute flaccid paralysis. India, 2007-2009. Emerg. Infect. Dis. 18, 1833–1840. doi: 10.3201/eid1811.111457
Rao, D. C., Ananda Babu, M., Raghavendra, A., Dhananjaya, D., Kumar, S., and Maiya, P. P. (2013). Non-polio enteroviruses and their association with acute diarrhea in children in India. Infect. Genet. Evol. 17, 153–161. doi: 10.1016/j.meegid.2013.04.011
Sadeuh-Mba, S. A., Bessaud, M., Massenet, D., Joffret, M. L., Endegue, M. C., Njouom, R., et al. (2013). High frequency and diversity of species C enteroviruses in Cameroon and neighboring countries. J. Clin. Microbiol. 51, 759–770.
Sadeuh-Mba, S. A., Kavunga-Membo, H., Joffret, M. L., Yogolelo, R., Endegue-Zanga, M. C., Bessaud, M., et al. (2019). Genetic landscape and macro-evolution of co-circulating Coxsackieviruses A and Vaccine-derived polioviruses in the democratic republic of Congo, 2008-2013. PLoS Negl. Trop. Dis. 13:e0007335. doi: 10.1371/journal.pntd.0007335
Song, Y., Zhang, Y., Fan, Q., Cui, H., Yan, D., Zhu, S., et al. (2017). Phylogenetic Characterizations of highly mutated EV-B106 recombinants showing extensive genetic exchanges with other EV-B in Xinjiang. China. Sci. Rep. 7:43080. doi: 10.1038/srep43080
Sun, Q., Zhang, Y., Zhu, S., Tian, H., Huang, G., Cui, H., et al. (2013). Transmission of human enterovirus 85 recombinants containing new unknown serotype HEV-B donor sequences in Xinjiang Uighur autonomous region. China. PLoS One 8:e55480. doi: 10.1371/journal.pone.0055480
Tang, J., Li, Q., Tian, B., Zhang, J., Li, K., Ding, Z., et al. (2016). Complete Genome Analysis of an Enterovirus EV-B83 Isolated in China. Sci. Rep. 6:29432. doi: 10.1038/srep29432
Tang, J., Tao, Z., Ding, Z., Zhang, Y., Zhang, J., Tian, B., et al. (2014). Complete genome characterization of a novel enterovirus type EV-B106 isolated in China, 2012. Sci. Rep. 4:4255. doi: 10.1038/srep04255
Tao, Z., Cui, N., Lin, X., Wang, H., Song, L., Li, Y., et al. (2014a). Complete genome sequence analysis of an enterovirus 75 isolate from China. Arch. Virol. 159, 2809–2813. doi: 10.1007/s00705-014-2014-2016
Tao, Z., Yuan, Q., Lin, X., Wang, S., Liu, Y., Ji, F., et al. (2014b). Molecular characterization of enteroviruses including a new type EV-C99 isolated from Xinjiang students in Shandong. China in 2011. Sci. Rep. 4:6564. doi: 10.1038/srep06564
Tao, Z., Cui, N., Xu, A., Liu, Y., Song, L., Li, Y., et al. (2013). Isolation and genomic characterization of three enterovirus 90 strains in Shandong. China. Arch. Virol. 158, 479–483. doi: 10.1007/s00705-012-1517-1512
Tao, Z., Cui, N., Xu, A., Wang, H., Song, L., Li, Y., et al. (2010). Genomic characterization of an enterovirus 97 strain isolated in Shandong. China. Virus Genes 41, 158–164. doi: 10.1007/s11262-010-0496-491
Walker, P. J., Siddell, S. G., Lefkowitz, E. J., Mushegian, A. R., Dempsey, D. M., Dutilh, B. E., et al. (2019). Changes to virus taxonomy and the international code of virus classification and nomenclature ratified by the international committee on taxonomy of viruses (2019). Arch. Virol. 164, 2417–2429. doi: 10.1007/s00705-019-04306-w
Wang, Q., Ji, F., Wang, S., Lin, X., Tao, Z., and Xu, A. (2019). Complete genome characterization of three enterovirus C96 isolates in China. Arch. Virol. 164, 2183–2186. doi: 10.1007/s00705-019-04291-4290
Wu, W., Xu, W. B., Chen, L., Chen, H. L., Liu, Q., Wang, D. L., et al. (2013). Molecular identification and analysis of human enteroviruses isolated from healthy children in Shenzhen. China from 2010 to 2011. PLoS One 8:e64889. doi: 10.1371/journal.pone.0064889
Xiang, Z., Xie, Z., Wang, Z., Ren, L., Xiao, Y., Li, L., et al. (2013). Human enterovirus genotype C104. China. Emerg. Infect. Dis. 19, 689–691. doi: 10.3201/eid1904.121435
Xiao, J., Zhang, Y., Hong, M., Han, Z., Zhang, M., Song, Y., et al. (2020). Phylogenetic characteristics and molecular epidemiological analysis of novel enterovirus EV-B83 isolated from Tibet. China. Sci. Rep. 10:6630. doi: 10.1038/s41598-020-63691-63694
Xu, A., Tao, Z., Lin, X., Liu, Y., Zhang, Y., Song, L., et al. (2011a). The complete genome sequence of an enterovirus 76 isolate in China reveals a recombination event. Arch. Virol. 156, 1685–1689. doi: 10.1007/s00705-011-1067-z
Xu, A., Tao, Z., Wang, H., Zhang, Y., Song, L., Li, Y., et al. (2011b). The complete genome analysis of two enterovirus 96 strains isolated in China in 2005 and 2009. Virus Genes 42, 323–330. doi: 10.1007/s11262-011-0584-x
Zhang, J., Liu, H., Huang, X., Sun, H., Yang, Z., and Ma, S. (2019). Complete genome sequence of a human enterovirus C99 strain isolated from a healthy child in Yunnan. China, in 2013. Microbiol. Resour. Announc. 8, e01489-1418. doi: 10.1128/MRA.01489-1418
Zhang, Y., Hong, M., Sun, Q., Zhu, S., Tsewang, Li, X., et al. (2014). Molecular typing and characterization of a new serotype of human enterovirus (EV-B111) identified in China. Virus Res. 183, 75–80. doi: 10.1016/j.virusres.2014.01.002
Zhang, Y., Sun, Q., Cui, H., Yan, D., Fan, Q., Song, Y., et al. (2016). Circulation of multiple serotypes of highly divergent enterovirus C in the Xinjiang Uighur autonomous region of China. Sci. Rep. 6:33595. doi: 10.1038/srep33595
Zhang, Y., Tan, X. J., Wang, H. Y., Yan, D. M., Zhu, S. L., Wang, D. Y., et al. (2009). An outbreak of hand, foot, and mouth disease associated with subgenotype C4 of human enterovirus 71 in Shandong. China. J. Clin. Virol. 44, 262–267. doi: 10.1016/j.jcv.2009.02.002
Zhang, Y., Zhu, Z., Yang, W., Ren, J., Tan, X., Wang, Y., et al. (2010). An emerging recombinant human enterovirus 71 responsible for the 2008 outbreak of hand foot and mouth disease in Fuyang city of China. Virol. J. 7:94. doi: 10.1186/1743-422X-7-94
Keywords: enterovirus B75, seroprevalence, temperature sensitivity, Tibet, phylogenetic analyses
Citation: Zhang K, Hong M, Zhang Y, Han Z, Xiao J, Lu H, Song Y, Yan D, Wang D, Zhu S, Xu W and Wu G (2021) Molecular Epidemiological, Serological, and Pathogenic Analysis of EV-B75 Associated With Acute Flaccid Paralysis Cases in Tibet, China. Front. Microbiol. 11:632552. doi: 10.3389/fmicb.2020.632552
Received: 23 November 2020; Accepted: 22 December 2020;
Published: 13 January 2021.
Edited by:
Lihua Xiao, South China Agricultural University, ChinaReviewed by:
Zexin Tao, Shandong Center for Disease Control and Prevention, ChinaJing Lu, University of Electronic Science and Technology of China, China
Copyright © 2021 Zhang, Hong, Zhang, Han, Xiao, Lu, Song, Yan, Wang, Zhu, Xu and Wu. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Yong Zhang, eW9uZ3poYW5nNzVAc2luYS5jb20=; Guizhen Wu, d2d6Y2RjQGhvdG1haWwuY29t