PMCCPMCCPMCC

Search tips
Search criteria 

Advanced

 
Logo of tropmedLink to Publisher's site
 
Am J Trop Med Hyg. 2016 April 6; 94(4): 786–793.
PMCID: PMC4824219

Enzootic Arbovirus Surveillance in Forest Habitat and Phylogenetic Characterization of Novel Isolates of Gamboa Virus in Panama

Abstract

Landscape changes occurring in Panama, a country whose geographic location and climate have historically supported arbovirus transmission, prompted the hypothesis that arbovirus prevalence increases with degradation of tropical forest habitats. Investigations at four variably degraded sites revealed a diverse array of potential mosquito vectors, several of which are known vectors of arbovirus pathogens. Overall, 675 pools consisting of 25,787 mosquitoes and representing 29 species from nine genera (collected at ground and canopy height across all habitats) were screened for cytopathic viruses on Vero cells. We detected four isolates of Gamboa virus (family: Bunyaviridae; genus: Orthobunyavirus) from pools of Aedeomyia squamipennis captured at canopy level in November 2012. Phylogenetic characterization of complete genome sequences shows the new isolates to be closely related to each other with strong evidence of reassortment among the M segment of Panamanian Gamboa isolates and several other viruses of this group. At the site yielding viruses, Soberanía National Park in central Panama, 18 mosquito species were identified, and the predominant taxa included A. squamipennis, Coquillettidia nigricans, and Mansonia titillans.

Introduction

The geographical location of Panama at the isthmus of North and South America, with boat traffic from around the world through the Panama Canal, creates a potential harbor for a diverse range of vector-borne pathogens, for example, dengue virus (DENV) and Leishmania. Previous work in Panama shows that vector diversity is high with some species capable of transmitting pathogens that can infect humans.13 Mosquito-borne viruses such as yellow fever virus (YFV), Venezuelan equine encephalitis virus (VEEV), and eastern equine encephalitis virus (EEEV) have caused fatal outbreaks in the past, particularly around the Canal region.48 The worldwide spread of DENV has also significantly affected Panama, and recent cases of the globally emerging chikungunya virus have been detected for the first time in this country.912

Panama has a tropical climate and a historically forested landmass with rich biodiversity. However, locales in Panama are undergoing extensive changes in landscape due to increased tourism, immigration, and population growth. Deforestation and encroachment of tropical forest changes the ecosystem's functionality and microclimate, including changes in temperature, sunlight, rainfall patterns, and species composition.13 Such anthropogenic modifications of a habitat can bring wild and domestic animals and humans and their respective pathogens into closer contact and lead to the emergence of infectious diseases/zoonotic pathogens.14 Carrera and others (2013) suggest that human cases of EEEV in Latin America may result from ecological changes bringing humans into increased contact with enzootic transmission cycles.15 Trypanosoma cruzi infection rate of the Chagas disease vector, Rhodnius pallescens, in Panama was found to be higher in deforested habitats and forest fragments compared with contiguous forests.16

It is important to consider that human disturbance of mosquito-breeding habitats and feeding behavior may influence maintenance and transmission of arboviruses. Herein, we performed an arbovirus-mosquito surveillance study in Panama, at four forest sites differing in the level of environmental perturbation, under the supposition that emerging infectious pathogens are circulating in sylvatic cycles and that human modification of native forests affects their prevalence in mosquitoes, possibly providing more opportunities for spillover into new host–vector systems. More explicitly, we tested the hypothesis that arboviral infection prevalence and diversity in mosquitoes might increase with further degradation of forest habitats. We describe the screening of over 25,000 mosquitoes and characterize four novel strains of Gamboa virus (GAMV) from the acrodendrophilic and avian-feeding mosquito, Aedeomyia squamipennis. In addition, at the main (relatively pristine) study site, we compared mosquito species diversity and abundance at two vertical strata.

Materials and Methods

Mosquito collections.

Mosquitoes were trapped over a 3-year period (2009–2012; details given in Table 1) at four tropical-forested sites of central Panama that varied in anthropogenic disturbance and original habitat quality. At first, two degraded sites Achiote and Las Pavas, comprising patches of second-growth forest representative of intermediate level of disturbance (forest cover > 35% and < 55%, respectively) and a more pristine site (Barro Colorado Island) represented by old-growth forest (forest cover > 55%) were assessed as part of a collaborative project between the Smithsonian Tropical Research Institute (STRI) and the U.S. Environmental Protection Agency. A fourth key site (Soberanía National Park) was later included, to involve collections at a pristine (forest cover > 55%), non-island site during November and December 2012. Locations of the four sites are mapped in Figure 1 .

Figure 1.
Map indicating study sites in central Panama.
Table 1
Collections of mosquitoes tested from four Panamanian forest sites

Mosquitoes were collected using CDC light traps (John W. Hock Company, Gainesville, FL), baited with CO2 from a dry ice source. Traps were set overnight at six points (spaced over 200 m apart) along transects situated away from main paths. Traps were placed at ground (1.5 m) or canopy (> 30 m) height, alternating each night. Geographic coordinates (values of latitude and longitude) for each collecting point were recorded using a handheld Global Positioning System unit (Garmin International, Olathe, KS) in WGS84 datum and imported into ArcView Geographic Information System (GIS) software (ESRI, Redlands, CA) to calculate the percentage of forest cover under a GIS scheme. Updated GIS data sets were obtained from Naos molecular biology and evolution laboratories at STRI, in Panama City. At dawn, captured mosquitoes were separated from other insect fauna and identified to species level using a chill table and taxonomic keys by Pecor and others, Wilkerson and others, and Sallum and Forattini.1719

Mosquitoes were pooled according to collection date, species, trap height, and location (up to 40 individuals per pool), and stored at −80°C before transfer for further analysis. Pools were sent to the Arbovirus Laboratory, Wadsworth Center (New York State Department of Health, Slingerlands, NY) for screening by cell culture.

Virus isolation and identification.

At the Wadsworth Center, pools were homogenized in mosquito diluent (phosphate-buffered saline, 20% fetal bovine serum, and antibiotics), using a TissueLyser (Qiagen, Valencia, CA). They were then clarified by centrifugation at 12,000 rpm for 4 minutes, and the supernatant was removed, respun, and then screened for cytopathic effect (CPE) on Vero (African green monkey kidney) cell culture. In brief, 100 μL of each homogenate supernatant was used to inoculate a sub-confluent monolayer of Vero cells on 6-well plates (Costar, Corning, NY) and maintained at 37°C with 5% CO2. Cells were monitored daily, for up to 10 days, and samples showing CPE were transferred to fresh plates for confirmation, before virus was harvested and stored at −80°C.

Virus isolates were identified either using indirect immunofluorescent antibody assays (IFA) (with mouse hyperimmune ascetic fluid provided by the Centers for Disease Control and Prevention, against a variety of arbovirus groups including alphavirus, flavivirus, and eight bunyavirus groups) or via reverse transcription polymerase chain reaction following RNA extraction using QIAamp Viral RNA Mini kit (Qiagen) according to the manufacturer's instructions, using primers specific for certain orthobunyavirus serogroups, generic alphaviruses, or flaviviruses.2023

Sequencing and phylogenetic analyses.

CPE-positive samples were submitted to the University of Texas Medical Branch (UTMB) Next Generation Sequencing Core. Virus was amplified and concentrated as previously described24; RNA was extracted and Illumina sequencing performed to determine the complete genome sequences for these cytopathic viruses. In brief, viral RNA was fragmented by incubation at 94°C for 8 minutes in 19.5 mL of fragmentation buffer (Illumina Inc., San Diego, CA). First and second strand synthesis, adapter ligation, and amplification of the library were performed using an Illumina TruSeq RNA Sample Preparation kit, under conditions prescribed by the manufacturer. Cluster formation of the library DNA templates was performed using the TruSeq PE Cluster Kit (Illumina Inc.) and the Illumina cBot workstation. Paired-end 50-base sequencing by synthesis was performed using TruSeq SBS kit v3 (Illumina Inc.) on an Illumina HiSeq 1000 using protocols defined by the manufacturer. Cluster density per lane was 820–940 K clusters/mm2 and post filter reads ranged from 148 to 218 million per lane. Base call conversion to sequence reads was performed using CASAVA-1.8.2 (Illumina Inc.). Reads were filtered for quality and adapter sequences were removed, then viral contigs were assembled de novo using AbySS software.25 Assembled contigs were checked using bowtie2 to align reads to the contigs, followed by visualization using the integrative genomics viewer.26,27

For the orthobunyavirus isolates, sequences from the small (S) (around 700 bp), medium (M) (around 4,500 bp), and large (L) (around 6,800 bp) segments of the genome were aligned using Seaview (version 4.5) and neighbor-joining (NJ) and maximum likelihood (ML) trees were constructed using a HKY85 substitution model, and 1,000 bootstrap replicates in Paup (version 4.0) software (Sinauer Associates, Sunderland, MA).28,29 Isolate sequences were compared with sequences of viruses within the wider Orthobunyavirus genus. We also focused specifically on Gamboa serogroupings (obtaining similar viruses using BLAST, http://blast.ncbi.nlm.nih.gov/Blast.cgi) and created a genetic distance matrix (HKY85) in Paup.28

Statistics.

Species richness (S) is the total number of species at the site. Mosquito species diversity was calculated using the formula for Shannon's H (where Pi is the proportion of species at each site), and evenness of species (EH) by dividing H by ln(S).

equation image

Total counts of mosquitoes, as per species, were recorded and compared with vertical strata at the main (relatively pristine) study site, Soberanía National Park. Mann–Whitney U tests were used to compare species diversity of canopy and ground level collections, using a 5% significance level with critical values of test statistic U based on Scanopy and Sground. Arboviral infection prevalence in mosquitoes was calculated for cytopathic viruses using maximum likelihood estimates (MLE) of pooled samples, for example, infection rate per 1,000 mosquitoes, using the calculator available at http://www.cdc.gov/westnile/resourcepages/mosqSurvSoft.html and compared across sites and vertical strata.

Results

Mosquito collection.

A total of 25,787 mosquitoes representing > 60% of the specimens trapped at four sites (i.e., two fairly pristine forest habitats at Soberanía National Park and Barro Colorado Island, and two more degraded habitats, Las Pavas and Achiote) were screened for arboviruses. The collection comprised 29 species from nine genera, with the most frequent taxa captured overall being Culex coronator (20.7%), A. squamipennis (14.4%), and Coquillettidia venezuelensis (8.3%). The number of mosquitoes and pools tested at various sites and elevations are shown in Table 1. At the main site discussed here, Soberanía, we trapped 4,139 mosquitoes (later screened as 139 pools) with a species richness (S) of 18 (representing seven genera). Species diversity entropy in this site was H = 2.061 and evenness EH = 0.713. Comparing across vertical strata at Soberanía, there was no significant difference between ground and canopy level, that is, species richness and overall diversity were similar in each stratum (Sground = 16, Scanopy = 14; Hground = 1.79, Hcanopy = 1.81; Mann–Whitney U = 75, P > 0.05). However, significantly more individuals were captured at canopy level (average/night ground 38.2, canopy 65.4; adjusting for an effort of 45 trap nights at ground level and 37 at canopy level; P < 0.05) than at understory. Furthermore, the predominant species in each stratum differed: Cx. nigripalpus, Psorophora cingulata, and Cx. conservator were only found at ground level, while Mansonia flaveola and Sabethes spp. were only detected at canopy level. The most frequently trapped species at Soberanía were A. squamipennis (31%, N = 1,275), the majority of which were detected at canopy level (93%); Cq. nigricans (21%, N = 838; of which 78% were trapped at ground level), and M. titillans (17%, N = 719; 71% collected at ground level). The fourth most frequently detected mosquito, Anopheles triannulatus, was detected 67% of the time at canopy level, whereas the majority (85%) of the fifth commonest species, An. punctimacula, were detected at ground level.

Virus isolation and identification.

In total, across all Panamanian sites, 675 mosquito pools (of 29 different species) were screened for arboviruses. Four mosquito pools showed positive CPE seven days postinoculation. All four isolates were detected from pools of A. squamipennis captured at canopy level of Soberanía in November or December 2012 and were shown by IFA to belong to the complex Gamboa (family Bunyaviridae; enveloped, tripartite, negative-sense RNA virus). Infection prevalence for GAMV among all 3,677 A. squamipennis screened was MLE = 1.11 (per 1,000 mosquitoes) (95% CI: 0.36, 2.67), and MLE = 3.3 (1.1, 9.1) at Soberanía site alone (4/29 pools; N = 1,275 individuals).

Sequencing and phylogenetic analyses.

We sequenced the full-length genomes of all four isolates of GAMV (GenBank accession no's: KT950259 (S), KT950267 (M), and KT950263 (L) (GAM118); KT950260, KT950268, and KT950264 (GAM122); KT950261, KT950269, and KT950265 (GAM130); and KT950262, KT950270, and KT950266 (GAM131). HKY85 genetic distance matrices (Table 2) show the relatedness among these and other closely related GAMV strains.30 When comparing our new strains to other virus groups in the Orthobunyavirus genus, tree topologies varied depending on the genome segment used for analysis. Their phylogenetic relationships based on the S and L segments (Figure 2A and andCC ) show that strains GAM118 and GAM130, and GAM122 and GAM131, form evolutionary pairings (genetic distances between them being 0.000 and 0.007 respectively for S segment, and 0.002 and 0.01 for L segment) and have a more recent common ancestry with Alajuela and Calchaqui virus as outgroups of this complex. On the basis of the M segment, GAM118 and GAM130 cluster together along with Gamboa strain GML435718 (accession number: KM272175; isolated in Colon, Panama, 1986), with genetic distances of 0.002 between GAM118 and GAM130, and 0.008 between GAM118 or GAM130 with GML435718). Furthermore, GAM131 groups with the prototype Gamboa strain MARU10962 (accession number: KM272181; isolated in Panama, 1962), showing a distance matrix score of 0.015, while GAM122 groups with Alajuela virus MARU11079 (accession number: KM272187; isolated in Panama, 1963) with a distance matrix score of 0.012, with these pairings showing evolutionary distance from each other (Figure 2B). Both NJ and ML trees showed similar topologies, thus we show only the former here. The Gamboa group of viruses is evolutionary distinct from other viruses in the genus (Figure 2), but based on the S segment, the Gamboa complex sits nearest to Group C orthobunyaviruses; on the L segment to the California serogroup; and on the M segment, to both California and Bunyamwera serogroups.

Figure 2.
Phylogeny tree showing the phylogenetic relationships of the (A) S segment, (B) M segment, and (C) L segment of Gamboa virus isolates and selected Orthobunyavirus sequences. On the basis of complete segment, open reading frame nucleotide sequences were ...
Table 2
Distance matrices of Gamboa group viruses, using a HKY85 model, based on the M segment, L segment, or S segment

Discussion

On the basis of previous disease outbreaks, virus isolations and conjecture regarding its geographic location, Panama would appear to be a “hot spot” for vector-borne pathogen activity. Yellow fever and malaria have historically caused considerable morbidity and mortality in the Canal region, and more recently, cases and isolates of both EEEV and VEEV have been recognized in eastern and central Panama.6,15,31 However, upon screening over 25,000 mosquitoes from four sites of central Panama, relatively few virus isolates were obtained in this study. Despite not detecting EEEV, VEEV, and sylvatic YFV from this area of the country, we make advancements in understanding enzootic vector distributions of GAMV. Our study also provides an opportunity to examine potential arbovirus vectors involved in disease emergence, and via specific examination of the Soberanía National Park, information on mosquito diversity at a novel habitat is reported.

We did not compare trends of mosquito species diversity-richness and abundance across our four sites for two main reasons: 1) Sampling efforts were different between Soberanía National Park and the other sites (15 days of November–December in Soberanía versus 96–125 days from any month of rainy season over the course of 3 years at the other three sites) and 2) due to logistics, the mosquitoes tested for virus from other sites were only subsets of larger collections, therefore Soberanía was the only site for which the total number of specimens collected were screened. For these reasons, and also due to data ownership, data on mosquito diversity-richness and abundance at Barro Colorado Island, Las Pavas, and Achiote will be the topic of a separate publication. Of the wide range of mosquitoes collected overall, many species have previously been reported as arbovirus vectors, for example, Culex (Melanoconion) spp. are associated with the enzootic cycle of VEEV and group C viruses, Cx. nigripalpus is a vector of the flavivirus St. Louis encephalitis, and YFV has been isolated from Aedes serratus.6,32,33 Frequently captured, Cx. declarator, M. titillans, and Cq. venezuelensis are also reported as potential arbovirus bridge vectors.34 In general, there are indications that the predominant mosquito species would vary according to habitat quality. A. squamipennis and M. titillans were most frequently observed in pristine habitats; conversely in a degraded forest habitat, species such as Cx. coronator and Cx. declarator were more common. These and other mosquitoes found at disturbed forest sites (where lower species richness overall was observed) are known vectors of medically important arboviruses elsewhere,3436 and we might have expected a higher virus prevalence at these disturbed sites. Future studies will have to test the hypothesis concerning the impact of forest disturbance on vector-species diversity by characterizing habitat degradation in more depth and using less biased trapping methods such as larval mosquito collections. Similar differences in the pattern of mosquito species dominance were found between vertical strata, indicating that Cx. nigripalpus, Ps. cingulata, and Cx. conservator are mainly understory species while M. flaveola, Sabethes spp., An. triannulatus, and A. squamipennis were most commonly collected at canopy level. Studies in pristine forest habitats elsewhere (e.g., Jones and others, 2004) have similarly found that certain mosquito species prevail more frequently according to a vertical niche.37 The potential implications of these differences for arbovirus transmission are unknown at this point, but they may include, for example, canopy mosquitoes such as An. triannulatus and A. squamipennis shifting from a canopy feeding tendency in pristine site to a ground feeding tendency in degraded habitat, thus increasing the transmission risk for some pathogens in human-altered forest ecosystems.

Our GAMV isolates were obtained from mosquitoes at a relativity pristine forest site, where the virus vector A. squamipennis was a predominant mosquito species, particularly at canopy level. The only known vector of GAMV, A. squamipennis is associated with ground water and abundant vegetation; it has been shown to have avian feeding preferences, and to be a vector for bird parasites such as avian Plasmodia.3,38 Transovarial transmission of GAMV in A. squamipennis has been demonstrated (Mmnimum field infection rate = 5.1/1,000 mosquitoes),39 and therefore could be an important maintenance mechanism of the virus in Panama.

GAMV was first described by Calisher and others (1981) using complement-fixation techniques, suggesting that there are two antigenic complexes and at least eight serotypes.39 Distribution of GAMV is noted to include Panama, Honduras, Argentina, Surinam, and Ecuador.40,41 There is pathology in mice with death by day 5–6 with an estimated viremia of 6.4 log10/mL, yet these viruses are not known to impact human or veterinary health.42 Hemagglutination inhibition studies (of prototype strain, MARU10962) suggested a close relationship to Capim virus42; however, genetic sequences for many bunyaviruses including Capim complex viruses are not available for comparison.

Figure 2 presents three phylogenies of the genus Orthobunyavirus and shows the distinct relationship of GAMV within the genus. Particularly contrasting the GAMV prototypes, Alajuela and Calchaqui viruses, recently placed on GenBank,30 we generate a more detailed evolutionary comparison for this group of viruses. We highlight the fact that Alajuela and Calchaqui viruses are variants within the complex of GAMVs. It is interesting that the phylogenetic relationship differs depending on which segment of the genome is used for comparison, with relationships in the M segment suggestive of a previous reassortment event. Particularly noteworthy is that in the M segment phylogeny, GAM122, GAM131, Alajuela, and the prototype Gamboa MARU10962 strains appear to be reassortants, whereas GAM118, GAM130, and other Gamboa strains show the same cluster relationship as in the S and L segments. There is also evidence of reassortment between Calchaqui and the Gamboa group in the M segment, compared with the outgrouping of Calchaqui in the L and S segments. We expect that further viruses exist that could complement the relationships.

We screened all mosquito pools on mammalian cell lines, to focus on detection of disease agents that might directly affect humans; however, the use of other cell lines such as C6/36 can lead to the detection of additional pathogens. Testing of the same Panamanian mosquito pools for insect-specific viruses, some of which may have important effects on vector-pathogen ecology (e.g., modifying the vector competence of mosquitoes for more serious human arboviral pathogens) will be described in a future publication.

For future studies, alternative trapping methods (especially baits that may attract different vectors) and further sylvatic regions could be assessed. In particular, when virus vectors are identified, mosquito-host interactions would be determined (e.g., by examining host DNA in blood-meals of engorged vectors to ascertain feeding patterns). This can indicate ecological cycles and vertebrates potentially exposed to arboviruses, or acting as reservoir hosts.43 For example, Diallo and others (2013) suggest Aedes taylori to have a role in sylvatic YFV and chikungunya virus transmission in Senegal after confirming feeding on monkeys.44

Despite no significant differences in the diversity of mosquito species found at canopy versus ground strata of Soberanía forests, we would not have detected GAMV without canopy sampling, and so highlight the importance of surveillance within all micro niches. Some daytime collections were attempted during our study and showed the potential for the capture of Culex, Sabethes, and Haemagogus species, the latter being important yellow fever vectors,45 therefore we recommend that future sampling also considers daytime vector species. Similarly, seasonality and climate are likely important factors affecting the population dynamics of both mosquito species and arboviruses. Soberanía collections were made in November and December, toward the end of the rainy season, although other sites were sampled at variable times of year and in multiple years. Importantly, we highlight the need to consider the ecosystem as a whole, assessing seasonality, and including the response of multiple vector types (e.g., mosquitoes, sandflies, biting midges, and ticks) and their preferred host species, to degradation in habitat. We hypothesize that opportunity for arbovirus emergence could vary seasonally according to habitat quality, with different vectors gaining prominence in disease transmission throughout the year and habitat quality gradients.

ACKNOWLEDGMENTS

We kindly acknowledge the Smithsonian Tropical Research Institute in Panama and the U.S. Environmental Protection Agency for enabling the extension of the original project idea comparing mosquito species diversity between differentially modified habitats. Jose R. Rovira provided outstanding entomological skills. We highly appreciate the field skills and support of Omar Lopez (INDICASAT-AIP), Gaudenia Mendoza (University of Panama), and David Walton. The Centers for Disease Control and Prevention (CDC) kindly provided MHIAF grouping fluids for key virus groups that were supportive in virus identification. Joseph Maffei of NYSDOH was supportive in IFA assays and screening. We are grateful to Steven Widen, Tom Wood, and Jill Thompson of UTMB Next Generation Sequencing core for virus sequencing, and Naomi Forrester of UTMB for guidance in phylogenetics and manuscript review. We are thankful to Octavio Smith from CATHALAC for generating the map in Figure 1. Finally, we thank Panama's National Environmental Authority (ANAM) for supporting scientific collecting and exporting of mosquitoes.

Footnotes

Financial support: Funding for this project was provided via the Robert E. Shope International Fellowship Award granted by the American Society of Tropical Medicine and Hygiene to Gillian Eastwood in 2012. Gillian Eastwood and Albert J. Auguste are also funded by the James W. McLaughlin Endowment fund. Jose R. Loaiza is funded by the National Research Investigator Board (SNI) of the Secretariat for Science, Technology and Innovation of Panama (SENACYT). The U.S. Environmental Protection Agency project “Mosquito Species Diversity and Landscape Change DW33-92296801-0” awarded to Montira Pongsiri covered the costs of mosquito sampling at Las Pavas, Achiote, and Barro Colorado Island, and some of the laboratory supplies and reagents used at NYSDOH.

Authors' addresses: Gillian Eastwood and Laura D. Kramer, Wadsworth Center, New York State Department of Health, Slingerlands, NY, E-mail: moc.liamtoh@g2llig and vog.yn.htlaeh@remark.arual. Jose R. Loaiza, Instituto de Investigaciones Científicas y Servicios de Alta Tecnología (INDICASAT), Ciudad del Saber, Panama, E-mail: ap.gro.tasacidni@aziaolj. Montira J. Pongsiri, U.S. Environmental Protection Agency, Washington, DC, E-mail: vog.ape@aritnom.irisgnop. Oris I. Sanjur, Molecular Biology and Evolution Laboratories, Smithsonian Tropical Research Institute, Ancon, Panama, E-mail: ude.is@orujnas. James E. Pecor, Department of Entomology, Walter Reed Army Institute of Research, Washington, DC, E-mail: ude.is@jrocep. Albert J. Auguste, Department of Pathology, University of Texas Medical Branch, Galveston, TX, E-mail: ude.bmtu@sugua1ja.

References

1. Estripeaut D, Aramburú MG, Sáez-Llorens X, Thompson HA, Dasch GA, Paddock CD, Zaki S, Eremeeva ME. Rocky Mountain spotted fever, Panama. Emerg Infect Dis. 2007;13:1763–1765. [PMC free article] [PubMed]
2. Galindo P, Carpenter SJ, Trapido H. Ecological observations on forest mosquitoes of an endemic yellow fever area in Panama. Am J Trop Med Hyg. 1951;31:98–137. [PubMed]
3. Loaiza JR, Bermingham E, Scott ME, Rovira JR, Conn JE. Species composition and distribution of adult Anopheles (Diptera: Culicidae) in Panama. J Med Entomol. 2008;45:841–851. [PMC free article] [PubMed]
4. Galindo P, Srihongse S. Evidence of recent jungle yellow-fever activity in eastern Panama. Bull World Health Organ. 1967;36:151–161. [PubMed]
5. Galindo P, Grayson MA. Culex (Melanoconion) aikenii: natural vector in Panama of endemic Venezuelan encephalitis. Science. 1971;172:594–595. [PubMed]
6. Quiroz E, Aguilar PV, Cisneros J, Tesh RB, Weaver SC. Venezuelan equine encephalitis in Panama: fatal endemic disease and genetic diversity of etiologic viral strains. PLoS Negl Trop Dis. 2009;3:e472. [PMC free article] [PubMed]
7. Dietz WH, Jr, Galindo P, Johnson KM. Eastern equine encephalomyelitis in Panama: the epidemiology of the 1973 epizootic. Am J Trop Med Hyg. 1980;29:133–140. [PubMed]
8. Aguilar PV, Estrada-Franco JG, Navarro-Lopez R, Ferro C, Haddow AD, Weaver SC. Endemic Venezuelan equine encephalitis in the Americas: hidden under the dengue umbrella. Future Virol. 2011;6:721–740. [PMC free article] [PubMed]
9. Armien B, Suaya JA, Quiroz E, Sah BK, Bayard V, Marchena L, Campos C, Shepard DS. Clinical characteristics and national economic cost of the 2005 dengue epidemic in Panama. Am J Trop Med Hyg. 2008;79:364–371. [PubMed]
10. Bayard V, Quiroz E, Mójica J. Re-emergence of dengue in Panama [in Spanish] Rev Med Panama. 1996;21:85–92. [PubMed]
11. Quiroz E, Ortega M, Guzmán MG, Vázquez S, Pelegrino JL, Campos C, Bayard V, Vázquez M, Kourí G. Dengue in Panama, 1993 [in Spanish] Rev Cubana Med Trop. 1997;49:86–93. [PubMed]
12. Díaz Y, Carrera JP, Cerezo L, Arauz D, Guerra I, Cisneros J, Armién B, Botello AM, Araúz AB, Gonzalez V, López Y, Moreno L, López-Vergés S, Moreno BA. Chikungunya virus infection: first detection of imported and autochthonous cases in Panama. Am J Trop Med Hyg. 2015;92:482–485. [PMC free article] [PubMed]
13. IUFRO (International Union of Forest Research Organizations) Thompson ID, editor. Chapter 2: forest biodiversity, carbon and other ecosystem services: relationships and impacts of deforestation and forest degradation. Understanding Relationships between Biodiversity, Carbon, Forests and People: The Key to Achieving REDD+ Objectives. 2012. http://www.iufro.org/science/gfep/biodiv-forman-redd-panel/report/download-chapter Available at. Accessed November 8, 2014.
14. Daszak P, Cunningham AA, Hyatt AD. Emerging infectious diseases of wildlife: threats to biodiversity and human health. Science. 2000;287:443–449. [PubMed]
15. Carrera JP, Forrester N, Wang E, Vittor AY, Haddow AD, López-Vergés S, Abadia I, Castaño E, Sosa N, Báez C, Estripeaut D, Díaz Y, Beltrán D, Cisneros J, Cedeño HG, Travassos da Rosa AP, Hernandez H, Martinez-Torres AO, Tesh RB, Weaver SC. Eastern equine encephalitis in Latin America. N Engl J Med. 2013;369:732–744. [PMC free article] [PubMed]
16. Gottdenker NL, Chaves LF, Calzada JE, Saldaña A, Carroll CR. Host life history strategy, species diversity, and habitat influence Trypanosoma cruzi vector infection in changing landscapes. PLoS Negl Trop Dis. 2012;6:e1884. [PMC free article] [PubMed]
17. Wilkerson RC, Strickman D, Litwak TR. Illustrated key to the female anopheline mosquitoes of Central America and Mexico. J Am Mosq Control Assoc. 1990;6:7–34. [PubMed]
18. Pecor JE, Mallampalli VL, Harbach RE, Peyton EL. Catalog and illustrated review of the subgenus Melanoconion of Culex (Diptera: Culicidae) Contrib Am Entomol Inst. 1992;27:1–228.
19. Sallum MA, Forattini OP. Revision of the Spissipes section of Culex (Melanoconion) (Diptera: Culicidae) J Am Mosq Control Assoc. 1996;12:517–600. [PubMed]
20. Dunn EF, Pritlove DC, Elliott RM. The S RNA genome segments of Batai, Cache Valley, Guaroa, Kairi, Lumbo, Main Drain and Northway bunyaviruses: sequence determination and analysis. J Gen Virol. 1994;75:597–608. [PubMed]
21. Kuno G, Mitchell CJ, Chang GJ, Smith GC. Detecting bunyaviruses of the Bunyamwera and California serogroups by a PCR technique. J Clin Microbiol. 1996;34:1184–1188. [PMC free article] [PubMed]
22. O'Guinn ML, Lee JS, Kondig JP, Fernandez R, Carbajal F. Field detection of eastern equine encephalitis virus in the Amazon Basin region of Peru using reverse transcription-polymerase chain reaction adapted for field identification of arthropod-borne pathogens. Am J Trop Med Hyg. 2004;70:164–171. [PubMed]
23. Sanchez-Seco MP, Rosario D, Domingo C, Hernandez L, Valdes K, Guzman MG, Tenorio A. Generic RT-nested-PCR for detection of flaviviruses using degenerated primers and internal control followed by sequencing for specific identification. J Virol Methods. 2005;126:101–109. [PubMed]
24. Auguste AJ, Kaelber JT, Fokam EB, Guzman H, Carrington CVF, Erasmus JH, Kamgang B, Popov VL, Jakana J, Liu X, Wood TG, Widen SG, Vasilakis N, Tesh RB, Chiu W, Weaver SC. A newly isolated reovirus has the simplest genomic and structural organization of any reovirus. J Virol. 2015;89:676–687. [PMC free article] [PubMed]
25. Simpson JT, Wong K, Jackman SD, Schein JE, Jones SJ, Birol I. ABySS: a parallel assembler for short read sequence data. Genome Res. 2009;19:1117–1123. [PubMed]
26. Langmead B, Salzberg SL. Fast gapped-read alignment with Bowtie 2. Nat Methods. 2012;9:357–359. [PMC free article] [PubMed]
27. Robinson JT, Thorvaldsdottir H, Winckler W, Guttman M, Lander ES, Getz G, Mesirov JP. Integrative genomics viewer. Nat Biotechnol. 2011;29:24–26. [PMC free article] [PubMed]
28. Gouy M, Guindon S, Gascuel O. SeaView version 4: a multiplatform graphical user interface for sequence alignment and phylogenetic tree building. Mol Biol Evol. 2010;27:221–224. [PubMed]
29. Swofford DL. PAUP*: Phylogenetic Analysis Using Parsimony (*and Other Methods). Version 4. Sunderland, MA: Sinauer Associates, Inc.; 2002.
30. Nunes MR, Chiang JO, de Lima CP, Martins LC, Aragao Dias A, Cardoso JF, Silva SP, Da Silva DE, Oliveira LF, Vasconcelos JM, Ferreira JP, Travassos da Rosa AP, Guzman H, Tesh RB, Vasconcelos PF. New genome sequences of Gamboa viruses (family Bunyaviridae, genus Orthobunyavirus) isolated in Panama and Argentina. Genome Announc. 2014;2:e00940–e14. [PMC free article] [PubMed]
31. Monath TP. Yellow fever: Victor, Victoria? Conqueror, conquest? Epidemics and research in the last forty years and prospects for the future. Am J Trop Med Hyg. 1991;45:1–43. [PubMed]
32. Day JF, Curtis GA. Blood feeding and oviposition by Culex nigripalpus (Diptera: Culicidae) before, during, and after a widespread St. Louis encephalitis virus epidemic in Florida. J Med Entomol. 1999;36:176–181. [PubMed]
33. Cardoso JC, de Almeida MAB, dos Santos E, da Fonseca DF, Sallum MAM, Noll CA, Monteiro HAO, Cruz ACR, Carvalho VL, Pinto EV, Castro FC, Nunes Neto JP, Segura MNO, Vasconcelos PFC. Yellow fever virus in Haemagogus leucelaenus and Aedes serratus mosquitoes, southern Brazil, 2008. Emerg Infect Dis. 2010;16:1918–1924. [PMC free article] [PubMed]
34. Parra-Henao G, Suarez L. Mosquitoes (Diptera: Culiciadae) as potential vectors of arboviruses in the Uraba region, northwest of Colombia [in Spanish] Biomedica. 2012;32:252–262. [PubMed]
35. Confalonieri UEC, Costa Neto C. Diversity of mosquito vectors (Diptera: culicidae) in Caxiuana, Para, Brazil. Interdiscip Perspect Infect Dis. 2012;2012:741273. [PMC free article] [PubMed]
36. Moulis RA, Russell JD, Lewandowski HB, Jr, Thompson PS, Heusel JL. Culex coronator in coastal Georgia and South Carolina. J Am Mosq Control Assoc. 2008;24:588–590. [PubMed]
37. Jones JW, Turell MJ, Sardelis MR, Watts DM, Coleman RE, Fernandez R, Carbajal F, Pecor JE, Calampa C, Klein TA. Seasonal distribution, biology, and human attraction patterns of culicine mosquitoes (Diptera: Culicidae) in a forest near Puerto Almendras, Iquitos, Peru. J Med Entomol. 2004;41:349–360. [PubMed]
38. Loaiza JR, Miller MJ. Seasonal pattern of avian Plasmodium-infected mosquitoes and implications for parasite transmission in central Panamá Parasitol Res. 2013;112:3743–3751. [PubMed]
39. Dutary BE, Petersen JL, Peralta PH, Tesh RB. Transovarial transmission of Gamboa virus in a tropical mosquito, Aedeomyia squamipennis. Am J Trop Med Hyg. 1989;40:108–113. [PubMed]
40. Calisher CH, Lazuick JS, Justines G, Francy DB, Monath TP, Gutierrez E, Sabattini MS, Bowen GS, Jakob WL. Viruses isolated from Aedeomyia squamipennis mosquitoes collected in Panama, Ecuador, and Argentina: establishment of the Gamboa serogroup. Am J Trop Med Hyg. 1981;30:219–223. [PubMed]
41. Calisher CH, Lazuick JS, Sudia WD. Brus Laguna virus, a Gamboa bunyavirus from Aedeomyia squamipennis collected in Honduras. Am J Trop Med Hyg. 1988;39:406–408. [PubMed]
42. Justines G. Catalogue of arboviruses: Gamboa (No. 240) Am J Trop Med Hyg. 1970;19:1153–1154. [PubMed]
43. Navia-Gine WG, Loaiza JR, Miller MJ. Mosquito-host interactions during and after an outbreak of equine viral encephalitis in eastern Panama. PLoS One. 2013;8:e81788. [PMC free article] [PubMed]
44. Diallo D, Chen R, Diagne CT, Ba Y, Dia I, Sall AA, Weaver SC, Diallo M. Bloodfeeding patterns of sylvatic arbovirus vectors in southeastern Senegal. Trans R Soc Trop Med Hyg. 2013;107:200–203. [PMC free article] [PubMed]
45. Barrett AD, Higgs S. Yellow fever: a disease that has yet to be conquered. Annu Rev Entomol. 2007;52:209–229. [PubMed]

Articles from The American Journal of Tropical Medicine and Hygiene are provided here courtesy of The American Society of Tropical Medicine and Hygiene