Dengue virus (DENV) transmission from humans to mosquitoes is a poorly documented, but critical component of DENV epidemiology. Magnitude of viremia is the primary determinant of successful human-to-mosquito DENV transmission. People with the same level of viremia, however, can vary in their infectiousness to mosquitoes as a function of other factors that remain to be elucidated. Here, we report on a field-based study in the city of Iquitos, Peru, where we conducted direct mosquito feedings on people naturally infected with DENV and that experienced mild illness. We also enrolled people naturally infected with Zika virus (ZIKV) after the introduction of ZIKV in Iquitos during the study period. Of the 54 study participants involved in direct mosquito feedings, 43 were infected with DENV-2, two with DENV-3, and nine with ZIKV. Our analysis excluded participants whose viremia was detectable at enrollment but undetectable at the time of mosquito feeding, which was the case for all participants with DENV-3 and ZIKV infections. We analyzed the probability of onward transmission during 50 feeding events involving 27 participants infected with DENV-2 based on the presence of infectious virus in mosquito saliva 7-16 days post blood meal. Transmission probability was positively associated with the level of viremia and duration of extrinsic incubation in the mosquito. In addition, transmission probability was influenced by the day of illness in a non-monotonic fashion; i.e., transmission probability increased until 2 days after symptom onset and decreased thereafter. We conclude that mildly ill DENV-infected humans with similar levels of viremia during the first two days after symptom onset will be most infectious to mosquitoes on the second day of their illness. Quantifying variation within and between people in their contribution to DENV transmission is essential to better understand the biological determinants of human infectiousness, parametrize epidemiological models, and improve disease surveillance and prevention strategies.

A novel orthobunyavirus, Bellavista virus, was isolated from Culex (Melanoconion) portesi mosquitoes in the Bellavista neighborhood of Iquitos, Peru, in 2009. The assembled segment L, M, and S sequences of strain PRD0552 are 6,950, 4,469, and 1,256 bases in length, respectively, comprising complete protein-coding sequences and partial terminal untranslated sequences.

Recent years have seen rising incidence of dengue and large outbreaks of Zika and chikungunya, which are all caused by viruses transmitted by Aedes aegypti mosquitoes. In most settings, the primary intervention against Aedes-transmitted viruses is vector control, such as indoor, ultra-low volume (ULV) spraying. Targeted indoor residual spraying (TIRS) has the potential to more effectively impact Aedes-borne diseases, but its implementation requires careful planning and evaluation. The optimal time to deploy these interventions and their relative epidemiological effects are, however, not well understood. We used an agent-based model of dengue virus transmission calibrated to data from Iquitos, Peru to assess the epidemiological effects of these interventions under differing strategies for deploying them. Specifically, we compared strategies where spray application was initiated when incidence rose above a threshold based on incidence in recent years to strategies where spraying occurred at the same time(s) each year. In the absence of spraying, the model predicted 361,000 infections [inter-quartile range (IQR): 347,000-383,000] in the period 2000-2010. The ULV strategy with the fewest median infections was spraying twice yearly, in March and October, which led to a median of 172,000 infections [IQR: 158,000-183,000], a 52% reduction from baseline. Compared to spraying once yearly in September, the best threshold-based strategy utilizing ULV had fewer median infections (254,000 vs. 261,000), but required more spraying (351 vs. 274 days). For TIRS, the best strategy was threshold-based, which led to the fewest infections of all strategies tested (9,900; [IQR: 8,720-11,400], a 94% reduction), and required fewer days spraying than the equivalent ULV strategy (280). Although spraying twice each year is likely to avert the most infections, our results indicate that a threshold-based strategy can become an alternative to better balance the translation of spraying effort into impact, particularly if used with a residual insecticide.

This study describes the evolution of knockdown resistance (kdr) haplotypes in Aedes aegypti in response to pyrethroid insecticide use over the course of 18 years in Iquitos, Peru. Based on the duration and intensiveness of sampling (~10,000 samples), this is the most thorough study of kdr population genetics in Ae. aegypti to date within a city. We provide evidence for the direct connection between programmatic citywide pyrethroid spraying and the increase in frequency of specific kdr haplotypes by identifying two evolutionary events in the population. The relatively high selection coefficients, even under infrequent insecticide pressure, emphasize how quickly Ae. aegypti populations can evolve. In our examination of the literature on mosquitoes and other insect pests, we could find no cases where a pest evolved so quickly to so few exposures to low or nonresidual insecticide applications. The observed rapid increase in frequency of resistance alleles might have been aided by the incomplete dominance of resistance-conferring alleles over corresponding susceptibility alleles. In addition to dramatic temporal shifts, spatial suppression experiments reveal that genetic heterogeneity existed not only at the citywide scale, but also on a very fine scale within the city.

Using a large, passive, clinic-based surveillance program in Iquitos, Peru, we characterized the prevalence of rickettsial infections among undifferentiated febrile cases and obtained evidence of pathogen transmission in potential domestic reservoir contacts and their ectoparasites. Blood specimens from humans and animals were assayed for spotted fever group rickettsiae (SFGR) and typhus group rickettsiae (TGR) by ELISA and/or PCR; ectoparasites were screened by PCR. Logistic regression was used to determine associations between patient history, demographic characteristics of participants and symptoms, clinical findings and outcome of rickettsial infection. Of the 2,054 enrolled participants, almost 2% showed evidence of seroconversion or a 4-fold rise in antibody titers specific for rickettsiae between acute and convalescent blood samples. Of 190 fleas (Ctenocephalides felis) and 60 ticks (Rhipicephalus sanguineus) tested, 185 (97.4%) and 3 (5%), respectively, were positive for Rickettsia spp. Candidatus Rickettsia asemboensis was identified in 100% and 33% of the fleas and ticks tested, respectively. Collectively, our serologic data indicates that human pathogenic SFGR are present in the Peruvian Amazon and pose a significant risk of infection to individuals exposed to wild, domestic and peri-domestic animals and their ectoparasites.

During the last 50 years, the geographic range of the mosquito Aedes aegypti has increased dramatically, in parallel with a sharp increase in the disease burden from the viruses it transmits, including Zika, chikungunya, and dengue. There is a growing consensus that vector control is essential to prevent Aedes-borne diseases, even as effective vaccines become available. What remains unclear is how effective vector control is across broad operational scales because the data and the analytical tools necessary to isolate the effect of vector-oriented interventions have not been available. We developed a statistical framework to model Ae. aegypti abundance over space and time and applied it to explore the impact of citywide vector control conducted by the Ministry of Health (MoH) in Iquitos, Peru, over a 12-year period. Citywide interventions involved multiple rounds of intradomicile insecticide space spray over large portions of urban Iquitos (up to 40% of all residences) in response to dengue outbreaks. Our model captured significant levels of spatial, temporal, and spatio-temporal variation in Ae. aegypti abundance within and between years and across the city. We estimated the shape of the relationship between the coverage of neighborhood-level vector control and reductions in female Ae. aegypti abundance; i.e., the dose-response curve. The dose-response curve, with its associated uncertainties, can be used to gauge the necessary spraying effort required to achieve a desired effect and is a critical tool currently absent from vector control programs. We found that with complete neighborhood coverage MoH intra-domicile space spray would decrease Ae. aegypti abundance on average by 67% in the treated neighborhood. Our framework can be directly translated to other interventions in other locations with geolocated mosquito abundance data. Results from our analysis can be used to inform future vector-control applications in Ae. aegypti endemic areas globally.

The mosquito Aedes aegypti is the vector of a number of medically-important viruses, including dengue virus, yellow fever virus, chikungunya virus, and Zika virus, and as such vector control is a key approach to managing the diseases they cause. Understanding the impact of vector control on these diseases is aided by first understanding its impact on Ae. aegypti population dynamics. A number of detail-rich models have been developed to couple the dynamics of the immature and adult stages of Ae. aegypti. The numerous assumptions of these models enable them to realistically characterize impacts of mosquito control, but they also constrain the ability of such models to reproduce empirical patterns that do not conform to the models' behavior. In contrast, statistical models afford sufficient flexibility to extract nuanced signals from noisy data, yet they have limited ability to make predictions about impacts of mosquito control on disease caused by pathogens that the mosquitoes transmit without extensive data on mosquitoes and disease. Here, we demonstrate how the differing strengths of mechanistic realism and statistical flexibility can be fused into a single model. Our analysis utilizes data from 176,352 household-level Ae. aegypti aspirator collections conducted during 1999-2011 in Iquitos, Peru. The key step in our approach is to calibrate a single parameter of the model to spatio-temporal abundance patterns predicted by a generalized additive model (GAM). In effect, this calibrated parameter absorbs residual variation in the abundance time-series not captured by other features of the mechanistic model. We then used this calibrated parameter and the literature-derived parameters in the agent-based model to explore Ae. aegypti population dynamics and the impact of insecticide spraying to kill adult mosquitoes. The baseline abundance predicted by the agent-based model closely matched that predicted by the GAM. Following spraying, the agent-based model predicted that mosquito abundance rebounds within about two months, commensurate with recent experimental data from Iquitos. Our approach was able to accurately reproduce abundance patterns in Iquitos and produce a realistic response to adulticide spraying, while retaining sufficient flexibility to be applied across a range of settings.

Ticks are obligatory hematophagous ectoparasites that transmit pathogens among various vertebrates, including humans. The microbial and viral communities of ticks, including pathogenic microorganisms, are known to be highly diverse. However, the factors driving this diversity are not well understood. The tropical horse tick, Dermacentor nitens, is distributed throughout the Americas and it is recognized as a natural vector of Babesia caballi and Theileria equi, the causal agents of equine piroplasmosis. In this study, we characterized the bacterial and viral communities associated with partially fed Dermacentor nitens females collected using a passive survey on horses from field sites representing three distinct geographical areas in the country of Colombia (Bolivar, Antioquia, and Cordoba). RNA-seq and sequencing of the V3 and V4 hypervariable regions of the 16S rRNA gene were performed using the Illumina-Miseq platform (Illumina, San Diego, CA, USA). A total of 356 operational taxonomic units (OTUs) were identified, in which the presumed endosymbiont, Francisellaceae/Francisella spp., was predominantly found. Nine contigs corresponding to six different viruses were identified in three viral families: Chuviridae, Rhabdoviridae, and Flaviviridae. Differences in the relative abundance of the microbial composition among the geographical regions were found to be independent of the presence of Francisella-like endosymbiont (FLE). The most prevalent bacteria found in each region were Corynebacterium in Bolivar, Staphylococcus in Antioquia, and Pseudomonas in Cordoba. Rickettsia-like endosymbionts, mainly recognized as the etiological agent of rickettsioses in Colombia, were detected in the Cordoba samples. Metatranscriptomics revealed 13 contigs containing FLE genes, suggesting a trend of regional differences. These findings suggest regional distinctions among the ticks and their bacterial compositions.

Enzootic strains of Venezuelan equine encephalitis virus (VEEV) have been isolated from febrile patients in the Peruvian Amazon Basin at low but consistent levels since the early 1990s. Through a clinic-based febrile surveillance program, we detected an outbreak of VEEV infections in Iquitos, Peru, in the first half of 2006. The majority of these patients resided within urban areas of Iquitos, with no report of recent travel outside the city. To characterize the risk factors for VEEV infection within the city, an antibody prevalence study was carried out in a geographically stratified sample of urban areas of Iquitos. Additionally, entomological surveys were conducted to determine if previously incriminated vectors of enzootic VEEV were present within the city. We found that greater than 23% of Iquitos residents carried neutralizing antibodies against VEEV, with significant associations between increased antibody prevalence and age, occupation, mosquito net use, and overnight travel. Furthermore, potential vector mosquitoes were widely distributed across the city. Our results suggest that while VEEV infection is more common in rural areas, transmission also occurs within urban areas of Iquitos, and that further studies are warranted to identify the precise vectors and reservoirs involved in urban VEEV transmission.

In the Americas, as in much of the rest of the world, the dengue virus vector Aedes aegypti is found in close association with human habitations, often leading to high population densities of mosquitoes in urban settings. In the Peruvian Amazon, this vector has been expanding to rural communities over the last 10-15 years, but to date, the population genetic structure of Ae. aegypti in this region has not been characterized. To investigate the relationship between Ae. aegypti gene flow and human transportation networks, we characterized mosquito population structure using a panel of 8 microsatellite markers and linked results to various potential mechanisms for long-distance dispersal. Adult and immature Ae. aegypti (>20 individuals per site) were collected from Iquitos city and from six neighboring riverine communities, i.e., Nauta, Indiana, Mazan, Barrio Florida, Tamshiaco, and Aucayo. FST statistics indicate significant, but low to moderate differentiation for the majority of study site pairs. Population structure of Ae. aegypti is not correlated with the geographic distance between towns, suggesting that human transportation networks provide a reasonable explanation for the high levels of population mixing. Our results indicate that Ae. aegypti gene flow among sub-populations is greatest between locations with heavy boat traffic, such as Iquitos-Tamshiaco and Iquitos-Indiana-Mazan, and lowest between locations with little or no boat/road traffic between them such as Barrio Florida-Iquitos. Bayesian clustering analysis showed ancestral admixture among three genetic clusters; no single cluster was exclusive to any site. Our results are consistent with the hypothesis that human transportation networks, particularly riverways, are responsible for the geographic spread of Ae. aegypti in the Peruvian Amazon. Our findings are applicable to other regions of the world characterized by networks of urban islands connected by fluvial transport routes.