Intra-host growth kinetics of dengue virus in the mosquito Aedes aegypti
Authors:
Mario Novelo aff001; Matthew D. Hall aff001; Damie Pak aff003; Paul R. Young aff004; Edward C. Holmes aff005; Elizabeth A. McGraw aff001
Authors place of work:
School of Biological Sciences, Monash University, Melbourne, Victoria, Australia
aff001; Center for Infectious Disease Dynamics, Department of Entomology, The Huck Institutes of the Life Sciences, The Pennsylvania State University, University Park, Pennsylvania, United States of America
aff002; Center for Infectious Disease Dynamics, Department of Biology, The Huck Institutes of the Life Sciences, The Pennsylvania State University, University Park, Pennsylvania, United States of America
aff003; Australian Infectious Diseases Research Centre, School of Chemistry and Molecular Biosciences, The University of Queensland, St Lucia, Queensland, Australia
aff004; Marie Bashir Institute for Infectious Diseases and Biosecurity, Charles Perkins Centre, School of Life and Environmental Sciences and Sydney Medical School, The University of Sydney, New South Wales, Australia
aff005
Published in the journal:
Intra-host growth kinetics of dengue virus in the mosquito Aedes aegypti. PLoS Pathog 15(12): e32767. doi:10.1371/journal.ppat.1008218
Category:
Research Article
doi:
https://doi.org/10.1371/journal.ppat.1008218
Summary
Dengue virus (DENV) transmission by mosquitoes is a time-dependent process that begins with the consumption of an infectious blood-meal. DENV infection then proceeds stepwise through the mosquito from the midgut to the carcass, and ultimately to the salivary glands, where it is secreted into saliva and then transmitted anew on a subsequent bite. We examined viral kinetics in tissues of the Aedes aegypti mosquito over a finely graded time course, and as per previous studies, found that initial viral dose and serotype strain diversity control infectivity. We also found that a threshold level of virus is required to establish body-wide infections and that replication kinetics in the early and intermediate tissues do not predict those of the salivary glands. Our findings have implications for mosquito GMO design, modeling the contribution of transmission to vector competence and the role of mosquito kinetics in the overall DENV epidemiological landscape.
Keywords:
Blood – Viral load – Viral replication – Dengue virus – Dose prediction methods – salivary glands – Mosquitoes – Aedes aegypti
Introduction
Dengue is the most prevalent arboviral disease globally [1], with estimates of more than 390 million infections per year [2] and over half of the world’s population at risk [3]. Although mortality is generally low, the morbidity caused by dengue disease is associated with a substantial socioeconomic burden [4]. Dengue virus (DENV), transmitted to humans through the bite of the Aedes aegypti mosquito [5], is spreading in large part due to the expanding geographic range of the vector [6]. Previously confined to Africa, Aedes aegypti is now present in tropical and temperate regions worldwide, and its spread is assisted by climate change, globalization and ineffective vector control programs [7].
DENV infection in the mosquito is a complex and dynamic process [8]. The virus must circumvent multiple tissue barriers, including the midgut and salivary glands, and infect a range of intermediate tissues in a stepwise fashion [9]. Following ingestion, DENV enters the midgut epithelial cells, where it replicates. Subsequently, the virus disseminates and infects secondary tissues, including hemocytes, fat body and reproductive tissue, ultimately reaching the salivary glands [10]. The midgut is thought to represent the primary barrier to the process of infection [9], capable of preventing many mosquitoes from reaching the stage of disseminated infections [11]. The rate of this progression dictates the extrinsic incubation period (EIP), or the delay before a mosquito can infect another human on a subsequent bite [12]. The EIP plays an important role in shaping transmission rates [13], with longer time windows reducing the number of opportunities for pathogen transmission over a mosquito’s lifetime.
The kinetics of infection are equally influenced by genetic variation in the virus. Dengue fever is caused by four different serotypes of virus (DENV 1–4) [14] that are 65–70% similar at the DNA sequence level across their ~11-kb genomes [15], while also exhibiting a high average within-serotype diversity of ~3% at the amino acid level [16]. Comparisons between strains within single serotypes, have revealed variation in both infection rates and EIP in mosquitoes [17, 18], most likely due to differences in viral replication rates [11]. In humans, genetic variation within and between serotypes also determines relative viral fitness (replication rate), epidemic potential and virulence [19, 20], with particular strains linked to more severe clinical manifestations. The intra-host/vector diversity of DENV can also play a role in transmission, such as variants with a replicative advantage can spread more rapidly overall, eventually displacing those with lower fitness [21]. For example, an uncharacteristically large outbreak of dengue in Cairns, Australia in 2008/2009 was attributed to the very short EIP of the DENV-3 strain in the mosquito [22].
Surprisingly, little is known about DENV kinetics in the mosquito and how the virus interacts with individual tissues during infection. Mosquitoes have evolved both systemic and tissue specific antiviral mechanisms to limit viral replication [23]; thus, viral kinetics may differ between tissues. Given that the virus moves in a stepwise fashion, selective pressures in an initial tissue might therefore have cascading effects on the viral kinetics in downstream tissues [24]. Additionally, different tissue types may offer a diversity of cell types and cellular niches that may vary in their capacity to support DENV replication. For example, the midgut epithelium is a dynamic niche, with cells being shed frequently [25], whereas other cell types in intermediate tissues may be more stable sites of virus production.
A better understanding of intra vector kinetics will assist with developing optimal intervention points for genetic modification, improve our model of vectorial capacity [13] by adding components of dose and serotype, etc. [26] and, expand our understanding of DENV epidemiology, and vector-virus interactions such as peaks in viral load and latency periods in the host. Herein we compared midgut, carcass and salivary gland loads for 4 strains of DENV representing each of the 4 serotypes, fed at either high or low infectious doses, daily over a period of 3 weeks. In so doing, we were able to assess how the factors of dose and serotype strain diversity define susceptibility, EIP and transmissibility but also how tissue specific differences and the inter relationships between tissue kinetics drive transmissibility.
Results
We orally challenged inbred wild-type Ae. aegypti mosquitoes with two infectious doses (108 and 105 DENV copies/ml), representative of plasma viremia ranges in humans [27], of the 4 DENV serotypes. In each vector competence experiment, mosquitoes were blood fed with a strain representing each of the 4 DENV serotypes at the 2 infectious doses. Ten individuals were then collected daily for 20 days, starting from the first day post-infection (DPI), to measure DENV infection status and infection kinetics in 3 key tissues—midgut (MG), carcass (CA) and salivary glands (SG)—per mosquito.
The proportion of DENV-infected mosquitoes is affected by both infectious dose and DPI
Infection status was determined by analyzing the proportion of individuals positive for DENV in each tissue at each DPI, based on the presence of a positive qRT-PCR. Infection prevalence (Fig 1) in the MG, CA and SG in all the serotypes was significantly influenced by the infectious dose and DPI (Table 1). In only one case was there a significant interaction between these 2 factors (Table 1, DENV-3 SG, GLM, F = 5.23, p = 0.02), indicating that, in general, dose and age of infection act independently to shape infection status.
Mosquito populations fed with higher doses had higher infection rates across tissues and strains (Fig 1). The DENV-1 and -2 strains exhibited a higher infection rate in all tissues and at both infectious doses compared to the DENV-3 and -4 strains. DENV-4 exhibited especially low infection rates in all tissues at the low infectious dose. The significance of DPI demonstrates that infection status is heavily time-dependent. For the DENV-1 and -2 strains, in the high infectious dose, infection rates were constantly high throughout 20 DPI in contrast to the low infectious dose, where there is noticeable difference between the infection rates between early and late DPI. For the DENV-3 strain, both infectious doses peaked in infection rate at ~12 days in all tissues. In the DENV-4 strain, infection levels are weaker than the other strains, which may contribute to the greater variation in the effect of time across tissues and doses.
EIP is affected by the serotype strain and initial infectious dose
If we take SG detection of viral load as a proxy for presence of virus in the saliva, the susceptibility data (Fig 1) can also be used to shed light on EIP. EIP can be measured in multiple ways, including the day of first arrival of virus in the tissue/saliva or the point at which 50% of the population has infected tissue/saliva. For all strains and infectious doses, viral load can be detected in the SG already by day 2. In many cases, especially for low doses, viral load appears in the SG and then detection is lost for several days before returning and becoming stable (Fig 1). If we calculate an EIP50 for populations where infection continues for more than two days, we see that infectious dose and serotype strain diversity are predictors; for the strains representing DENV 1–4 at low-dose, they are 5, 6, 9 and 15 days, respectively, while for the high-dose they are 1, 3, 7 and 7 days.
Intra-host DENV kinetics are influenced by initial infectious dose and strain
To quantify the kinetics of DENV infection, we assessed whether the cumulative change in DENV load in each tissue followed a trend that was either a sigmoidal or linear over time. Thus, we first fit a 3-parameter logistic model to DENV loads in each tissue, dose and serotype combination, which modeled the upper asymptote or the maximum DENV load, the slope or growth rate at the midpoint, and the midpoint or infection age which is halfway between the lower and upper asymptotes (i.e. ED50). We then estimated separate trends for each treatment combination. From the candidate models, the best fitting model for each treatment combination was selected using Akaike’s information criterion (AIC, S4 Table). For MG and CA tissues the best fitting model was the 3-parameter logistic model, while for the SG tissue we identified a linear relationship between DPI and DENV load. We used the slope and maximum DENV load parameters to compare the behavior of the different serotypes at different infectious doses and across tissues and times post-infection (Fig 2).
Some treatments do not show substantial DENV replication and others exhibit subpopulation structure of DENV load in infected mosquitoes
For a number of treatments DENV loads did not rise to very high levels (Fig 2). For example, in the low-dose categories neither the DENV-3 nor -4 strains reached a load higher than log 103 in any tissue as measured by population mean. In these cases, although we were able to fit a dose response curve (DRC) to the data, the parameters of max load and growth rate are not informative (S1 Table). Therefore, the subsequent discussion of growth parameter estimates focuses only on those treatments with substantial DENV loads and replication (Tables 2 and 3). For the remaining treatments we aimed to fit a single DRC to the different treatment combinations (infectious dose, tissue and strain). In several treatments, however, we identified potential subpopulation structure, with a proportion of the mosquitoes exhibiting high loads and the rest exhibiting very low loads. Those treatments were the low infectious dose in the MG and CA for strains of serotype DENV-1 and -2 and the high infectious dose in the MG and CA for strains of serotypes DENV-3 and -4 (Fig 2). After noticing subpopulation structure in our data, we applied Hartigan’s dip test for all datasets. For any with a significant result, we then used a bimodality coefficient test to partition the data. We found that using a coefficient of bimodality of >0.75 effectively split our datasets into two clear subpopulations (S2 Table) using the following criteria: (1) for the y-axis, we selected the lowest load in the histogram located between the two highest peaks (S1 Fig); and (2) for the x-axis, we identified the DPI where the two populations began to diverge.
DENV load is always higher in the midgut than other tissues for high infectious dose
In all serotypes, maximum DENV loads were consistently higher in the MG than the CA. Because we fit linear models to the SG data, we examined the average DENV load at the last DPI (20 days) in lieu of a maximum load (S2 Fig). For all serotypes and high infectious doses, the SG load at this point was lower than the CA loads. This strong pattern suggests a “trickledown” effect for DENV load across the sequentially infected tissues (Fig 3). In contrast, low infectious doses resulted in low DENV loads in all tissues and for all strains, and DENV-3 and -4 could not be evaluated due to their lack of successful infection and failure of the selected model to run. Additionally, we ran generalized linearized models on DENV load in the MG and CA (for high dose only) as a predictor for the SG DENV load. We found no significant predictive ability regardless of the serotype and tissue (S5 Table).
Tissue-specific DENV kinetics vary by strain in the high infectious dose
In the MG, the maximum loads of DENV-2 (12.2 log c/ng) and -3 (11.5 log c/ng) are similarly high and greater than the strains representing serotypes 1 (9.11 log c/ng) and 4 (9.74 log c/ng) (Table 2). In the CA, the DENV-1 (8.99 log c/ng) and -2 (9.53 log c/ng) strains have the highest loads and are both significantly different from DENV-3 (8.85 log c/ng) and -4 (6.44 log c/ng) strains (Table 3). In the SG, the average final load (DPI 20) for the DENV-2 strain is not different from the DENV-1 and -3 strains but is greater than the DENV-4 (S2 Fig). Although we have examined only one representative strain from each serotype, these patterns suggest that DENV-2 may be better at replicating in mosquitoes and that the other serotypes do not show consistent patterns of dominance over one another across tissues. For growth rate, no one strain appeared to dominate over the others across multiple tissues (Tables 3 and 4). For example, the DENV-3 strain has the highest growth rate (change in viral copies/day) in the MG (0.66 log DENV copies per DPI), whereas DENV-2 strain has the lowest (0.31 log DENV copies per DPI). In the CA, DENV-4 has the highest growth rate (1.36 log DENV copies per DPI), whereas DENV-1 has the lowest (0.21 log DENV copies per DPI). Interestingly, in the SG all strains have distinctly different growth rates, ranked in the following order: DENV-1, 2, 3, and 4 (Table 4). This suggests that the growth rates in early process/upstream tissues are not predictive of rates in the final tissue controlling transmissibility.
The DENV-2 strain has a higher maximum DENV load than the DENV-1 strain in the low infectious dose
For the low infectious dose, we only considered mosquitoes infected with either DENV-1 or -2 strains, as the remaining strains did not exhibit “successful” infections meaning that even though infection rate was high, the actual DENV load was substantially low. Here, we saw that the DENV-2 strain had a significantly higher max DENV load in the MG and a higher average final load in the SG than the DENV-1 strain. For growth rate, there were no differences between the DENV-1 and DENV-2 strains in the MG or CA, but in SG, the DENV-1 strain was significantly higher (Table 5).
Discussion
To explore the role of infectious dose, serotype and tissue in viral infection kinetics we sampled DENV loads in populations of infected mosquitoes over numerous, sequential time-points of the mosquito’s lifespan relevant to transmission in the field. We reveal that the kinetics of DENV infection in the midgut, carcass and salivary glands of the mosquito Aedes aegypti are strikingly different among the strains selected for this study, and that these differences are also driven by the initial infectious dose. Specifically, we showed that (1) initial infectious dose dictates infection frequency, with the DENV-1 and -2 strains analyzed here infecting a greater proportion of mosquitoes than the DENV-3 and -4 strains; (2) for DENV-3 and -4 strains only, the high infectious dose led to strong infection, with lower doses producing lower viral replication; (3) viral growth rates in particular tissues did not predict the total DENV load in that same tissue, with some strains having higher growth rate but lower total DENV load; (4) in the salivary glands all strains have independent growth rates not predicted by growth rates in previous tissues; and (5) for every tissue and infectious dose, the DENV-2 strain reached a higher viral load according to both the average final load at 20 DPI and max DENV load.
The initial infectious dose has been shown to influence the likelihood of a mosquito becoming infected following an infectious blood meal and whether there is ultimately dissemination to the salivary glands after a blood meal [28, 29]. The amount of virus in human blood required to infect ~50% of a mosquito population varies between serotypes, with DENV-1 and -2 requiring ~10-fold less than DENV-3 and -4. This same rank order is observed here. However, strain variation within each serotype is substantial [30], so further comparisons across multiple strains within each serotype are clearly needed to draw any conclusions about whether the strains we have selected are predictive of the entire serotype.
We were also able to identify a potential threshold of virus needed in the blood meal to initiate successful infections and assure likely transmissibility via the salivary glands. A similar threshold effect has been observed in Drosophila melanogaster populations challenged with systemic bacterial infections [31].The study identified two possible outcomes of infection; (1) where hosts died with high bacterial burden or (2) where hosts survived with chronic infection and low pathogen load. Bacterial load and the time it took the host to establish an immune response were found to determine the trajectory of infection. Similarly, we saw either high infection (leading to transmissibility) or low/loss of infection. Interestingly, for our study, the threshold effect was driven by the initial infectious dose and by serotype. DENV-1 and -2 strains required less virus to establish an infection, whereas, DENV-3 and -4 strains required substantially more. The failure to initiate a successful infection could be due to poor replication in the midgut or the midgut preventing dissemination [10, 11, 32]. This pattern is not unique to the midgut, however, with low and inefficient replication also in disseminated tissues, such as the carcass. The data may simply be the result of low population sizes of virus in combination with the large number of vacant cellular niches. During infection DENV selectively annexes and manipulates host metabolism and machinery to increase viral translation and replication [33–35], while also enhancing replication of mosquito cells to increase niche availability in the mosquito [36]. The more virions that infect a given cell, the more efficient this process becomes and, thus, more viruses are produced. Initial viral load combined with inherent replicative differences between strains is likely to produce the pattern we see in our data, that the DENV-1 and 2 strains require a lower dose to produce successful infections.
Interestingly, higher viral growth rates in tissues did not necessarily lead to higher DENV loads. For example, for the high infectious dose of DENV-1, the midgut had the highest max DENV load compared to the carcass, however its growth rate was smaller. This relationship was seen for the other treatments as well. It is possible that the migration of viruses between tissues via the tracheal system or hemolymph [37] or directly from nearby tissues is impacting the total viral counts. In support of tracheal or hemolymph travel is the initial spike in infection rate in salivary glands followed by a reduction to zero and then a second front of infection several days later. It is possible that the latter wave represents the larger migration of virus from bodily tissues that then establishes permanent infection. Given how the circulatory system of the mosquito works, viral particles could be carried from tissues that are not in immediate proximity to the salivary glands. This process could directly affect growth rate in upstream tissues. The salivary glands also represent a special case of the entry and exit balance, with virus being excreted into the saliva and hence lost from tissue measures. Under a more natural setting additional factors would likely affect kinetics including, exposure to non-infectious blood meals prior to consumption of an infectious [38], consumption of sequential dengue infectious meals [39] or interactions between DENV and other viruses present in blood meals or the mosquito body [40].
The disconnect between kinetics in intermediate tissues and those of the salivary glands suggests that there are either stochastic processes driving kinetics in early tissues that do not predict how many viruses make it to the salivary glands, or that factors affecting replication rate in the salivary glands are somewhat independent [23]. For example, previous research has demonstrated that the strength and efficacy of immunity varies across tissues [23] and in response to different serotypes [41]. It is also possible that DENV is modulating the vector immune response in a strain- or tissue-dependent manner [42]. Differences between the salivary glands and other tissues may also be due to the number of cellular niches, and/or unique replication rates in those niches. For example, in contrast to the midgut epithelium, DENV tropism in the salivary glands appears to be heterogenous within the tissue; this might be due to the distribution of cellular receptors [43], with a possible higher concentration of viral entry receptors in the lateral distal lobes of the tissue [44]. In future studies, it would be useful to quantify the production rate of virus using negative strand strain-specific PCR methods [45] in the midgut and in the salivary glands relative to the availability of cellular niches.
The differences in viral kinetics between the DENV strains analyzed here, particularly in the salivary glands, might explain disparities in transmission potential in the field. A range of studies have linked different genotypes to overall dominance in terms of circulation in human populations both at local and global scales [46–49]. Additionally, DENV variants with a replicative advantage in both the host and vector can spread more rapidly, eventually displacing those with lower fitness [21]. Our findings indicate that 2 of the 4 DENV strains studied here show greater capacity to replicate in the mosquito. In some cases, these differences may have consequences for public health, as some specific variants of DENV have been linked to more severe clinical manifestations than others [50–53]. These findings also have implications for the mosquito’s contribution to the DENV epidemiological landscape, with the understanding that aspects of human immunity, including cross-reactivity and non-specific serotype responses also shape the nature of circulating variants [54]. Lastly, our work provides better insight into ideal design of vector competence experiments and the selection of gene candidates for the engineering of virus refractory mosquitoes. Specifically, our findings support a greater focus on salivary gland or saliva-based measures of infection and genes that either completely reduce viral loads in the midgut or that control viral replication in the endpoint tissue, salivary glands.
Materials and methods
Ae. aegypti rearing and virus preparation
Approximately 8000 eggs of an Ae. aegypti inbred line collected in Townsville, Australia 13 generations previously were vacuum hatched and divided into individual 30 × 40 × 8-cm plastic trays containing 1 liter of reverse osmosis (RO) autoclaved water. Larvae (~150) were placed in individual trays containing 3 liters of RO autoclaved water, supplemented with common fish food (Tetramin®, Melle, Germany). Pupae were collected and placed into breeding cages, containing approximately 450 mosquitoes each. Adults were provided 10% sucrose ad libitum. All mosquitoes were reared in a controlled environment at 26°C, 75% relative humidity and a 12-hr light/dark cycle.
The DENV serotypes/strains used for this experiment are listed in Table 6. The virus was propagated in cell culture, as described previously [55]. Ae. albopictus C6/36 cells were grown at 26°C in RPMI 1640 medium (Invitrogen, Carlsbad, CA) supplemented with 10% fetal bovine serum (FBS), 1× Glutamax (Invitrogen) and HEPES buffer. Cells were first allowed to form monolayers of around 60–80% confluence in T175 flasks (Sigma Aldrich, St. Louis, MO), and then inoculated with DENV and maintained in RPMI medium supplemented with 2% FBS. After 7 days post-inoculation, live virus was harvested, titrated via absolute quantification PCR and adjusted to two final viral loads of 105 and 108 DENV copies per ml. Single-use aliquots were stored at −80°C for subsequent plaque assay titration.
Mosquito infections
Prior to infection, female mosquitoes were sorted and placed in 1-liter plastic cups with a density of ~150 individuals. Sucrose was removed from the mosquitoes 24 hrs, prior to oral infection. Double-chamber glass feeders were covered with pig intestine previously immersed in a 10% sucrose solution. Water heated to 37°C was circulated in the outer chamber of the feeders, and a 1:1 mix of defibrinated sheep blood and the previously titrated DENV virus was placed inside the feeder. The mosquitoes were split into 2 feeding groups, each allowed to feed for ~2 hrs. One group was fed with a DENV infectious blood meal concentration of 108 and the other with 105 DENV copies/ml. After 24 hrs, blood fed mosquitoes were identified by visual inspection and separated into cups of 10 individuals.
Mosquito dissection and RNA extraction
Blood fed Ae. aegypti mosquitoes were collected daily, 10 per feeding group, over 20 days. Mosquitoes were anesthetized by chilling and dissected on a chill plate in a drop of sterile phosphate buffer saline (PBS). Midgut, salivary glands and carcass tissues were placed in separated 1.5-ml microcentrifuge tubes (Sarstedt, Nümbrecht, Germany) containing 200 μl of TRIzol reagent (Invitrogen, Carlsbad, CA, USA) and 2-mm glass beads. Samples were homogenized and frozen at −80°C. RNA extraction was performed using the TRIzol RNA extraction protocol according to the manufacturer’s instructions, and RNA was eluted in 25 μl of nucleic acid-free water. Samples were DNase treated (Life Technologies, Carlsbad, CA, USA) according to the manufacturer instructions [56]. Total RNA was determined with a NanoDrop lite spectrophotometer (Thermo Scientific, Waltham, MA).
DENV absolute quantification via qRT-qPCR
Quantitative real-time PCR (qRT-PCR) reactions for DENV detections were performed with 4× TaqMan fast virus 1-step master mix (Roche Applied Science, Switzerland), PCR-grade water, 10 μM of DENV primers and probe and 2.5 μl of RNA in a final volume of 10 μl. Reactions were run in microseal 96 microplates (Life Technologies Life, Technologies, Carlsbad, CA) covered with optically clear film. LightCycler 480 (Roche Applied Science, Switzerland) thermal cycling conditions were 50°C for 5 min for reverse transcription, 95°C for 10 s for RT inactivation/denaturation followed by 50 amplification cycles of 95°C for 3 s, 60°C 30 s and 72°C for 1 s.
Standard curves were generated from triplicate samples on each plate spanning the range of 10 to 108 copies/reaction of DENV fragment copies. DENV standards were constructed as described elsewhere [57], and they contained the 3’UTR of the DENV genome. The limit of detection was set at 10 copies in this experiment; water was used as a negative control and standards were run in triplicates. The concentration of DENV genome copies in each sample was extrapolated from the standard curve as DENV copies per nanogram of total RNA.
The primer sequences used for the detection of DENV, as previously used [58], were F: 5′-AAGGACTAGAGGTTAGAGGAGACCC-3′, R: 5′CGTTCTGTGCCTGGAATGATG-3′ and P: 5′-HEX-AACAGCATATTGACGCTGGGAGAGACCAGA-BHQ1-3′.
Survival assay
A second set of Townsville mosquitoes were reared in the same conditions as above. Mosquitoes were blood fed with the two viral concentrations as above (108 and 105 DENV copies/ml). A total of 80 fed mosquitoes per concentration were separated into clear 200-ml cups, each containing 10 individuals. They were monitored daily to assess death until all had expired. Neither DENV serotypes nor infectious doses had an effect on survival compared to a blood-only fed control (S3 Fig, S3 Table).
Data analysis
Data analysis was carried out in R v 3.6.0 (http://www.r-project.org/). For statistical analysis of the infection frequency, quasibinomial GLM was used to correct for overdispersion of data, and Tukey for contrasts was used for post hoc comparisons. For the kinetics of DENV virus, first we fit a 3-parameter logistic DRC model to DENV loads in each tissue, dose and serotype combination, which modeled the upper asymptote or the maximum DENV load, the slope or growth rate at the midpoint, and the midpoint or infection age which is halfway between the lower and upper asymptotes (i.e. ED50). We then estimated separate trends for each treatment combination. From the candidate models, the best fitting model for each treatment combination was selected using Akaike’s information criterion (AIC). For MG and CA tissues the best fitting model was the 3-parameter logistic model, while for the SG tissue we identified a linear relationship between DPI and DENV load. DRC parameters and pairwise comparisons were estimated with the DRC package [59]. The best model was selected based on the AIC (yielding the lowest) values. In general, for the salivary glands, the best fit was a linear model; thus, they were subsequently analyzed in this fashion. For the average DENV load at DPI 20, significant differences were based on Tukey post hoc comparison following ANOVAs. All DENV loads were reported on a log scale given the value range.
Supporting information
S1 Table [docx]
DRC estimates for unsuccessful infections by tissue, infectious dose and DENV strain.
S2 Table [docx]
Multimodal distribution analysis.
S3 Table [docx]
Survival curve statistical analysis.
S4 Table [docx]
AIC values for the DRC model.
S5 Table [docx]
Results of the GLM fitted to salivary glands DENV load.
S1 Fig [docx]
Density plots.
S2 Fig [docx]
Average salivary gland DENV load at the last DPI (20 days).
S3 Fig [docx]
Survival curves for all treatments.
Zdroje
1. Murray N.E., Quam M.B., and Wilder-Smith A., Epidemiology of dengue: past, present and future prospects. Clin Epidemiol, 2013. 5: p. 299–309. 10.2147/CLEP.S34440 doi: 10.2147/CLEP.S34440 23990732
2. Bhatt S., Gething P.W., Brady O.J., Messina J.P., Farlow A.W., Moyes C.L., Drake J.M., Brownstein J.S., Hoen A.G., Sankoh O., Myers M.F., George D.B., Jaenisch T., Wint G.R., Simmons C.P., Scott T.W., Farrar J.J., Hay S.I., The global distribution and burden of dengue. Nature, 2013. 496(7446): p. 504–507. 10.1038/nature12060 doi: 10.1038/nature12060 23563266
3. Brady O.J., Gething P.W., Bhatt S., Messina J.P., Brownstein J.S., Hoen A.G., Moyes C.L., Farlow A.W., Scott T.W., Hay S.I., Refining the global spatial limits of dengue virus transmission by evidence-based consensus. PLOS Negl Trop Dis, 2012. 6(8): p. e1760. 10.1371/journal.pntd.0001760
4. Gubler D.J., Dengue, urbanization and globalization: the unholy trinity of the 21(st) century. Trop Med Health, 2011. 39(4 Suppl): p. 3–11. 10.2149/tmh.2011-S05 doi: 10.2149/tmh.2011-S05 22500131
5. Simmons C.P., Farrar J.J., Nguyen V.V., Willis B., Dengue. New Eng J Med, 2012. 366(15): p. 1423–1432. 10.1056/NEJMra1110265 doi: 10.1056/NEJMra1110265 22494122
6. Powell J.R. and Tabachnick W.J., History of domestication and spread of Aedes aegypti—A Review. Mem Instit Oswaldo Cruz, 2013. 108(Suppl 1): p. 11–17. 10.1590/0074-0276130395
7. Colón-González F.J., Fezzi C., Lake I.R., Hunter P.R., The Effects of Weather and Climate Change on Dengue. PLOS Negl Trop Dis, 2013. 7(11): p. e2503. 10.1371/journal.pntd.0002503 doi: 10.1371/journal.pntd.0002503 24244765
8. Lequime S., Fontaine A., Ar Gouilh A., Moltini-Conclois I., Lambrechts L., Genetic drift, purifying selection and vector genotype shape dengue virus intra-host genetic diversity in mosquitoes. PLOS Genet, 2016. 12(6): p. e1006111. 10.1371/journal.pgen.1006111 doi: 10.1371/journal.pgen.1006111 27304978
9. Franz A.W.E., Kantor A.M., Passarelli A.L., Clem R.J., Tissue barriers to arbovirus infection in mosquitoes. Viruses, 2015. 7(7): p. 3741–3767. 10.3390/v7072795 doi: 10.3390/v7072795 26184281
10. Khoo C.C., Doty J.B., Held N.L., Olson K.E., Franz A.W., Isolation of midgut escape mutants of two American genotype dengue 2 viruses from Aedes aegypti. Virol J, 2013. 10(1): p. 1–11. 10.1186/1743-422X-10-257
11. Cox J., Brown H.E., and Rico-Hesse R., Variation in vector competence for dengue viruses does not depend on mosquito midgut binding affinity. PLOS Negl Trop Dis, 2011. 5(5): p. e1172. 10.1371/journal.pntd.0001172 doi: 10.1371/journal.pntd.0001172 21610852
12. Black W.C., Bennett K.E., Gorrochotegui-Escalante N., Barillas-Mury C.V., Fernandez-Salas I., de Lourdes Munoz M., Farfan-Ale J.A., Olson K.E., Beaty B.J., Flavivirus susceptibility in Aedes aegypti. Arch Med Res, 2002. 33(4): p. 379–388. doi: 10.1016/s0188-4409(02)00373-9 12234528
13. Macdonald G., The epidemiology and control of malaria. 1957, Amen House, Warwick Square, London E.C.4.: Oxford University Press. xiv + 201 + xl + 11 pp.
14. Deen J.L., Harris E., Willis B., Balmaseda A., Hammond S.N., Rocha C., Dung N.M., Hung N.T., Hien T.T., Farrar J.J., The WHO dengue classification and case definitions: time for a reassessment. The Lancet, 2006. 368(9530): p. 170–173.
15. Perera R. and Kuhn R.J., Structural proteomics of dengue virus. Curr Op Microbiol, 2008. 11(4): p. 369–377. 10.1016/j.mib.2008.06.004
16. Azhar E.I., Hashem A.M., El-Kafrawy S.A., Abol-Ela S., Abd-Alla A.M.M., Sohrab S.S., Farraj S.A., Othman N.A., Ben-Helaby H.G., Shshi A., Madani T.A., Jamjoom G., Complete genome sequencing and phylogenetic analysis of dengue type 1 virus isolated from Jeddah, Saudi Arabia. Virol J, 2015. 12: p. 1. 10.1186/s12985-014-0235-7 doi: 10.1186/s12985-014-0235-7 25591713
17. Anderson J.R. and Rico-Hesse R., Aedes aegypti vectorial capacity is determined by the infecting genotype of dengue virus. Am J Trop Med Hyg, 2006. 75.
18. Armstrong P.M. and Rico-Hesse R., Efficiency of dengue serotype 2 virus strains to infect and disseminate in Aedes aegypti. Am J Trop Med Hyg, 2003. 68.
19. Lambrechts L., Fansiri T., Pongsi A., Thaisomboonsuk B., Klungthong C., Richardson J.H., Ponlawat A., Jarman R.G., Scott T.W. Dengue-1 virus clade replacement in Thailand associated with enhanced mosquito transmission. J Virol, 2012. 86(3): p. 1853–1861. 10.1128/JVI.06458-11 doi: 10.1128/JVI.06458-11 22130539
20. Ty Hang V.T., Holmes E.C., Duong V., Nguyen T.Q., Tran T.H., Quail M., Churcher C., Parkhill J., Cardosa J., Farrar J., Willis B., Lennon N.J., Birren B.W., Buchy P., Henn M.R., Simmons C.P., Emergence of the Asian 1 genotype of dengue virus serotype 2 in Vietnam: in vivo fitness advantage and lineage replacement in South-East Asia. PLOS Negl Trop Dis, 2010. 4(7):e757. 10.1371/journal.pntd.0000757 doi: 10.1371/journal.pntd.0000757 20651932
21. Guzman M.G. and Harris E., Dengue. The Lancet, 2015. 385(9966): p. 453–465. 10.1016/S0140-6736(14)60572-9
22. Ritchie S.A., Pyke A.T., Hall-Mendelin S., Day A., Mores C.N., Christofferson R.C., Gubler D.J., Bennett S.N., van den Hurk A.F., An explosive epidemic of DENV-3 in Cairns, Australia. PLOS One, 2013. 8(7): p. e68137. 10.1371/journal.pone.0068137 doi: 10.1371/journal.pone.0068137 23874522
23. Cheng G., Liu Y., Wang P., Xiao X, Mosquito defense strategies against viral infection. Trends Parasitol, 2016. 32(3): p. 177–186. 10.1016/j.pt.2015.09.009 doi: 10.1016/j.pt.2015.09.009 26626596
24. Hall M. D., Bento G., and Ebert D. (2017). The evolutionary consequenses of stepwise infection processes. Tr Ecol Evol, 32(8) 612–623. 10.10.16/j.tree.2017.05.009
25. Zieler H., Zieler H., Garon C. F., Fischer E. R., Shahabuddin M. A tubular network associated with the brush-border surface of the Aedes aegypti midgut: implications for pathogen transmission by mosquitoes. J Exp Biol, 2000. 203(Pt 10): p. 1599–611. 10769222
26. Christofferson R.C. and Mores C.N., Estimating the magnitude and direction of altered arbovirus transmission due to viral phenotype. PLOS ONE, 2011. 6(1): p. e16298. doi: 10.1371/journal.pone.0016298 21298018
27. Duong Veasna., Lambrechts Louis., Paul Richard E., Ly Sowath., Lay Rath Srey., Long Kanya C., Huy Rekol., Tarantola Arnaud., Scott Thomas W., Sakuntabhai Anavaj., Buchy, Philippe, Asymptomatic humans transmit dengue virus to mosquitoes. PNAS USA 2015. 112(47): p. 14688–14693. 10.1073/pnas.1508114112 doi: 10.1073/pnas.1508114112 26553981
28. Gubler D.J., Nalim S., Tan R., Saipan H., Sulianti Saroso J., Variation in susceptibility to oral infection with dengue viruses among geographic strains of Aedes aegypti. Am J Trop Med Hyg, 1979. 28(6): p. 1045–1052. doi: 10.4269/ajtmh.1979.28.1045 507282
29. Rosen L., Rosenbloom L.E., Gubler D.J., Lien J.C., Chaniotis B.N., Comparative susceptibility of mosquito species and strains to oral and parenteral infection with dengue and Japanese encephalitis viruses. Am J Trop Med Hyg, 1985. 34(3): p. 603–615. doi: 10.4269/ajtmh.1985.34.603 2860816
30. Fontaine A., Lequime S., Moltini-Conclois I., Jiolle D., Leparc-Goffart I., Reiner R.C., Lambrechts L., Epidemiological significance of dengue virus genetic variation in mosquito infection dynamics. PLOS Pathog, 2018. 14(7): p. e1007187. 10.1371/journal.ppat.1007187 doi: 10.1371/journal.ppat.1007187 30005085
31. Duneau D., Ferdy J. B., Revah J., Kondolf H., Ortiz G. A., Lazzaro B. P., Buchon N. Stochastic variation in the initial phase of bacterial infection predicts the probability of survival in D. melanogaster. Elife, 2017. 6. 10.7554/eLife.28298
32. Taracena M.L., Bottino-Rojas V., Talyuli O.A.C., Walter-Nuno A.B., Oliveira J.H.M., Anglero-Rodriguez Y., Wells M.B., Dimopoulos G., Regulation of midgut cell proliferation impacts Aedes aegypti susceptibility to dengue virus. PLOS Negl Trop Dis, 2018. 12(5): p. e0006498. 10.1371/journal.pntd.0006498 doi: 10.1371/journal.pntd.0006498 29782512
33. Chotiwan N., Andre B.G., Sanchez-Vargas I., Islam M.N., Grabowski J.M., Hopf-Jannasch A., Gough E., Nakayasu E., Blair C.D., Belisle J.T., Hill C.A., Kuhn R.J., Perera R., Dynamic remodeling of lipids coincides with dengue virus replication in the midgut of Aedes aegypti mosquitoes. PLOS Pathog, 2018. 14(2): p. e1006853. 10.1371/journal.ppat.1006853 doi: 10.1371/journal.ppat.1006853 29447265
34. Choy M.M., Sessions O.M., Gubler D.J., Ooi E.E., Production of infectious dengue virus in Aedes aegypti is dependent on the ubiquitin proteasome pathway. PLOS Negl Trop Dis, 2015. 9(11): p. e0004227. 10.1371/journal.pntd.0004227 doi: 10.1371/journal.pntd.0004227 26566123
35. Reid D.W., Campos R.K., Child J.R., Zheng T., Chan K.W.K., Bradrick S.S., Vasudevan S.G., Garcia-Blanco M.A., Nicchitta C.V., Dengue virus selectively annexes endoplasmic reticulum-associated translation machinery as a strategy for co-opting host cell protein synthesis. J Virol, 2018. 92(7): p. e01766–17. 10.1128/JVI.01766-17 doi: 10.1128/JVI.01766-17 29321322
36. Helt A.-M. and Harris E., S-phase-dependent enhancement of dengue virus 2 replication in mosquito cells, but not in human cells. J Virol, 2005. 79(21): p. 13218–13230. doi: 10.1128/JVI.79.21.13218-13230.2005 16227245
37. Forrester N.L., Guerbois M., Seymour R.L., Spratt H., Weaver S.C., Vector-borne transmission iImposes a severe bottleneck on an RNA virus population. PLOS Pathog, 2012. 8(9): p. e1002897. 10.1371/journal.ppat.1002897 doi: 10.1371/journal.ppat.1002897 23028310
38. Zhu Y., Zhu, Yibin., Zhang, Rudian., Zhang, Bei., Zhao, Tongyan., Wang, Penghua., Liang, Guodong., Cheng, Gong. Blood meal acquisition enhances arbovirus replication in mosquitoes through activation of the GABAergic system. Nat Comm, 2017. 8(1): p. 1262. 10.1038/s41467-017-01244-6
39. Muturi E.J., Buckner E., and Bara J., Superinfection interference between dengue-2 and dengue-4 viruses in Aedes aegypti mosquitoes. Trop Med Int Health, 2017. 22(4): p. 399–406. 10.1111/tmi.12846 doi: 10.1111/tmi.12846 28150899
40. Le Coupanec A., Le Coupanec A., Tchankouo-Nguetcheu S., Roux P., Khun H., Huerre M., Morales-Vargas R., Enguehard M., Lavillette D., Misse D., Choumet V. Co-infection of mosquitoes with chikungunya and dengue viruses reveals modulation of the replication of both viruses in midguts and salivary glands of Aedes aegypti mosquitoes. Int J Mol Sci, 2017. 18(8). 10.3390/ijms18081708
41. Smartt C.T., Shin D., and Alto B.W., Dengue serotype-specific immune response in Aedes aegypti and Aedes albopictus. Mem Instit Oswaldo Cruz, 2017. 112(12): p. 829–837. 10.1590/0074-02760170182
42. Sim S. and Dimopoulos G., Dengue virus inhibits immune responses in Aedes aegypti cells. PLOS ONE, 2010. 5(5): p. e10678. 10.1371/journal.pone.0010678 doi: 10.1371/journal.pone.0010678 20502529
43. Cao-Lormeau V.-M., Dengue viruses binding proteins from Aedes aegypti and Aedes polynesiensis salivary glands. Virol J, 2009. 6(1): p. 35. 10.1186/1743-422X-6-35
44. Salazar M.I., Richardson J.H., Sanchez-Vargas I., Olson K.E., Beaty B.J., Dengue virus type 2: replication and tropisms in orally infected Aedes aegypti mosquitoes. BMC Microbiol, 2007. 7: p. 9. doi: 10.1186/1471-2180-7-9 17263893
45. Raquin V. and Lambrechts L., Dengue virus replicates and accumulates in Aedes aegypti salivary glands. Virol, 2017. 507: p. 75–81.
46. Dumre S.P., Bhandari R., Shakya G., Shrestha S.K., Cherif M.S., Ghimire P., Klungthong C., Yoon I.K., Hirayama K., Na-Bangchang K., Dengue virus serotypes 1 and 2 responsible for major dengue outbreaks in Nepal: clinical, laboratory, and epidemiological features. Am J Trop Med Hyg, 2017. 97(4): p. 1062–1069. 10.4269/ajtmh.17-0221 doi: 10.4269/ajtmh.17-0221 29031282
47. Shrivastava S., Tiraki D., Diwan A., Lalwani S.K., Modak M., Mishra A.C., Arankalle V.A., Co-circulation of all the four dengue virus serotypes and detection of a novel clade of DENV-4 (genotype I) virus in Pune, India during 2016 season. PLOS ONE, 2018. 13(2): p. e0192672. 10.1371/journal.pone.0192672 doi: 10.1371/journal.pone.0192672 29470509
48. Liebman K.A., Stoddard S.T., Morrison A.C., Rocha C., Minnick S., Sihuincha M., Russell K.L., Olson J.G., Blair P.J., Watts D.M., Kochel T., Scott T.W., Spatial dimensions of dengue virus transmission across interepidemic and epidemic periods in Iquitos, Peru (1999–2003). PLOS Negl Trop Dis, 2012. 6(2): p. e1472. 10.1371/journal.pntd.0001472 doi: 10.1371/journal.pntd.0001472 22363822
49. Sasmono R.T., Taurel A.F., Prayitno A., Sitompul H., Yohan B., Hayati R.F., Bouckenooghe A., Hadinegoro S.R., Nealon J., Dengue virus serotype distribution based on serological evidence in pediatric urban population in Indonesia. PLOS Negl Trop Dis, 2018. 12(6): p. e0006616. 10.1371/journal.pntd.0006616 doi: 10.1371/journal.pntd.0006616 29953438
50. Balmaseda A., Hammond S.N., Perez L., Tellez Y., Saborio S.I., Mercado J.C., Cuadra R., Rocha J., Perez M.A., Rocha C., Harris E., Serotype-specific differences in clinical manifestations of dengue. Am J Trop Med Hyg, 2006. 74: p. 449–456. 16525106
51. OhAinle M., Balmaseda A., Macalalad A.R., Tellez Y., Zody M.C., Saborio S., Nunez A., Lennon N.J., Birren B.W., Gordon A., Henn M.R., Harris E., Dynamics of dengue disease severity determined by the interplay between viral genetics and serotype-specific immunity. Sci Transl Med, 2011. 3(114): p. 114ra128. 10.1126/scitranslmed.3003084
52. Vaughn D.W., Green S., Kalayanarooj K., Innis B.L., Nimmannitya S., Suntayakorn S., Endy T.P., Raengsakulrach B., Rothman A.L., Ennis F.A., Nisalak A., Dengue viremia titer, antibody response pattern, and virus serotype correlate with disease severity. J Infect Dis, 2000. 181(1): p. 2–9. doi: 10.1086/315215 10608744
53. Halsey E.S., Marks M.A., Gotuzzo E., Fiestas V., Suarez L., Vargas J., Aguayo N., Madrid C., Vimos C., Kochel T.J., Laguna-Torres V.A., Correlation of serotype-specific dengue virus infection with clinical manifestations. PLOS Negl Trop Dis, 2012. 6(5): p. e1638. 10.1371/journal.pntd.0001638 doi: 10.1371/journal.pntd.0001638 22563516
54. R Reich Nicholas G., Shrestha Sourya., King Aaron A., Rohani Pejman., Lessler Justin., Kalayanarooj Siripen., Yoon In-Kyu., Gibbons Robert V., Burke Donald S., Cummings Derek A. T. Interactions between serotypes of dengue highlight epidemiological impact of cross-immunity. J Royal Soc Interface, 2013. 10(86). 10.1098/rsif.2013.0414 doi: 10.1098/rsif.2013.0414 23825116
55. Frentiu F.D., Robinson J., Young P.R. McGraw E.A., O'Neill S.L, Wolbachia-mediated resistance to dengue virus infection and death at the cellular level. PLOS ONE, 2010. 5(10): p. e13398. 10.1371/journal.pone.0013398 doi: 10.1371/journal.pone.0013398 20976219
56. Amuzu H.E., Simmons C.P., and McGraw E.A., Effect of repeat human blood feeding on Wolbachia density and dengue virus infection in Aedes aegypti. Parasit Vectors, 2015. 8: p. 246. 10.1186/s13071-015-0853-y doi: 10.1186/s13071-015-0853-y 25903749
57. Ye Y.H., Carrasco A.M., Frentiu F.D., Chenoweth S.F., Beebe N.W., van den Hurk A.F., Simmons C.P., O'Neill S.L., McGraw E.A., Wolbachia reduces the transmission potential of dengue-infected Aedes aegypti. PLOS Negl Trop Dis, 2015. 9(6): p. e0003894. 10.1371/journal.pntd.0003894 doi: 10.1371/journal.pntd.0003894 26115104
58. Ye Y.H., Ng T. S., Frentiu F. D., Walker T., van den Hurk A. F., O'Neill S. L., Beebe N. W., McGraw E. A., Comparative susceptibility of mosquito populations in North Queensland, Australia to oral infection with dengue virus. Am J Trop Med Hyg, 2014. 90(3): p. 422–30. 10.4269/ajtmh.13-0186 doi: 10.4269/ajtmh.13-0186 24420782
59. Ritz C., Baty F., Streibig J.C., Gerhard D., Dose-response analysis using R. PLOS ONE, 2016. 10(12): p. e0146021. 10.1371/journal.pone.0146021
Štítky
Hygiena a epidemiologie Infekční lékařství LaboratořČlánek vyšel v časopise
PLOS Pathogens
2019 Číslo 12
- Farmakovigilanční studie perorálních antivirotik indikovaných v léčbě COVID-19
- 10 bodů k očkování proti COVID-19: stanovisko České společnosti alergologie a klinické imunologie ČLS JEP
Nejčtenější v tomto čísle
- Coxiella burnetii Type 4B Secretion System-dependent manipulation of endolysosomal maturation is required for bacterial growth
- IL-22 produced by type 3 innate lymphoid cells (ILC3s) reduces the mortality of type 2 diabetes mellitus (T2DM) mice infected with Mycobacterium tuberculosis
- The pandemic Escherichia coli sequence type 131 strain is acquired even in the absence of antibiotic exposure
- A role of hypoxia-inducible factor 1 alpha in Mouse Gammaherpesvirus 68 (MHV68) lytic replication and reactivation from latency