Next Article in Journal
Conserved Oligomeric Golgi (COG) Complex Proteins Facilitate Orthopoxvirus Entry, Fusion and Spread
Previous Article in Journal
Replication of a Dog-Origin H6N1 Influenza Virus in Cell Culture and Mice
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Viruses
Volume 12
Issue 7
10.3390/v12070706
viruses-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Recombination Analysis of Non-Poliovirus Members of the Enterovirus C Species: Restriction of Recombination Events to Members of the Same 3DPol Cluster

by
Lieke Brouwer
1,*,
Kimberley S.M. Benschop
2,
Dung Nguyen
3,
Everlyn Kamau
3,
Dasja Pajkrt
4,
Peter Simmonds
3,† and
Katja C. Wolthers
1,†
1
Department of Medical Microbiology, Academic Medical Center (AMC), Amsterdam University Medical Centers (Amsterdam UMC), 1105AZ Amsterdam, The Netherlands
2
Centre for Infectious Disease Control, National Institute for Public Health and the Environment (RIVM), 3721MA Bilthoven, The Netherlands
3
Nuffield Department of Medicine, University of Oxford, Oxford OX1 3SY, UK
4
Department of Pediatric Infectious Diseases, Academic Medical Center (AMC), Amsterdam University Medical Centers (Amsterdam UMC), 1105AZ Amsterdam, The Netherlands
*
Author to whom correspondence should be addressed.
†
These authors contributed equally to this work.
Viruses 2020, 12(7), 706; https://0-doi-org.brum.beds.ac.uk/10.3390/v12070706
Submission received: 17 June 2020 / Accepted: 23 June 2020 / Published: 30 June 2020
(This article belongs to the Section Animal Viruses)
Browse Figures
Versions Notes

Abstract

:
Enteroviruses (EVs) are highly prevalent viruses worldwide. Recombination is known to occur frequently in EVs belonging to species Enterovirus A, Enterovirus B, and Enterovirus C. Although many recombinant vaccine-derived poliovirus (VDPV) strains have been reported, our knowledge on recombination in non-polio EVs in the species Enterovirus C is limited. Here, we combined a dataset consisting of 11 newly generated full-length Enterovirus C sequences and 180 publicly available sequences to study recombination dynamics in non-polio EVs. To identify recombination patterns, maximum likelihood phylogenetic trees of different genomic regions were constructed, and segregation analyses were performed. Recombination was observed between members of the same 3DPol cluster, but was rarely observed between members of different clusters. We hypothesize that this restriction may have arisen through their different compartmentalization in respiratory and enteric tracts related to differences in cellular tropisms so that the opportunity to recombine may not be available.

1. Introduction

Enteroviruses (EVs) are single stranded positive-sense RNA viruses in the family of Picornaviridae. Upon infection, EVs can cause a wide variety of symptoms and disease outcomes, ranging from mild respiratory or gastro-intestinal symptoms to meningitis, encephalitis, and acute flaccid paralysis [1]. All currently known EVs belong to one of 15 species (Enterovirus A-Enterovirus L and Rhinovirus A-Rhinovirus C), with Enterovirus A–Enterovirus D and all three rhinovirus species known to infect humans [2]. The EV genome is approximately 7500 nucleotides (nts) in length. It encodes a polyprotein consisting of a P1, P2, and P3 region, flanked by untranslated regions (UTRs) at both the 5′ and 3′ end. The P1 region is further cleaved into 4 structural proteins that make up the viral capsid (VP1-VP4), while the P2 and P3 region are further cleaved into non-structural proteins (2A - 2C and 3A - 3D, respectively) required for virus replication.
Recombination is often observed within members of species Enterovirus A, Enterovirus B, and Enterovirus C, but recombination events are rare in Enterovirus D [3]. Recombination events in EVs are almost exclusively detected at the edges of the structural P1 region, or within the non-structural 5′UTR, P2 or P3 region [4,5]. Though recombination has been shown to occur across the whole genome, not all recombinants may be viable, and there would thus be selection of recombinants with breakpoints in these specific locations [6].
Recombination is usually demonstrated by showing phylogeny violations; dissimilar clustering of viral strains between evolutionary trees constructed from nt sequences of different regions of the genome [7,8,9,10,11,12,13,14]. Recombinant forms (RFs) can then be assigned based on clustering of the strains in the nonstructural region [9,10,12,13,14], and breakpoints can be detected by performing bootscanning [7,8,11]. Among members of Enterovirus A, recombination is common, although EV-A71 shows lower recombination rates compared to other types in this species [13,15,16]. Recombination rates in Enterovirus B are high [4] and a vast number of recombinant forms have been identified in this species [9,10,12]. Among members of Enterovirus C, three genetic subgroups have been described (Ca, Cb, and Cc) between which recombination appears to be restricted, although there is frequent recombination within these subgroups [7]. In the current study, for clarity, we will refer to these subgroups as subspecies C1, C2, and C3, respectively. Recombinant forms of the three serotypes of poliovirus in subspecies C3—and specifically their live attenuated vaccine strains—are frequently reported, and can contain fragments of coxsackievirus A11 (CVA11), CVA13, CVA17, and CVA20 strains, all members of subspecies C3 [17,18,19]. Outbreaks of poliomyelitis caused by these vaccine derived polioviruses (VDPVs) have been reported in several countries where the live attenuated vaccine is used [20,21,22]. Furthermore, it has been shown that in vitro constructed recombinants of PV and other subspecies C3 strains are viable in cell culture [23]. Though PVs were highly endemic worldwide, the other Enterovirus C types are mainly endemic in African countries whereas they are rarely detected in Europe, Asia, and North America [24,25,26,27,28,29,30,31,32]. As a result, only a small number of sequences are available for members of Enterovirus C, compared to Enterovirus A or Enterovirus B, and few recombination analyses have been conducted on this species.
In the current study, we used full length sequences obtained from Enterovirus C positive samples from a previous study [24] and publicly available Enterovirus C sequences to study recombination dynamics within this species.

2. Materials and Methods

2.1. Sample Processing and Full Length Sequencing

Stool samples positive for Enterovirus C (n = 41) by VP1 sequencing from a previous study [24] were selected for full length sequencing. For viral RNA extraction, the QIAamp Viral RNA Mini Kit (Qiagen, Hilden, Germany) was used according to manufacturer’s instructions. For all samples, 4 µL of extracted RNA was added to 1 µL of random primer (N9). The mixture was incubated at 65 °C for 5 min. 5 µL of Superscript First Strand mastermix (2 µL 5X first strand buffer, 0.75 µL water, 1.25 µL 0.1 M dNTP, 0.5 µL 0.1 M DTT, 0.5 µL Superscript III RT (Thermo Fischer, Waltham, USA)) was added and the samples were incubated at 42 °C for 60 min. 5 µL of Sequenase Mix 1 (1 µL 5x Sequanase buffer (Thermo Fischer), 3.85 µL water, 0.15 µL Sequenase Enzyme (Thermo Fischer)) was added and the samples were incubated at 37 °C for 8 min. 0.6 µL of Sequenase Mix 2 (0.45 µL Sequenase dilution buffer, 0.15 µL Sequenase enzyme) was added and the samples were incubated at 37 °C for 8 min. The obtained double stranded cDNA was purified using AMPure XP beads (Beckman Coulter, Brea, USA) according to manufacturers’ instructions. The DNA concentration was measured by Qubit (Thermo Fischer) (according to manufacturer’s instructions), and all samples were normalized to 0.2 ng/µL. For tagmentation and library preparation, the Nextera XT DNA Library Preparation Kit (Illumina, San Diego, USA) was used according to manufacturers’ instructions, and samples were purified using AMPure XP beads. The Qubit cDNA kit was used to quantify the libraries, and the size distribution was determined using the Agilent 2200 TapeStation Bioanalyzer. Libraries were normalized to 1.6 ng/µL (4 nM) with Tris-Cl 10mM pH 8.5 with 0.1% Tween. The normalized libraries were incubated at 95 °C for 5 min while shaking at 1800 rpm. 5 µL of each library was transferred to a single low bind tube. 15 µL of the DNA pool was added to 15 µL of 0.2M NaOH to denature the DNA into single strands. The denatured DNA pool was added to pre-chilled Hybridization Buffer (HT1) to obtain 600 µL of 16pM DNA. 6 µL of 12.5 pM Phix was added to 594 µL of the denatured pool for a 1% PhiX concentration. Sequencing was performed using the MiSeq Illumina cartridge and sequencer. The raw data were processed using GenomeDetective version 1.111 [33]. The obtained sequences were deposited in the DDBJ/GenBank/EMBL database under accession numbers MN914196-MN914206.

2.2. Database Construction

A database was constructed comprising all our study sequences that included at least the entire VP1 through 3DPol genomic regions (nt positions 2457–7334 based on reference strain CVA21 Kuykendall, accession number AF546702), and all unique non-polio Enterovirus C sequences of at least 7000 nts in length available in the DDBJ/GenBank/EMBL database (April 2019). Strains that were reported as clones and one EV-C99 strain (MH144606) that was found to contain sequencing artifacts were excluded. An Echovirus 1 strain (AF029859) was added as an outgroup for the construction of phylogenetic trees. A multiple sequence alignment was made using Mafft version 7.429 [34]. Information on year of detection, country of detection, and sample type (e.g., stool or respiratory sample) was extracted from GenBank for all sequences. If information on a sequence was not available in GenBank, it was retrieved from the original article.

2.3. Detecting Phylogenetic Violation

Phylogenetic trees of the VP1 (nt positions 2457–3350, reference genome AF546702), 2C (nt positions 4089–5072), and 3DPol (nt positions 5952–7334) genomic region were constructed using the maximum likelihood (ML) method with the general time reversible (GTR) + Gamma substitution model in RAxML version 8.2.10 [35]. The trees were visualized using the ggtree package version 1.16.4 [36] in R version 3.6 [37]. Log-likelihood values of the three topologies were evaluated for each of the three alignments (VP1, 2C and 3DPol) using the Shimodaira-Hasegawa method [38] implemented in RaxML version 8.2.10 [35].
Segregation analyses were performed using the TreeOrder scan function in SSE version 1.3 [35] using a window of 300 nts and steps of 24 nts. Segregation analyses were performed on the complete dataset for segregation into subspecies, and on each of the subspecies C1, C2, and C3 for segregation into types.

2.4. Plotting p-Distance and Determining Recombinant Forms

Pairwise p-distances between the VP1 and 3DPol regions of all sequences were extracted using SSE version 1.3. Histograms of the 3DPol pairwise p-distances were constructed in R version 3.6. Plots mapping VP1 and 3DPol p-distance for each sequence pair were constructed for each subspecies (C1, C2 and C3).
Clusters in the 3DPol genomic region were defined based on the pairwise p-distances in this region. Recombinant forms were assigned using a sequence p-distance threshold of 0.06–0.12 in the 3DPol region, determined through analysis of the distance distribution; these thresholds are comparable to those estimated in previous studies of EV-A71 and E30 [12,13].

3. Results

Sequences including the full VP1 through 3DPol sequence were obtained from 11 of 41 study strains and were included for further analyses. Based on their clustering in the phylogenetic tree of the VP1 region (Figure 1A), the study strains were typed as EV-C99 (03-1000, 04-4491), CVA13 (04-4444, 04-4517, 04-1438, 04-1450), CVA20 (03-4166, 03-4107, 04-1378, 03-4209), and CVA22 (03-4101). A total of 180 non-polio Enterovirus C sequences of more than 7000 base pairs were extracted from GenBank (April 2019). The complete database contained 191 sequences. Phylogenetic trees were constructed from the VP1, 2C, and 3DPol region. The Shimodaira-Hasegawa tests showed that VP1, 2C, and 3DPol alignments were significantly better explained by their respective topologies than by those of other genomic regions (p < 0.01), consistent with the occurrence of recombination events.

3.1. Segregation of Sequences in the VP1 Genomic Region

Based on their clustering in the VP1 coding region, Enterovirus C strains can be divided into three subspecies, which we named here C1 (CVA-1, -19,-22 and EV-C104, -C105, -C109, -C116, -C117 and -C118), C2 (CVA-21, -24 and EV-C96 and –C99) and C3 (CVA11, -13, -17, -20 and EV-C102) for clarity (previously named Ca, Cb, and Cc, respectively [7]). In addition to the clustering in these subspecies, all strains segregated into type-specific clusters in the VP1 coding region, except for the CVA22 strains and one EV-C116 strain, which clustered together (Figure 1A).

3.2. Segregation of Sequences in the 2C Genomic Region

In the phylogenetic tree of the 2C coding region (Figure 1B), six CVA13 strains clustered with subspecies C2, and two EV-C96 strains clustered within subspecies C1 (all reported before [7]), while all other strains remained in their subspecies-specific clusters. Segregation into type-specific clusters was violated for all types, except for EV-C105, EV-C109, EV-C118, EV-C104, and EV-C117 (in subgroup C1).

3.3. Segregation of Sequences in the 3DPol Genomic Region

Clusters were defined based on the pairwise p-distances in the 3DPol coding region in the complete dataset, which showed a bimodal distribution (Figure 2A). The two peaks represented intra- and inter cluster distances. A cut-off p-distance of 0.21 was set in between the two local maxima, separating the sequences into four distinct clusters, as shown in the phylogenetic tree (Figure 1C). This approach has been applied previously in defining different VP1-based clades in parechoviruses [39,40]. Strains belonging to subspecies C2 and C3 grouped together in cluster IV, while strains belonging to types CVA1, CVA19, CVA22, and EV-C116 in subspecies C1 and two EV-C96 strains formed cluster I. Strains belonging to types EV-C104 and EV-C117 formed cluster II while strains in EV-C105, EV-C109, and EV-C118 formed cluster III. The definition of these clusters was further supported by mapping the pairwise VP1 and 3DPol p-distances (Figure S1). Violation of type-specific clustering occurred across the phylogeny, although there were two clusters of only CVA21 and only CVA24 sequences, respectively.

3.4. Segregation of Sequences Across the Genome

Segregation of the sequences into subspecies was close to perfect in the structural P1 region of the genome. Segregation scores remained high throughout the P2 region but were lower in the 5′UTR and P3 region (Figure 3). Segregation into type-specific clusters in subspecies C1 decreased instantly at the start of the P2 region. In subspecies C2, segregation remained high throughout the P2 region but decreased in the P3 region. In subspecies C3, segregation vastly decreased at the start of the P2 region and dropped even further in the P3 region (Figure 3). These results are in line with the clustering of sequences in the phylogenetic trees constructed for the VP1, 2C, and 3DPol coding regions in the P1, P2, and P3 regions, respectively.

3.5. Recombinant Forms

P-distance cut-offs for defining recombinant forms (RFs) were based on the p-distance distribution in the 3DPol region (Figure 2B–E). The cut-off for cluster I, II, and III was set at a p-distance of 0.06. For cluster IV the cut-off was set at 0.12. A total of 89 RFs were identified in the current dataset. Most of the recombinant forms were represented by single strains, while 26 RFs each comprised multiple strains. The three RFs that occurred most frequently were in EV-C104/EV-C117, CVA21, and CVA24 strains (12, 14, and 25 strains with the same RF, respectively). One RF was shared by three CVA20 study strains (04-1378, 03-4166, 03-4107), and another was shared by two EV-C99 study strains (03-1000, 04-4491).

3.6. Clustering Associated with Site of Origin

The viral strains were categorized as derived from gastro-intestinal, respiratory, or ocular samples, based on information on sample type (e.g., stool sample or throat swab) available in GenBank. All of the strains with known sample type information in 3DPol clusters I and IV were derived from gastro-intestinal samples, except for a cluster of CVA21 strains, which were mainly of respiratory origin, and a cluster of CVA24 strains, which were mainly of ocular origin. These CVA21 and CVA24 clusters stayed monophyletic throughout the genome, except for one CVA-24 strain from Uganda, which was closest to a Nigerian CVA20 strain in its 3DPol region. All of the strains in 3DPol clusters II and III with sample type information were of respiratory origin, except for an EV-C105 strain, which was isolated from a stool sample [41] (Figure S2).

4. Discussion

In the current study, we used a set of full-length sequences extracted from the DDBJ/GenBank/EMBL database and newly generated sequences from study samples to investigate recombination dynamics within species Enterovirus C.
In line with other studies [8,17], we see high levels of phylogenetic violation in trees based on different regions of the genome, indicating frequent recombination within Enterovirus C. This is further supported by the low extent of segregation in the non-coding and non-structural regions of the genome. The sharp decreases in segregation at the borders of the P1 region are in line with previous reports and suggests the existence of favored recombination sites at these locations [4,7,42]. The decrease in segregation in the beginning of the P3 region implies another recombination hotspot at this location for subspecies C2 and C3. It has been shown before that recombination hotspots coincide with changes in sequence diversity; mainly the sequence diversity in the structural P1 region is higher than the diversity in the nonstructural 5′UTR, P2, P3, and 3′UTR regions [43].
Recombination seems to be restricted to occur within, but not between, several clusters, here named I, II, III, and IV. This is in line with reported Enterovirus C recombinants, which usually consist of parental genomes of types in the same cluster [5], although cluster I-II recombinants have been reported [41]. EV-C96 is the only type in our dataset that showed the ability to recombine with strains in two different clusters (cluster I and IV).
The distinct clusters within which recombination occurs may be separated due to differential tissue tropism. While viruses in subspecies C1, cluster I (C1/I) and viruses in C2/IV and C3/IV are known to infect the gastro-intestinal tract, viral strains in C1/II and C1/III are isolated almost exclusively from respiratory samples, a distinction that has been made previously for viruses in subspecies C1 [44]. This difference in tropism would imply that the groups of viruses never encounter each another, precluding recombination. The inability of the groups with a similar tissue tropism, such as C1/I and C2/IV or C3/IV, to recombine, could result from a higher sequence distance between their genomes, possibly leaving them unable to form viable viral strains upon recombination. Alternatively, these groups of viruses might infect different cell types within their target organs, thus remaining physically separated from one another.
Within cluster IV, there are two noticeable subclusters consisting of a single RF; one containing several CVA21 sequences (C2/IV), and one containing CVA24 sequences (C2/IV). The CVA24 subcluster contains almost exclusively strains isolated from conjunctival fluid in conjunctivitis patients, except for three gastro-intestinal strains from India. As has been suggested previously, these subclusters may contain a single RF due to their tissue tropism [7].
Although the genetic variety in the 3DPol region within each cluster is limited compared to the structural region of the genome, there are a vast number of RFs present, which is in line with previous reports [4,9]. The presence of different RFs in our study strains, despite some of the strains being closely related in the VP1 region, further illustrates the high rate of recombination.
The main limitation of our study is the amount of full-length sequences available, as typing is usually performed using only the (partial) VP1 sequence. As such, the available data may be biased towards specific types, geographical locations and time points originating from laboratories that conduct full length sequencing. Addition of more full-length sequences in the future that more accurately represent the diversity of circulating strains, may lead to different results and conclusions. It is therefore important that analyses on EV evolution and recombination are performed regularly and that full-length sequencing is encouraged [8].
In conclusion, using a dataset composed of Enterovirus C sequences, we have shown frequent recombination throughout the non-structural region of Enterovirus C strains, while the structural region seems to not be affected by recombination events, a phenomenon commonly seen for other enterovirus species. These recombination events are restricted to occur within, but not between, distinct clusters. We hypothesize that these clusters have arisen due to differences in tropism. As recombination plays an important role in EV evolution, understanding the recombination patterns of these viruses is essential in characterizing EV strains in outbreaks and to gain a deeper understanding of EV circulation. Furthermore, it provides us with valuable information for future studies on phenotypical characteristics and pathogenesis of EVs.

Supplementary Materials

The following are available online at https://0-www-mdpi-com.brum.beds.ac.uk/1999-4915/12/7/706/s1, Figure S1: Pairwise VP1 and 3DPol distances, Figure S2: Clustering by site of origin.

Author Contributions

Conceptualization, L.B., K.C.W., D.P. and P.S.; methodology, L.B., E.K., D.N., P.S., K.S.M.B.; writing—original draft preparation, L.B.; writing—review and editing, L.B.; All authors have read and agreed to the published version of the manuscript.

Funding

For the current study, L.B. was supported by grants from the AMC Young Talent Fund, Ketel 1 Study Fund, FEMS and KNAW. K.W. and D.P. were funded by EU H2020 ITN grant nr. 812673 (OrganoVIR).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Suresh, S.; Rawlinson, W.D.; Andrews, P.I.; Stelzer-Braid, S. Global epidemiology of nonpolio enteroviruses causing severe neurological complications: A systematic review and meta-analysis. Rev. Med. Virol. 2019, e2082. [Google Scholar] [CrossRef] [PubMed]
  2. Kong, K.L.; Lau, J.S.Y.; Goh, S.M.; Wilson, H.L.; Catton, M.; Korman, T.M. Myocarditis caused by human parechovirus in adult. Emerg. Infect Dis. 2017, 23, 1571–1573. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  3. Nikolaidis, M.; Mimouli, K.; Kyriakopoulou, Z.; Tsimpidis, M.; Tsakogiannis, D.; Markoulatos, P.; Amoutzias, G.D. Large-scale genomic analysis reveals recurrent patterns of intertypic recombination in human enteroviruses. Virology 2019, 526, 72–80. [Google Scholar] [CrossRef] [PubMed]
  4. Simmonds, P.; Welch, J. Frequency and dynamics of recombination within different species of human enteroviruses. J. Virol. 2006, 80, 483–493. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  5. Kyriakopoulou, Z.; Pliaka, V.; Amoutzias, G.D.; Markoulatos, P. Recombination among human non-polio enteroviruses: Implications for epidemiology and evolution. Virus Genes 2015, 50, 177–188. [Google Scholar] [CrossRef]
  6. Runckel, C.; Westesson, O.; Andino, R.; DeRisi, J.L. Identification and manipulation of the molecular determinants influencing poliovirus recombination. PLoS Pathog. 2013, 9, e1003164. [Google Scholar] [CrossRef]
  7. Smura, T.; Blomqvist, S.; Vuorinen, T.; Ivanova, O.; Samoilovich, E.; Al-Hello, H.; Savolainen-Kopra, C.; Hovi, T.; Roivainen, M. Recombination in the evolution of enterovirus C species sub-group that contains types CVA-21, CVA-24, EV-C95, EV-C96 and EV-C99. PLoS ONE 2014, 9, e94579. [Google Scholar] [CrossRef]
  8. Bessaud, M.; Joffret, M.L.; Holmblat, B.; Razafindratsimandresy, R.; Delpeyroux, F. Genetic relationship between cocirculating Human enteroviruses species C. PLoS ONE 2011, 6, e24823. [Google Scholar] [CrossRef]
  9. McWilliam Leitch, E.C.; Cabrerizo, M.; Cardosa, J.; Harvala, H.; Ivanova, O.E.; Kroes, A.C.M.; Lukashev, A.; Muir, P.; Odoom, J.; Roivainen, M.; et al. Evolutionary dynamics and temporal/geographical correlates of recombination in the human enterovirus echovirus types 9, 11, and 30. J. Virol. 2010, 84, 9292–9300. [Google Scholar] [CrossRef] [Green Version]
  10. Cabrerizo, M.; Trallero, G.; Simmonds, P. Recombination and evolutionary dynamics of human echovirus 6. J. Med. Virol. 2014, 86, 857–864. [Google Scholar] [CrossRef]
  11. Lukashev, A.N.; Lashkevich, V.A.; Ivanova, O.E.; Koroleva, G.A.; Hinkkanen, A.E.; Ilonen, J. Recombination in circulating human enterovirus B: Independent evolution of structural and non-structural genome regions. J. Gen. Virol. 2005, 86 Pt 12, 3281–3290. [Google Scholar] [CrossRef]
  12. McWilliam Leitch, E.C.; Bendig, J.; Cabrerizo, M.; Cardosa, J.; Hyypia, T.; Ivanova, O.E.; Kelly, A.; Kroes, A.C.M.; Lukashev, A.; MacAdam, A.; et al. Transmission networks and population turnover of echovirus 30. J. Virol. 2009, 83, 2109–2118. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  13. McWilliam Leitch, E.C.; Cabrerizo, M.; Cardosa, J.; Harvala, H.; Ivanova, O.E.; Koike, S.; Kroes, A.C.M.; Lukashev, A.; Perera, D.; Roivainen, M.; et al. The association of recombination events in the founding and emergence of subgenogroup evolutionary lineages of human enterovirus 71. J. Virol. 2012, 86, 2676–2685. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  14. Puenpa, J.; Vongpunsawad, S.; Osterback, R.; Waris, M.; Eriksson, E.; Albert, J.; Kroes, A.C.M.; Lukashev, A.; Perera, D.; Roivainen, M. Molecular epidemiology and the evolution of human coxsackievirus A6. J. Gen. Virol. 2016, 97, 3225–3231. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  15. Lukashev, A.N.; Shumilina, E.Y.; Belalov, I.S.; Ivanova, O.E.; Eremeeva, T.P.; Reznik, V.I.; Trotsenko, O.E.; Drexler, I.F.; Drostenet, C. Recombination strategies and evolutionary dynamics of the Human enterovirus A global gene pool. J. Gen. Virol. 2014, 95 Pt 4, 868–873. [Google Scholar] [CrossRef]
  16. van der Sanden, S.; van Eek, J.; Martin, D.P.; van der Avoort, H.; Vennema, H.; Koopmans, M. Detection of recombination breakpoints in the genomes of human enterovirus 71 strains isolated in the Netherlands in epidemic and non-epidemic years, 1963–2010. Infect Genet Evol. 2011, 11, 886–894. [Google Scholar] [CrossRef]
  17. Adeniji, J.A.; Faleye, T.O. Enterovirus, C. strains circulating in Nigeria and their contribution to the emergence of recombinant circulating vaccine-derived polioviruses. Arch. Virol. 2015, 160, 675–683. [Google Scholar] [CrossRef]
  18. Arita, M.; Zhu, S.L.; Yoshida, H.; Yoneyama, T.; Miyamura, T.; Shimizu, H. A Sabin 3-derived poliovirus recombinant contained a sequence homologous with indigenous human enterovirus species C in the viral polymerase coding region. J. Virol. 2005, 79, 12650–12657. [Google Scholar] [CrossRef] [Green Version]
  19. Combelas, N.; Holmblat, B.; Joffret, M.L.; Colbere-Garapin, F.; Delpeyroux, F. Recombination between poliovirus and coxsackie A viruses of species C: A model of viral genetic plasticity and emergence. Viruses 2011, 3, 1460–1484. [Google Scholar] [CrossRef] [Green Version]
  20. Mbaeyi, C.; Alleman, M.M.; Ehrhardt, D.; Wiesen, E.; Burns, C.C.; Liu, H.; Ewetola, R.; Seakamela, L.; Mdodo, R.; Ndoutabe, M.; et al. Update on vaccine-derived poliovirus outbreaks—Democratic Republic of the Congo and Horn of Africa, 2017–2018. MMWR Morb. Mortal. Wkly. Rep. 2019, 68, 225–230. [Google Scholar] [CrossRef]
  21. Mbaeyi, C.; Wadood, Z.M.; Moran, T.; Mjourn, A.F.; Stehling-Ariza, T.; Nikulin, J.; Al Safadi, M.; Iber, J.; Zomahoun, L.; Abourshaid, N.; et al. Strategic Response to an Outbreak of Circulating Vaccine-Derived Poliovirus Type 2—Syria, 2017–2018. MMWR Morb. Mortal. Wkly Rep. 2018, 67, 690–694. [Google Scholar] [CrossRef] [PubMed]
  22. Adamu, U.S.; Archer, W.R.; Braka, F.; Damisa, E.; Siddique, A.; Baig, S.; Higgins, J.; Etapelong Sume, G.; Banda, R.; Kipkoech Korir, C.; et al. Progress Toward Poliomyelitis Eradication—Nigeria, January 2018–May 2019. MMWR Morb. Mortal. Wkly. Rep. 2019, 68, 642–646. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  23. Bessaud, M.; Joffret, M.L.; Blondel, B.; Delpeyroux, F. Exchanges of genomic domains between poliovirus and other cocirculating species C enteroviruses reveal a high degree of plasticity. Sci. Rep. 2016, 6, 38831. [Google Scholar] [CrossRef] [PubMed]
  24. Brouwer, L.; Van Der Sanden, S.M.G.; Calis, J.C.J.; Bruning, A.H.L.; Wang, S.; Wildenbeest, J.G.; Rebers, S.P.H.; Phiri, K.S.; Westerhius, B.M.; Van Hensbroek, M.B.; et al. High frequency of Polio-like Enterovirus C strains with differential clustering of CVA-13 and EV-C99 subgenotypes in a cohort of Malawian children. Arch. Virol. 2018, 163, 2645–2653. [Google Scholar] [CrossRef] [Green Version]
  25. Cabrerizo, M.; Trallero, G.; Echevarria, J.E.; Moreno-Docon, A.; Pena, M.J.; Perez-Ruiz, M.; Avellón, A. Molecular characterization of enteroviruses associated with neurological infections in Spain, 2008. J. Med. Virol. 2013, 85, 1975–1977. [Google Scholar] [CrossRef]
  26. Bubba, L.; Martinelli, M.; Pellegrinelli, L.; Primache, V.; Tanzi, E.; Pariani, E.; Sandro, B. A 4-year Study on Epidemiologic and Molecular Characteristics of Human Parechoviruses and Enteroviruses Circulating in Children Younger Than 5 Years in Northern Italy. Pediatr. Infect. Dis. J. 2017, 36, 13–19. [Google Scholar] [CrossRef]
  27. Cristanziano, V.D.; Bottcher, S.; Diedrich, S.; Timmen-Wego, M.; Knops, E.; Lubke, N.; Kaiser, R.; Pfister, H.; Kabore, Y.; D’Alfonso, R. Detection and characterization of enteroviruses and parechoviruses in healthy people living in the South of Cote d’Ivoire. J. Clin. Virol. 2015, 71, 40–43. [Google Scholar] [CrossRef]
  28. Abedi, G.R.; Watson, J.T.; Nix, W.A.; Oberste, M.S.; Gerber, S.I. Enterovirus and Parechovirus Surveillance—United States, 2014–2016. MMWR Morb. Mortal Wkly. Rep. 2018, 67, 515–518. [Google Scholar] [CrossRef]
  29. Abedi, G.R.; Watson, J.T.; Pham, H.; Nix, W.A.; Oberste, M.S.; Gerber, S.I. Enterovirus and human parechovirus surveillance—United States, 2009–2013. MMWR Morb. Mortal. Wkly. Rep. 2015, 64, 940–943. [Google Scholar] [CrossRef]
  30. Jeong, E.J.; Lee, J.H.; Kim, M.S.; Bae, G.R.; Jung, C.; Lee, K.; Choi, S.-M.; Kim, D.-K.; Lee, D.-S.; Kim, W.-D.; et al. Molecular characterization of enteroviruses detected in Gyeong-Ju and Po-Hang provinces of Korea in 2003. Arch. Virol. 2010, 155, 1707–1712. [Google Scholar] [CrossRef]
  31. Kowada, K.; Takeuchi, K.; Hirano, E.; Toho, M.; Sada, K. Development of a multiplex real-time PCR assay for detection of human enteric viruses other than norovirus using samples collected from gastroenteritis patients in Fukui Prefecture, Japan. J. Med. Virol. 2018, 90, 67–75. [Google Scholar] [CrossRef] [PubMed]
  32. Kumar, A.; Shukla, D.; Srivastava, S.; Idris, M.Z.; Dhole, T.N. High frequency of enterovirus serotype circulation in a densely populated area of India. J. Infect. Dev. Ctries. 2013, 7, 475–483. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  33. Vilsker, M.; Moosa, Y.; Nooij, S.; Fonseca, V.; Ghysens, Y.; Dumon, K.; Pauwels, R.; Alcantra, L.C.; Eynden, E.V.; Vandamme, A.-M.; et al. Genome Detective: An automated system for virus identification from high-throughput sequencing data. Bioinformatics 2019, 35, 871–873. [Google Scholar] [CrossRef] [Green Version]
  34. Katoh, K.; Standley, D.M. MAFFT multiple sequence alignment software version 7: Improvements in performance and usability. Mol. Biol Evol. 2013, 30, 772–780. [Google Scholar] [CrossRef] [Green Version]
  35. Stamatakis, A. RAxML version 8: A tool for phylogenetic analysis and post-analysis of large phylogenies. Bioinformatics 2014, 30, 1312–1313. [Google Scholar] [CrossRef] [PubMed]
  36. Yu, G.; Lam, T.T.; Zhu, H.; Guan, Y. Two methods for mapping and visualizing associated data on phylogeny using Ggtree. Mol. Biol. Evol. 2018, 35, 3041–3043. [Google Scholar] [CrossRef]
  37. Team RC. R: A Language and Environment for Statistical Computing. Vienna, Austria, 2019. Available online: https://www.R-project.org/ (accessed on 30 April 2020).
  38. Shimodaira, H.; Hasegawa, M. Multiple comparisons of log-likelihoods with applications to phylogenetic inference. Mol. Biol. Evol. 1999, 16, 114. [Google Scholar] [CrossRef] [Green Version]
  39. Benschop, K.S.; de Vries, M.; Minnaar, R.P.; Stanway, G.; van der Hoek, L.; Wolthers, K.C.; Simmond, P. Comprehensive full-length sequence analyses of human parechoviruses: Diversity and recombination. J. Gen. Virol. 2010, 91 Pt 1, 145–154. [Google Scholar] [CrossRef]
  40. Calvert, J.; Chieochansin, T.; Benschop, K.S.; McWilliam Leitch, E.C.; Drexler, J.F.; Grywna, K.; da Costa Ribeiro, H., Jr.; Drosten, C.; Harvala, H.; Poovorawan, Y.; et al. Recombination dynamics of human parechoviruses: Investigation of type-specific differences in frequency and epidemiological correlates. J. Gen Virol. 2010, 91 Pt 5, 1229–1238. [Google Scholar] [CrossRef]
  41. Lukashev, A.N.; Drexler, J.F.; Kotova, V.O.; Amjaga, E.N.; Reznik, V.I.; Gmyl, A.P.; da Costa Riberiro, H.J.; Drosten, C.; Harvala, H.; Poovorawan, Y.; et al. Novel Serotypes 105 and 116 are members of distinct subgroups of human enterovirus C. J. Gen Virol. 2012, 93 Pt 11, 2357–2362. [Google Scholar] [CrossRef]
  42. Simmonds, P. Recombination and selection in the evolution of picornaviruses and other Mammalian positive-stranded RNA viruses. J. Virol. 2006, 80, 11124–11140. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  43. Simmonds, P. SSE: A nucleotide and amino acid sequence analysis platform. BMC Res. Notes. 2012, 5, 50. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  44. Tokarz, R.; Haq, S.; Sameroff, S.; Howie, S.R.; Lipkin, W.I. Genomic analysis of coxsackieviruses A1, A19, A22, enteroviruses 113 and 104: Viruses representing two clades with distinct tropism within enterovirus C. J. Gen. Virol. 2013, 94 Pt 9, 1995–2004. [Google Scholar] [CrossRef]
Viruses 12 00706 g001 550
Figure 1. Phylogenetic trees (Maximum Likelihood method with GTR-Gamma nt substitution model) of the VP1 (A), 2C (B) and 3DPol (C) genomic regions containing all study strains and Enterovirus C strains retrieved from GenBank. GenBank extracted strains are named as type/accession number/two-letter country code of country of detection/ year of detection. Study strains are named with their studynumber including the year of isolation (P02-XXXX, 2002; P03-XXXX, 2003; P04-XXXX, 2004). The subspecies and subgroups (C1, C2, C3) are marked in the VP1 tree and strain names are colored accordingly. Tips are colored according to type. The identified clusters (I, II, III, IV) are marked in the 3DPol tree. Bootstrap values >70% are shown. The length of the branch of the outgroup (Echovirus 1/AF029859/1997/US, branch shown with dotted line) was adjusted to fit the figure margins.
Figure 1. Phylogenetic trees (Maximum Likelihood method with GTR-Gamma nt substitution model) of the VP1 (A), 2C (B) and 3DPol (C) genomic regions containing all study strains and Enterovirus C strains retrieved from GenBank. GenBank extracted strains are named as type/accession number/two-letter country code of country of detection/ year of detection. Study strains are named with their studynumber including the year of isolation (P02-XXXX, 2002; P03-XXXX, 2003; P04-XXXX, 2004). The subspecies and subgroups (C1, C2, C3) are marked in the VP1 tree and strain names are colored accordingly. Tips are colored according to type. The identified clusters (I, II, III, IV) are marked in the 3DPol tree. Bootstrap values >70% are shown. The length of the branch of the outgroup (Echovirus 1/AF029859/1997/US, branch shown with dotted line) was adjusted to fit the figure margins.
Viruses 12 00706 g001
Viruses 12 00706 g002 550
Figure 2. Histograms of pairwise p-distances in the 3DPol region of (A) All sequences and (B,C,D,E) Sequences within cluster I, II, III, and IV, respectively.
Figure 2. Histograms of pairwise p-distances in the 3DPol region of (A) All sequences and (B,C,D,E) Sequences within cluster I, II, III, and IV, respectively.
Viruses 12 00706 g002
Viruses 12 00706 g003 550
Figure 3. Segregation plots of all sequences, showing the segregation into subspecies (in black) and of sequences belonging to subspecies C1, C2, and C3 (green, red, and blue, respectively) showing the segregation into types. The position in the genome is shown on the x-axis. The dispersion value is shown on the right x-axis. The percentage of segregation (inverse of the dispersion) is shown on the right y-axis.
Figure 3. Segregation plots of all sequences, showing the segregation into subspecies (in black) and of sequences belonging to subspecies C1, C2, and C3 (green, red, and blue, respectively) showing the segregation into types. The position in the genome is shown on the x-axis. The dispersion value is shown on the right x-axis. The percentage of segregation (inverse of the dispersion) is shown on the right y-axis.
Viruses 12 00706 g003

Share and Cite

MDPI and ACS Style

Brouwer, L.; Benschop, K.S.M.; Nguyen, D.; Kamau, E.; Pajkrt, D.; Simmonds, P.; Wolthers, K.C. Recombination Analysis of Non-Poliovirus Members of the Enterovirus C Species: Restriction of Recombination Events to Members of the Same 3DPol Cluster. Viruses 2020, 12, 706. https://0-doi-org.brum.beds.ac.uk/10.3390/v12070706

AMA Style

Brouwer L, Benschop KSM, Nguyen D, Kamau E, Pajkrt D, Simmonds P, Wolthers KC. Recombination Analysis of Non-Poliovirus Members of the Enterovirus C Species: Restriction of Recombination Events to Members of the Same 3DPol Cluster. Viruses. 2020; 12(7):706. https://0-doi-org.brum.beds.ac.uk/10.3390/v12070706

Chicago/Turabian Style

Brouwer, Lieke, Kimberley S.M. Benschop, Dung Nguyen, Everlyn Kamau, Dasja Pajkrt, Peter Simmonds, and Katja C. Wolthers. 2020. "Recombination Analysis of Non-Poliovirus Members of the Enterovirus C Species: Restriction of Recombination Events to Members of the Same 3DPol Cluster" Viruses 12, no. 7: 706. https://0-doi-org.brum.beds.ac.uk/10.3390/v12070706

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop