Abstract
Dengue viruses (DENVs) cause the greatest public health burden globally among the arthropod-borne viruses. DENV transmission risk has also expanded from tropical to subtropical regions due to the increasing range of its principal mosquito vector, Aedes aegypti. Focal out-breaks of dengue fever (dengue) in the state of Florida (FL) in the USA have increased since 2009. However, little is known about the competence of Ae. aegypti populations across different regions of FL to transmit DENVs. To understand the effects of DENV genotype and serotype variations on vector susceptibility and transmission potential in FL, we orally infected a colony of Ae. aegypti (Orlando/ORL) with laboratory and wild type DENV-1 through -4. Wild type DENVs were more infectious to and transmissible by ORL mosquitoes. We used these same DENVs to examine natural Ae. aegypti populations to determine whether spatial distributions correlated with differential vector competence. Vector competence across all DENV serotypes was greater for mosquitoes from areas with the highest dengue incidence in south FL compared to north FL. Vector competence for wild type DENVs was significantly higher, revealing that transmission risk is influenced by virus/vector combinations. These data support a targeted mosquito-plus-pathogen screening approach to more accurately estimate DENV transmission risk.
1. Introduction
Dengue viruses (DENVs) belong to the genus Flavivirus in the family Flaviviridae. There are four genetically related dengue viruses termed DENV-1, DENV-2, DENV-3, and DENV-4. These arthropod-borne viruses (arboviruses) are transmitted to humans through the bite of an infected Aedes species mosquito. The primary vector of these viruses is Aedes aegypti (Ae. aegypti). DENVs are the causative agents of dengue fever (“dengue”) and a severe form of dengue termed dengue hemorrhagic fever. Dengue manifests a wide spectrum of clinical signs that range from subclinical (inapparent) disease to illnesses accompanied by fever, aches and pains (eye pain, and muscles, joint, and bone aches and pains), nausea, vomiting, and a rash [1]. Severe dengue can result in internal bleeding, shock, and death; bleeding can occur from the nose or gums, and blood can be present in vomit or stool. According to the World Health Organization, dengue is the most critical mosquito-borne disease in the world [1]. Between 1990 and 2013, symptomatic infections from DENV quadrupled, and nearly half of the world’s population is currently at-risk for DENV infection since the principal mosquito vector has extended its home range and now thrives in countries wherein it was previously absent [2-4]. Thus, DENVs have by now caused outbreaks of dengue in over 120 countries with an estimated 400,000 infections and 10,000 deaths globally each year [1, 4]. Despite a lengthy historical record of dengue outbreaks, there is still no effective treatment for this disease, and only one licensed vaccine is available. However, as some seronegative recipients have developed clinical signs reminiscent of antibody dependent enhancement of DENV infection leading to severe dengue, this vaccine is limited to those with proof of a previous dengue infection [5]. Due to the lack of preventative interventions, the current best practices for DENV prevention rely heavily on mosquito vector control, although this approach is insufficient to completely prevent DENV transmission [6-11].
Though Ae. aegypti are generally considered susceptible to infection with and capable of transmitting DENVs (i.e., competent vectors), field-caught and laboratory bred Ae. aegypti exhibit considerable variation in vector competence for DENV. This variation is influenced by virus serotype and strain, as well as mosquito genetics, mosquito metagenome, and environmental factors [12-22]. As such, the overwhelming number of estimates of Ae. aegypti DENV competency derived from DENV-2 laboratory strain New Guinea C (NGC) [12]) and of mosquito laboratory colony infections lack the granularity and biological complexity to inform the vector competency for other DENV serotypes and strains in other mosquito populations. This narrow insight limits our understanding of potential serotype- and strain-specific interactions that could influence real-world scenarios. As such, it is prudent to test the vector competence of local mosquitoes from a geographic region of interest to evaluate the focally relevant risk for DENV transmission. In Thailand, Ae. aegypti from areas with high human dengue incidence were shown to have higher vector competence than Ae. aegypti from regions of low dengue incidence [23]. However, in Cuba, Ae. aegypti sampled from areas categorized as low- and high-risk based on long-term trends in local human dengue case data demonstrated no significant differences in vector competence for a local circulating strain of DENV-1, with transmission rates ranging from 10 – 25% [24]. Therefore, to understand the risk that dengue poses to a community it is important to prioritize the comprehensive examination of vector competence of local vectors collected across spatially risk-partitioned zones.
Having previously broken the cycle of dengue from the United States (USA) through insecticide use and source reduction, the state of Florida (FL) is vulnerable to the establishment of all four serotypes of DENV in mosquito populations leading to subsequent local transmission in human populations [6]. Ae. aegypti mosquito populations have resurged following widespread displacement by Ae. albopictus in the mid-1990’s and are now present throughout much of the state [25]. In addition, hundreds of travel-associated dengue cases are recorded each year in FL, and these correlate with sporadic local transmission, acting as seeding events [26]. Caribbean countries, particularly Haiti, are a key source of imported DENV in FL [27, 28]. Nearly all local transmission of DENV in the state has occurred in its southern half with the majority of cases in the south-ernmost counties (Co.) [26, 29]. There were 28 locally acquired DENV-1 cases in Monroe Co. in 2009 and 65 cases in 2010, as well as 29 cases in Martin Co. in 2013 and 71 cases in Monroe and MiamiDade Co.’s in 2020 [27, 30-33]. To date, vector competence studies for DENV in FL have been completed for Ae. aegypti from two counties, Monroe (Key West) and Indian River (Vero Beach), and these studies used DENV-1 strains from either Key West or Puerto Rico (PR), or a DENV-2 strain from PR, respectively [15, 29, 34]. In spite of sporadic local transmission of DENV in FL during the past twelve years [26], little proactive work has been done to quantify and characterize vector competence of Ae. aegypti mosquito populations from different regions of the state.
It has become increasingly clear that insecticide resistance can influence vector competence. Pyrethroids are commonly used insecsticides for controlling both nuisance mosquitoes and disease vectors globally, and pyrethroid resistance in Ae. aegypti is a significant issue in FL [35]. Pyrethroid use is widespread globally and insecticide resistance is prevalent but variable throughout the state of Florida [36]. There are two common nonsynonymous mutations to the voltage-gated sodium channel transmembrane protein gene in Ae. aegypti resulting from a 1016 valine to isoleucine and/or a 1534 phenylalanine to cysteine change; these mutations are known as knockdown resistance (kdr), and they reduce pyrethroid activity through target-site insensitivity. It was shown that a pyrethroid-resistant Ae. aegypti strain with kdr mutations selected from mosquitoes collected in Monroe Co., FL had a higher DENV-1 dissemination rate compared to its unselected, susceptible counterpart [37]. A study of a laboratory colony of Ae. aegypti from Orlando, FL (ORL) with both pyrethroid-susceptible and resistant phenotypes found that the pyrethroid-resistant mosquitoes had higher disseminated Zika virus titers [38]. Although the mechanism by which pyrethroid resistance-conferring loci influence vector competence is unclear, the pyrethroid resistance status of field-derived mosquitoes should also be considered in estimating vector competence and arbovirus transmission risk.
Herein, we examined two facets of vector competence, susceptibility (midgut infection) and transmission potential (virus detection in saliva) to understand if Ae. aegypti from various DENV risk zones (defined based on local DENV case counts, low: north, medium: central, and high: south) in FL (Figure 1) have varying levels of competence. We assessed competency across DENV 1-4, using four laboratory (L) stock viruses (DENV-1 L, DENV-2 L, DENV-3 L, DENV-4 L), and two low-passage wild type (WT) strains (DENV-1 WT and DENV-4 WT). We also measured and compared infection rates and virus titers in mosquito tissues (infection intensities) between pyrethroid-susceptible ORL mosquitoes and Ae. aegypti populations from the aforementioned risk zones. These data were used to understand how intra- and inter-serotype variation may affect vector competence and consequently DENV transmission risk in the state.
Aedes aegypti collection locations for four field-derived colonies and one laboratory colony in Florida, USA., with reported kdr genotype frequency distributions shown as part-of-whole graphs.
2. Materials and Methods
2.1. Mammalian Cell Culture and Virus Propagation
Mammalian cell culture and virus propagation methods for this experimental design were conducted similarly to previously published methods [17]. Briefly, for the preparation of virus stocks, Vero E6 cells grown as monolayers were seeded onto the 75 cm2 growing surfaces of filter-cap T75 cell culture flasks and incubated until they attained 80% confluency. The cells were then individually inoculated with DENV-1 through DENV-4 laboratory stock viruses (DENV-1 L, DENV-2 L, DENV-3 L, DENV-4 L) obtained from the American Type Culture Collection (ATCC.org) as well as a DENV-1 wild-type (DENV-1 WT) isolate from a child in Haiti (strain: Haiti/1207/2014: KT279761.2) and a DENV-4 wild-type isolate also from a child in Haiti (strain: Haiti/0075/2015: MK514144.1) (Table A1). The virus stocks were collected between 7- and 9-days post-inoculation when approximately 50% of the cells showed virus-specific cytopathic effects (CPE), and they were maintained and tittered as described in full previously [17].
2.2. Risk Zone Classification
We chose to test Ae. aegypti from four counties to represent present day low, medium, or high risk areas for dengue virus transmission in FL. St. Johns Co. in north FL (Figure 1) was classified as “low risk” because no dengue-positive mosquito pools nor local human cases have been reported from this county since dengue cases returned in FL in 2009 after over fifty years without local cases [39]. Manatee Co. in southwest FL was classified as “medium risk” bescause there was a report of dengue-positive Ae. aegypti (DENV-4) in 2016 and 2017, but no previous local dengue transmission was detected since 2009 [40]. Miami-Dade Co. in south FL was classified as “high risk” because there have been reports of both dengue-positive mosquito pools (unpublished) and human dengue cases in multiple years since 2009 [30, 33]. Collier Co. was added as an additional “low risk” county (with no dengue-positive mosquito pools nor human cases) to act as a “low risk” south FL comparator for Miami-Dade Co., since both counties are further south in the sub-tropical zone, and reflect climates with warmer temperatures and higher relative humidity on average than St. Johns Co. in northeastern FL [41, 42].
2.3. Field-caught Ae. aegypti
We obtained field-caught G0 Ae. aegypti eggs from: Anastasia Mosquito Control District in St. Johns Co., Collier Mosquito Control District (Collier Co.) and Miami-Dade Mosquito Control District (Miami-Dade Co.). We also received Manatee Co. Ae. aegypti eggs (F7) from Dr. Derrick Mathias at the Florida Medical Entomology Laboratory at the University of Florida. Upon receipt, the eggs were hatched, and adult mosquitoes were cold anesthetized and morphologically verified as Ae. aegypti by expert identification. F2 generation Ae. aegypti from Miami-Dade (MD) and St. Johns Co’s (SJ) were used for the following studies. The F7 generation of Ae. aegypti from Manatee Co. (MAN) was used. The Collier Co. (COL) comparator group was obtained from the field in 2018 and was reared in the laboratory since 2018 (exact generation unknown). The field colonies (homogenized pooled legs from F1 generation of MD and SJ, F7 generation of MAN and unknown generation of COL) were screened via RT-PCR for the presence of Cell Fusing Agent virus (CFAV) (Table A2), a potentially confounding insect-specific virus (ISV); all mosquito pools tested negative.
Ae. aegypti Orlando strain (ORL) was obtained as adults from the United States Department of Agriculture’s Center for Medical, Agricultural and Veterinary Entomology (USDA-CMAVE) in Gainesville, Florida and then established as a colony in our laboratory. Mosquitoes were reared in growth chambers set to 28°C with 80% relative humidity and a neutral photoperiod regimen (12 h light/12 h dark).
2.4. Mosquito Rearing, Infection and Tissue/Specimen Preparations
Four to seven days post-eclosion, female mosquitoes were separated and used in per os DENV infections. These females were artificially fed with a 2:2:1 human hematocrit type O+ blood (hematocrit) (Lifesouth Community Blood Centers, Gainesville, FL): virus stock: heat-inactivated human serum (HIHS). After one hour, blood-engorged mosquitoes were kept for further analysis and all other mosquitoes were discarded. We maintained the mosquitoes for 14 days on a 10% sucrose solution. On day 14, we cold anesthetized the mosquitoes, removed their legs and wings, and manually inserted their proboscis into a capillary tube (Drummond, Broomall, PA) filled with 3 µL of 1:1 hematocrit: HIHS [17]. Each proboscis was withdrawn after 45 minutes or after 2 µL of hematocrit was imbibed. The remaining contents of the capillary tube were expelled into microcentrifuge tubes. Next, the midgut of each mosquito was dissected and individually stored in microcentrifuge tubes in DMEM, and we recorded which mosquitoes had blood in their bodies from the salivation assay. These experiments were repeated in triplicate for each of the six virus groups across all four field-caught mosquito strains and the ORL laboratory strain.
2.5. Mosquito Midgut and Saliva Analysis
BHK-21 cells (Baby hamster kidney fibroblast cells, ATCC® CCL-10™) (a generous gift from the Dimopoulos laboratory at Johns Hopkins University) were seeded onto 24-well plates at a density of 5×104/well and incubated at 37°C with 5% CO2 until nearly confluent. Each mosquito midgut was homogenized using a Bullet Blender (Nextadvance, NY) adjusted to speed setting 8 for 3 minutes and serially diluted 10-fold before adding 100 µL of each dilution series to the 24-well plate. We used 0.8% w/v methyl cellulose media for our overlay and incubated the plates for six days (except for DENV-4 that was incubated for 5 days due to quicker development of plaques). The plates were fixed in a 1:1 methanol/acetone solution with 1% w/v crystal violet for an hour, after which plaques were counted, and virus titer expressed as plaque forming units per mL (PFU/ mL).
The DENV-1 laboratory strain (DENV-1 L) did not clearly form plaques on BHK-21 or Vero E6 cells after numerous attempts. Cytopathic effects could be visualized under a phase contrast inverted microscope when the virus was inoculated onto Vero E6 cells. Therefore, we seeded 48-well plates with Vero E6 cells at a density of 3×104/ well, and upon reaching 70% cell confluency we individually inoculated each well with 50 µL of midgut homogenate from ORL mosquitoes exposed to DENV-1 L. The cultures were monitored for a total of 27 days wherein media was replenished, and inoculum re-passaged as needed or when cells became overly confluent; this was done to ensure that CPE could be detected even if starting virus concentration was low. The readout for positive midgut samples was a dichotomous qualitative outcome of “infected” or “not infected” based on presence or absence of CPE compared to the mock-inoculated (with DMEM) wells grown in parallel conditions throughout the twenty-seven days.
RNA was purified from saliva samples from ORL mosquitoes that had a DENV-positive midgut, either by plaque assay (DENV-1 WT, DENV-2 L, DENV-3 L, DENV-4 WT and DENV-4 L) or by cell culture (DENV-1 L) using a QIAamp Viral RNA Mini Kit (Qiagen, Valencia, CA, USA). Each RNA sample was subjected to rtRT-qPCR using a dual-target pan-DENV system (Table A2) (Superscript III, Invitrogen) and run as technical duplicates on a BioRad CFX96 Touch Real-Time PCR Detection System. A standard curve of each virus was run for each assay to determine the virus-specific regression analysis for estimating PFUe / mL.
2.6. Statistical Analyses
All statistical analyses were performed using SPSS statistical software (IBM Corp. Released 2020. IBM SPSS Statistics for Windows, Version 27.0. Armonk, NY: IBM Corp). Data for ORL and COL exposed to DENV-4 WT and DENV-4 L are reanalyzed from Stephenson et al., 2021 [44] to provide a complete analysis and comparison between all virus groups and mosquito populations [17]. In all tables and figures, groups bearing the same alphabetical letter were determined to be statistically similar in their values (infection rates or infection intensities), whereas groups bearing different alphabetical letters were determined to have statistical differences in their values.
To examine infection rates, a multilevel logistic regression was used with one fixed effect predictor; the binary response variable, midgut infection status, “non-infected = 0” or “infected = 1”, and the fixed variable included each strain of mosquitoes exposed to one of the six virus groups (e.g., St. Johns Ae. aegypti exposed to DENV-1 WT). Sample sizes are listed in Table A3. We evaluated the potential for “batch effect” by including experimental replicate number as a random effect variable. The random effect predictor was significant (P = 0.027) (as such, it was incorporated as a factor in the full model). We report odds ratios (ORs) and P values for each pairwise combination of IRs. Significant statistical significance (α = 0.05) is denoted via groups assigned different letters.
The transmission potentials (TPs) of field mosquitoes were analyzed via binary logistic regression, with sample sizes listed in Table A3. The binary saliva infection status was “non-infected = 0” or “infected = 1”. There was not a significant batch effect for saliva samples, so we moved forward with a univariate analysis of saliva-positivity rates from field mosquitoes exposed to one of six different DENVs. We report ORs and P values for each pairwise combination of TPs. Statistical significance (α = 0.05) is denoted via groups assigned different letters.
Average titers in log10 PFU / mL (transformed for normality) were compared using a mixed methods ANOVA between each virus group and field mosquito population for both midgut and saliva samples. We included batch effect as the random effect variable to control for replicate effects. The fixed effect was a combination of virus strain and mosquito strain. Both midgut and saliva residuals for log10 titer passed the Shapiro-Wilk test for normality and the outcome variable passed Levene’s test for equality of variances.
2.7. kdr Genotyping
A previously published qPCR assay for the detection of two kdr single nucleotide polymorphisms (SNPs) coding for amino acids located at positions 1016 and 1534 of the NAV gene was adapted for this work (Table A3) [36, 43, 44]. Forty eggs from each field mosquito colony were homogenized [COL (generation unknown), MAN (F8), MD (F3) and SJ (F3)]. Fifteen individual mosquito midgut samples (COL and MD) or 15 individual whole bodies (MAN and SJ) were also homogenized to compensate for desiccation of some eggs that likely resulted in no signal for the assay. Final sample sizes were n = 31 for COL, n = 30 for MAN, n = 29 for MD and n = 31 for SJ. ORL was used as the pyrethroid-susceptible control and Puerto Rico (PR) Ae. aegypti (a kind gift from the USDA-CMAVE), as the pyrethroid-resistant control for this assay. We also included a heterozygote control that contained one ORL and one PR mosquito. The Applied Biosystems™ SYBR™ Select Master Mix for CFX was used with cycling conditions of 50°C for 2 minutes, 95°C for 2 minutes and 45 cycles of 95°C for 15 seconds and 60°C for thirty seconds. There was a final melt curve with a ramp from 60°C to 95°C with reads every 0.3°C and genotypes were determined by melt curve analysis. For the 1016 SNP, amplicons from a pyrethroid-susceptible mosquito had a melting temperature at 85.6°C, amplicons from resistant mosquitoes had a melting temperature of 77°C, and heterozygotes had a peak at both temperatures. For the 1534 SNP assay, amplicons from a susceptible mosquito had a melting temperature of 79.3°C and amplicons from a resistant mosquito had a melting temperature of 84.4°C. Heterozygotes had two peaks at both previously listed temperatures. We used genotype information from each mosquito per population to understand the distribution of genotype frequencies across each colony and visualized these results using pie charts.
3. Results
3.1. Serotype-specific vector competence of the Ae. aegypti (ORL) colony
We developed a vector competence baseline for Ae. aegypti ORL mosquitoes to refine the community’s perception for the utility of this colony in such studies. The ORL colony used herein have been continuously reared at the USDA (Gainesville, FL) since 1952, is pyrethroid susceptible and is not infected with cell-fusing agent virus (CFAV), an insect-specific virus that has been shown to influence arbovirus infections of mosquito cells [22]. We performed pairwise comparisons of midgut infection rates (IRs) to evaluate susceptibility to infection (Figure 2 and Table 1). IRs represent the average number of DENV-infected midguts out of the total number of exposed mosquitoes per group, across three biological replicates. In pairwise comparisons of DENV-exposed ORL mosquitoes, midgut IR was highest for DENV-1 WT with 97%, followed by DENV-4 WT (61%), DENV-2 L (28%), DENV-4 L (21%), DENV-3 L (20%) and DENV-1 L (10%). The odds of midgut infection were significantly higher for ORL mosquitoes exposed to DENV-1 WT compared to DENV-1 L (Odds ratio (OR): 219.9, P value (P) < 0.001), DENV-2 L (OR: 56.2, P < 0.001), DENV-3 L (OR: 120.4, P < 0.001), DENV-4 WT (OR: 14.5, P < 0.001) and DENV-4 L (OR: 79.5, P < 0.001). There were also higher odds of midgut infection for DENV-4 WT compared to DENV-1 L (OR: 15.1, P < 0.001), DENV-2 L (OR: 3.9, P = 0.001), DENV-3 L (OR: 8.3, P < 0.001) and DENV-4 L (OR: 5.5, P < 0.001). Lastly, DENV-2 L had higher odds of midgut infection than DENV-1 L (OR: 3.8, P = 0.012).
Pairwise comparisons DENV infection rates and transmission rates of Florida Aedes aegypti populations reveal virus-specific and mosquito population-dependent differences
Infection rates and transmission rates of Floridan Aedes aegypti mosquitoes for DENV-1 through DENV-4 are highest for south FL (COL and MD) and laboratory (ORL) populations with wild type virus strains (DENV-1 WT and DENV-4 WT). Capital alphabetical letters denote statistical significance between infection rates and lowercase letters represent statistical significance between transmission rates. Bars with different letters are significantly different from one another with a P - value < 0.05 via multilevel logistic regression (infection rates) or binary logistic regression (transmission rates).
To assess the transmission potential, we performed pairwise comparisons of virus titers in saliva, i.e., transmission rates (TR), which is the total number of mosquitoes with DENV detected in their saliva out of the total number of virus-exposed mosquitoes (Figure 2 and Table 1). We identified significantly higher TRs for DENV-1 WT compared to DENV-2 L (OR: 10.6, P = 0.004), DENV-3 L (OR: 11.4, P = 0.029), DENV-4 WT (OR: 6.0, P < 0.001) and DENV-4 L (OR: 7.3, P = 0.018).
3.2. Wild type DENV infections have elevated midgut and saliva infection intensities in Ae. aegypti ORL mosquitoes
We compared infection intensities among ORL infected with DENVs and saw higher titers for wild-type viruses (Figure 3 and Table 2). We observed higher mean midgut titers for mosquitoes infected with DENV-1 WT compared to DENV-2 L (P = 0.044). The average titer for midgut tissues infected with DENV-4 WT was significantly higher than titers for DENV-1 WT (P = 0.049), DENV-2 L (P < 0.001) and DENV-4 L (P = 0.038). The average saliva infection intensity of ORL for DENV-1 WT was significantly higher than DENV-1 L (P = 0.024) and DENV-2 L (P = 0.035). Similarly, the saliva infection intensity for DENV-4 WT was significantly higher than DENV-1 L and DENV-2 L (P = 0.034 and 0.042).
Average DENV titer (log10 PFU / mL) across Floridian Aedes aegypti reveals higher titers for the laboratory colony (ORL) as well as increased midgut titers for DENV-4 wild type
Infection intensities for Floridan Aedes aegypti are elevated for wild type DENV strains (DENV-1 WT and DENV-4 WT) compared to laboratory virus strains. Capital alphabetical letters denote statistical significance between midgut infection intensities and lowercase letters represent statistical significance between saliva infection intensities. Groups with different letters are significantly different from one another with a P - value < 0.05 via mixed methods ANOVA.
3.3. Field-derived Ae. aegypti from south FL have comparatively greater vector competence
We compared IRs for each DENV risk group: Ae. aegypti from “low risk” St. Johns Co. (SJ) in north FL, “low risk” Collier Co. (COL) in south FL, “medium risk” Manatee Co. (MAN), and “high risk” Miami-Dade Co. (MD). There were significantly higher odds of midgut infection for MD exposed to DENV-4 WT compared to both MAN (OR: 3.0, P = 0.024) and SJ (OR: 6.4, P < 0.001). There were also significantly higher odds (OR: 2.8, P = 0.035) of midgut infection for MD Ae. aegypti exposed to DENV-4 L compared to SJ. When comparing the two south FL mosquito groups, COL and MD had similar midgut IRs across all virus pairings.
Transmission potential was higher for MD Ae. aegypti exposed to DENV-1 WT (42%) compared to MAN (17%) (OR: 4.8, P = 0.012). There was also significantly higher DENV-4 WT transmission potential for both MD (44%) and MAN (27%) compared to SJ (6%). MD were 8.2 times more likely to have DENV-4 WT detected in their saliva compared to SJ, and MAN were 5.9 times more likely to have DENV-4 WT detected in their saliva compared to SJ. There were no significant differences among the TRs of laboratory viruses between SJ, MAN or MD. When comparing the two south FL mosquito groups, COL had higher DENV-4 L transmission potential (OR: 15.0, P = 0.03) compared to MD but the remaining comparisons were similar.
3.4. Higher infection and transmission rates for wild type DENVs in Floridan field Ae. aegypti
Overall, we found higher IRs and TRs for the wild type viruses compared to laboratory viruses. For example, COL had higher odds of DENV-4 WT midgut infection compared to every other virus group for COL (DENV-1 WT, OR: 6.7, P < 0.001; DENV-1 L, OR: 80.3, P < 0.001; DENV-2 L, OR: 3.7, P = 0.008; DENV-3 L, OR: 23.4, P < 0.001; DENV-4 L, OR: 9.4, P < 0.001). MD had a significantly higher odds of DENV-1 WT transmission potential compared to DENV-2 L (OR: 4.8, P = 0.012), DENV-3 L (OR: 5.07, P = 0.038) and DENV-4 L (OR: 32.3, P < 0.001).
3.5. DENV-4 wild type attained highest titers in field Ae. aegypti midgut samples but not in saliva
There were no significant differences in infection intensities for SJ Ae. aegypti. MAN had significantly higher average DENV-4 WT midgut titers, compared to DENV-1 WT (P = 0.002) and DENV-4 L (P = 0.034), but DENV-4 L reached higher average saliva titers compared to DENV-4 WT (P = 0.05). MD had the highest average DENV-4 WT midgut titers in comparison to DENV-1 WT (P = 0.008), DENV-2 L (P = 0.008) and DENV-3 L (P < 0.001). COL midgut infection intensity was higher for DENV-4 WT than both DENV-1 WT (P = 0.045) and DENV-2 L (P < 0.001). Even with higher DENV-4 WT midgut titers observed across most of the mosquito populations, DENV-4 WT saliva titers were comparable or lower than other virus groups for field Ae. aegypti.
3.6. Higher infection intensity for laboratory Ae. aegypti (Orlando) compared to Floridan field Ae. aegypti
ORL reached the highest average midgut titers compared to COL, MAN, MD and SJ, except for DENV-2 L titers for MAN. ORL had significantly higher midgut titers for DENV-1 WT compared to SJ (P = 0.002), MD (P < 0.001) and MAN (P < 0.001). Additionally, ORL had significantly higher midgut titers for DENV-3 L compared to MD (P < 0.001), and for DENV-4 WT compared to SJ (P < 0.001) and MD (P < 0.001). ORL also had higher DENV-4 L midgut titers compared to SJ (P = 0.039). Finally, ORL had higher DENV-4 WT titers in saliva compared to COL (P = 0.05) and MAN (P = 0.005).
3.7. Knockdown resistance allele profiles in field-derived Ae. aegypti insversely trends with vector competence for DENV
The frequency distributions of kdr alleles are presented in Figure 1 and Table A4. SJ had the highest resistance profile, with 100% of tested mosquitoes bearing both resistance alleles (IICC). MD had the second highest percentage of resistant mosquitoes (38% IICC), but had wide variability across kdr genotypes, and did not have any fully pyrethroid susceptible mosquitoes. COL had 16% IICC mosquitoes, and 3.2% of mosquitoes carried the second most resistant genotype (IIFC). COL had the most variability in genotypes of the four groups tested and had a similar distribution to MD. Lastly, MAN had 13% IICC mosquitoes, 50% VICC (heterozygote for the 1016 allele and homozygous resistant for the 1534 allele), and 37% VVCC (homozygous susceptible for the 1016 allele and homozygous resistant for the 1534 allele). The field mosquito populations with higher pyrethroid resistance, based on their kdr genotypes, had the lowest measured vector competence. ORL is 100% susceptible at both loci and had the highest vector competence measures, followed by MD with 38% of mosquitoes bearing resistance genotypes. However, the SJ population had 100% resistance at these loci and was found to have the lowest competence for DENVs.
4. Discussion
We provide evidence that field Ae. aegypti populations from different DENV risk zones in FL vary in their vector competence profiles between the four DENV serotypes, and between WT and laboratory virus strains. We observed a gradient of increasing competence from “low risk” SJ through “high risk” MD. Overall, MD had the highest IRs, second only overall to the ORL colony, and highest susceptibility and transmission potential for both WT strains compared to the laboratory viruses. In general, the highest IRs, TRs, and midgut titers occurred with either of the WT viruses (DENV-1 WT and DENV-4 WT) compared to the laboratory viruses (DENV-1 L, DENV-2 L, DENV-3 L and DENV-4 L), even though all, except DENV-1 L (105 PFU / mL), had comparable blood meal titers (106 PFU / mL). The DENV-1 L strain would not propagate to higher titers than 105 PFU / mL regardless of cell line or duration of culturing in our hands. It is known that infectious dose is a contributing factor to the establishment of a midgut infection [45], but it is evident that the differences we identified in susceptibility occurred regardless of infectious dose standardization.
The four DENV laboratory stocks used in this study were previously found to share the same putative receptors in Ae. aegypti midguts [46]. Little is currently known about whether polymorphisms in these midgut receptors and associated proteins influence epithelial cell infections for field and laboratory DENVs. In humans, there is evidence of altered susceptibility of cell types by DENV-2 strains for both laboratory and field isolates [47]. In mosquitoes, similar DENV-2 midgut binding affinity was observed across various field isolates that had altered vector competence outcomes in Ae. aegypti [48]. Mutations that arise in arbovirus genomes can enhance virus dissemination and vector competence [49, 50]. Previous work identified several DENV-2 mutants that bear amino acid changes that enhance virus dissemination out of the midgut [49]. Genetic differences among our virus strains, particularly to the envelope gene, could alter their midgut binding and cell tropism. At the nucleotide and amino acid level, DENV-2 and DENV-4 are the most divergent among the four serotypes, with DENV-1 and DENV-3 being most similar [51-53]. In particular, DENV-4 WT and DENV-4 L share 93% nucleotide identity and 98% amino acid identity but the DENV-4 L sequence has an additional 15 nucleotides in its 3’ UTR not seen in contemporary strains within genotype IIb [40]. We suspect that certain genetic changes arising from continuous cell culture of the DENV laboratory stocks could be responsible for their lower fitness in Ae. aegypti. Now that serotype and strain-specific trends in vector competence have been established for both laboratory and field Ae. aegypti, as well as across laboratory and field strains of DENV, future work is needed to uncover the driving factors responsible for the phenotypic differences.
Several studies previously classified ORL mosquitoes as “refractory” compared to the Rockefeller strain of Ae. aegypti, but we suspect that this misclassification is likely a result of suboptimal vector-virus combinations that were used, especially since a labmaintained, historical DENV-2 NGC strain was primarily tested [54, 55]. An additional factor at play is that colonies referred to as “ORL” in the literature but from different labs have had different evolutionary histories post sub-colonization from the original ORL colony, which undoubtedly influences the degree of genetic similarity among colonies. In our study, ORL generally had higher susceptibility and competence measures, particularly for the wild type DENV strains, compared to field Ae. aegypti. The highest DENV titers were reported from ORL midgut samples compared to every other mosquito strain, except for DENV-2 titers for MAN. Additionally, DENV titers from ORL midgut and saliva samples were higher than several other mosquito strains for the various virus combinations. In the published literature, ORL has been reported to have a median DENV-2 midgut PFU of zero [56] and has also been found to have a DENV-2 midgut prevalence of zero in comparison to the laboratory colony known as Rockefeller [55]. Those studies had comparable or higher DENV-2 infectious doses (106 or 107 PFU) to our study, yet we measured an average and median DENV-2 concentration around 102 PFU / midgut with an infection rate of 28% and transmission rate of 4%. Due to higher DENV competence of ORL, compared to FL field Ae. aegypti, we conclude that it would not serve well as a colony to use to ascertain risk of DENV transmission in a real-world setting but that it could serve as a useful model to study in vivo vector-virus interactions when high infection rates of WT virus are needed. It appears that ORL has been previously discounted as a refractory colony, but our results show vector competence for ORL to be more nuanced and highly virus strain-dependent.
Our work highlights the potential system for DENV-4 persistence in the mosquito population in FL. We previously identified and fully sequenced a DENV-4 strain with high identity to our DENV-4 WT from Haiti in Ae. aegypti from Manatee Co. in 2016 and 2017, in the absence of a human case across both years [40]. It is possible that there were unreported asymptomatic dengue cases [57], but this cannot be known with certainty unless paired with serology studies. However, DENV-3 maintenance by vertical transmission in Ae. albopictus from Brazil occurred in a time period without autoch-thonous transmission [58]. Additionally, experimental infection of field-acquired Ae. aegypti has shown efficient vertical and venereal transmission of DENV-2 during second and third egg laying cycles, especially when virus dissemination was high [59]. MAN had significantly higher DENV-1 WT and DENV-4 WT IRs compared to DENV-2 L, DENV-3 L and DENV-4 L, but did not have significantly higher TRs than those same groups. The high midgut IR and infection intensity of DENV-4 WT in MAN, and an unremarkable TR or saliva infection intensity suggest that vertical transmission could have been favored to maintain this virus in nature. In support of this hypothesis, our research group recently showed that an in-fectious clone of the DENV-4 sequence found in Manatee Co. Ae aegypti is capable of being both horizontally and vertically transmitted [60], and we also found DENV-4 WT to be vertically transmitted at a higher rate than DENV-4 L in ORL [17].
We noted that the presence of both the 1016 and 1534 pyrethroid-resistant kdr SNPs were associated with the most refractory, lowest DENV susceptible, phenotypes of the field lines. This trend is opposite to early reports showing insecticide resistance increases virus dissemination and competence [37, 38]. Those laboratory studies differ from ours in that one study actively selected for phenotypic pyrethroid resistance in a field Ae. aegypti colony and then evaluated DENV-1 dissemination in mosquito bodies of F13 and F20 generations [37]. The other study backcrossed ORL to obtain both a pyrethroid-susceptible and a pyrethroid-resistant colony, in which they evaluated ZIVK competence [38]. Our study used lower generation field-derived mosquitoes (F2 and F7) and did not select for phenotypic pyrethroid resistance. Without direct assessment of resistance phenotype, we cannot say definitively that possessing kdr alleles conferred insecticide resistance or that resistance is correlated with higher vector competence. However, it is still important to consider that in nature, mosquitoes from these collection sites are likely exposed to pyrethroids from vector control mediated by a combination of publicly funded control programs and private pesticide applicators. Considering, within mosquito population structuring has been observed previously [33], additional studies are needed to comparatively measure other components of vectorial capacity, such as differences in host-seeking behavior or locomotor activity, longevity, and vector competence across multiple mosquito populations between regions, e.g., St. Johns Co. vs. Miami-Dade Co.
We show that DENV serotype-specific variation in vector competence across four field Ae. aegypti populations from FL. Vector competence was highest for the two south FL mosquito groups, from Miami-Dade and Collier Co’s and for the low passage wild type strains (DENV-1 WT and DENV-4 WT). These observations emphasize the importance of choosing geographically and medically relevant mosquito and virus strains in vector competence studies to estimate more accurately DENV risk. Moreover, these data offer practical information to pinpoint areas and human populations at increased DENV risk. We expect that allocating resources proactively towards highly targeted vector control efforts in FL that have Ae. aegypti with high DENV transmission potential can prevent future or stop dengue outbreaks in FL.
Data Availability
All data are available in the manuscript
Author Contributions
C.J.S, J.A.L, R.R.D conceived the study design. C.J.S., H.C. conceived and planned experiments. C.J.S., H.C. C.M.W. and N.M.N-M. conducted experiments and interpreted results. J.A.L, R.R.D and D.K.M. provided oversight and supervision for experiments, data analysis and interpretation. R.R.D. and D.K.M. provided resources for this study. R.R.D. provided funding for the study. C.J.S. wrote the manuscript with editing and revision support from H.C., C.M.W., N.M.N-M., D.K.M., J.A.L. and R.R.D.
Funding
This research was supported in part by the United States Centers for Disease Control (CDC) Grant 1U01CK000510-03: Southeastern Regional Center of Excellence in Vector-Borne Diseases: The Gateway Program. The CDC did not have a role in the design of the study, the collection, analysis, or interpretation of data, nor in writing the manuscript. Support was also provided by the University of Florida Emerging Pathogens Institute, the University of Florida Preeminence Initiative through the UF College of Veterinary Medicine and the University of Florida Graduate Student Funding Award.
Institutional Review Board Statement
Not applicable.
Informed Consent Statement
Not applicable.
Data Availability Statement
Data can be made available upon request.
Conflicts of Interest
The authors declare no conflict of interest.
Acknowledgments
We thank the United States Department of Agriculture Center for Medical, Agricultural and Veterinary Entomology in Gainesville, Florida for providing Ae. aegypti (Orlando strain) and to Dr. Keira Lucas and the Collier Mosquito Control District, Dr. Derrick Mathias, Dr. Eva Buckner, and the Manatee Mosquito Control District, Dr. Rudy Xue and the Anastasia Mosquito Control District, and Chalmers Vasquez and the Miami-Dade Mosquito Control District for providing field-collected Ae. aegypti eggs. We also thank Drs. Laurence Thirion and Remi Charrel from Unite des Virus Emergents, Aix Marseille University, France for providing primers and probes for the pan-dengue RT-PCR assay. C.M.W. is a contracted employee of the U.S. Government, and the views expressed in this manuscript are those of the authors and do not necessarily reflect the official policy or position of the Department of the Navy, Department of Defense, nor the U.S. Government, this work was prepared as part of her official duties. Title 17, U.S.C., 105 provides that copyright protection under this title is not available for any work of the U.S. Government. Title 17 U.S.C., 101 defines a U.S. Government work as a work prepared by a military Service member or employee of the U.S. Government as part of that person’s official duties.
Appendix A
Dengue virus (DENV) strains and sources used in experiments
Primer and probe sequences for pan-DENV, CFAV and kdr allele detection in Aedes aegypti.
Sample sizes per experimental group of Aedes aegypti exposed to dengue viruses
Knockdown resistance (kdr) genotype percentages for Florida Aedes aegypti field populations (St. Johns, Manatee, Miami-Dade, Collier) and a laboratory colony (Orlando; ORL)
Footnotes
Citation: Stephenson, C.; Coatsworth, H.; Waits, C.; Nazario-Maldonado, N.; Mathias, D.; Dinglasan, R.; Lednicky, J. Geographic partitioning of dengue virus transmission risk in Florida. Viruses 2021, 13, x. https://doi.org/10.3390/xxxxx