Characterization of Enteroviruses from Non-Human Primates in Cameroon Revealed Virus Types Widespread in Humans along with Candidate New Types and Species

Enteroviruses (EVs) infecting African Non-Human Primates (NHP) are still poorly documented. This study was designed to characterize the genetic diversity of EVs among captive and wild NHP in Cameroon and to compare this diversity with that found in humans. Stool specimens were collected in April 2008 in NHP housed in sanctuaries in Yaounde and neighborhoods. Moreover, stool specimens collected from wild NHP from June 2006 to October 2008 in the southern rain forest of Cameroon were considered. RNAs purified directly from stool samples were screened for EVs using a sensitive RT-nested PCR targeting the VP1 capsid coding gene whose nucleotide sequence was used for molecular typing. Captive chimpanzees (Pan troglodytes) and gorillas (Gorilla gorilla) were primarily infected by EV types already reported in humans in Cameroon and elsewhere: Coxsackievirus A13 and A24, Echovirus 15 and 29, and EV-B82. Moreover EV-A119, a novel virus type recently described in humans in central and west Africa, was also found in a captive Chimpanzee. EV-A76, which is a widespread virus in humans, was identified in wild chimpanzees, thus suggesting its adaptation and parallel circulation in human and NHP populations in Cameroon. Interestingly, some EVs harbored by wild NHP were genetically distinct from all existing types and were thus assigned as new types. One chimpanzee-derived virus was tentatively assigned as EV-J121 in the EV-J species. In addition, two EVs from wild monkeys provisionally registered as EV-122 and EV-123 were found to belong to a candidate new species. Overall, this study indicates that the genetic diversity of EVs among NHP is more important than previously known and could be the source of future new emerging human viral diseases.


Introduction
Enteroviruses (EVs), in the Picornaviridae family, form a diversified genus infecting numerous mammalian species including humans, Non-Human Primates (NHP), sheeps, cows and pigs [1,2]. Most EV infections are asymptomatic, but they can cause a wide range of diseases that can be severe and occasionally fatal in humans and animals [1,3,4]. Human EVs can cause symptoms ranging from mild febrile illness to severe forms, such as common cold, gastroenteritis, upper and lower respiratory diseases, acute hemorrhagic conjunctivitis, aseptic meningitis, myocarditis, encephalitis, myelitis and acute flaccid paralysis (AFP) [5]. An efficient surveillance system of AFP has been set up in the world thanks to the global poliomyelitis eradication program. Apart from AFP, other enteroviral diseases syndromes are rarely taken into consideration in the developing world including in sub-Saharan Africa despite the evidence of extensive EV circulation [6][7][8][9].
EVs are small non-enveloped viruses having a capsid with icosahedral symetry. Their genome is made up of a single polyadenylated positive RNA strand of about 7.5 kb, which is covalently linked to a small viral protein VPg at the 59 terminus. The single open-reading frame (ORF) is flanked by two untranslated regions (59UTR and 39UTR). The large polyprotein translated from the genomic RNA strand is cleaved into four structural proteins (VP1-VP4) and the non-structural proteins.
The Executive Committee of the International Committee on the Taxonomy of Viruses (ICTV) approved on February 2013 to rename the EV species without their host names. The Enterovirus Even though EVs are ubiquitous in humans, it is noteworthy that the worldwide distribution of EV types and species is not homogenous. Some human EV types and variants are more prevalent or geographically restricted in certain parts of the world. For instance, some EV-A such as EV-A71, Coxsackievirus A6 (CV-A6) and CV-A16 are involved in large epidemics primarily in Asia. EV-B species are seemingly predominant in temperate countries [10][11][12] even though the application of molecular characterization directly on clinical samples (without cell culture techniques) has recently uncovered several new EV-C types [13][14][15][16][17]. High rates of EV-C types have been repeatedly reported in tropical countries [7,8,18,19]. Our recent studies revealed an extensive circulation with a tremendous genetic diversity of EVs in humans living in Cameroon as well as in the Central African Republic. A remarkable number of human EVs belonged to phylogenetic lineages seemingly specific to sub-Saharan Africa: a new EV-A71 genogroup E, EV-C105, EV-B93, EV-D94, EV-D111, EV-D120 as well as several variants of the types CV-A13 and E-7 [8,14,18,20,21].
In contrast to human EVs which are extensively studied, simian EVs (NHP EVs) have been hardly considered. Most NHP EVs were primarily isolated in the 1950s to 1970s from primary cell cultures or clinical specimens originating from captive or wildcaught primates used in biomedical research [22][23][24][25][26][27][28][29]. The molecular characterization of 19 historical simian virus types thought to be EVs revealed that only 13 of them were actually members of the Enterovirus genus whereas the others defined a novel picornavirus genus termed Sapelovirus [30][31][32]. Seven of the confirmed simian EV types shared the same species with human-derived EV types: BaEV, SV19, SV26, SV35, SV43 and SV46 belonged to the EV-A species while SA5 belonged to the EV-B species. Six other simian EV types were more distantly related to those of human origin and were thus classified as separate species, SEV-A (A2plaque, SV4, and SV28) and EV-J (SV6, N125 and N203) respectively [30,33].
As other central African countries, Cameroon houses a wide range of Old World primates' species including apes (chimpanzees and gorillas) and at least 28 species of Old World monkeys. A recent study of EVs infecting wild apes in Cameroon identified 3 virus types including EV-A76 and EV-D111 [34] which also circulate among humans in central Africa [8,18,21]. However, the circulation and genetic diversity of EVs infecting African NHP is still poorly documented.
We further investigated the genetic diversity of EV types and species in the stool specimens originating from captive and wildliving apes and monkeys in Cameroon. We found virus types known to infect humans primarily in captive chimpanzees and gorillas but also in wild chimpanzees. In addition, some EVs from wild monkeys and chimpanzees were genetically distinct from all existing virus types, thus indicating that the overall picornavirus landscape is much more diversified than previously suspected.

Samples
Stool specimens of NHP were collected with the permission of the Cameroonian Ministries of Health, Environment and Forestry, and Research. Stool specimens of captive apes and monkeys were collected in April 2008, in a single day, at the Mvog-Beti zoo in Yaounde and the apes' sanctuary in Mfou (a rural village located at about 40 Km from Yaounde). Captive NHP were primarily wild-born orphaned animals living mainly in Southern rain forest areas of Cameroon. These animals have been brought to the zoos or sanctuaries, generally after their mother had been killed by hunters. About 5 g of stool specimens collected from freshly deposited stools were transported in refrigerated containers and kept frozen at 220uC until use. Further stool specimens originating from wild NHP were collected between June 2006 and October 2008 at forest sites located in the southern part of Cameroon. These specimens were part of a larger stool collection constituted for SIV screening [35,36]. About 2 g of stool specimen was collected from freshly deposited dung in a cryotube. These tubes were kept frozen in nitrogen containers for about 3 weeks and subsequently transported to the laboratory for storage at 2 20uC.
The preparation and clarification of stool suspensions at 20% (Wt/Vol) were performed as described elsewhere [37].

Fecal nucleic acids extraction
The fecal RNAs were extracted and purified from 140 mL of clarified stool suspensions using the QIAamp Viral RNA Mini Kit according to the manufacturer's instructions (Qiagen, France). RNAs were eluted in 60 mL of nuclease-free water.
Fecal DNAs were extracted from 220 mg of EV-positive stool samples using the QIAamp DNA Stool Mini kit (Qiagen, France), according to the manufacturer's instructions. DNAs were eluted with nuclease-free water in a final volume of 100 mL DNA solution.
Total RNAs and DNAs were stored at 280uC until use.
All RNA samples containing EVs detected with the VP1-specific amplification, including weak signals, insufficient for amplicon sequencing, were also subjected to RT-PCR targeting 59UTR and

Author Summary
Enteroviruses (EVs) are among the most common viruses infecting humans. The circulation and diversity of EVs infecting humans have been extensively studied in most parts of the world but those infecting Non-Human Primates (NHP) are still poorly documented. As most African countries, Cameroon houses a wide range of Old World primates' species including apes (chimpanzees and gorillas) and at least 28 species of Old World monkeys. In this study, we investigated the diversity of EV types and species from the stool specimens obtained from captive and wild-living apes and monkeys in Cameroon. We found that captive chimpanzees and gorillas were primarily infected by EV types having a worldwide distribution in humans, thus confirming that close contacts can promote cross-species transmission of EVs between humans and NHP. We also found proposed new EV types infecting wild NHP that were genetically distinct from all known EVs. Given the potential cross-species transmission of at least some EV types and the increasing contacts between human and NHP, EVs infecting NHP could serve as a source of new emerging human viral pathogens with unpredictable symptomatology.
3D pol regions. In the 59UTR region, a 450 base pairs DNA fragment was amplified with the primer pair UG52/UC53 [40]. RT-PCR products were then used as templates in semi-nested PCR using the PCR forward primer UG52 and a new reverse primer (AS3: 59-ACGGACACCCAAAGTAGTCGGT-39, nt 559 to 538 relative to SV46), designed from a highly conserved region in the alignment of all database available 59UTR sequences of simian EVs and representative human EVs. A portion of the 3D pol region was amplified with primers and reaction conditions previously reported [38,41].
In order to refine virus type identification, complete VP1 sequence was determined for a set of viruses selected as representatives of all lineages revealed by the analysis of partial VP1 sequences. Overlapping amplicons encompassing the VP3-2A region were amplified by PCR with a combination of consensus degenerate primers (AMTH, 224, 011 and GDCL) [39,41,42] and semi-nested PCR carried out with strain-specific forward and reverse primers.

Species determinations and search for prospective human-derived biological contaminations
For all stool specimens positive for EV detection and/or isolation, primate species were determined by mitochondrial DNA (mtDNA) analysis as described elsewhere [35,43,44].
In addition, in order to verify that EVs detected in NHP samples did not result from prospective contaminations by humanderived biological materials, fecal DNA from all EV-positive samples was screened for human DNA by using the primers pair FOXP2-humF/FOXP2-commonR [45]. This sensitive PCR can specifically amplify a portion of the forkhead box P2 (FOXP2) gene of the human DNA in a mixture containing DNA from other vertebrates, including chimpanzees and gorillas. Human DNA was extracted from a stool sample and used as positive control to validate that the sensitivity of the amplification was as previously reported [45].

Sequencing of amplicons and sequence analysis
Specific PCR products were either directly purified with QIAquick PCR Purification kit or gel-isolated and purified using the QIAquick Gel Extraction kit (Qiagen, France) following the manufacturer's protocol. Amplicons were subjected to direct sequencing using the BigDye terminator v3.1 kit (Applied Biosystems) and an ABI Prism 3140 automated sequencer (Applied Biosystems). They were sequenced in both directions using amplification primers.
Partial or complete VP1 capsid coding gene of each virus was pairwise compared with the prototype sequences of all known human and simian EV prototype strains retrieved from databases (http://www.picornaviridae.com/sequences/sequences.htm). Field isolates were considered to belong to a given type when they featured a nt identity $75% and an aa identity $88% with the prototype strain of this type [46].
Full-length VP1 sequences of viruses that were noticeably divergent from any known EVs (nt identity,75%), were submitted to the Picornaviridae Study Group of the International Committee for the Taxonomy of Viruses for assignment as a new type.
Multiple sequence alignments were performed with CLC Main Workbench 5.7.2 software (CLC bio, Aarhus, Denmark). Phylograms were inferred both by the distance and maximum likelihood methods implemented in MEGA, version 5.10 [47] and PhyML 3.0 [48] respectively. Distance based phylogenetic trees were reconstructed by the neighbor-joining method with the Kimura two-parameter method for computing evolutionary distances [49]. All alignment gaps were removed from the analysis for each Table 1. Global description of the primates samples and overall PCR, isolation and VP1 sequencing results. , Overall 12 VP1 amplicons showed very low signal insufficient for sequencing reactions. However, 11 of the corresponding viruses could be sequenced in the 59UTR region while the remaining one was refractory to 59UTR amplification but showed a divergent 3D pol sequence (GenBank nu KF648605 and KF648607).  Nucleotide sequences generated in this study were submitted to databases (accession numbers KF614472-KF614507 for partial 59UTR sequences, KF541625-KF541633 for partial VP1 sequences, KF541634-KF541649 and KF700264-KF700265 for full-length VP1 sequences, and KF648600-KF648622 for partial 3D pol sequences).

Cell lines and virus isolation
Stool samples from wild NHP that were positive for EV molecular detection and all samples from captive NHP were tested for virus isolation by cell cultures. A volume of 0.2 mL of chloroform-treated and clarified stool suspensions were inoculated onto monolayered Human rhabdomyosarcoma (RD), human larynx epidermoid carcinoma (HEp-2c), African Green monkey kidney (Vero), and rhesus monkey Kidney (LLC-MK2) cells. RD, HEp-2c and Vero cells were maintained at 36uC in D-MEM medium while LLC-MK2 cells were maintained at 36uC in the 199 medium. Both media were supplemented with 2% foetal calf serum and 2 mM L-glutamine. Infected tubes were microscopically checked for 7 days to detect the appearance of cytopathogenic effects (CPE). The supernatant of the cultures which had remained negative were frozen, thawed, and inoculated to fresh cell cultures. Then newly inoculated cultures were checked for the next 7 days. Virus isolates harvested from cell cultures were analyzed using the molecular approaches described above.
To increase the virus detection rate and to further characterize the viruses detected, all stool samples from captive NHP and stool samples from wild NHP that were EV-positive in molecular screening were inoculated to the cultures of RD, HEp-2c, Vero and LLC-MK2 cell lines. Virus isolation produced 7 isolates propagated on one or two of the cell lines used (Tables 1 and 2). All these isolates originated from captive NHP, five of which were negative using RT-snPCR targeting both 59UTR and VP1 genomic regions. The final rate of EV detection/isolation among captive NHP was thus 20.2% (20/99) ( Table 1).
In this study, the PCR based screening for human DNA was performed on all EV-positive samples regardless of the virus type identified. In contrast to the positive control (a human stool sample), none of the EV-positive stool samples, including those containing virus types previously reported in humans, showed a PCR-based evidence of contamination by human-derived tissues.

Molecular typing
Partial or complete VP1 sequence could be successfully determined for 27 EVs detected and/or isolated (Table 1). VP1 sequences determined from isolates Z057 and Z108-1 were identical to those obtained directly from viral RNA extracted from corresponding fecal samples. Overall VP1 amplicons from 12 samples showed very low signal insufficient for sequencing reactions. However, the viruses from 11 of these samples could be sequenced in the 59UTR region while the one from the remaining sample was refractory to 59UTR amplification and was shown to have a divergent 3D pol sequence (GenBank nu KF648607) (Table  S1). Apes derived enteroviruses. Of the 21 apes-derived EVs identified, 13 from chimpanzees and 5 from gorillas belonged to virus types previously reported in humans. These known human EVs belonged to the species EV-A (5 EV-A76 and 2 EV-A71), EV-B (2 E-15, 1 E-29 and 1 EV-B82), and EV-C (6 CV-A13 and 1 CV-A24) ( Table 2). Apart from E-15 and EV-B82, all these EV types have been previously identified among humans in our recent study in Cameroun [8]. For the first time, EV-A71, which is known to have a high potential to cause severe disease outbreak in humans [51,52], was found in captive NHP-derived samples. Interestingly EV-A76, a known human EV type, was found in wild chimpanzees, in accordance with a recent study reporting EV-A76 in a wild chimpanzee in Cameroon [34].
Besides the EV types well known in humans, recently discovered virus types as well as proposed new virus types were identified in chimpanzees. The full-length VP1 sequence of the virus Z046 (captive chimpanzee) enabled its unambiguous assignment as an EV-A119, an EV type recently discovered in a human in Cameroon [6]. The virus CHE20 displayed sequence similarities lower than 75% and 85% with all databases available field and prototype strains at nucleotide (nt) and aminoacid (aa) levels respectively ( Table 2). The Picornaviridae Study Group (PSG) of the International Committee for the Taxonomy of Viruses proposed the assignment of the virus CHE20 (wild chimpanzee) into a new type, designated EV-J121 (http://WWW.picornastudygroup.com/ types/enterovirus/ev_types.htm). The partial VP1 sequence of the unique virus from a wild gorilla GOJ01 was somewhat related to the species EV-A. However, it was quite divergent from all known human and simian EVs. Further analysis is required to draw a conclusion about the type assignment of this potential new EV-A.
Two samples contained a mixture of 2 viruses: CVA-13 and a potentially divergent EV (which could be sequenced only in the 3D pol region) originating from a captive chimpanzee sample (Z108); EV-A71 and EV-B82 originating from a captive gorilla sample (Z098) ( Table 2 and Table S1).
Monkeys-derived enteroviruses. Among the 6 monkey EVs that could be typed, the virus Z088 from a captive mandrill was an E-29, a virus type circulating among humans in Cameroon [8] and elsewhere. The five other monkeys' viruses were clearly distinct from EV types isolated in humans. The virus Z057, originating from a captive mandrill, was a simian picornavirus 17 whose prototype strain SA5 was originally recovered from a vervet monkey, native to Africa. All remaining 4 viruses from wild monkeys did not reach the sequence similarity thresholds (75% nt and 85% aa identities) allowing their classification into existing virus types. Indeed, related viruses RCMH5 and RCMH6 on one hand, RCMH3 and RCMH7 on the other hand, were distinct to all prototype strains and were thus registered as proposed new EV types, designated EV-122 and EV-123 respectively ( Table 2).

Phylogenetic relationships of the NHP enteroviruses
Since phylograms inferred using Maximum Likelihood (ML) and Neighbor Joining (NJ) algorithms displayed very similar patterns, only NJ trees were considered.
Some NHP viruses featured close genetic relationships with EVs which have been already shown to circulate in this region. This was exemplified by chimpanzees-derived EV-A76 which displayed a close relationship with another chimpanzee strain LM1677 recently detected in a Cameroonian wild chimpanzee [34]. Both viruses formed together a chimpanzee-specific group (bootstrap value of 100%) that was somewhat distinct to human-derived EV-A76 including those from Cameroon (CMR) and the Democratic Republic of Congo (COD) (Figure 1). This group also included other closely related EV-A76 (CHB2, CHB5 and CHB8) detected in chimpanzees' samples (not shown in figure 1 because only partial VP1 sequences were available for these viruses -see table 2). In the same manner, EV-A119 from a captive chimpanzee grouped reliably with the homotypic prototype recently identified in a human in Cameroon [6] (Figure 1). Two other EV-A119 fulllength VP1 sequences from a previous study originated from wild chimpanzees and a gorilla stool samples collected in Cameroon in 2010 [20] (Figure 1). Gorillas-derived EV-A71 fell into EV-A71 genogroup E, recently identified in sub-Saharan Africa [8,18,53] ( Figure 1). Similar clustering pattern was displayed by chimpanzees' CV-A13 ( Figure 2). All six partial and complete CV-A13 VP1 sequences defined a unique phylogenetic sub-cluster (98.5 to 100.0% nt identity between each other) gathering viruses originating from the same enclosure in the sanctuary of Mfou (data not shown). These viruses fell reliably within the CV-A13 cluster C gathering isolates recovered from humans in Cameroon and Madagascar from 2002 to 2009 [8,18,54] (Figure 2).
In contrast, a number of EVs from NHP did not show a clear relatedness with viruses already known to circulate in central and west Africa or elsewhere. For instance, the gorilla-derived CV-A24 was somewhat distinct from isolates originating from Cameroon and did not display peculiar relatedness with any of the known African lineages of CV-A24 [8] (Figure 2). Similarly, the chimpanzee and mandrill-derived E-29 were highly similar to each other but differed significantly from their counterparts originating from humans in Cameroon (Figure 3).
There were no E-15 and SA5 isolate identified so far in humans and NHP in the considered region. The E-15 isolate Z097 from a captive gorilla was quite distinct from homotypic strains infecting humans (Figure 3). The mandrill's simian picornavirus 17 Z057 was distinct from the wild chimpanzee EV-B110 from Cameroon and clearly belonged to the SA5 type originally described in a vervet monkey, native to Africa (Figure 3). This provided further support to the fact that at least some EV types could naturally infect different NHP species.
This study also identified hitherto unknown virus types that were phylogenetically distinct from all known human and animal EVs. In accordance with pairwise comparison with prototype strains, RCMH3/RCMH7 and RCMH5/RCMH6 from wild monkeys did not fit into any of the existing EV types and species. They defined two proposed new EV types: EV-122 and EV-123 ( Figure 4). With only 55.6 to 57.9% (nt) and 52.2 to 53.2% (aa) VP1 sequence identities to the closest EV type SV6 [30], EV-122 and EV-123 clearly defined a candidate new EV species. This proposed new EV species was interspersed between EV-A and EV-J (Figure 4). Compared to other EV species, EV-A, EV-J and the candidate new species constituted a separate branch supported Figure 2. Phylogenetic relationships among Enterovirus C strains from human and non-human primates. The NJ tree is based on the alignment of the full-length VP1 sequences. Studied viruses from apes are highlighted in bold red with host names indicated in brackets while human-derived isolates from Cameroon are specified in bold blue. Prototype strains are indicated by triangles (m). For clarity, type and lineage-specific clusters containing exclusively human isolates have been collapsed. The scale is shown at the bottom as substitutions per site. Viral isolates belonging to enterovirus species commented in the text are gathered in grey-shaded boxes. doi:10.1371/journal.pntd.0003052.g002 by a bootstrap value of 99% (Figure 4). The wild chimpanzee virus CHE20 was strikingly divergent compared to EVs originating from humans and chimpanzees. It grouped reliably with previously described types SV6, EV-J103 and EV-J108 (Figure 4). Within the EV-J cluster, the virus CHE20 defined a new type designated EV-J121. It was more closely related to EV-J103 POo-1 (nt & aa identity of 62.2% and 66.4%, respectively - Table 3) than to the other EV-J types (nt & aa identity ,58.7 & 58.8%, respectively). However, the similarity scores between EV-J121 CHE20 and EV-J103 POo-1 VP1 sequences were very heterogeneous along the whole VP1 genomic region ( Figure S1). In the 59 two thirds of the VP1 region (nt 1-591, according to CHE20 numbering), sequences featured high nt and aa identities of 65.8% and 66.4%, respectively. In addition these CHE20 and POo-1 sequences gathered in the same branch of a phylogenetic tree supported by a high bootstrap value (100% - Figure S1). In contrast sequences of the last one third of the VP1 region (nt 592-867) showed low nt and aa identities (54.3% and 55.4%, respectively) and fell in separate branches ( Figure S1). These findings suggested that EV-J121 VP1 region may have resulted from a recombination event between two distinct viruses.

Comparative analysis of non-structural genomic sequences of EVs from NHP and humans
To analyze the genetic relationships between non-structural genomic regions of EV from NHP, sequences of the 59UTR and 3D pol regions of most studied viruses were compared with those of human origin from Cameroon ( Figure 5).
Interestingly, the 59UTR and 3D pol sequences of the gorillasderived EV-A71 featured a close relatedness with the human EV-A71 isolate C08-146 belonging to the African genogroup that was recently reported in Cameroon [8] (Figure 5). This indicates that both human and gorilla EV-A71 have evolved from a recent common ancestor without apparent recombination events involving the other considered EV strains.
All EV-A76 sequences determined in this study were closely related to those of previously reported human and simian EV-A76 from Cameroon ( Figure 5). In addition, most EV-A76 sequences and in particular those of the 59UTR and 3D pol regions were closely related to those of the simian EV-A119 Z046 and to a  certain extent to those of human-derived EV-A89, EV-A90 and EV-A91. This suggested that these EV-A types have evolved through circulation in both humans and NHP and may be through recombination among each others.
All other 59UTR and 3D pol sequences of NHP EVs clustered apart from those of homotypic viruses of human origin. This result corroborated those of the VP1 region and argued in favor of a relative specificity of NHP-derived EVs compared to their homotypic human-derived isolates ( Figure 5). Similar results were found when the 59UTR and 3D pol sequences of NHP-derived CVA-13 and CVA-24 were compared to those of homotypic strains previously reported among humans in Cameroon [8] (data not shown).
Concerning the 3D pol region in particular, most of the sequences grouped into two main clusters. The cluster I (bootstrap value of 100%) only contained EV-A sequences while the cluster II (bootstrap value of 98%) contained only sequences of EV-B. This clustering pattern is in accordance with previous observations suggesting that in a majority of cases, recombination events in the non-structural part of the genome involve virus types belonging to the same species. Nonetheless, sequences of types BA13 (EV-A species), SA5 and EV-B110 (EV-B species) were found not to belong to these two clusters and were more closely related to sequences from other species. Interestingly, all these EV-A and -B sequences falling outside the two EV-A and EV-B-specific 3D pol clusters were of NHP origin. In particular, the 3D pol (and the 59UTR) sequences of some simian SA5 and EV-B110 (EV-B species) isolates grouped consistently with those of the novel NHPderived EV-122 and EV-123, suggesting evolution through interspecies recombination.

Discussion
The main goal of this study was to characterize the genetic diversity of EVs in both captive and wild NHP in Cameroon and to compare this diversity with that found in humans. Our findings revealed the co-circulation of some EVs in humans and NHP and uncovered proposed new EV types and species currently circulating in wild apes and monkeys.
All but one of the 16 EVs identified among captive NHP belonged to virus types that circulate among humans (Table 1), thus confirming that the NHP species are highly sensitive to human EVs and are at risk of infections with these pathogens. These results are similar to those documenting the presence of a wide range of human EVs in synanthropic NHP living in close association with humans in Bangladesh, 100% of the detected picornaviruses belonging to virus types reported in humans [55]. In contrast, in captive monkeys living in Bangladesh zoos only 15% of picornaviruses belonged to virus types reported in humans [56]. This last result may appear conflicting with our findings; however it could be explained by differences in captivity conditions for NHP in Bangladesh zoos and Cameroon sanctuaries. In Bangladesh, NHP were kept out of the way of humans in cages limiting physical contacts whereas in Cameroon most captive NHP were living in vast enclosures, with more easy interactions with employees and the public [56].
Among EVs of human origin detected in captive settings in Cameroon were Echoviruses and CV-A which have been shown to be worldwide distributed in humans. In particular CV-A13 and CV-A24 have been shown to have high endemicity and diversity among humans in central Africa [8,18]. Human to NHP transmission may have occurred as previously reported in another context, in Bangladesh. However, a relative sequence divergence was featured by VP1 gene of NHP derived E-29 and CV-A24 compared to their homotypic human counterparts from Cameroon. This suggests that some EV types known in humans may have been co-circulating and evolving independently for a while in NHP. However, this observation may be biased by the fact that some human and NHP EVs considered for this comparison originated from distant regions in Cameroon. At this stage, it cannot be excluded that captive NHP EVs belonging to virus types circulating in humans originated from wild NHP.
Despite the fact that EV-A76 was described recently [57], it is well distributed in humans with high rates reported in southeastern Asia and central Africa [55,[57][58][59]. Strikingly, both central Africa and southeastern Asia, where EV-A76 is highly prevalent in humans, are homes to several endemic NHP species. Strengthening previous observations [34], our study identified EV-A76 in wild chimpanzees. Even though their human origin cannot be wiped out at this stage, it is likely these EV-A76 were indigenous to wild chimpanzees since they originated from stool samples collected in the remote jungle with minimal human contacts.
The identification of EV-A119 and EV-A71 genogroup E in captive apes and their apparent association to sub-Saharan Africa [6,8,9,18] suggest their co-circulation in both human and NHP in this region. An apparent geographic association to sub-Saharan Africa has also been observed for EV-D111 which was originally found in a wild chimpanzee from Cameroon [34] and was subsequently shown to be distributed among humans in Central Africa [8,18,21]. These data suggest that NHP, especially apes, could serve as reservoirs of potential new EV types in humans.
Major discoveries of this study were that wild NHP were infected in most cases by typical simian EVs and that still unknown EV types and species were circulating in wild apes and monkeys. This indicated that either wild NHP are living out of the way of humans, with minimum physical contacts in Cameroon forests or, although less likely, that some EV types, species or genotypes are poorly adapted to humans. In addition, it is likely that the abundance and diversity of NHP EVs have been underestimated in this study, some EVs being refractory to PCR amplification using generic primers. In particular, wild gorilla showed a remarkably low EV detection rate at 0.25% (1/403). In accordance with this hypothesis, 5 samples of captive NHP that were PCR-negative in initial screening were shown to contain EVs that were successfully isolated in human cell cultures (RD and/or HEp-2c) (Table 1). Moreover, 12 samples yielding very low VP1 amplification signal could not be clearly identified. However, 7 of them showed 59UTR sequences closely related to those of known human or NHP-derived EVs (data not shown). The remaining 5 samples contained viruses that showed quite divergent 59UTR sequences (from a captive Cercopithecus neglectus and 3 wild chimpanzees) and 3D pol sequences (from a wild chimpanzee), suggesting unknown EVs or picornaviruses (data not shown). The fact that some studied viruses were refractory to viral isolation or PCR amplification in one of the target genomic regions provided further support to the fact that our approach did not take exhaustively into account the actual diversity of EVs in NHP. However, as other studies, this one has demonstrated the existence of known and hitherto unknown EV types and variants in NHP [20,34]. It is likely that hitherto unknown EVs remain to be discovered among NHP, in particular in the wild fauna. More powerful molecular approaches (like high throughput sequencing) will be helpful for an exhaustive assessment of the genetic diversity of EVs infecting NHP.
Although wild NHP were infected in most cases by typical simian EVs, it cannot be excluded that at times, these simian EVs could spread from NHP to humans living in adjacent habitats. Such transmission could happen directly from wild NHP to humans or could transit through the NHP populations kept in zoos and sanctuaries. These scenarios may have happened for a number of EVs discovered in humans and wild apes recently: EV-A76, EV-D111 and EV-A119 [6,8,18,21,57]. Wild NHP could play the role of reservoir and/or source of future emerging EVs with unpredictable symptomatology in humans. It is important to document the picornavirus landscape among NHP, especially in chimpanzees which are our closest primate relatives, in order to better predict the risks of their emergence and their pandemic potential in humans.
In Cameroon, wild chimpanzees have been shown to be infected with EVs related to those circulating in humans, suggesting possible cross-species transmission of EVs among primates (34). In the present study we showed that captive NHP of this country can be infected with well-known human EVs, as well as with typical simian EVs or EVs known to infect both humans and NHP. This supports the idea that zoos and sanctuaries may favor contacts and EV infection among humans and NHP. The finding of candidate new EV types and species as well as other potentially more divergent viruses revealed that the diversity of picornaviruses in NHP is broader than previously suspected. This argues for studying further EVs circulating in wild NHP and humans living in sympatry in order to estimate the exposure and potential dissemination of new EVs among humans in the forested regions in Cameroon. Figure S1 Comparative analysis of the full-length VP1 sequences of SV6, EV-J103, EV-J108, and the newly sequenced EV-J121 types. The similarity plots were generated using SimPlot version 3.5.1 with 100 nucleotide (nt) windows, 20 nt increments, and the Kimura 2-parameter method with a transition-transversion ratio of 8.0. Partial VP1 based trees depicting the relationships between the strains considered are presented at the bottom of the corresponding target portion.

(TIF)
Table S1 Summary of the detection and sequencing in the 59UTR, VP1 and 3D pol regions of enteroviruses among wild and captive non-human primates. (DOCX)