Arbovirus is an ecological term applied to viruses transmitted by arthropods, particularly mosquitoes and ticks 1. Most arboviruses are maintained in enzootic cycles and are transmitted between mosquitos and vertebrates. Humans become infected and may develop diseases ranging from subclinical or mild infections to systemic febrile illness, hemorrhagic fever, and meningoencephalitis 1. Examples of arboviral diseases include yellow fever, Zika, and dengue fever, which have reemerged or increased their geographic range causing public health emergencies and urging governments to enhance their surveillance programs and research 2. Around 534 arboviruses have been described, and at least 135 of them are of medical concern 1,3. Almost all arboviruses have RNA genomes and are classified into nine families and nine genera 1,4; arbovirus species of current medical importance in South America belong to Flavivirus (Flaviviridae), Alphavirus (Togaviridae), and Orthobunyavirus (Perybunyaviridae) genera 4.
In Ecuador, three flaviviruses (yellow fever virus, YFV, dengue virus, DENV, and Zika virus, ZIKV) and two alphaviruses (Venezuelan equine encephalitis virus, VEEV, and Chikungunya virus, CHIKV) have caused epidemic outbreaks 5-9. YFV has two well-recognized transmission cycles: urban and sylvatic 10. Recognized urban epidemics occurred in 1842, 1853-1856, 1867-1869, 1877-1878 and 1880 5. Urban YF was eradicated in Ecuador in 1919 11, but sylvatic YF has been documented since 1949 in Amazonian basin provinces 12,13. Vaccination campaigns against YFV have been conducted in Ecuador since 1944, particularly in the Amazonian and coastal provinces. DENV caused the first epidemic in 1988 6 and currently, it is the most commonly detected arbovirus 14. ZIKV and CHIKV are more recent introductions; the first epidemics by these viruses occurred in 2014 and 2016, respectively. VEEV has caused two major epidemics, one in 1969 and the other in 1972; as usual for this virus, epidemic outbreaks were preceded by an epizootic event in equines 7,15. According to the Centers for Disease Control and Prevention (CDC), the 1969 outbreak caused 20,000 cases and 200 fatalities in Ecuador 16. Subsequently, a major outbreak of VEEV was reported in Guatemala, Central America, northern México, and Texas in the United States of America. The introduction of the epidemic subtype from Ecuador to Guatemala was speculated, but no convincing mechanism was identified 17. Later, it was shown that the different outbreaks originated from a residual live virus in incompletely formalin-inactivated vaccines 18,19.
Other arboviruses of medical concern have also been reported in Ecuador. They include the Saint Louis encephalitis virus (SLEV), the West Nile virus (WNV), and the Ilheus virus (ILHV) belonging to the Flavivirus genus; the eastern equine encephalitis virus (EEEV), the western equine encephalitis virus (WEEV), and the Mayaro virus (MAYV) from the Alphavirus genus, and the Oropouche virus (OROV) belonging to the Orthobunyavirus genus 20-24. SLEV, EEEV, and WEEV were isolated from sentinel hamsters and mosquitoes between 1974 and 1978 20. Then, SLEV as well as ILHV, OROV, MAYV, and VEEV were reported after the detection of specific antibodies (Abs) in febrile humans sampled in the Ecuadorian Amazonia in 1997 and between 2001 and 2004 23-25. ILHV and OROV were later isolated from febrile patients of the Ecuadorian Coast and the Amazonia, respectively 22,24. WNV and EEEV have been reported by the detection of specific Abs in horses of the Ecuadorean Coast: WNV from healthy equines sampled between 2007 and 2009 and EEEV from two equines showing neurological signs sampled in 2013 21,26. Eight other enzootic arboviruses have been isolated in Ecuador and named according to the local geographic region where they were found, but little is known about their medical importance 20. One example is the “Playas virus”, of which two strains (75V3066 and 75V938) are classified as Cache Valley virus and one (75V5758) as Maguari virus, both belonging to the genus Orthobunyavirus and associated with human disease 27.
Serosurveys are a key tool for determining the proportion of individuals exposed and of those with immunity to infectious agents, as well as for the detection of subclinical infections 28. In Ecuador, little is known about these aspects of arboviruses. The surveillance activities of the Ministerio de Salud Pública - MSP (Ministry of Public Health) are focused on the diagnosis of DENV, YFV, CHIKV, and ZIKV, mainly using nucleic acid detection; other arboviruses are surveyed using generic primers, a method with intrinsic limitations that may lead to missing arboviruses of medical concern. In this study, we aimed to investigate the proportion of three Ecuadorian populations with IgG Abs (past exposure/immunity) against selected flaviviruses and alphaviruses. This population included conscripts of a military fort of the Ecuadorian Amazonia recently vaccinated against YFV with the “YFV-17D” vaccine. We also determined the percentage of IgM Abs (recent exposure) against DENV and WNV in febrile and healthy individuals to look for clinical and subclinical infections. Finally, we emphasize the importance of widening and strengthening arbovirus surveillance in Ecuador.
Materials and methods
Ethical statement
This project was approved by the institutional review board of the Instituto Nacional de Higiene y Medicina Tropical “Dr. Leopoldo Izquieta Pérez” (INHMTLI) now Instituto Nacional de Investigación en Salud Pública. Human samples were collected under verbal informed consent. Domestic animals were sampled given owner consent and supervision. Hamsters were used as sentinel species in the field and suckling mice “Balb/C4” were required for inoculation of hamster tissue and mosquito pools. Hamsters and mice were grown and housed at the laboratory animal facility of INHMTLI, which has adequate infrastructure to guarantee animal welfare and dedicated personnel for cage cleaning and animal feeding. After exposure/inoculation, the animals were euthanized using chloroform as an inhaled anesthetic. All these procedures were performed in strict accordance with the institutional guidelines for the care and use of laboratory animals. We followed the ethical guidelines of the Declaration of Helsinki and the “Public Health Service Policy on Humane Care and Use of Laboratory Animals” published by the National Institute of Health of the United States of America.
Sampling and sampling location
This study comprised serum samples of humans and equines for arbovirus seroprevalence analysis, serum samples of hamsters for arbovirus seroconversion analysis, mosquito pools, and hamster tissues for viral isolation. Sampling was performed from 2009 to 2011 at three locations in Ecuador, two in the coastal region (Santa Elena Province-Manglaralto and Los Ríos Province-Vinces) and one in the Amazonia basin (Pastaza Province- Puyo) (figure 1). Manglaralto and Vinces were chosen for the presence of marshes and wetlands that contribute to the presence of mosquitoes. Puyo is in the Amazon rainforest, where many arboviruses circulate in their natural cycles. Hamsters were exposed as sentinel species at each sampling location for one week. During the study period, each location was visited six times, more or less every three months. The number of human samples collected by location is shown in figure 1. Detailed information on the number of human samples, sentinel hamsters, equines, and mosquitoes collected by location and date is shown in table S1.
The procedure for sampling humans and domestic animals is described by Beaty, et al. 29. The human population sampled was expected to have been vaccinated against YFV in different immunization campaigns. In one location (Puyo), we sampled conscripts of the Jungle Brigade 17 Pastaza (Amazonas Military Fort), who had all recently been vaccinated with YFV- 17D. Two different groups of conscripts were sampled three times during the year of their stay in the military fort to evaluate IgG Abs against YFV and their seroprevalence to other flaviviruses.
Mosquitoes were trapped with CDC light traps supplemented with CO230 and then frozen in dry ice and transported to the laboratory.
Laboratory analysis
Seroprevalence/seroconversion to arbovirus was evaluated mainly through the analysis of IgG Abs. For humans, each sample was analyzed against a battery of arboviruses comprising DENV, YFV, SLEV, WNV, EEEV, WEEV, and VEEV. Only monotypic results were considered positive for the corresponding virus. IgG Abs against DENV and WNV were evaluated by commercial ELISA methods of PANBIO (Windsor, Queensland, Australia, specificity 81.6 - 97.2% and sensitivity 73.9 - 85.1%) and FOCUS (Cypress, CA, USA, specificity 97% and sensitivity 36%), respectively. IgG Abs against YFV, SLEV, EEEV, WEEV, and VEEV were evaluated with in-house ELISAs following Johnson, et al.31. Sucrose-acetone-extracted suckling mouse brain antigens and positive serum against these viruses were donated by the CDC’s National Center of Infectious Diseases, Division of Vector-Borne Infectious Diseases (Fort Collins, CO, United States of America). IgM Abs against WNV and DENV were also analyzed in individuals who were febrile at the time of sampling or reported to be febrile one or two weeks before sampling (table S2). Since arboviruses also cause subclinical infections, IgM Abs against these viruses were further analyzed in a subset of nonfebrile individuals using commercial ELISAs of PANBIO (specificity: 85.4 - 98.9%; sensitivity: 46.6 - 64.7%) and FOCUS (specificity: 100%; sensitivity: 100%), as described above. Positive and negative serum samples were included for the respective ELISA test. The IBM SPSS Statistics package for Windows, version XX (IBM Corp., Armonk, N.Y., USA), was used for frequency analysis of monotypic samples according to individual age.
Equine and hamster IgG Abs were evaluated with the hemagglutination inhibition (HI) test on kaolin-adsorbed sera according to the technique of Clarke, et al.32, adapted to microtiter plates and following the procedures described in the document Técnicas de laboratorio para el diagnóstico y la caracterización de los virus del dengue33. The test was conducted for EEEV, WEEV, and VEEV, and nine hamsters were also evaluated for YFV. We used the same antigens donated by the CDC. A shortage of antigens did not allow testing other arboviruses or more hamsters for YFV. An HI titer of 1:20 or higher inhibiting hemagglutination produced by four antigen units was considered positive. Hamsters 5-7 (table S1) were analyzed before and after their exposure at the sampling site 34. After exposure, hamsters were euthanized and their tissues (brain, heart, lungs, kidneys, liver, spleen, and skeletal muscle) were processed for IC inoculation in suckling mice Balb/C4 and analyzed by RT-PCR.
Mosquitoes were identified, classified, and grouped by species. Spearman correlation was carried out using SPSS to test the relationship between the number of mosquitoes captured by species and location and the proportion of human samples showing IgG Abs reactive against the arboviruses tested. This nonparametric test was chosen given that the proportion of samples having Abs was calculated on categorical data (yes/no). Mosquitoes pooled by species in groups of 50 were crushed in a sterile environment and processed for intracerebral (IC) inoculation in suckling mice Balb/C4 and RT-PCR analysis.
RT-PCR for flaviviruses and alphaviruses was conducted on samples from febrile humans, hamster tissues, and mosquito pools according to Ayers, et al.35 and Sánchez-Seco, et al.36, respectively.
Results
Analysis of IgG Abs against selected arboviruses was performed on 1,842 human sera. We found Abs for all the viruses investigated: DENV, YFV, SLEV, WNV, EEEV, WEEV, and VEEV. The percentages of Abs against different arboviruses in 1,263 samples (without including conscripts) are shown in Table 1. The flaviviruses showed the highest seroprevalence with WNV (63.61%), followed by YFV (62.67%), SLEV (52.29%), and DENV (20.63%); the alphaviruses showed the lowest prevalence, particularly WEEV (4.67%). Some samples reacted with only one of the viruses tested (monotypic samples): WNV: 58, SLEV: 20, EEEV: 54, VEEV: 64, and WEEV: 18. The frequency of monotypic samples according to the age of the patient is shown in figures S1-S5. Individuals between 2 and 82 years of age showed monotypic reactions.
Virus | Percentage | (Positivies/analyzed) |
---|---|---|
DENV | 20,63 | (260/1,260) |
YFV | 62,67 | (769/1,227) |
SLEV | 51,29 | (635/1,238)* |
WNV | 63,61 | (778/1,223)* |
EEEV | 8,63 | (109/1,263)* |
WEEV | 4,67 | (59/1,263)* |
VEEV | 7,60 | (96/1,263)* |
DENV: Dengue virus, YFV: Yellow fever Virus, SLEV: Saint Louis encephalitis Virus, WNV: West Nile encephalitis Virus, EEEV: Eastern equine encephalitis virus, WEEV: Western equine encephalitis virus, VEEV: Venezuelan equine encephalitis virus
* There were monotipic results as follows: WNV 69 (34 Vinces, 21 Manglaralto, and 14 Puyo); SLEV 20 (5 Vinces, 2 Manglaralto, 13 Puyo); EEEV (11 Vinces, 39 Manglaralto, 4 Puyo); WEEV (3 Vinces, 10 Manglaralto, 5 Puyo); VEEV (25 Vinces, 23 Manglaralto, 16 Puyo).
The results of IgG analysis by location are shown in figure 2. A total of 474 samples were studied in Vinces; the highest seroprevalence was found for WNV, 63.59% (276/434), and the lowest was for WEEV, 4.85% (23/474). In Manglaralto, we screened 496 samples. The highest seroprevalence was found for WNV, 87.7% (435/496), and the lowest for WEEV, 5.84% (29/496). In Puyo, we analyzed 872 samples comprising 293 samples from the civilian population and 579 from conscripts. For the civilian population, the highest seroprevalence was found for YFV, 52.75% (153/290), and the lowest for WEEV, 2.38% (7/293). Puyo is one of the Amazonian provinces subject to regular YFV vaccination.
The percentage of conscripts with IgG Abs against arbovirus is shown in figure 3a and b. Two groups (group 1: 61 individuals and group 2: 67 individuals) recently vaccinated against YFV were sampled three times. For group 1, the percentage of response to the vaccine was 85.05% and between 23.15% and 37.84% showed Abs against DENV, WNV, and SLEV. In group 2, 81.31% showed Abs against YFV, and between 5.90% and 42.74% showed Abs against DENV, WNV, and SLEV.
A total of 102 individuals were registered as febrile and analyzed for IgM Abs against WNV and DENV, and the detailed results are shown in Table S2. Four individuals from Puyo (30.76%-4/13) showed IgM Abs reactive against DENV, one in the civilian population, and three in the second group of conscripts (two were positive in the second sampling and the other in the third sampling). Two samples showed monotypic IgM Abs-reactive against WNV: one from a 19-year-old female in Vinces (1.66%-1/60) and another from a 53- year old woman in Manglaralto (5.55% -1/18).
A total of 1,343 samples of nonfebrile individuals (including conscripts) were analyzed for IgM Abs against DENV of which 5.06% (18/356) in Vinces, 2.42% (12/496) in Manglaralto, and 7.94% (39/491) in Puyo showed Abs. A total of 893 samples of nonfebrile individuals (including conscripts) analyzed for IgM Abs against WNV demonstrated 0.20% (1/496) positivity from Manglaralto and 0.20% (1/491) positivity from Puyo.
The HI results are shown in tables S3 and S4; 149 equines were analyzed. Hemagglutination-inhibiting Abs were found for all the antigens evaluated at serum dilutions of 1:20-1:1280. Counting the monotypic results and the highest dilution showing complete inhibition of hemagglutination in heterotypic reactions, 1.34% (2/149) of equines demonstrated Abs to EEEV, 5.37% (8/149) to WEEV, and 6.04% (9/149) to VEEV across all sampling sites. Abs to EEEV were found in Vinces, Abs to WEEV were detected in Vinces and Manglaralto, and Abs against VEEV were detected in all the locations sampled. Hemagglutination-inhibiting Abs at higher serum dilutions were found for WEEV in one animal in Vinces (1:1280) (table S3).
A total of 84 hamsters were used as sentinels. After exposure in the field, 22.62% of them demonstrated Abs to EEEV (18.51% -5/27 in Vinces, 15.63% -5/32 in Manglaralto, and 25.71% -9/35 in Puyo), 20.24% showed Abs to VEEV (7.41% -2/27 in Vinces, 31.25% -10/32 in Manglaralto, and 20% -7/35 in Puyo), and 2.38% seroconverted to YFV (22.22% -2/9 in Puyo) (table S4). Hemagglutination-inhibiting Abs at higher serum dilutions were found for YFV in one hamster in Puyo (1:640) and for EEEV in two hamsters in Manglaralto (1:1280). None of the hamsters became sick.
A total of 73,505 mosquitoes were captured. The number by locality and species is described in table S5. The most common species was Culex (Mel) spp. (N=63.057). The most abundant by location were Culex (Mel) spp. in Vinces (N=62.736), Aedes (Stg) aegypti in Manglaralto (N=912), and Coquillettidia venezuelensis in Puyo (N=160). Correlation results are shown in table S6. There was a very strong negative correlation (ρ=-1, p≥0.1) between Psorophora ferox and samples with Abs against DENV, SLEV, WNV, EEV, WEEV, and between C. venezuelensis or A. albimanus and samples showing Abs for YFV (ρ=-1 p≥0.1). A strong positive correlation was found between Mansonia pseudotitilians/indubitans or Anopheles (Nys) triannulatus and the proportion of individuals showing Abs against VEEV (ρ=1, p≥0.1). There was a positive correlation between A. aegypti and DENV, SLEV, WNV, EEV, and WEEV, but it was not significant. A list of mosquito species identified in this study and the arboviruses they harbor according to viral isolation studies elsewhere is shown in table S7.
RT-PCR analysis of 128 mosquito pools and 119 hamster tissues did not lead to the amplification of Flavivirus or Alphavirus RNA.
We inoculated 807 suckling Balb/C mice, 576 with hamster tissue and 231 with suspensions of polled mosquitoes. No virus was isolated.
Discussion
Serosurvey of arboviruses is constrained by Ab cross-reactivity because of antigenic similarities between members of the same virus family cocirculating in the same region. The plaque reduction neutralization test (PRNT) is the gold standard for the specific detection of Abs, but it is difficult to implement in resource-limited countries. ELISA and HIA tests are valuable alternatives when monotypic results are accounted for. Using these techniques, we report the percentage of individuals with IgG Abs reactive against DENV, YFV, WNV, SLEV, EEEV, WEEV, and VEEV and with IgM Abs reactive against DENV and WNV. We found monotypic results for WNV, SLEV, and each of the encephalitis viruses evaluated while the frequency was found in individuals of different ages. We also demonstrated the seroconversion of naive sentinel hamsters to the EEEV and VEEV in all sampling locations and for YFV in Amazonia. Mosquito genera recognized as vectors of arboviruses were found in all locations studied. We discuss our findings in the context of arboviruses whose activity has been confirmed in Ecuador and with the current surveillance schemes.
Monotypic IgG Abs against WNV were found in Manglaralto-Santa Elena Province and Vinces-Los Ríos Province (N=58: 34 Vinces, 24 Manglaralto). IgM Abs against this virus were also found in two IgG-positive samples, one from Manglaralto and one from Vinces. The result in Vinces is consistent with the finding of specific IgG Abs against this virus in equines sampled during 2007-2009 by Coello 21; thus, it is probable that WNV infection occurred in humans in this region during these years. The confirmation of human infection would require the use of PNRT; however, these findings warrant further investigation of the activity of this virus in humans, vectors, and potential reservoirs. Serosurvey of birds in Guayaquil (a coastal city around 100 km from Vinces) did not find evidence of WNV infection at this location during 2011 37. Analysis of SLEV showed 20 monotypic results (2 in Manglaralto, 13 in Puyo, and 5 in Vinces). From 1974 to 1978, SLEV was isolated on the Ecuadorian Coast from Culex nigripalpus mosquitoes (isolated 76 V-1177 and 24684) collected in Huertas Negras-Tenguel/Guayas and Puerto Nuevo- Vinces/Los Ríos and sentinel hamsters (isolated 78V-5682) exposed in Playas/Guayas (20). From 2001 to 2004, specific IgM Abs against this virus were found in febrile patients sampled from Ecuadorian Amazonia 25.
A higher percentage of Abs was found against YFV, which was probably due to YFV vaccination history. Overall, the civilian population showed ~60% Abs against YFV, with Puyo showing a higher percentage of civilian people with YFV Abs (52.75%). In conscripts, the percentage of Abs was between 81.31% and 85.05%, which is not optimal but agrees with vaccine efficacy percentages 38,39. Two hamsters seroconverted to YFV in Puyo, which added to the cases of sylvatic yellow fever occurring in Ecuadorian Amazonia showing the active circulation of the virus in this region and the presence of susceptible people. According to information provided by the MSP, mass preventive YFV vaccination in Ecuador was performed in 1944, 1975, 1997, 2005, and 2016, and data about routine immunization of infants have been publicly available since 2007 40. A worldwide analysis of YFV vaccination from 1970 to 2016 41 shows a cumulative vaccination coverage of ~ 60% in the Ecuadorian Amazonia and of ~ 30% in the entire country and recommends reaching the 80% threshold to prevent or control outbreaks.
Regarding Alphavirus, we found IgG Abs reactive against VEEV, EEEV, and WEEV in humans and equines and Abs to VEEV and EEEV in naive hamsters. All these viruses were isolated in Ecuador by Calisher, et al. between 1974 and 1978 20. VEEV comprises six antigenic subtypes (I-VI) of which I-AB and I-C are associated with epizootic/epidemic activity in equines and humans while other serotypes circulate in natural enzootic cycles 42. The epidemic subtype IB (now I-AB) and the enzootic 1-D have been reported in Ecuador; IB, in particular, was reported from a horse showing neurological disease in 1944 and during the outbreak of 1969 7. The isolates from 1969 were antigenically similar to those causing outbreaks in Central and North America 42 and likely originated from incompletely formalin-inactivated vaccines 18,19. IgM Abs against VEEV have been detected in Ecuadorian individuals since 1958, then in 1960, and from febrile patients sampled during 2000-2004 7. The emergence of epizootic/epidemic strains from enzootic VEEV has also been shown and justifies surveillance attention 43,44. For EEEV, specific IgM Abs were detected in two equines showing neurological signs from Chongón, Guayas Province, in 2013 26. EEEV circulating in South America is now known as Madariaga virus45 and has been linked to neurologic disease in humans and horses 46. WEEV has not been reported to cause epizootics or epidemics in Ecuador.
As regards mosquitoes, most of the species/genera identified in this study are known vectors of arboviruses elsewhere (table S7). For example, some Culex species are vectors of WNV, SLEV, and VEEV, and the presence of these mosquitoes in all locations is consistent with the finding of Abs against the viruses they transmit. A. aegypti, the main vector of DENV, was found only in Manglaralto, although Abs against dengue were found in all locations. The absence of this mosquito species in the other locations studied could be due to sampling bias; in Manglaralto, traps were placed near human settlements and in Puyo and Vinces they were placed in more forested areas. Correlation data between species and the proportion of human individuals showing IgG Abs against arbovirus should be interpreted with caution because few locations were analyzed (N=3). Anopheles nuneztovari has not been reported in Ecuador before, but this species is difficult to identify and will need further confirmation. Anopheles are not usual vectors for arboviruses, but the finding of unreported species highlights the importance of carrying out more studies on mosquitoes.
By 2009-2012, arboviral surveillance in Ecuador comprised DENV and YFV through ELISA and nucleic acid detection. After the arrival of ZIKV and CHIKV and the confirmation of OROV, the surveillance was extended to these viruses. Following our results, a percentage of negative samples are currently analyzed for VEEV, EEEV, and MAYV by RT-PCR. Given the emergence and reemergence of arboviral diseases worldwide and their continuous spread, we strongly recommend widening the analysis and surveillance of other arboviruses, especially those isolated previously or detected in febrile patients in the country. These must include the arbovirus surveyed here but also the ILHV and the Playas virus. This last virus, isolated in Ecuador in 1974-1978, is genetically similar to the Cache Valley virus and the Maguari virus, which have been reported to cause human disease. Likewise, we emphasize the need for updated studies and surveillance of mosquito ecology and their vectorial capacity for arbovirus transmission.
We also urge for the establishment of a set of recommendations given by the Global Arbovirus Group of Experts (GAGE) 47, which include strengthening laboratory capacity through the implementation of multiplex diagnostic tests reliable for differential diagnosis, access to viral isolation, sequencing, and genotyping, and encouraging the development of research programs to study the mechanism of arboviral evolution, emergence, and dispersal. Implementation of viral isolation would allow phenotypic characterization and availability of reagents for serological studies, which, in its turn, will reduce the dependence on international reference centers for local diagnosis 48. A constraint of this research was the lack of a local method for the specific analysis of Abs, such as PNRT, which requires cell culture, a procedure used intermittently by the MSP surveillance laboratories. Sequencing and genotyping could now be approached through genomic technologies, whole genome sequencing, and the simultaneous analysis of several samples making it useful for rapid outbreak investigations, viral phylogeography, and evolutionary studies 49-51. Improving epidemiological surveillance and implementation of the recommendations provided by the GAGE could be achieved locally and regionally by promoting collaboration between public laboratories and academia, as recommended by Miranda, et al.52; this approach would be particularly useful for developing countries where infrastructure and technical and economic resources are scarce.