INTRODUCTION
More than 400 species have been described for the genus Anopheles, with approximately 40 regarded as vectors of interest,1 mostly of Plasmodium parasites. Anopheline mosquitoes are primarily known for transmitting malaria,2 but they also have the potential to transmit viruses.3,4 In general, they are highly mobile and thrive by using human activities and movement to disperse around the globe.5,6 Although the most known and well-studied species are the African Anopheles gambiae and the Indo-Iranian Anopheles stephensi, due to the malaria-associated socioeconomic and health burdens they cause in those regions,2 less well studied anopheline species predominate in other areas of the planet with the potential to spread different pathogens.
Anopheles albimanus is the main anopheline inhabiting northern South America, Central America (reviewed in detail in reference 7), and the Caribbean islands.8 Its broad geographical distribution may be aided by the species’ ability to survive in both fresh9,10 and brackish water.11,12 Although it remains incompletely understood, An. albimanus has been described as a much more zoophilic, crepuscular, and exophagic mosquito with higher adaptive capabilities than other anopheline species,13–16 all of which affect the species’ success in transmitting Plasmodium. However, host availability and environmental conditions appear to influence its host choice17,18 and resting behavior.19 The flexible behavior of this species may be facilitating its spread into the southern United States, lower latitudes (more temperate areas) in South America, and urban and peri-urban settlements where it encounters human hosts. Its expansion may also be helped by climate change, which is broadening the species’ geographical habitat range.8 Despite the work that has been done to characterize its ecology and behavior, little is known about the capacity of the species to harbor and transmit classic and emerging tropical mosquito-borne viruses.
Many arboviruses produce similar disease symptoms in humans that include fever, headache, rash, diarrhea, and joint pain, which can last for months. Because treatment is not specific to the etiological agent and neither are many clinical surveillance and diagnostic tools, the prevalence of emerging viruses can be misdiagnosed and hence underestimated in areas with more common viral outbreaks such as Chikungunya or dengue viruses, as initially occurred with the Zika virus epidemic.20 However, despite displaying similar clinical symptoms, viruses may differ in intrinsic replication rates or use different cellular receptors in the mosquito to achieve a successful human-to-mosquito-to-human viral transmission route. Arboviruses rely on the rapid infection of a mosquito after feeding on an infectious host, and must penetrate and overcome multiple tissue and immune barriers to propagate throughout the body and reach the salivary glands.21 The virus must also replicate in the salivary glands efficiently to later infect a naive vertebrate host through salivation during a second bloodmeal. To date, only one arbovirus is known to be primarily transmitted through the bite of Anopheles mosquitoes in the field (O’nyong-nyong virus22; Anopheles funestus and An. gambiae), but other Anopheles species have been shown to be capable vectors of alphaviruses in the laboratory.23,24
Here, we ask whether An. albimanus is a competent vector of arboviruses. We orally challenged adult females with infectious bloodmeals containing one of the following togaviruses (genus Alphavirus): Mayaro (MAYV; -D and -L genotypes), Chikungunya (CHIKV; Asian lineage), O’nyong-nyong (ONNV), and Sindbis virus (SINV), as well as the flavivirus dengue virus serotype 2 (DENV-2; Asian lineage), due to its prevalence and health impact worldwide. Collectively, these togaviruses produce human disease that spans the planet (Figure 1A): MAYV in Central/South America, ONNV in Africa, CHIKV across the tropics, and SINV in colder climates. The tested viruses have differences in both virion structure (T = 3 [flavivirus] versus T = 4 [alphavirus] icosahedral geometry, lack of M structural protein in alphaviruses) and envelope proteins that create variation in their capacity for cellular entry and invasion, as well as differences in genomic structure (flaviviruses do not have a poly-A 3′ end) (Figure 1B). Following viral challenge, we monitored mosquitoes for arbovirus infection, dissemination of virus beyond the midgut and throughout the body, and secretion of virus in saliva. We report for the first time that An. albimanus can become infected with and transmit multiple alphaviruses—including MAYV, a human pathogen that is already spreading within this mosquito’s geographic range in the Americas.
RESULTS
Growth kinetics of arboviruses in vitro.
Prior to in vivo infection experiments, we assessed the in vitro growth of each virus over time. Each virus was propagated in Vero (alphaviruses) or C6/36 (flavivirus) cells. Alphaviruses typically replicate more quickly than flaviviruses, which can be seen in the rapid and severe cytopathic effects they produce in invaded cells of vertebrate origin.25 To ensure virus for subsequent experiments was harvested during its replicative stage, we collected samples for alphaviruses across different time points up to 48 hours postinfection (hpi), when the presence of abundant cellular death and abnormal media pH were visibly evident. Infectious titers were assessed by plaque-forming assay (PFA; SINV) or focus-forming assay (FFA; the remaining viruses). All alphavirus titers peaked at 24 hpi with viral titers > 107 focus-forming units (ffu) (or plaque-forming units; pfu) per milliliter (Figure 1C), after which they declined only slightly until the final time point (48 hpi). Note that CHIKV’s viral titer (3 × 108 ffu/mL) was assessed only at 24 hpi due to biosafety level 3 (BSL-3) laboratory constraints. Conversely, DENV-2 infections were sampled at days 0, 3, and 7, and we found viral titers in supernatant were highest at day 7 (1 × 107 ffu/mL).
Both MAYV genotypes infect, replicate, and are transmitted efficiently by An. albimanus.
For MAYV mosquito challenges, we used two different strains of the virus (BeAn343102, genotype -D; BeAr505411, genotype -L). Results are reported as 1) infection rate (IR), the proportion of challenged mosquitoes with infected midguts, 2) dissemination rate (DR) and efficiency (DE), the proportion of infected (DR) or challenged (DE) mosquitoes with infected bodies, and 3) transmission rate (TR) and efficiency (TE), the proportion of infected (TR) or challenged (TE) mosquitoes with infected saliva. These ratios were calculated using the number of samples deemed infection positive by viral titers measured from the midgut, carcass (rest of the body), and saliva, respectively, and are reported with numerical subscripts to indicate sampling day. Viral titers are reported as means at individual time points (shown as Mx, Cx, or Sx[in], where “M” represents midgut, “C” represents carcass, “S” represents saliva samples, “x” denotes the day of collection (days postinfection; dpi), and “in” represents only infected subsets. Average titers across all time points are reported without a time point subscript.
For both MAYV genotypes, we found that An. albimanus was highly susceptible to infection, dissemination, and transmission. MAYV-D successfully established infections in the midgut and disseminated to the rest of the body in nearly all mosquitoes at all surveyed time points (Figure 2A, Table 1; IR: 97.6%; DR: 97.6%; DE: 95.3%). Although infection prevalence was consistent across time, disseminated viral titers rose until 10 dpi (Figure 2A; C7 versus C10: U = 250, P = 0.034), when both infection (M10: 1.3 × 106 ffu/mL) and dissemination (C10: 3.5 × 107 ffu/mL) titers were highest. The infection and dissemination patterns of MAYV-L were similar to those of MAYV-D, with the virus infecting and disseminating through mosquitoes at high rates (Figure 2B, Table 1; IR: 94.3%; DR: 98.8%; DE: 93.2%). As with MAYV-D, MAYV-L titers also peaked at 10 dpi both in midguts (M10: 7.6 × 106 ffu/mL; M7 versus M10: U = 125, P < 0.0001) and carcasses (C10: 2.5 × 106 ffu/mL; C7 versus C10: U = 119, P < 0.0001). Titers then decayed significantly from this peak (M10 versus M14: U = 151, P = 0.0004; C10 versus C14: U = 191.5, P = 0.007). Despite high prevalence at all time points, we observed a higher variation in infection and dissemination for MAYV-D compared with MAYV-L (i.e., M14: F = 3.23, P = 0.0026; C14: F = 4.84, P < 0.0001). Transmission trends differed slightly between the two genotypes. Unlike its infection and dissemination rates, which were high and steady, MAYV-D transmission efficiency (Figure 2A) increased over time (Table 1; TE7: 20%; TE10: 32%; TE14: 48.4%), although viral titers detected in those infectious mosquitoes remained constant (S7in: 1.3 × 103 ffu/mL; S14in: 1.25 × 103 ffu/mL). Note that MAYV-D was not sampled at 21 dpi due to high mortality that differed from controls (Figure 2C,Supplemental Figure 1; 96.9%, χ2 = 43.72, degree of freedom [df] = 1, P < 0.0001). For MAYV-L, virus was present in the saliva samples of about one quarter of mosquitoes at most time points (Figure 2B, Table 1; TE7: 26.9%; TE10: 23%; TE14: 22.2%) but neither prevalence nor viral titers (S7in: 1.76 × 102 ffu/mL; S14in: 3.11 × 102 ffu/mL) increased significantly over time. At 21 dpi, no mosquitoes were able to transmit the virus (i.e., TE21: 0%), which can be due to viral clearance from the body. Indeed, lower titers were also detected in midgut and carcass at 21 dpi (Figure 2B; MAYV-L: M14 versus M21: U = 19.5, P < 0.0001; C14 versus C21: U = 8.5, P < 0.0001). However, these low titers may reflect selection bias, because mortality was also very high in MAYV-L-challenged mosquitoes at 21 dpi (Figure 2C,Supplemental Figure 1; 93.75%, χ2 = 15.59, df = 1, P < 0.0001).
Parameters describing infections in An. albimanus
Virus and collection time | % IR (x/N) | % DR (y/x) | % DE (y/N) | % TR (z/y) | % TE (z/N) |
---|---|---|---|---|---|
MAYV-D | |||||
7 dpi | 100 (30/30) | 96.7 (29/30) | 96.7 (29/30) | 20.6 (6/29) | 20 (6/30) |
10 dpi | 96 (24/25) | 100 (24/24) | 96 (24/25) | 33.3 (8/24) | 32 (8/25) |
14 dpi | 96.8 (30/31) | 96.7 (29/30) | 93.5 (29/31) | 51.8 (15/29) | 48.4 (15/31) |
21 dpi | – | – | – | – | – |
MAYV-L | |||||
7 dpi | 92.3 (24/26) | 95.8 (23/24) | 88.5 (23/26) | 30.4 (7/23) | 26.9 (7/26) |
10 dpi | 96.2 (25/26) | 96 (24/25) | 92.3 (24/26) | 25 (6/24) | 23 (6/26) |
14 dpi | 96.3 (26/27) | 104 (27/26) | 100 (27/27) | 22.2 (6/27) | 22.2 (6/27) |
21 dpi | 88.9 (8/9) | 100 (8/8) | 88.9 (8/9) | 0 (0/8) | 0 (0/9) |
ONNV | |||||
7 dpi | 100 (26/26) | 96.2 (25/26) | 96.2 (25/26) | 8 (2/25) | 7.7 (2/26) |
10 dpi | 96.2 (25/26) | 100 (25/25) | 96.2 (25/26) | 4 (1/25) | 3.8 (1/26) |
14 dpi | 96.3 (26/27) | 100 (26/26) | 96.3 (26/27) | 3.8 (1/26) | 3.7 (1/27) |
21 dpi | 100 (16/16) | 100 (16/16) | 100 (16/16) | 31.3 (5/16) | 31.3 (5/16) |
SINV | |||||
7 dpi | 100 (24/24) | 95.8 (23/24) | 95.8 (23/24) | 4.3 (1/23) | 4.2 (1/24) |
10 dpi | 100 (27/27) | 92.6 (25/27) | 92.6 (25/27) | 4 (1/25) | 3.7 (1/27) |
14 dpi | 94.7 (18/19) | 100 (18/18) | 94.7 (18/19) | 11.1 (2/18) | 10.5 (2/19) |
21 dpi | 93.8 (15/16) | 100 (15/15) | 93.8 (15/16) | 26.7 (4/15) | 25 (4/16) |
CHIKV | |||||
7 dpi | 33 (10/30) | 20 (2/10) | 6.7 (2/30) | 0 (0/2) | 0 (0/30) |
10 dpi | 10 (3/30) | 0 (0/3) | 0 (0/30) | – | 0 (0/30) |
14 dpi | 10 (3/30) | 0 (0/3) | 0 (0/30) | – | 0 (0/30) |
21 dpi | – | – | – | – | – |
DENV-2 | |||||
7 dpi | 0 (0/17) | – | 0 (0/17) | – | 0 (0/17) |
10 dpi | 0 (0/14) | – | 0 (0/14) | – | 0 (0/14) |
14 dpi | 0 (0/13) | – | 0 (0/13) | – | 0 (0/13) |
21 dpi | 0 (0/8) | – | 0 (0/8) | – | 0 (0/8) |
DE = dissemination efficiency; dpi = days post infection; DR = dissemination rate; IR = infection rate; N = total number of challenged mosquitoes; TE = transmission efficiency; TR = transmission rate; x = number with virus present in midgut; y = number with virus present in carcass; z = number with virus present in saliva. IR, DR, DE, TR, and TE for all the assessed viruses are reported.
An. albimanus is not a competent vector of CHIKV or DENV-2.
Whereas MAYV was able to successfully infect and transmit through An. albimanus mosquitoes, this was not the case for the other assayed viruses endemic in its native range. Following challenge with CHIKV, An. albimanus was able to become infected at moderate levels; we detected viral presence in 33% of the mosquitoes’ midguts at 7 dpi (Figure 3A; M7in: 2 × 104 ffu/mL), but that decreased to 10% at 10 and 14 dpi (M10in: 3.67 × 104 ffu/mL; M14in: 1.44 × 104 ffu/mL). We found that CHIKV did not efficiently escape the midgut and disseminate, likely due to a midgut escape barrier. Only two carcass samples were CHIKV positive at 7 dpi (6.7%, 2/30), and 0% of infections were disseminated at both 10 and 14 dpi (Figure 3A, Table 1). None of these carcass-infected mosquitoes progressed to infected saliva. It therefore appears that CHIKV infection, at least with the H20235 virus strain and An. albimanus colony strain we tested, can only be established at the midgut level and rarely disseminates. CHIKV-associated mosquito mortality was not significantly different compared with the blood-fed-only controls (Figure 3C,Supplemental Figure 1; 45.9%, χ2 = 0.82, df = 1, P = 0.365).
Besides American alphaviruses, we also assessed An. albimanus’s vector competence for DENV-2 due to its socioeconomic importance worldwide and because it is endemic in An. albimanus’s habitat. Consistent with previous findings in anophelines, here we report that An. albimanus is refractory to DENV-2 infection. Despite challenging the mosquitoes with relatively high titers (5 × 106 ffu/mL) and detecting positive sampling at 0 dpi (which depicts mosquito intake of an infectious bloodmeal; Figure 3B), none of the collected individuals (Figure 3B, Table 1; 0/61) from 7 to 21 dpi was found to carry infectious virus in midgut, carcass, or saliva. These results led us to question whether anophelines possess a flavivirus-specific midgut barrier, which would restrict viral replication in the midgut upon the ingestion of an infective bloodmeal, or instead they possess a body-wide infection or replication barrier against flaviviruses that prevents replication in all tissues. To test this, we injected infectious DENV-2 into the hemolymph of An. albimanus, thereby bypassing the hypothetical midgut barrier. To test anophelines and flaviviruses more broadly than this one species and virus, we additionally tested DENV in An. gambiae, An. stephensi, and Anopheles quadrimaculatus while also testing Zika virus (ZIKV) in An. albimanus. After 3 days of infection, none of the mosquitoes presented infectious particles in their bodies as assessed by FFA (Supplemental Figure 2), indicating that the injected virus was not able to replicate within the mosquito and that the tested anophelines are completely refractory to human flavivirus infections.
Alphaviruses present outside the Americas can infect, disseminate, and be transmitted by An. albimanus.
We also asked whether An. albimanus may be a suitable vector of ONNV and SINV. These alphaviruses are not known to have caused outbreaks in the Americas yet, but they could emerge in currently unaffected areas due to globalization, travel, and climate change—just as ZIKV spread to new continents.26,27 We found infectious ONNV virions (p5′dsONNic) were able to both infect and disseminate from the midgut in nearly all challenged An. albimanus when they were fed at 1 × 107 ffu/mL (Figure 4A, Table 1; IR: 97.9%; DR: 98.9%; DE: 96.8%), showing that the species is susceptible to the virus. Unlike the pattern observed for MAYV, the highest ONNV infection intensity in the midgut was detected at 7 dpi (Figure 4A; M7 = 4 × 105 ffu/mL; M7 versus M10: U = 90.5, P < 0.0001), which then dropped slightly and remained stable throughout the remaining time points (M10 versus M21: U = 167, P = 0.29). Dissemination viral titers were similar to MAYV, peaking at 10 dpi (C10 = 2.5 × 106 ffu/mL; C7 versus C10: U = 148.5, P = 0.0004) followed by a stable plateau (C10 versus C21: U = 155.5, P = 0.18). Despite high viral prevalence in the midgut and body, transmission rates were very low early in ONNV infections (Figure 4A; TE/TR7–14: 4–7%) until an abrupt increase at 21 dpi (TE/TR21: 31.3%; S14 versus S21: U = 157, P = 0.014). Our data thus show that the virus can efficiently invade the salivary glands only late in the course of infections. Mortality was not significantly different in ONNV-positive mosquitoes compared with their blood-fed-only counterparts (Figure 4C,Supplemental Figure 1; 70.4%, χ2 = 2.06, df = 1, P = 0.151).
An. albimanus was also susceptible to being infected and transmitting SINV, and infection dynamics were broadly similar to those of ONNV. Almost all mosquitoes presented infected midguts (Figure 4B, Table 1; IR: 97.7%), with dissemination to the rest of their bodies also occurring at high frequencies (Figure 4B, Table 1; DR: 96.4%; DE: 94.2%). SINV titers were stable over time—only increasing in carcasses (i.e., dissemination) from 7 to 14 dpi (U = 115, P = 0.005). This may indicate a slower replication rate within the mosquito tissues that delays and reduces its vectorial capacity. Early in the course of infections, few saliva samples presented infectious SINV (S7–10: 4%; S14: 11%) and those that did had low titers. Similar to ONNV, viral loads in saliva rose significantly at the latest time point (TE21: 25%; S21in: 3.7 × 102 pfu/mL), again representing a slower transition from midguts to the salivary glands.
DISCUSSION
Here, we demonstrate that An. albimanus, the most common Anopheles mosquito in Mesoamerica and the Caribbean,8,16,28 is a competent laboratory vector for a broad range of alphaviruses but refractory to flavivirus infection. Although arboviral spread through anophelines has received little research attention, we found this species was readily infected and transmitted three of five tested viruses—suggesting that it (and possibly other anophelines) may be susceptible to infection with a substantial number of viruses.3 Although the tested alphaviruses (except CHIKV) were able to establish disseminated infections and be transmitted by this species, they presented different patterns of infection and transmission.
Notably, an alphavirus that is endemic to the Americas (MAYV)—and that therefore has the potential to infect An. albimanus in the wild—had the highest transmission efficiency of the tested viruses. Both MAYV genotypes (-D and -L; Figure 2A and B, respectively) are able to escape the mosquito’s tissue barriers,29 replicating freely within the body of the mosquito and reaching the salivary gland lumen. We found MAYV infections progressed the most rapidly, as indicated by higher titers in the midguts and body, as well as earlier and greater prevalence in saliva. Saliva samples tested positive for MAYV at a substantial rate (20–48%) at all surveyed time points except 21 dpi, when most mosquitoes that potentially carried high viremia may have suffered from MAYV-associated mortality (Figure 2C,Supplemental Figure 1). This increased mortality likely affects their vectorial capacity in the wild, because mosquitoes may die prematurely, or because they may be infectious but not vigorous enough to seek out and bite hosts. However, it remains unclear whether (and how often) An. albimanus becomes infected with MAYV in the wild.
Our data also suggest there may be differences in how the two MAYV subgroups infect the tested colony of An. albimanus. For example, we detected higher variance in MAYV-D infection titers compared with those of MAYV-L infections (Figure 2A versus Figure 2B). Although both show a short extrinsic incubation period30 that represents heightened vectorial capacity, MAYV-D’s transmission increased over time whereas MAYV-L remained consistent at around 25%, although this pattern was not statistically significant. Similar variation between MAYV genotypes in infection was previously observed in An. gambiae.20 These differences suggest there may be interesting phenotypical variation in infection that could be explored in future studies.
To date, few MAYV outbreaks have been detected in humans, with most infections occurring in the Amazon Basin.31 However, recent studies indicate that MAYV’s true incidence may be underestimated32 because its symptoms resemble those of other pathogens in the area (DENV, CHIKV). Epidemiological concern is rising due to MAYV’s ability to colonize urban areas and its potential to be transmitted through lesser-known vectors of arboviruses in Central and South America.31,33 With the geographical range expansion of An. albimanus into lower latitudes of South America as well as northward into the United States, coupled with the high viral susceptibility and transmission observed in a laboratory colony, there is vectorial potential for this species to drive such outbreaks.
Not all alphaviruses endemic to the Americas were successful in infecting An. albimanus. Most notably, CHIKV had poor prevalence in challenged mosquitoes. Differences in transmission capacity can be explained by dynamic and complex virus–vector interactions that affect the vectorial capacity of the mosquito species,34 and which may trace back to variation in host-to-pathogen genotype interactions (GxG). GxG interactions govern outcomes in many systems, including arboviral infections35 and symbiosis36 in insects. In some cases,26,27 the combination of a specific insect genotype and a specific pathogen genotype strongly shapes infection status and progression. Although Anopheles are not typically vectors for CHIKV,37 some species including An. stephensi can harbor and transmit it.24 In these rare cases, GxG interactions may be key to the ability of the mosquito to transmit the virus. Consistent with what has been observed for most anophelines, An. albimanus were able to harbor the infection following a viral bloodmeal but were refractory to dissemination and transmission of CHIKV (Figure 3A). One hypothesis is that An. albimanus has coevolved with this strain of CHIKV (collected on the Caribbean island of St. Martin) and the GxG interaction between mosquito and virus is not ideal.35,38 Of course, our findings do not exclude the possibility that CHIKV may be transmitted through An. albimanus mosquitoes in some circumstances, because these results were obtained using a single laboratory colony. It remains possible that wild strains from other locations may present a more suitable combination for the virus to transmit—by differing either in their nuclear genotype35 or microbiome composition.39
We found that An. albimanus can also carry and transmit the non-American alphaviruses ONNV and SINV (Figure 4A and B). These two viruses affect different geographical and climate regions and are primarily transmitted by other mosquito species.40,41 ONNV is closely related to CHIKV and MAYV (all belong to the Semliki Forest antigenic complex) and causes outbreaks in humans in Central Africa. In contrast, SINV is a distant relative. It belongs to the Western equine encephalitis complex and is the causal agent of Pogosta disease, mostly in northern Europe42 and southern Africa,43 with sporadic cases detected in Australia.44 Despite their biological differences, both ONNV and SINV are able to infect An. albimanus, albeit with slower transmission kinetics compared with MAYV. More specifically, although infected early at the midgut level, very few ONNV- or SINV-challenged mosquitoes were saliva positive at the earliest time point tested (7 dpi). Rather, their highest transmission was detected at 21 dpi for both viruses. This longer time of incubation (and the mortality observed over time) suggests that, even though viral transmission rates are somewhat elevated at that time, mosquito survival is decreasing rapidly due to aging and infection status and should not necessarily pose a higher risk of transmission epidemiologically. Although neither ONNV nor SINV has been detected in the Americas yet, globalization and human travel45 pose a risk for these viruses to spread there, especially ONNV, whose primary vectors in Africa include the anopheline species An. gambiae and An. funestus46,47 whereas SINV is mainly transmitted by culicine species.48
We found multiple anopheline species were refractory to infection with the tested flaviviruses DENV-2 and ZIKV, even when virus was injected directly into the hemolymph. The inability to harbor these flaviviruses could trace back to a variety of different mechanisms, for example, lack of a particular replication machinery component49 or absence of a specific cellular receptor or factor required for infection success.50–52 However, there are reports of insect-specific flavivirus infections in anophelines,53,54 raising questions about the underlying mechanism and host restriction features of different flaviviruses, and whether some specific human flavivirus infections do occur in anophelines given permissive GxG interactions.
In short, this study tested the capacity of An. albimanus to be an important vector of alphaviruses in the Americas. Our data show that, although An. albimanus is unlikely to drive CHIKV and flavivirus infections, other alphaviruses (especially MAYV) can infect and be transmitted through this mosquito species very efficiently. Our results highlight both the importance for Anopheles mosquitoes to be recognized as potential vectors of arboviruses as well as the knowledge gaps that still need to be filled on the genus’s vectorial capacities worldwide.
MATERIALS AND METHODS
Mosquitoes.
The An. albimanus mosquito colony (STECLA strain, MRA-126) was kept and reared at the Millennium Science Complex insectary (Center for Infectious Disease Dynamics, The Pennsylvania State University) at a continuous 27 ± 1°C, 85% relative humidity, 12:12-hour light:dark cycle. Larvae were fed koi pellets (Tetra Pond Koi Vibrance; Tetra, Melle, Germany) from hatching to pupation. Adult mosquitoes were reared in 30 × 30 × 30-cm metal cages and provided 10% sucrose solution ad libitum, as well as fed weekly with anonymous human blood (Biological Specialty, Colmar, PA) for reproduction and colony maintenance using a membrane feeder.
Cells.
African green monkey kidney cells (Vero; CCL-81) and Aedes albopictus larval cells (C6/36; CCL-126) (ATCC, Manassas, VA) were cultured in complete media consisting of Dulbecco’s modified Eagle’s medium or RPMI-1640, respectively, complemented with 10% fetal bovine serum (FBS) and 1% penicillin/streptomycin—all reagents purchased from Gibco, Thermo Fisher Scientific (Waltham, MA). For each passage, cells were detached by trypsinization (0.25% trypsin; Corning Inc., Corning, NY) and diluted in fresh complete media or plated for experiments.
Viruses.
A total of six different viruses was used for experimental infections (Figure 1A). Two strains of MAYV were used: BeAn343102 (BEI Resources, Manassas, VA) is a genotype D strain (MAYV-D) isolated in May 1978 from a monkey in Para, Brazil, and BeAr505411 (BEI Resources) is a genotype L strain (MAYV-L) also isolated in Para, Brazil in March 1991 from Haemagogus janthinomys mosquitoes. The full-length ONNV and SINV infectious clones (p5′dsONNic and p5′dsMRE16ic) derive from the Uganda SG-650 strain of ONNV55 and wild-type MRE16 strain of SINV isolated from Malaysia,56 respectively. Both infectious clones were obtained on filter paper and transfected, and virions were collected prior to passaging in Vero cells for the experiments. The DENV-2 ET300 strain was isolated from a human patient in Timor-Leste in 2000 (GenBank accession number EF440433.1). Lastly, the CHIKV H20235 strain (NR-49901; BEI Resources) was isolated from a human in St. Martin in 2013. All work with this strain (from cell culture to mosquito infections) was performed in the Eva J. Pell Laboratory for Advanced Biological Research BSL-3 facility at The Pennsylvania State University.
All alphaviruses were passed in African green monkey kidney (Vero) cells at 37°C in a humidified 5% CO2 incubator, whereas DENV-2 was passed in Ae. albopictus RNA interference-deficient C6/36 cells at 28°C. Viruses were allowed to infect cells at a multiplicity of infection of 0.1 for 1 hour and then removed and replaced with media containing 2% FBS. Virus-infected supernatant was aliquoted at different time points (typically 24 hpi for alphaviruses and 7 dpi for DENV-2) and stored at −80°C until further titration or use for mosquito infections. Viral stock titers were obtained using FFAs (ffu/mL) or PFAs (pfu/mL), as described below.
Vector competence assays.
To determine the vectorial competence of An. albimanus, adult females were orally challenged with an infected bloodmeal containing a high-titer dose of one of the five togaviruses or DENV-2. Specifically, 6- to 8-day-old non-blood-fed females were allowed to feed on human blood for 1 hour through a synthetic membrane at the bottom of a glass feeder jacketed with 37°C water and containing either 1 × 107 ffu/mL (alphaviruses) or 5 × 106 ffu/mL (DENV-2) of the stocks obtained above (Figure 1C, gray boxes). Fully engorged mosquitoes were sorted from non-fed ones and split evenly in cups for each collection time point. A small subset of mosquitoes was collected at 0 dpi to confirm that the viral intake was infectious and similar across samples.
Infection rate, dissemination and transmission rates (DR and TR), as well as dissemination and transmission efficiency (DE and TE) were assessed at 7, 10, 14, and 21 dpi. IR was measured as the rate of mosquitoes with infected midguts among the total number of mosquitoes. DR and DE were measured as the rate of mosquitoes with infected carcasses among the mosquitoes with infected midguts or over the total assessed samples, respectively. TR and TE were measured as the rate of mosquitoes with infectious saliva among the positive bodies or the total number of assessed mosquitoes, respectively.
At all time points, mosquitoes were anesthetized using triethylamine (Sigma, St. Louis, MO) before individual forced salivation. Saliva was collected by placing the female’s proboscis into a pipette tip containing 20 µL of a 50% sucrose, 50% FBS solution, as previously described,57 for 30 minutes. Solution was then released into a tube filled with 100 µL of mosquito diluent (20% heat-inactivated FBS, 50 µg/mL penicillin/streptomycin, 50 µg/mL gentamicin, and 2.5 µg/mL fungizone in Dulbecco’s phosphate-buffered saline [PBS]) and placed on ice. Each female’s midgut was dissected and placed in a 2-mL tube containing 300 µL of mosquito diluent. The rest of the body (carcass) was also collected in an identical tube. Tissue samples were homogenized at 30 Hz for 2 min using a TissueLyser II (QIAGEN, Hilden, Germany). All samples were stored at −80°C until viral testing.
Intrathoracic injections.
Anopheles gambiae, An. stephensi, An. quadrimaculatus, and An. albimanus females were briefly anesthetized on a chill block (BioQuip Products, Compton, CA) cooled to 4°C and DENV-2 and ZIKV stocks were injected intrathoracically under a microscope using a pulled glass capillary with a manual microinjector (Nanoject II, Drummond Scientific, Broomall, PA) to ensure uniformity of dosage. Sixty-nine microliters of diluted virus stock (∼70 DENV-2 pfu) were delivered intrathoracically into each female. After injection, mosquitoes were maintained under standard housing conditions of 27°C with 80% relative humidity and 12:12-hour light:dark cycle and fed 10% sucrose solution ad libitum.
Focus-forming assay.
The presence of infectious particles of all viruses except SINV in saliva, midguts, and carcasses was tested by FFAs in Vero (alphaviruses) or C6/36 (DENV-2) cells. Cells were counted using a hemacytometer (Hausser Scientific, Horsham, PA) and plated in complete media the day before infection to achieve 80–90% confluency (Vero: 3 × 104 cells per well; C6/36: 2 × 105 cells per well) in 96-well plates. The following day, media were removed from all wells, and cells were incubated for 1 hour with 30 µL of 10-fold dilutions (10−1 to 10−4) of each homogenized tissue sample in FBS-free media. Saliva samples were not diluted, due to their lower titers. Viral media were removed from the wells after 1 hour, replaced with 100 µL of overlay (final 0.8% methylcellulose [or CMC] in complete media), and incubated at 37°C for 24 hours or 28°C for 3 days, depending on the cell culture used. Cells were then fixed using 4% paraformaldehyde in PBS (Sigma-Aldrich) for 20 min and permeabilized with 0.2% Triton-X in PBS for another 20 minutes. Samples were washed two or three times with cold 1× PBS after both fixation and permeabilization steps. Viral antigens in infected cells were labeled overnight using mouse monoclonal anti-CHIKV E2 envelope glycoprotein clone CHK-48 (for all alphaviruses except SINV; α-CHK-48, BEI Resources) or mouse monoclonal anti-flavivirus clone D1-4G2-4-15 (for DENV-2; BEI Resources) diluted 1:500 in PBS. The next day, cells were again washed thoroughly with cold PBS to remove unbound primary antibody. Bound primary antibody was then labeled for 1 hour at room temperature using an Alexa 488 goat anti-mouse IgG secondary antibody (Invitrogen, OR Waltham, MA) at a 1:750 dilution in PBS, which was then rinsed off with reverse osmosis water before evaluation by fluorescence microscopy. Green fluorescence was observed using a fluorescein isothiocyanate filter on an Olympus BX41 microscope with a UPlanFI 4× objective. Foci were counted by eye in the appropriate dilution (10–100 foci) and viral titers were backcalculated to ffu/mL.
Plaque-forming assay.
The α-CHK-48 antibody used in FFAs does not cross-react with SINV, which is evolutionarily the most distantly related to CHIKV.58 Thus, we elected to assess mosquito SINV infections by traditional PFAs, which are comparable to FFAs because both detect the presence of infectious viral particles in a sample using a cell-based method.
Mosquito samples were tested for SINV infectious particles by plaque assay on Vero cells with minimal modifications.59 The day before infection, cells were counted as described above and plated (5 × 106 cells per well) in six-well plates. For saliva infections, media were removed and replaced for 100 µL of undiluted sample. For midguts or carcasses, 10-fold dilutions using FBS-free media were performed and 180 µL of each dilution (in most instances, 10−2 to 10−4) was used for cellular infection. Inoculated plates were placed in a 37°C incubator with 5% CO2 for 1 hour for viral entry to occur. Then, virus-containing media were removed and replaced with 1.5 mL of an agar overlay (equal parts complete media and 1.2% agarose) and placed back into the incubator. After 2 days, 1.5 mL of a second agar overlay (identical to the first agar overlay but containing a 1% final concentration of neutral red [Amresco, Solon, OH] to allow for cellular staining) was added to each well and plates were incubated. The following day, agar discs were removed, and samples were treated with 4% formaldehyde to inactivate any remaining virus. Each stained well was rinsed thoroughly with water and set aside to dry well. Wells that produced 10–100 plaques were used to ensure accurate counts, and viral titers for mosquito saliva, midguts, and carcasses were calculated in pfu/mL. When samples produced too many plaques to count, additional plaque assays were performed with extra 10-fold dilutions.
Statistical analysis and figure generation.
Differences in viral titer between midgut, carcass, and saliva samples were assessed by two-tailed nonparametric Mann-Whitney U tests due to the nonnormality of the samples. In a few cases where samples were normally distributed, we used a parametric Welch’s t test instead. Survival data for each virus were compared pairwise to a blood-fed-only control sample using a log-rank Mantel-Cox test. All P values that were below 0.05 (P < 0.05) were considered significant. All data were initially plotted and analyzed using Prism software version 9.2.0 (283) (GraphPad, San Diego, CA). Final figures were assembled using Adobe Illustrator 2021 (25.4.1; Adobe, San Jose, CA).
Supplemental Materials
ACKNOWLEDGMENTS
ONNV and SINV infectious clones were kindly provided by Brian Foy (Colorado State University, Ft. Collins, CO). The DENV-2 ET300 strain was provided by the McGraw lab (Penn State University, University Park, PA). Mosquitoes and antibodies were provided by BEI Resources. We thank Kaylee Montanari and Amelia Romo from the Rasgon lab, as well as personnel from the Eva J. Pell BSL-3 laboratory for Advanced Biological Research at Penn State, for their continual technical support.
REFERENCES
- 2.↑
World Health Organization , 2021. World Malaria Report. Geneva, Switzerland: World Health Organization.
- 3.↑
Brault AC, Foy BD, Myles KM, Kelly CLH, Higgs S, Weaver SC, Olson KE, Miller BR, Powers AM, 2004. Infection patterns of O’nyong nyong virus in the malaria-transmitting mosquito, Anopheles gambiae. Insect Mol Biol 13: 625–635.
- 4.↑
Nanfack Minkeu F, Vernick KD, 2018. A systematic review of the natural virome of Anopheles mosquitoes. Viruses 10: 222.
- 5.↑
Service MW, 1997. Mosquito (Diptera: Culicidae) dispersal—the long and short of it. J Med Entomol 34: 579–588.
- 6.↑
Ravinet M, Faria R, Butlin RK, Galindo J, Bierne N, Rafajlović M, Noor MAF, Mehlig B, Westram AM, 2017. Interpreting the genomic landscape of speciation: a road map for finding barriers to gene flow. J Evol Biol 30: 1450–1477.
- 8.↑
Fuller DO, Ahumada ML, Quiñones ML, Herrera S, Beier JC, 2012. Near-present and future distribution of Anopheles albimanus in Mesoamerica and the Caribbean Basin modeled with climate and topographic data. Int J Health Geogr 11: 13.
- 9.↑
Rodríguez AD, Rodríguez MH, Meza RA, Hernández JE, Rejmankova E, Savage HM, Roberts DR, Pope KO, Legters L, 1993. Dynamics of population densities and vegetation associations of Anopheles albimanus larvae in a coastal area of southern Chiapas, Mexico. J Am Mosq Control Assoc 9: 46–58.
- 10.↑
Savage HM, Rejmankova E, Arredondo-Jimenez JI, Roberts DR, Rodriguez MH, 1990. Limnological and botanical characterization of larval habitats for two primary malarial vectors, Anopheles albimanus and Anopheles pseudopunctipennis, in coastal areas of Chiapas State, Mexico. J Am Mosq Control Assoc 6: 612–620.
- 11.↑
Hernandez JE, Epstein LD, Rodriguez MH, Rodriguez AD, Rejmankova E, Roberts DR, 1997. Use of generalized regression tree models to characterize vegetation favoring Anopheles albimanus breeding. J Am Mosq Control Assoc 13: 28–34.
- 12.↑
Vázquez-Martínez MG, Rodríguez MH, Arredondo-Jiménez JI, Méndez-Sánchez JD, Bond-Compeán JG, Gold-Morgan M, 2002. Cyanobacteria associated with Anopheles albimanus (Diptera: Culicidae) larval habitats in southern Mexico. J Med Entomol 39: 825–832.
- 13.↑
Hobbs JH, Sexton JD, St Jean Y, Jacques JR, 1986. The biting and resting behavior of Anopheles albimanus in northern Haiti. J Am Mosq Control Assoc 2: 150–153.
- 14.↑
Mekuria Y, Tidwell MA, Williams DC, Mandeville JD, 1990. Bionomic studies of the Anopheles mosquitoes of Dajabon, Dominican Republic. J Am Mosq Control Assoc 6: 651–657.
- 15.↑
Roberts DR, Manguin S, Rejmankova E, Andre R, Harbach RE, Vanzie E, Hakre S, Polanco J, 2002. Spatial distribution of adult Anopheles darlingi and Anopheles albimanus in relation to riparian habitats in Belize, Central America. J Vector Ecol 27: 21–30.
- 16.↑
Sinka ME, Rubio-Palis Y, Manguin S, Patil AP, Temperley WH, Gething PW, Van Boeckel T, Kabaria CW, Harbach RE, Hay SI, 2010. The dominant Anopheles vectors of human malaria in the Americas: occurrence data, distribution maps and bionomic précis. Parasit Vectors 3: 72.
- 17.↑
Loyola EG, González-Cerón L, Rodríguez MH, Arredondo-Jiménez JI, Bennett S, Bown DN, 1993. Anopheles albimanus (Diptera: Culicidae) host selection patterns in three ecological areas of the coastal plains of Chiapas, southern Mexico. J Med Entomol 30: 518–523.
- 18.↑
Solarte Y, Hurtado C, Gonzalez R, Alexander B, 1996. Man-biting activity of Anopheles (Nyssorhynchus) albimanus and An. (Kerteszia) neivai (Diptera: Culicidae) in the Pacific lowlands of Colombia. Mem Inst Oswaldo Cruz 91: 141–146.
- 19.↑
Bown DN, Rodriguez MH, Arredondo-Jimenez JI, Loyola EG, Rodriguez MC, 1993. Intradomiciliary behavior of Anopheles albimanus on the coastal plain of southern Mexico: implications for malaria control. J Am Mosq Control Assoc 9: 321–324.
- 20.↑
Oidtman RJ, España G, Perkins TA, 2021. Co-circulation and misdiagnosis led to underestimation of the 2015–2017 Zika epidemic in the Americas. PLoS Negl Trop Dis 15: e0009208.
- 21.↑
Franz AW, Kantor AM, Passarelli AL, Clem RJ, 2015. Tissue barriers to arbovirus infection in mosquitoes. Viruses 7: 3741–3767.
- 22.↑
Rezza G, Chen R, Weaver SC, 2017. O’nyong-nyong fever: a neglected mosquito-borne viral disease. Pathog Glob Health 111: 271–275.
- 23.↑
Brustolin M, Pujhari S, Henderson CA, Rasgon JL, 2018. Anopheles mosquitoes may drive invasion and transmission of Mayaro virus across geographically diverse regions. PLoS Negl Trop Dis 12: e0006895.
- 24.↑
Yadav P, Gokhale MD, Barde PV, Singh DK, Mishra AC, Mourya DT, 2003. Experimental transmission of Chikungunya virus by Anopheles stephensi mosquitoes. Acta Virol 47: 45–47.
- 25.↑
Akhrymuk I, Lukash T, Frolov I, Frolova EI, 2019. Novel mutations in nsP2 abolish Chikungunya virus-induced transcriptional shutoff and make the virus less cytopathic without affecting its replication rates. J Virol 93: e02062-18.
- 26.↑
Duffy MR et al., 2009. Zika virus outbreak on Yap Island, Federated States of Micronesia. N Engl J Med 360: 2536–2543.
- 27.↑
Zanluca C, de Melo VCA, Mosimann ALP, dos Santos GIV, dos Santos CND, Luz K, 2015. First report of autochthonous transmission of Zika virus in Brazil. Mem Inst Oswaldo Cruz 110: 569–572.
- 28.↑
Loaiza JR, Bermingham E, Scott ME, Rovira JR, Conn JE, 2008. Species composition and distribution of adult Anopheles (Diptera: Culicidae) in Panama. J Med Entomol 45: 841–851.
- 29.↑
Hardy JL, Houk EJ, Kramer LD, Reeves WC, 1983. Intrinsic factors affecting vector competence of mosquitoes for arboviruses. Annu Rev Entomol 28: 229–262.
- 30.↑
Reisen WK, 1989. Estimation of vectorial capacity, Palm Springs CA, November 16, 1988. Bull Soc Vector Ecol 14: 39–70.
- 31.↑
Caicedo E-Y, Charniga K, Rueda A, Dorigatti I, Mendez Y, Hamlet A, Carrera J-P, Cucunubá ZM, 2021. The epidemiology of Mayaro virus in the Americas: a systematic review and key parameter estimates for outbreak modelling. PLoS Negl Trop Dis 15: e0009418.
- 32.↑
Brunini S, França DDS, Silva JB, Silva LN, Silva FPA, Spadoni M, Rezza G, 2017. High frequency of Mayaro virus IgM among febrile patients, central Brazil. Emerg Infect Dis 23: 1025–1026.
- 33.↑
Esposito DLA, Fonseca B, 2017. Will Mayaro virus be responsible for the next outbreak of an arthropod-borne virus in Brazil? Braz J Infect Dis 21: 540–544.
- 34.↑
Lambrechts L, Fellous S, Koella JC, 2006. Coevolutionary interactions between host and parasite genotypes. Trends Parasitol 22: 12–16.
- 35.↑
Lambrechts L, Chevillon C, Albright RG, Thaisomboonsuk B, Richardson JH, Jarman RG, Scott TW, 2009. Genetic specificity and potential for local adaptation between dengue viruses and mosquito vectors. BMC Evol Biol 9: 160.
- 36.↑
Smee MR, Raines SA, Ferrari J, 2021. Genetic identity and genotype × genotype interactions between symbionts outweigh species level effects in an insect microbiome. ISME J 15: 2537–2546.
- 37.↑
Vanlandingham DL, Tsetsarkin K, Klingler KA, Hong C, McElroy KL, Lehane MJ, Higgs S, 2006. Determinants of vector specificity of O’nyong nyong and Chikungunya viruses in Anopheles and Aedes mosquitoes. Am J Trop Med Hyg 74: 663–669.
- 38.↑
Fansiri T, Fontaine A, Diancourt L, Caro V, Thaisomboonsuk B, Richardson JH, Jarman RG, Ponlawat A, Lambrechts L, 2013. Genetic mapping of specific interactions between Aedes aegypti mosquitoes and dengue viruses. PLoS Genet 9: e1003621.
- 39.↑
Cansado-Utrilla C, Zhao SY, McCall PJ, Coon KL, Hughes GL, 2021. The microbiome and mosquito vectorial capacity: rich potential for discovery and translation. Microbiome 9: 111.
- 40.↑
Powers AM, Brault AC, Tesh RB, Weaver SC, 2000. Re-emergence of Chikungunya and O’nyong-nyong viruses: evidence for distinct geographical lineages and distant evolutionary relationships. J Gen Virol 81: 471–479.
- 41.↑
Bergman A, Dahl E, Lundkvist Å, Hesson JC, 2020. Sindbis virus infection in non-blood-fed hibernating Culex pipiens mosquitoes in Sweden. Viruses 12: 1441.
- 42.↑
Skogh M, Espmark Å, 1982. Ockelbo disease: epidemic arthritis-exanthema syndrome in Sweden caused by Sindbis-virus like agent. Lancet 1: 795–796.
- 43.↑
McIntosh BM, Jupp PG, Dos Santos I, Meenehan GM, 1976. Epidemics of West Nile and Sindbis viruses in South Africa with Culex (Culex) univittatus Theobald as vector. S Afr J Sci 72: 295–300.
- 44.↑
Doherty KL, Bodey AS, Carew JS, 1969. Sinbis virus infection in Australia. Med J Aust 2: 1016–1017.
- 45.↑
Tappe D, Kapaun A, Emmerich P, de Mendonca Campos R, Cadar D, Günther S, Schmidt-Chanasit J, 2014. O’nyong-nyong virus infection imported to Europe from Kenya by a traveler. Emerg Infect Dis 20: 1766–1767.
- 46.↑
Williams MC, Woodall JP, Corbet PS, Gillett JD, 1965. O’nyong-nyong fever: an epidemic virus disease in East Africa. 8. Virus isolations from Anopheles mosquitoes. Trans R Soc Trop Med Hyg 59: 300–306.
- 47.↑
Johnson BK, Gichogo A, Gitau G, Patel N, Ademba G, Kirui R, Highton RB, Smith DH, 1981. Recovery of O’nyong-nyong virus from Anopheles funestus in western Kenya. Trans R Soc Trop Med Hyg 75: 239–241.
- 48.↑
Ling J, Smura T, Lundström JO, Pettersson JHO, Sironen T, Vapalahti O, Lundkvist Å, Hesson JC, 2019. Introduction and dispersal of Sindbis virus from Central Africa to Europe. J Virol 93: e00620-19.
- 49.↑
Slonchak A et al., 2020. Zika virus noncoding RNA suppresses apoptosis and is required for virus transmission by mosquitoes. Nat Commun 11: 2205.
- 51.↑
Clark LE et al., 2022. VLDLR and ApoER2 are receptors for multiple alphaviruses. Nature 602: 475–480.
- 52.↑
Haqshenas G, Terradas G, Paradkar PN, Duchemin J-B, McGraw EA, Doerig C, 2019. A role for the insulin receptor in the suppression of dengue virus and Zika virus in Wolbachia-infected mosquito cells. Cell Rep 26: 529–535.e3.
- 53.↑
Colmant AMG et al., 2017. A new clade of insect-specific flaviviruses from Australian Anopheles mosquitoes displays species-specific host restriction. mSphere 2: e00262-17.
- 54.↑
Guarido MM, Govender K, Riddin MA, Schrama M, Gorsich EE, Brooke BD, Almeida AP, Venter M, 2021. Detection of insect-specific flaviviruses in mosquitoes (Diptera: Culicidae) in northeastern regions of South Africa. Viruses 13: 2148.
- 55.↑
Haddow AJ, Davies CW, Walker AJ, 1960. O’nyong-nyong fever: an epidemic virus disease in East Africa 1. Introduction. Trans R Soc Trop Med Hyg 54: 517–522.
- 56.↑
Sammels LM, Lindsay MD, Poidinger M, Coelen RJ, Mackenzie JS, 1999. Geographic distribution and evolution of Sindbis virus in Australia. J Gen Virol 80: 739–748.
- 57.↑
Dodson BL, Rasgon JL, 2017. Vector competence of Anopheles and Culex mosquitoes for Zika virus. PeerJ 5: e3096.
- 58.↑
Forrester NL, Palacios G, Tesh RB, Savji N, Guzman H, Sherman M, Weaver SC, Lipkin WI, 2012. Genome-scale phylogeny of the alphavirus genus suggests a marine origin. J Virol 86: 2729–2738.
- 59.↑
Blow JA, Dohm DJ, Negley DL, Mores CN, 2004. Virus inactivation by nucleic acid extraction reagents. J Virol Methods 119: 195–198.