Next Article in Journal
Intercropping Rosemary (Rosmarinus officinalis) with Sweet Pepper (Capsicum annum) Reduces Major Pest Population Densities without Impacting Natural Enemy Populations
Next Article in Special Issue
Effect of Selection for Pyrethroid Resistance on Abiotic Stress Tolerance in Aedes aegypti from Merida, Yucatan, Mexico
Previous Article in Journal
Comparison of New Kairomone-Based Lures for Cydia pomonella (Lepidoptera: Tortricidae) in Italy and USA
Previous Article in Special Issue
Modernizing the Toolkit for Arthropod Bloodmeal Identification
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

The Genetic Basis for Salivary Gland Barriers to Arboviral Transmission

Center for Vector-borne Infectious Diseases, Department of Microbiology, Immunology, and Pathology, Colorado State University, Ft. Collins, CO 80523, USA
*
Author to whom correspondence should be addressed.
Submission received: 2 December 2020 / Revised: 8 January 2021 / Accepted: 11 January 2021 / Published: 15 January 2021

Abstract

:

Simple Summary

Mosquito salivary glands are infected with an arbovirus prior to entering saliva. The virus in the mosquito’s saliva can then infect a vertebrate host when the vector acquires its next bloodmeal. Salivary gland infection and escape barriers (SGIB and SGEB, respectively) have been identified that modulate arbovirus transmission. SGIB are manifest as the absence of virus in the salivary glands of mosquitoes having a disseminated infection, while SGEB are evident as an absence of virus in the saliva of mosquitoes even though salivary glands are infected. The interaction between Aedes aegypti and viruses is dynamic and related to the genetic background of each vector population and virus variant. There is little understanding of the genetic basis for SGIB and SGEB. Here, nested, half-sib families of Aedes aegypti were used to estimate genetic and environmental variances, in which daughters from an individual dam differ in their phenotype due to the sire. We found that SGIB has a strong genetic basis with dengue virus infections but not with Zika or chikungunya virus infections, in which all salivary glands became infected. SGEB has a moderate genetic basis for Zika or chikungunya infections but not with dengue virus infections.

Abstract

Arthropod-borne viruses (arboviruses) infect mosquito salivary glands and then escape to saliva prior to virus transmission. Arbovirus transmission from mosquitoes can be modulated by salivary gland infection barriers (SGIBs) and salivary gland escape barriers (SGEBs). We determined the influence of SGIBs and SGEBs by estimating the quantitative genetic contributions of Aedes aegypti half-sib families (Mapastepec, Mexico) infected with three dengue 2 (DENV2), two chikungunya (CHIKV), and two Zika (ZIKV) genotypes. We determined virus titer per salivary gland and saliva at seven days post-infection and virus prevalence in the half-sib population. CHIKV or ZIKV genotypes did not present SGIB, whereas DENV2 genotypes showed low rates of SGIB. However, virus titer and prevalence due to additive genetic factors in the half-sib family displayed a significant narrow-sense heritability (h2) for SGIB in two of the three DENV2 genotypes and one CHIKV and one ZIKV genotype. SGEBs were detected in all seven virus strains: 60–88% of DENV2 and 48–62% of CHIKV or ZIKV genotype infections. SGEB h2 was significant for all CHIKV or ZIKV genotypes but not for any of the DENV2 genotypes. SGIBs and SGEBs exhibited classical gene-by-gene interaction dynamics and are influenced by genetic factors in the mosquito and the virus.

1. Introduction

Vector competence (VC) for an arbovirus is affected by intrinsic (genetic) factors and mechanisms that control the ability of a vector mosquito to acquire, maintain, and transmit an arbovirus. There are several tissue barriers in the mosquito that the virus must overcome in order to establish a productive infection: the midgut infection barrier (MIB), midgut escape barrier (MEB), salivary gland infection barrier (SGIB), and salivary gland escape barrier (SGEB). In mosquito vectors of disease, the salivary gland is the organ leading to transmission of an arbovirus from an infected mosquito [1]. In all mosquitoes, the two salivary glands are located in the anterior portion of the thoracic hemocele. In females, each gland consists of two lateral and one medial cylindrical lobe (Figure 1A). The distal end of each lobe is called the acinus, while the proximal end is called the neck (Figure 1A). Each lobe consists of a basal lamina bounding a single layer of acinar cells (Figure 1B). These cells are distributed around a central salivary duct containing an apical cavity for saliva storage. Apical plasma membranes have microvillar projections that are highly convoluted to increase the surface area for discharge of saliva into the apical cavity (diverticulum) (Figure 1B). Following successful dissemination from the midgut, arboviruses typically infect hemocytes, fat body, neural tissue, and occasionally muscle tissue before infecting the salivary glands. Virus amplification in fat body and specifically in hemocytes is an important prerequisite before the virus is ready to infect the salivary glands [2]. Virus in the hemolymph penetrates the basal lamina of the salivary glands surrounding acinar cells to generate a salivary gland infection (Figure 1C) [3]. Factors that prevent virus penetration of the salivary gland basal lamina constitute SGIB. Several studies suggest that tracheoles or nerve tissues may allow virus entry into acinar cells [4,5], where viral replication occurs. Barriers that prevent infection and virus replication in the salivary glands also constitute SGIB (SGI, Figure 1C). Viruses are eventually released in the apical cavities where saliva is stored prior to its release into the hypopharynx and deposition of virus into a vertebrate during probing and feeding. Barriers that block release of mature virus in the apical cavities are referred to as SGE1 and SGE2 (Figure 1C). Virus replication in the salivary glands has been reported for DENV [6], CHIKV [3], and ZIKV [7].
The existence of SGEB has been definitively demonstrated for Japanese encephalitis virus (JEV) in Culex tritaeniorhynchus [8], Snowshoe hare virus in Aedes triseriatus [9], La Crosse virus (LACV) in Aedes hendersoni [10], and Sindbis virus (SINV) in Culex theileri [11]. An exceptionally strong SGEB was detected in Ae. hendersoni for LACV, wherein 65% of females had salivary gland infections by plaque assay, but only 5% actually transmitted infectious virus [12]. These results were confirmed by immunofluorescence assays (IFA), which showed high concentrations of LACV antigens in the salivary glands. More recently, SGEB have been reported for Rift Valley fever virus (RVFV) in Culex nigripalpus, Anopheles crucians, and Aedes infirmatus [13,14]. A recent transmission electron microscopy study [3] nicely illustrated CHIKV replication and storage in Ae. albopictus salivary glands and showed that CHIKV nucleocapsids bud from two different sites in acinar cells to become mature particles. The first site is at vesicular membranes located in the apical cavity of acinar cells (SGE1, Figure 1C), while the second occurs at cellular plasma membranes (SGE2, Figure 1C). Mature viral particles clearly accumulated in acinar cell junctions (Figure 1D) [3].
Furthermore, it has been shown that different viruses preferentially replicate in different lobes of the salivary glands [15]. Arbovirus infection of salivary glands typically begins in the distal lateral lobes [15,16,17,18,19]. Arboviruses such as DENV2 and CHIKV infect the proximal lateral and medial [20] lobes of Ae. aegypti, whereas SINV does not seem to infect the median lobes of Aedes albopictus or Ae. aegypti [15,16,20,21]. Electron micrographic studies confirm that clumps of crystalline flavi- and alphaviruses [3,22] occur in the salivary glands and saliva. It has been reported that the lateral and medial salivary gland lobes interact with monoclonal antibodies differently, suggesting there is a biochemical difference between these two lobe surfaces [23]. Others have observed that specific proteins are differentially expressed by different salivary gland lobes [24], which could account for the observed differences between lobes or between mosquitoes in susceptibility to virus infection [2,25].
To reach the salivary glands, arboviruses must overcome natural barriers, infect various tissues, exploit physical conduits, and counteract innate and cellular defenses [15,26,27,28]. Meanwhile, vectorial capacity involves extrinsic factors, including temperature, availability of vertebrate hosts, feeding behavior, population density, longevity, and predation [29]. The interaction between Ae. aegypti and a virus is dynamic and is related to the genetic background of each particular vector population and arbovirus variant [30,31,32,33]. It is possible that the existence of mosquito genetic variants is related to these functions, which could affect the degree of susceptibility in mosquitoes [34]. There is little understanding of whether there is a genetic basis for SGIB and SGEB or whether they are random and subject to a few or many of the “environmental” factors listed in the previous paragraph. One of the principal applications of quantitative genetics is to determine the relative contributions of genetic and environmental factors to variation and covariation in one or more quantitative phenotypes [35]. Basically, the phenotypic variance (Vp) of a population is viewed as the sum of separate genetic (Vg) and environmental (Ve) variances. The genetic variance is in turn made up of additive (Va) and non-additive (dominance/epistatic) components. Va and the associated narrow-sense heritability (h2 = Va/Vp) can vary from 0 to 1, and indicate the proportion of phenotypic variance explained by additive genetic effects [35].
Herein, we assume that if the variances in SGIB and SGEB among siblings or between half-sib families are due to one or all of the random causes listed above, including variation in the amount of saliva produced by a female, then Ve will be large, Va will be small, and h2 will approach zero. Alternatively, if the amount of virus collected in saliva is due primarily to genetic, rather than environmental, factors, then the variation in SGIB and SGEB will have a large h2, approaching one. Herein, we perform a half-sib analysis to estimate the quantitative genetic and environmental basis of SGIB and SGEB for three DENV2 genotypes (American, Cosmopolitan, and Asian genotype), two CHIKV genotypes (Asian and ECSA-IOL genotype) and two ZIKV genotypes (Asian and African genotype). Basically, we are testing if, given all the environmental effects listed above, a significant genetic effect can be detected.

2. Materials and Methods

2.1. Cells, Viruses, and Mosquitoes

C6/36 cells (ATCC CRL-1660), Vero (ATCC CCL-81), and LLC-MK2 cells (ATCC number CCL-7.1) were maintained in minimum essential medium (MEM) supplemented with 7% heat-inactivated (56 °C, 30 min) fetal bovine serum (FBS, from Atlas Biologicals, cat# F-0500-A, Fort Collins, CO, USA), 1% L-glutamine solution (Glutamine from Mediatech, Corning catalog# 25-005-Cl, Tewksbury, MA, USA), 1% non-essential amino acids solution (Nonessential amino acids from Mediatech 25-025-Cl), and 1% penicillin/streptomycin (Pen-Strep from Mediatech, catalog# 30-001-Cl). The cells were incubated at 37 °C, 5% CO2. Virus isolates used in this study include DENV2 Quintana Roo 94 (Qro94, American genotype; Accession# DQ341196), DENV2 Merida BC17 (Cosmopolitan genotype; Accession# AY449677), DENV2 Guerrero C932 (Asian II genotype; Accession# AY449678), CHIKV R99659 (Asian genotype; Accession# KJ451624.1), CHIKV LR2006 OPY-1 (La Réunion, ECSA-IOL genotype; Accession # DQ443544.2), ZIKV PRVABC59 (Asian genotype; Accession# KU501215), and ZIKV Dakar 41525 (African genotype; Accession# KU955591). Hereafter, viruses are designated in the text as DENV2 American, DENV2 Cosmopolitan, DENV2 Asian, CHIKV Asian, CHIKV ECSA-IOL, ZIKV Asian, and ZIKV African, respectively. DENV2 isolates were endemic strains from Mexico. The virus strains used in this study have different levels of vector competence: low for DENV2 American, moderate for DENV2 Asian, CHIKV Asian, and ZIKV Asian, and high for DENV2 Cosmopolitan and ZIKV African. The American, Cosmopolitan, and Asian DENV2 isolates used here were available in our lab and afforded us the ability to determine heritability in the context of three DENV2 genotypes from Mexico. DENV2 were propagated in C6/36 cells at a multiplicity of infection (MOI) of 0.01. ZIKV and CHIKV were propagated similarly in Vero cells. Supernatant was harvested (3–5 days post-infection (dpi) for CHIKV and ZIKV and 10–12 dpi for DENV), clarified by centrifugation at 4 °C, and aliquoted into single-use vials before freezing at −80 °C. Aedes aegypti strain Mapastepec (15°26′44.05′′ N, 92°54′9.74′′ W), a domestic mosquito from an urban region in from the state of Chiapas, Mexico, was originally collected in 2016. Adult mosquitoes were fed on raisins and water, and the females were allowed to feed on uninfected blood meals to stimulate oogenesis. Mosquitoes were maintained at 28 °C, 80% relative humidity with 16–8 h light-dark periods. Half-sib families for this study were created using generations F4, F5, and F6.

2.2. Intrathoracic Inoculation

Intrathoracic injections of DENV, ZIKV, and CHIKV were conducted to assess virus replication in salivary glands and other mosquito tissues outside the midgut. Intrathoracic inoculation of mosquitoes with arboviruses usually leads to more rapid, and higher proportions of, disseminated infections, compared to oral challenge, and mosquitoes become infected with a constant dose of virus [14,15,36,37]. We acknowledge that intrathoracic injections bypass the midgut, arguably the most important barrier to mosquito infection. However, past experience [15] has shown that only a few orally-infected mosquitoes in a family actually go on to develop infected salivary glands. These small numbers preclude an accurate quantitative genetic analysis of SGIB and SGEB. Notably, the peak concentration of virus detected in the salivary glands is often similar, regardless of which method is initially used [15]. The method for intrathoracic inoculation of mosquitoes is described in [38] with some modifications. Adult females (7 days post-eclosion) were infected by injecting 500 plaque forming units (PFU) in a volume of 69 nL of the virus stock using a Nanojet II (Drummond Scientific Company, Broomall, PA, USA). All injections were performed under a dissecting microscope using glass needles that were prepared using a vertical pipette puller (P-30, Sutter Instrument Co., Novato, CA, USA). Mosquitoes were left for 7 days at 28 °C and 80% relative humidity before samples were collected.

2.3. Saliva Collection and Salivary Gland Dissections

At 7 dpi, saliva from each of 10–15 daughters/dam were collected [39]. Briefly, females were chilled, and their wings and legs removed. Their proboscises were inserted into a 1 µL micropipette (microcaps®, Drummond Scientific Company, PA, USA) filled with immersion oil type B, and the mosquitoes were allowed to salivate into the oil at room temperature. After 5 min (probing) and 45–60 min, each individual capillary was observed under the stereoscope Olympus SZH10 at 2X magnification for the presence of saliva. Mosquitoes that did not visibly salivate were eliminated. To determine the relationship between saliva volume and saliva titer, categories were assigned 1 (≈4.7 nL), 2 (≈7 nL), 3 (≈10.8 nL), or 4 (≈13.5 nL), based on the total observed volume compared with a reference capillary (Figure 2) [40].
The proboscis was removed from the capillary, and oil containing the saliva was expelled under pressure into 1.5 mL tubes containing 300 μL Dulbecco’s Modified Eagle Medium (DMEM supplemented with 20% of heat inactivated FBS, 1% penicillin/streptomycin, 1% glutamine, 1% non-essential amino acids) and frozen immediately on dry ice. After collecting saliva, salivary glands from mosquitoes were dissected and collected in 500 μL DMEM medium (7% heat inactivated FBS, 1% penicillin/streptomycin, 1% glutamine, 1% non-essential amino acids) and frozen immediately on dry ice. Virus in saliva and salivary glands were detected and quantified by plaque assay.

2.4. Infectious Virus Titration by Plaque Assay

We used LLC-MK2 cells for DENV2 plaque assays and Vero cells for ZIKV and CHIKV plaque assays. Plaque assays were performed on confluent monolayer cells in 24-well plates, infected with 10-fold serial dilutions of saliva samples or salivary gland homogenate that had been sterile-filtered with Acrodisc syringe filters with Supor membrane 0.2 μm (Pall Life Sciences, East Hills, NY, USA). After 1 h, infected cells were overlaid with a 1% agar–nutrient mixture [agar solution (1 g/79 mL ddH2O): nutrient solution (10 mL of 10×/100 mL Media 199, 5% of heat inactivated FBS, 7.5% sodium bicarbonate 4 mL/100 mL), 2% DEAE-dextrose in Hanks balanced solution (1 mL/100 mL), 0.5 mL/100 mL MEM essential amino acids (15× solution), and 0.5 mL/100 mL MEM Vitamins (100× solution)]. After 4, 6, or 10 days incubation at 37 °C (CHIKV, ZIKV, and DENV2 respectively), cells were stained with 150 µL/well of 3 mg/mL MTT (3-[4,5-dimethylthiazol-2-yl]-2,5-diphenyltetrazolium bromide) solution and incubated for at least 4 h [41,42]. Viral titers were determined by counting plaques. Individual saliva or salivary gland titers are reported as plaque forming units per milliliter (PFU/mL). All virus negative (saliva or salivary gland) samples were tested twice to be certain they were really negative. The limit of detection for the plaque assays was 1 PFU/mL for CHIKV and ZIKV virus and 1–10 PFU/mL for DENV.

2.5. Quantitative Genetic Analyses

A nested half-sib design was used to provide estimates of genetic and environmental variances and covariances [35,43,44]. A single male was mated to several (usually 3–5) females. The females were allowed to bloodfeed and lay eggs. Eggs from each female were collected, and 5 days later, were reared separately, awaiting adult eclosion (Figure 3). All 7-day-old daughters from an individual dam were intrathoracically inoculated with 500 PFU of virus, and after 7 days, were tested for the presence or absence of virus and the amount of virus in the salivary glands and saliva (Figure 3). Half-sib families differ in their phenotype due to the sire. If the SGIB or SGEB have a genetic basis, then sharing the same sire is why half-sibs are similar for the SGIB or SGEB phenotype. Variance in phenotype between half-sib families is equal to the covariance in phenotype among half-sibs. A high variance between families should cause high covariance within families, because members of a family are alike and the heritability approaches one. Conversely, a low variance in phenotype between families will lead to low covariance within families, and the heritability approaches zero.
Restricted maximum likelihood (REML) was used to estimate quantitative genetic parameters [44,45]. REML does not require balanced designs, and missing cells have a minimal effect. Furthermore, REML allows a diversity of hypotheses to be tested. The computational limits that formerly made REML impractical for large datasets have mostly been overcome. For these reasons, REML is now the method of choice for estimating quantitative genetic parameters from most designs. We used PROC MIXED in SAS (SAS User’s Guide, Cary, NC: Statistical Analysis System Institute, Inc., 1987) to obtain REML estimates and hypothesis tests for half-sib designs [44,45]. The SAS code for performing REML analyses appears in Figure S1. The total phenotypic variance (σ2P) is a function of three observed sources of variance: the average variance among offspring of a single sire is σ2S, the average variance among full-sib families from dams mated to a common sire is σ2d, and the average variance within members of a full-sib family is σ2w. Additive genetic variance (Va) is calculated as:
4 σ2S
The variance among dams (Vd) bred to a common sire is calculated as:
Vd = σ2d − σ2S
The variance among full-sibs is (Vw):
Vw = σ2w − 2 σ2S
Summing Equations (1)–(3),
Vp = Va + Vd + Ve = 4 σ2S + (σ2d − σ2S) + (σ2w − 2 σ2S) = σ2S + σ2d + σ2w
Narrow-sense heritability (h2) is:
h2 = Va/Vp
See SAS lines 5–10 for SGIB, and lines 24–30 for SGEB in Figure S1. The null hypothesis that σ2S = 0 was analyzed with a likelihood ratio test. The 2 restricted log likelihood for the full model (Equation (4)) was subtracted from 2 restricted log likelihood for the Va = 0 mode
Vp = Vd + Ve
See SAS lines 12–23 for SGIB and lines 31–42 for SGEB in Supplement 1. This is a likelihood ratio test because the difference in log likelihoods equals the log of the ratio of the likelihoods. There is one degree of freedom (d. f.) because Equations (4) and (6) differ by one parameter.

2.6. Statistical Analyses

All statistical analyses were performed with GraphPad Prism (version 5.0, La Jolla, CA, USA). A two-tailed Fisher’s exact test was used to compare infection rates for SGIB and SGEB and determine any correlation among virus titers in saliva and salivary glands (Pearson correlation). Analysis of variance (one-way and two-way ANOVA) was used to determine the statistical significance, and CORREL function (Excel 2016, Microsoft, Redmond, DC, USA) was used for determining the correlation coefficient. Significance was defined as p < 0.05.

3. Results

3.1. Infection Rates

Our breeding design was established with a minimum number of sires (15–20), a minimum of 2–3 half-sib families (dam/sire) and a minimum of ten F1 females/half-sib families. Under this criterion, we obtained 24 half-sib families for DENV2 American genotype. We also established 21 half-sib families for DENV2 Asian and Cosmopolitan genotypes, ZIKV Asian and African genotypes, and 22 half-sib families for CHIKV genotypes. Table 1 lists the number of paired samples tested for each of the seven viral strains. The number of salivary glands that were not infected divided by the total number of salivary glands tested is the SGIB%. The number of infected salivary glands that failed to produce any virus in the saliva are also listed. This number divided by the number of infected salivary glands is the SGEB%. The percentage transmission is 100 − SGEB% (Table 1).
These calculations assume that viral dissemination occurred in each intrathoracically-inoculated mosquito based on data that shows 100% infection of whole mosquitoes (n = 20–30) by plaque assay for each viral strain at 7 dpi (Figure 4).
Dissemination rates after intrathoracic inoculation were typically between 95–100% [46]. DENV2 American virus titer was significantly lower than other DENV genotypes as well as CHIKV and ZIKV titers (p < 0.0001). Significant differences were also detected in virus titer between CHIKV genotypes (p < 0.0001). CHIKV ECSA-IOL virus titers were higher (9 × 105 ± 9 × 104 PFU/mL) than CHIKV Asian virus titers (1.7 × 105 ± 1.9 × 105 PFU/mL). We observed no significant differences in virus titer between the two ZIKV genotypes (p = 0.8931) or between DENV2 Asian and Cosmopolitan genotypes (p = 0.6971) (Figure 4).

3.2. Salivary Gland Infection Barrier (SGIB)

CHIKV or ZIKV genotypes did not present SGIB, whereas DENV2 genotypes showed low rates of SGIB (Table 1). Specifically, after 7 dpi, 5.5% (38/691), 0.75% (5/667), and 0.59% (4/681) of Ae. aegypti salivary glands were uninfected with the DENV2 American, Asian, and Cosmopolitan genotypes, respectively. All salivary glands, on the other hand, were infected with CHIKV or ZIKV. Average viral titer (log10 (PFU/mL)) per salivary gland was calculated for all seven virus genotypes. Means and 95% confidence intervals are shown in Figure 5A.
When comparing genotypes across the same virus family (DENV2, CHIKV, and ZIKV), genotype titer means for DENV2 and ZIKV were significantly different (p < 0.0001). Of all viruses tested, DENV2 American had the lowest virus titer (2.699 ± 0.58 log10 (PFU/mL)), while ZIKV African had the highest titer (5.723 ± 0.331 log10 (PFU/mL)). DENV2 Cosmopolitan had the highest virus titer (4.742 ± 0.23 log10 (PFU/mL)) between DENV2 genotypes, while DENV2 Asian (3.607 ± 0.0421 log10 (PFU/mL)) was the second lowest of all viruses tested. CHIKV ECSA-IOL had slightly higher titers (4.921 ± 0.206 log10 (PFU/mL)) as compared to CHIKV Asian (4.687 ± 0.711 log10 (PFU/mL)). CHIKV genotypes had the same titers as DENV2 Cosmopolitan (p = 0.3294). Both ZIKV genotypes (Asian, 5.244 ± 0.136 log10 (PFU/mL); African, 5.723 ± 0.331 log10 (PFU/mL)) exhibited the highest titers. SGIB did affect the quantity of virus present in the salivary gland. It is very clear that different half-sib families can, and do, vary greatly in the amount of virus in salivary glands (Figure 6).
When comparing the spread of sire means between DENV2 genotypes, the greatest variation in log10 (PFU/mL) means/sire occurred with DENV2 American (Figure 6A) infection, and the smallest spread of sire means occurred with DENV2 Cosmopolitan infection (Figure 6C). Variation among sire families with CHIKV ECSA-IOL (Figure 6D) was slightly greater than that among CHIKV Asian families (Figure 6E). Variation among sire families with ZIKV Asian and ZIKV African were also small (Figure 6F,G). Narrow-sense heritabilities (h2) for SGIB were calculated to characterize each viral infection (Table 2) once the observational components σ2S, σ2d, and σ2w and the causal components of variance (Vp, Vsire, Vdam, Ve) for log10 (PFU/mL)/salivary gland were obtained (Table S1).
SGIB exhibited strong genetic bases with DENV American and Asian genotypes, but not with DENV2 Cosmopolitan genotype. For DENV2 American genotype and DENV2 Asian genotype, both the estimate of the variance among sire means and the h2 were significantly different (for American: σ2S = 0.2700, p = 0.0028; h2 = 0.8103, p = 5.34 × 10−7; for Asian: σ2S = 0.1029, p = 0.0348; h2 = 0.5799, p = 0.0072). However, for DENV2 Cosmopolitan genotype, the variance among sire means were small (0.0082) and non-significant (p = 0.3238). The narrow-sense heritability h2 was 0.0685 (p = 0.3274).
There were large differences between the CHIKV and ZIKV genotypes. SGIB exhibited strong genetic bases with CHIKV ECSA-IOL and ZIKV African, but not CHIKV or ZIKV Asian genotypes. For CHIKV ECSA-IOL and ZIKV African, the variance among sire means and the h2 were significantly different (for CHIKV ECSA-IOL: σ2S = 0.0337, p = 0.005; h2 = 0.5174, p = 6.2 × 10−6; for ZIKV African: σ2S = 0.0762, p = 0.0133; h2 = 0.3106, p = 0.0009). For CHIKV Asian genotype, the variance among sire means and the h2 were not significantly different (σ2S = 0.0012, p = 0.4537; h2 = 0.0385, p = 0.5). For ZIKV Asian genotype, the variance among sire means was 0.008 and not significant (p = 0.0918), and h2 was 0.1526 and marginally significant (p = 0.0416). There was a positive correlation between the variance among sire means and the narrow-sense heritability (h2). As predicted by Equation (5) (h2 = Va/Vp), the greater the variance in sire means, the larger the heritability. In other words, when the amount of virus was low (e.g., Figure 6A,B, DENV2 American and Asian), estimates of heritability and variance among sire means tended to be large and significant. However, once viral titers reach a high enough level in the salivary gland, (DENV2 Cosmopolitan, CHIKV, and ZIKV genotypes in Figure 6C–G), genes in the mosquito may no longer have an effect, and heritability estimates were small.
Variation among sire families (Figure 6) suggested that viral titer in the salivary gland is negatively correlated with the magnitude of heritability. We performed a correlation analysis between the average log10 PFU/mL virus in the salivary glands and the heritability to evaluate this observation (Figure 7). Pearson’s correlation coefficient was large and negative, as expected, when ZIKV African was excluded (R = −0.819, p = 0.0461). If ZIKV African was included, Pearson correlation coefficient was −0.24150 (p = 0.6019).

3.3. Salivary Gland Escape Barrier

In contrast to SGIB, all seven virus strains showed SGEB (Table 1). SGEB% (uninfected saliva/infected salivary glands) rates for DENV2 genotypes ranged, from 88.5% (578/653), 72.5% (480/662), 59.7% (404/677) for American, Asian, and Cosmopolitan, respectively. CHIKV genotypes showed 47.7% (327/686) and 61.5% (452/735) for ECSA-IOL and Asian, respectively. Finally, we obtained 52.1% (408/783) for ZIKV Asian and 32.1% (200/623) for ZIKV African. All three DENV2 genotypes varied significantly, as did the two CHIKV genotypes (p < 0.0001). DENV2 Cosmopolitan and CHIKV Asian did not significantly differ in SGEB% (p = 0.9447), nor did CHIKV ECSA-IOL or ZIKV Asian (p = 0.7332).
Clearly SGEB provides a barrier to transmission in mosquitoes with infected salivary glands. Of all viruses tested, DENV2 American had the lowest transmission rate (11.5%), while ZIKV African had the highest (67.9%) (Table 1). As with SGIB, average viral titer (log10(PFU/mL)) per saliva were calculated for all seven viral genotypes, and are shown in Figure 5B. When comparing genotypes across the same virus family genotypes, titer means for DENV2, CHIKV, and ZIKV were significantly different (DENV2, p <0.0001; CHIKV, p = 0.0129; ZIKV, p < 0.0001). Of all viruses tested, DENV2 American had the lowest average virus titer (0.2442 ± 0.039 log10 (PFU/mL)), while ZIKV African had the highest titer (1.538 ± 0.1302 log10 (PFU/mL)) in saliva. DENV2 Cosmopolitan had the highest virus titer (0.6324 ± 0.071 log10 (PFU/mL)) of the three DENV2 genotypes tested, while DENV2 Asian (0.4038 ± 0.055 log10 (PFU/mL)) was the second lowest of all viruses tested. CHIKV ECSA-IOL had higher titers (0.922 ± 0.085 log10 (PFU/mL)) as compared to CHIKV Asian (0.6191 ± 0.0.079 log10 (PFU/mL)). CHIKV Asian had the same titer as DENV2 Cosmopolitan (p = 0.9013). ZIKV Asian and CHIKV ECSA-IOL both had similar titers (p = 0.7030) (Figure 5B). Moreover, the variance among sire means in log10 (PFU/mL)/saliva was not significantly different between all three DENV2 genotypes tested (American, p = 0.1552; Asian, p = 0.1945; and Cosmopolitan, p = 0.136) (Figure 8A–C). In contrast, the variance among sire means log10 (PFU/mL)/saliva for CHIKV and ZIKV genotypes were all significantly different (CHIKV ECSA-IOL, p = 0.0154; CHIKV Asian, p = 0.0242; ZIKV Asian, p = 0.0173; and ZIKV African, p = 0.0173) (Figure 8D–F).
Narrow-sense heritabilities (h2) for SGEB were also calculated to characterize each viral infection as SGIB. The observational components σ2S, σ2d, and σ2w and the causal components of variance (Vp, Vsire, Vdam, Ve) for log10 (PFU/mL)/saliva were obtained (Table S2). In contrast to SGIB, the heritabilities for SGEB (log10 (PFU/mL)/saliva) were not significant for any of the three DENV2 genotypes tested (American, h2 = 0.319, p = 0.127; Asian, h2 = 0.160, p = 0.185; Cosmopolitan, h2 = 0.2888, p = 0.1103). Except for ZIKV African, the heritabilities for SGEB were likewise significant (for CHIKVECSA-IOL, h2 = 0.461, p = 0.00041; for CHIKV Asian, h2 = 0.437, p = 7.7932 × 10−5; for ZIKV Asian, h2 = 0.3286, p = 0.000226; and for ZIKV African, h2 = 0.4050, p = 0.0415). Clearly, genetic factors in the mosquito do influence the amount of virus in the saliva for CHIKV and ZIKV Asian (Table 3).
In examining Figure 8, note the positive trend (R = 0.86) between the variance among sire means and the h2. The average log10 (PFU/mL) virus in the saliva was not correlated with heritability (R = 0.552, p = 0.2557). Once the virus has infected the salivary gland, SGEB appears to have a major effect on the presence or absence of virus in the saliva (Table 1). SGEB has a moderate genetic basis for Zika or chikungunya infections, but not with dengue virus infections.

3.4. Relationship between SGIB and SGEB

We tested for a statistical correlation between log10 (PFU/mL)/salivary gland and log10 (PFU/mL)/saliva to determine whether the amount of virus in the salivary glands affected the amount of virus in the saliva. We calculated Pearson’s correlation coefficients between the log10 (PFU/mL)/salivary gland and the average log10 (PFU/mL)/saliva across all individual mosquitoes (Table 4) for all viruses tested. These correlations were only calculated among mosquitoes with virus in the salivary gland. Three of the seven genotypes have a significant correlation, but the magnitude of the correlations was small (DENV2 Asian, p = 0.0487; CHIKV ECSA-IOL1, p = 0.0026; ZIKV African, p < 0.0001).
An additional correlation analysis was completed between the average log10 (PFU/mL)/salivary gland and the average log10 (PFU/mL)/saliva across all sire families for a genotype. When combined across all genotypes, a correlation of R = 0.6533 (p < 0.0001) was obtained.
However, significant correlation was only obtained for DENV2 American (p = 0.0095), DENV2 Asian (p = 0.0124) and ZIKV African (p = 0.0016). This suggests that the overall correlation is an artifact of combining results from different viral genotypes. In general, the titer of virus in the salivary glands was not strongly correlated with the titer of virus in the saliva.

3.5. Relationship between Saliva Volume, Virus Titer, and Prevalence

To evaluate whether we were examining the genetics of SGEB underlying the amount of virus in the saliva or the SGEB genetics underlying the amount of saliva produced by mosquito, an additional 100 mosquitoes were intrathoracially injected with 500 PFU/mL for each virus, and at the time of saliva collection. All samples were assigned categories (1+, 2+, 3+, etc.) based on observed volume (Figure 2), using the stereoscope Olympus SZH10 at 2X magnification. Mosquitoes that did not visibly salivate or samples that were too difficult to measure were eliminated. Techniques have been designed to quantify saliva [47,48], but they prevent subsequent work with collected virus, and these techniques require excessive standardization to correlate the amount of label protein with the volume of saliva. We are aware this method did not (or could not) accurately determine the amount of saliva collected from individual mosquitoes, but it gave a basis for correlating the volume of saliva to the virus titer.
If all mosquitoes produce the same concentration of virus, then we would expect the amount of virus to be proportional to the amount of saliva. The amount of live virus would also be a function of the volume of saliva released into the host, which is also a continuous variable. Alternatively, if the saliva contains clusters of viruses, then there should be no relationship between the amount of saliva and the amount of virus. It has been reported that saliva often contains clusters of viruses rather than a uniform distribution [3,22]. There is no significant difference between the amount of saliva and the virus titer (log10 (PFU/mL)) for all virus strains tested in this study (DENV2 genotypes: p = 0.8530 for American, p = 0.7188 for Asian, and p = 0.7532 for Cosmopolitan; CHIKV genotypes: p = 0.4236 for ECSA-IOL and p = 0.5910 for Asian; ZIKV genotypes: p = 0.7377 for Asian and p = 0.0963 for African). It is very clear that the amount of virus in saliva varies greatly in samples of the same volume (Figure 9). Interestingly, the prevalence of infected saliva was similar across all volumes of saliva collected (DENV2 genotypes: df = 0.1073, p = 0.9909 for American; df = 0.2005, p = 0.0.9775 for Asian; df = 0.7809, p = 0.8540 for Cosmopolitan; CHIKV genotypes: df = 1.923, p = 0.5885 for ECSA-IOL; df = 2.670, p = 0.6145 for Asian; ZIKV genotypes: df = 3.522, p = 0.3179 for Asian and df = 1.309, p = 0.8599 for African (chi squares values were only calculated without highest volume of saliva due to small number of samples)). This is true for all DENV, CHIKV, and ZIKV genotypes tested.
An additional correlation analysis was completed between the log10 (PFU/mL)/saliva and saliva volume for all viruses (Table 5). These correlations were only calculated in saliva samples positive for virus. The only significant correlation that was obtained was for ZIKV African, but the magnitude of the correlation was small (R = 0.1506, p = 0.0018). In general, saliva titers were not correlated with saliva volumes.

4. Discussion

Herein, we report the results of a half-sib quantitative genetic analysis of SGIB and SGEBs in a freshly colonized strain of Ae. aegypti from Mapastepec, Mexico. Half-sib analysis allows for an assessment of whether there is a genetic basis for a phenotypic trait that varies in expression in nature or in the laboratory. The critical question addressed by this study is whether there is a genetic basis for SGIB and SGEB in Ae. aegypti. The answer is affirmative, but with the qualification that the heritability of these traits in the mosquito is dependent upon the viral genotype, or more specifically, on the growth rate of the virus in the salivary gland. In addition, many of the significant h2 estimates are small, suggesting that most of the phenotypic variation is still unexplained, and that many uncontrolled factors in the environment contribute to the Vp exhibited in SGIB and SGEB. For each half-sib analysis, we assumed that the variance contributed by the virus should be minimal, because high-passaged virus in cell culture were used in this study. Usually there is abundant genetic variation present in Ae. aegypti field collections [25,49], and this genetic variation may affect the growth of the virus in the mosquito’s salivary gland. Examining Figure 6A,B, it is very clear that different half-sib families vary greatly in the amount of virus in salivary glands. The high h2 suggests that this is due to the segregation of genetic factors that affect viral reproduction among and within half-sib families. In contrast with this idea, a DENV2 genotype (Cosmopolitan) and a CHIKV genotype (Asian) (Figure 6C,E) exhibited obviously less variation in the amount of virus and higher average viral titers (Figure 5A), and the h2 was close to zero. Table 2 indicates significant h2 values for DENV2 American and Asian genotypes, for CHIKV ECSA-IOL genotype, and for ZIKV Asian genotype. Thus, it would appear that an SGIB following the introduction of virus into the hemocele has a minimal effect on salivary gland infection rates. However, SGIB did affect the quantity of virus present in the salivary gland.
A low percentage of mosquitoes injected with DENV2 genotypes remained uninfected, but all mosquitoes injected with CHIKV or ZIKV genotypes were infected (Table 1). The extrinsic incubation period varies among virus genotype and mosquito strain, and it is greatly dependent on environmental factors, such as temperature and effective dose of virus [50,51,52]. However, a minimum of one week after initial infection is required for DENV to be released and detected in saliva, indicating an absence of SGEB and SGIB [51]. There is minimal information on relative growth rates of different arboviruses following intrathoracic inoculation in mosquitoes. Since we inoculated all mosquitoes with the same number of PFU, and since we collected all mosquitoes at 7 dpi, our observations suggest, for example, that DENV2 American genotype grows more slowly than the other genotypes and that ZIKV African genotype grows faster (Figure 4). To formally test this hypothesis, we processed salivary glands and saliva from DENV2 American infected mosquitoes at a later extrinsic incubation period (14 dpi). Average virus titers in salivary glands at 14 dpi (2.6 × 103 PFU/mL) were significantly higher (p = 0.0003) than at 7 dpi (1.7 × 104 PFU/mL). However, the percentage of salivary glands infected and the percentage of saliva infected was similar, at 7 dpi and 14 dpi.
One could also vary the amount of virus injected for ZIKV and CHIKV and then assess genetic contributions throughout the course of infection. We assumed that if the variances in SGIB and SGEB are due to one or all of the random causes listed in the Introduction (such as virus genotype or limit of detection of virus by plaque assay), then Ve will be large, Va will be small, and h2 will approach zero. The distribution of phenotypes among the progeny of half-sib families will appear as in Figure 6C,E. Alternatively, if the amount of virus detected in salivary glands or from saliva is due, at least in part, to genetic factors, then the variation in SGIB and SGEB will have a significant h2. The distribution of phenotypes among the progeny of half-sib families appear as in Figure 6A,B,D,F,G. Once the virus has infected the salivary gland, SGEBs appear to have a major effect on the presence or absence of virus in the saliva (Table 1, Figure 8) and the concentration of virus in the saliva (Table 3). Saliva often contains clusters of virus rather than a uniform distribution [3,22]. Indeed, this explains the poor correlation between the amount of saliva and the amount of virus obtained for all virus strains tested in this study, as has been reported in Ae. albopictus infected with different CHIKV genotypes [40]. A similar argument could be made regarding the dissection of the salivary gland. It is likely that the dissected glands will be contaminated by adjacent or proximal tissues. However, if variation in the amount of virus in the dissected glands or saliva was attributable to random causes, then none of the h2 estimates should have been significant. However, in fact, four of the seven half-sib analyses of virus titer/salivary gland had a significant h2, and in three of the seven half-sib analyses of virus titer/saliva, h2 was significant. Thus, despite the random contribution of saliva volume and the potential presence of attached tissues during salivary gland dissections, significant heritabilities were still estimated. Further, it is unlikely that we may have underestimated h2 for SGEB as a result of problems associated with determining the volume of saliva produced by individual females and other uncontrolled factors. It is very clear that the amount of virus can vary greatly in similar saliva volumes (Figure 9).
The heritabilities for DENV2 SGEB were not significant, nor were the variances among sires (Figure 8, Table 3). In contrast, the heritabilities were significant, as were the variances among sires, when infected with either of the CHIKV genotypes or with ZIKV. Clearly, genetic factors in the mosquito do influence the abundance of CHIKV and ZIKV during salivation. The observed negative correlation between viral titers at the time of harvest and h2 (Figure 7) is consistent with a threshold model wherein once a viral population has reached a critical density, then genetic factors in the mosquito that influence viral survival and titer no longer have an effect. As the viral population approaches the threshold, the amount of variation among means of sire and h2 will approach zero (Figure 7). This pattern suggests that when the amount of a virus is low, genetic factors in the mosquito may influence viral titer in the salivary glands. However, once viral titers reach a high enough level in the salivary glands, the mosquito’s genetic contribution may no longer have any effect. Indeed, genetic factors associated with the innate immune response to virus infections are in play here, and the particular combination of genetic factors in a given mosquito determines how successful the virus is, so low or high titers reflect a genetic basis in the mosquito. Negatives values of Vd were obtained from DENV2 American genotype and CHIKV ECSA-IOL genotype for SGIB (Supplementary Table S1), and Vd values were negative for all six genotypes, except the DENV2 Cosmopolitan genotype, for SGEB (Supplementary Table S2). Mathematically, this occurs because the average variation virus titer (log10 PFU/mL) among dams within sire is less than the average variation among all sires–families. The genetic explanation for this involves either (1) dominance wherein variance among heterozygotes is not intermediate in value between parental types or (2) common environmental variance arising when progeny of a single female is all raised in the same larval container. This leads to two types of environmental variance: variance arising from a common environment (all larvae in a family raised in the same container), and variance arising from half-sibs raised under different, albeit unapparent, conditions. Because Vp = Va + Vd + Ve and h2 = Va/Vp, a negative dominance reduces the denominator and increases h2. Raising half-sib larvae in a common environment could reduce Ve and inflate the covariance among siblings, and this also could inflate h2. There is abundant evidence for gene by gene (G × G) interactions between arboviral genotypes and Ae. aegypti [53]. The genetics underlying this vector–pathogen association is still incompletely understood, but variation in vector competence appears to result, in part, from the interplay between two genomes. The findings reported here are consistent with G × G interactions. Mosquito innate immunity genes [54,55,56] impacting viral growth may only be able to act when the viral titer is low. Thus, viral genes that enable rapid growth may overcome innate immunity genes in the mosquito. In our study, variation in the mosquito arises from sampling directly from a field population, where potentially many innate immunity genes are segregating. It is clear that the different genotypes of the virus in any of the three species examined here generate a great deal of phenotypic variation in the virus propagation rates and overcome infection and transmission barriers in mosquitoes. In terms of infection rates, it is clear that SGEB provides the greatest heritable barrier to transmission in mosquitoes with infected salivary glands.

5. Conclusions

Herein, we report the results of a half-sib quantitative genetic analysis of SGIB and SGEBs in a freshly colonized strain of Ae. aegypti from Mapastepec, Mexico. The critical question addressed by this study is whether there is a genetic basis for SGIB and SGEB in Ae. aegypti. The answer is affirmative, but with the qualification that the heritability of these traits in the mosquito is dependent upon the viral genotype, or more specifically, on the growth rate of the virus in the salivary gland. This pattern suggests that when the amount of a virus is low, genetic factors in the mosquito may have an influence on viral titer in the salivary glands. However, once viral titers reach a high enough level in the salivary glands, the mosquito’s genetic contribution may no longer have any effect. Once the virus reaches a certain threshold, the effects of mosquito genes in determining the amount of virus declines, as does the variance in viral titers among sires. In addition, many of the significant h2 estimates are small, suggesting that most of the phenotypic variation is still unexplained, and that many uncontrolled factors in the environment contribute to the Vp exhibited in SGIB and SGEB.

Supplementary Materials

The following are available online at https://0-www-mdpi-com.brum.beds.ac.uk/2075-4450/12/1/73/s1, Figure S1: The SAS code for performing REML analyses, Table S1: Quantitative genetic analysis of SGIB in Aedes aegypti, Table S2. Quantitative genetic analysis of SGEB in Aedes aegypti.

Author Contributions

Conceptualization: W.C.B.IV, K.E.O., and I.S.-V.; methodology: I.S.-V.; formal analysis: all authors; investigation: all authors; data curation: all authors; writing—original draft preparation: all authors; writing—review and editing: K.E.O.; supervision: W.C.B.IV and K.E.O.; funding acquisition: W.C.B.IV and K.E.O. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by NIH R21-AI128374-02 “Assessing the genetic basis of arboviral salivary gland escape barriers in Aedes aegypti” to W.C.B.IV, K.E.O., and I.S.V.

Data Availability Statement

All data produced from this study are included in this published paper.

Acknowledgments

The authors thank Adeline Williams for editorial assistance. The authors would also like to thank Gabby Ramirez and Susi Bennett for technical support in rearing mosquitoes and processing samples.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript; or in the decision to publish the results.

References

  1. Kenney, J.L.; Brault, A.C. The role of environmental, virological and vector interactions in dictating biological trans-mission of arthropod-borne viruses by mosquitoes. Adv. Virus Res. 2014, 89, 39–83. [Google Scholar] [PubMed]
  2. Franz, A.W.; Kantor, A.M.; Passarelli, A.L.; Clem, R.J. Tissue Barriers to Arbovirus Infection in Mosquitoes. Viruses 2015, 7, 3741–3767. [Google Scholar] [CrossRef] [PubMed]
  3. Vega-Rúa, A.; Schmitt, C.; Bonne, I.; Krijnse-Locker, J.; Failloux, A.-B. Chikungunya Virus Replication in Salivary Glands of the Mosquito Aedes albopictus. Viruses 2015, 7, 5902–5907. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  4. Chamberlain, R.W.; Sudia, W.D. Mechanism of Transmission of Viruses by Mosquitoes. Annu. Rev. Èntomol. 1961, 6, 371–390. [Google Scholar] [CrossRef] [PubMed]
  5. Murphy, F.A.; Whitfield, S.G.; Sudia, W.D.; Chamberlain, R.W. Invertebrate Immunity. In Interactions of Vector with Verte-Brate Pathogenic Viruses; Academic Press: Cambridge, MA, USA, 1975. [Google Scholar]
  6. Raquin, V.; Lambrechts, L. Dengue virus replicates and accumulates in Aedes aegypti salivary glands. Virology 2017, 507, 75–81. [Google Scholar] [CrossRef]
  7. Guedes, D.R.; Paiva, M.H.; Donato, M.M.; Barbosa, P.P.; Krokovsky, L.; Rocha, S.W.D.S.; La Saraiva, K.; Crespo, M.M.; Rezende, T.M.; Wallau, G.L.; et al. Zika virus replication in the mosquito Culex quinquefasciatus in Brazil. Emerg. Microbes Infect. 2017, 6, 1–11. [Google Scholar] [CrossRef] [Green Version]
  8. Takahashi, M. Differential transmission efficiency for Japanese encephalitis virus among colonized strains of Culex tri-taeniorhynchus. Med. Entomol. Zool. 1982, 33, 325–333. [Google Scholar] [CrossRef] [Green Version]
  9. Beaty, B.J.; Holterman, M.; Tabachnick, W.; Shope, R.E.; Rozhon, E.J.; Bishop, D.H. Molecular basis of bunyavirus transmission by mosquitoes: Role of the middle-sized RNA segment. Science 1981, 211, 1433–1435. [Google Scholar] [CrossRef]
  10. Paulson, S.L.; Grimstad, P.R.; Craig, G.B. Midgut and salivary gland barriers to La Crosse virus dissemination in mosquitoes of the Aedes triseriatus group. Med. Vet. Èntomol. 1989, 3, 113–123. [Google Scholar] [CrossRef]
  11. Jupp, P.G. Culex theileri and Sindbis virus; salivary glands infection in relation to transmission. J. Am. Mosq. Control. Assoc. 1985, 1, 374–376. [Google Scholar]
  12. Grimstad, P.R.; Paulson, S.L.; Craig, G.B. Vector Competence of Aedes Hendersoni (Diptera: Culicidae) for La Crosse Virus and Evidence of a Salivary-Gland Escape Barrier1. J. Med. Èntomol. 1985, 22, 447–453. [Google Scholar] [CrossRef] [PubMed]
  13. Romoser, W.S.; Turell, M.J.; Lerdthusnee, K.; Neira, M.; Dohm, D.; Ludwig, G.; Wasieloski, L. Pathogenesis of Rift Valley fever virus in mosquitoes—Tracheal conduits & the basal lamina as an extra-cellular barrier. Infect. Dis. Nat. Mech. Viral Emerg. Persistence 2005, 89–100. [Google Scholar] [CrossRef]
  14. Turell, M.J.; Britch, S.C.; Aldridge, R.L.; Kline, D.L.; Boohene, C.; Linthicum, K.J. Potential for Mosquitoes (Diptera: Cu-licidae) From Florida to Transmit Rift Valley Fever Virus. J. Med. Entomol. 2013, 50, 1111–1117. [Google Scholar] [CrossRef] [PubMed]
  15. Salazar, M.I.; Richardson, J.H.; Sánchez-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, 9. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  16. Bowers, D.F.; Abell, B.A.; Brown, D.T. Replication and tissue tropism of the alphavirus Sindbis in the mosquito Aedes al-bopictus. Virology 1995, 212, 1–12. [Google Scholar] [CrossRef] [Green Version]
  17. Raquin, V.; Wannagat, M.; Zouache, K.; Legras-Lachuer, C.; Moro, C.V.; Mavingui, P. Detection of dengue group viruses by fluorescence in situ hybridization. Parasites Vectors 2012, 5, 243. [Google Scholar] [CrossRef] [Green Version]
  18. Janzen, H.G.; Rhodes, A.J.; Doane, F.W. Chikungunya virus in salivary glands of Aedes aegypti (L.): An electron microscope study. Can. J. Microbiol. 1970, 16, 581–586. [Google Scholar] [CrossRef]
  19. Gaidamovich, S.Y.; Khutoretskaya, N.V.; Lvova, A.I.; Sveshnikova, N.A. Immunofluorescent staining study of the salivary glands of mosquitoes infected with group A arboviruses. Intervirology 1973, 1, 193–200. [Google Scholar]
  20. Tchankouo-Nguetcheu, S.; Bourguet, E.; Lenormand, P.; Rousselle, J.-C.; Namane, A.; Choumet, V. Infection by chikungunya virus modulates the expression of several proteins in Aedes aegypti salivary glands. Parasites Vectors 2012, 5, 264. [Google Scholar] [CrossRef] [Green Version]
  21. Ciano, K.A.; Saredy, J.J.; Bowers, D.F. Heparan Sulfate Proteoglycan: An Arbovirus Attachment Factor Integral to Mosquito Salivary Gland Ducts. Viruses 2014, 6, 5182–5197. [Google Scholar] [CrossRef] [Green Version]
  22. Girard, Y.A.; Schneider, B.S.; McGee, C.E.; Wen, J.; Han, V.C.; Popov, V.; Mason, P.W.; Higgs, S. Salivary gland morphol-ogy and virus transmission during long-term cytopathologic West Nile virus infection in Culex mosquitoes. Am. J. Trop. Med. Hyg. 2007, 76, 118–128. [Google Scholar] [CrossRef] [Green Version]
  23. Barreau, C.; Conrad, J.; Fischer, E.; Lujan, H.D.; Vernick, K.D. Identification of surface molecules on salivary glands of the mosquito, Aedes aegypti, by a panel of monoclonal antibodies. Insect Biochem. Mol. Biol. 1999, 29, 515–526. [Google Scholar] [CrossRef]
  24. Juhn, J.; Naeem-Ullah, U.; Guedes, B.A.M.; Majid, A.; Coleman, J.; Pimenta, P.F.P.; Akram, W.; James, A.A.; Marinotti, O. Spatial mapping of gene expression in the salivary glands of the dengue vector mosquito, Aedes aegypti. Parasites Vectors 2011, 4, 1. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  25. Dickson, L.B.; Campbell, C.L.; Juneja, P.; Jiggins, F.M.; Sylla, M.; Black, W.C., 4th. Exon-Enriched Libraries Reveal Large Genic Differences Between Aedes aegypti from Senegal, West Africa, and Populations Outside Africa. G3 (Bethesda) 2017, 7, 571–582. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  26. Dubrulle, M.; Mousson, L.; Moutailler, S.; Vazeille, M.; Failloux, A.-B. Chikungunya Virus and Aedes Mosquitoes: Saliva Is Infectious as soon as Two Days after Oral Infection. PLoS ONE 2009, 4, e5895. [Google Scholar] [CrossRef] [PubMed]
  27. McElroy, K.L.; Girard, Y.A.; McGee, C.E.; Tsetsarkin, K.A.; Vanlandingham, D.L.; Higgs, S. Characterization of the anti-gen distribution and tissue tropisms of three phenotypically distinct yellow fever virus variants in orally infected Aedes ae-gypti mosquitoes. Vector Borne Zoonotic Dis. 2008, 8, 675–687. [Google Scholar] [CrossRef]
  28. Chouin-Carneiro, T.; Vega-Rua, A.; Vazeille, M.; Yebakima, A.; Girod, R.; Goindin, D.; Dupont-Rouzeyrol, M.; Louren-ço-de-Oliveira, R.; Failloux, A.B. Differential Susceptibilities of Aedes aegypti and Aedes albopictus from the Americas to Zika Virus. PLoS Negl. Trop. Dis. 2016, 10, e0004543. [Google Scholar] [CrossRef]
  29. Richardson, J.; Molina-Cruz, A.; Salazar, M.I.; Black, W., 4th. Quantitative analysis of dengue-2 virus RNA during the extrin-sic incubation period in individual Aedes aegypti. Am. J. Trop. Med. Hyg. 2006, 74, 132–141. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  30. Bennett, K.E.; Olson, K.E.; Muñoz Mde, L.; Fernandez-Salas, I.; Farfan-Ale, J.A.; Higgs, S.; Black, W.C., 4th; Beaty, B.J. Variation in vector competence for dengue 2 virus among 24 collections of Aedes aegypti from Mexico and the United States. Am. J. Trop. Med. Hyg. 2002, 67, 85–92. [Google Scholar] [CrossRef] [Green Version]
  31. Myles, K.M.; Pierro, D.J.; Olson, K.E. Comparison of the Transmission Potential of Two Genetically Distinct Sindbis Viruses after Oral Infection ofAedes aegypti (Diptera: Culicidae). J. Med. Èntomol. 2004, 41, 95–106. [Google Scholar] [CrossRef] [Green Version]
  32. Coffey, L.L.; Failloux, A.-B.; Weaver, S.C. Chikungunya Virus–Vector Interactions. Viruses 2014, 6, 4628–4663. [Google Scholar] [CrossRef] [PubMed]
  33. Khoo, C.C.H.; 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, 257. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  34. Tabachnick, W.J. Nature, Nurture and Evolution of Intra-Species Variation in Mosquito Arbovirus Transmission Competence. Int. J. Environ. Res. Public Health 2013, 10, 249–277. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  35. Falconer, D.S. Introduction to Quantitative Genetics; Longmans Green: Harlow, Essex, UK, 1996. [Google Scholar]
  36. Chen, W.-J.; Wei, H.L.; Hsu, E.L.; Chen, E.R. Vector competence of Aedes albopictus and Ae. aegypti (Diptera: Culicidae) to dengue 1 virus on Taiwan: Development of the virus in orally and parenterally infected mosquitoes. J. Med. Èntomol. 1993, 30, 524–530. [Google Scholar] [CrossRef] [PubMed]
  37. Styer, L.M.; Kent, K.A.; Albright, R.G.; Bennett, C.J.; Kramer, L.D.; Bernard, K.A. Mosquitoes Inoculate High Doses of West Nile Virus as They Probe and Feed on Live Hosts. PLoS Pathog. 2007, 3, 1262–1270. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  38. Gubler, D.J.; Rosen, L. A Simple Technique for Demonstrating Transmission of Dengue Virus by Mosquitoes without the Use of Vertebrate Hosts. Am. J. Trop. Med. Hyg. 1976, 25, 146–150. [Google Scholar] [CrossRef]
  39. Phillips, A.; Mossel, E.; Sanchez-Vargas, I.; Foy, B.; Olson, K. Alphavirus transducing system: Tools for visualizing infection in mosquito vectors. J. Vis. Exp. 2010, 45, 2363. [Google Scholar] [CrossRef] [Green Version]
  40. Sanchez-Vargas, I.; Harrington, L.C.; Black, W.C.; Olson, K.E. Analysis of Salivary Glands and Saliva from Aedes albopic-tus and Aedes aegypti Infected with Chikungunya Viruses. Insects 2019, 10, 39. [Google Scholar] [CrossRef] [Green Version]
  41. Sladowski, D.; Steer, S.J.; Clothier, R.H.; Balls, M. An improved MTT assay. J. Immunol. Methods 1993, 157, 203–207. [Google Scholar] [CrossRef]
  42. Takeuchi, H.; Baba, M.; Shigeta, S. An application of tetrazolium (MTT) colorimetric assay for the screening of anti-herpes simplex virus compounds. J. Virol. Methods 1991, 33, 61–71. [Google Scholar] [CrossRef]
  43. Becker, W.A. Manual of Quantitative Genetics; Academic Enterprises: Pullman, WA, USA, 1992; pp. 55–66. [Google Scholar]
  44. Lynch, M.; Walsh, B. Genetics and Analysis of Quantitative Traits; Sinauer Associates: Sunderland, MA, USA, 1998. [Google Scholar]
  45. Shaw, R.G. Maximum-Likelihood Approaches Applied to Quantitative Genetics of Natural-Populations. Evolution 1987, 41, 812–826. [Google Scholar]
  46. Mathur, G.; Sanchez-Vargas, I.; Alvarez, D.; Olson, K.E.; Marinotti, O.; James, A.A. Transgene-mediated suppression of dengue viruses in the salivary glands of the yellow fever mosquito, Aedes aegypti. Insect Mol. Biol. 2010, 19, 753–763. [Google Scholar] [CrossRef] [Green Version]
  47. Devine, T.L.; Venard, C.E.; Myser, W.C. Measurement of salivation by Aedes aegypti (L.) feeding on a living host. J. Insect Physiol. 1965, 11, 347–353. [Google Scholar] [CrossRef]
  48. Hurlbut, H.S. Mosquito Salivation and Virus Transmission. Am. J. Trop. Med. Hyg. 1966, 15, 989–993. [Google Scholar] [CrossRef]
  49. Evans, B.R.; Gloria-Soria, A.; Hou, L.; McBride, C.; Bonizzoni, M.; Zhao, H.; Powell, J.R. A Multipurpose, High-Throughput Single-Nucleotide Polymorphism Chip for the Dengue and Yellow Fever Mosquito, Aedes aegypti. G3-Genes Genomes Genet. 2015, 5, 711–718. [Google Scholar] [CrossRef] [Green Version]
  50. Black, W.C., 4th; Bennett, K.E.; Gorrochótegui-Escalante, N.; Barillas-Mury, C.V.; Fernández-Salas, I.; de Lourdes Muñoz, M.; Farfán-Alé, J.A.; Olson, K.E.; Beaty, B.J. Flavivirus susceptibility in Aedes aegypti. Arch. Med. Res. 2002, 33, 379–388. [Google Scholar] [CrossRef] [Green Version]
  51. Chan, M.; Johansson, M.A. The Incubation Periods of Dengue Viruses. PLoS ONE 2012, 7, e50972. [Google Scholar] [CrossRef] [PubMed]
  52. Winokur, O.C.; Main, B.J.; Nicholson, J.; Barker, C.M. Impact of temperature on the extrinsic incubation period of Zika virus in Aedes aegypti. PLoS Negl. Trop. Dis. 2020, 14, e0008047. [Google Scholar] [CrossRef] [Green Version]
  53. Lambrechts, L. Quantitative genetics of Aedes aegypti vector competence for dengue viruses: Towards a new paradigm? Trends Parasitol. 2011, 27, 111–114. [Google Scholar] [CrossRef] [Green Version]
  54. Kumar, A.; Sunil, S.; Sirisena, P.; Dubey, S.K.; Kumar, R.; Shrinet, J.; Sunil, S. Mosquito Innate Immunity. Insects 2018, 9, 95. [Google Scholar] [CrossRef] [Green Version]
  55. Prasad, A.N.; Brackney, D.E.; Ebel, G.D. The Role of Innate Immunity in Conditioning Mosquito Susceptibility to West Nile Virus. Viruses 2013, 5, 3142–3170. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  56. Pan, X.L.; Pike, A.; Joshi, D.; Bian, G.W.; McFadden, M.J.; Lu, P.; Liang, X.; Zhang, F.R.; Raikhel, A.S.; Xi, Z.Y. The bacte-rium Wolbachia exploits host innate immunity to establish a symbiotic relationship with the dengue vector mosquito Ae-des aegypti. ISME J. 2018, 12, 277–288. [Google Scholar] [CrossRef] [PubMed] [Green Version]
Figure 1. Structure of Aedes spp. salivary glands [3]. (A) Schematic representation of salivary glands of a mosquito at 2 days post-emergence. (B) Transverse section through a median lobe of the salivary gland showing the acinar cells and their relationship to the salivary duct. (C) An enlargement of B showing a cartoon of a transverse-section of salivary gland indicating mature virus penetrating the basal lamina to establish an infection in the acinar cells (A,C) and then penetrating the plasma membrane to fill the apical cavity and/or accumulate in the acinar cell junctions. Pink arrows trace arbovirus entry and movement in individual acinar cells leading to the virus’ eventual release through the apical cavity (yellow lines). Salivary gland escape (SGE)1 and SGE2 refer to potential barriers (entry into the apical cavity and acinar cell junctions, respectively) for virus escape from salivary glands (D) Immunofluorescence of salivary gland lobe infected with DENV-2 Merida (Cosmopolitan) at 7 dpi shows that mature viral particles clearly accumulated in acinar cell junctions (red arrow). (E) Uninfected salivary glands showed no immunofluorescence signal. N: nuclei. SD: salivary duct. AC: apical cavity for saliva storage. BL: basal lamina. Cyt: Cytoplasm. CM: Cell membrane. DLL: distal lateral lobe, ML: medial lobe. * Represents a virus cluster (1C).
Figure 1. Structure of Aedes spp. salivary glands [3]. (A) Schematic representation of salivary glands of a mosquito at 2 days post-emergence. (B) Transverse section through a median lobe of the salivary gland showing the acinar cells and their relationship to the salivary duct. (C) An enlargement of B showing a cartoon of a transverse-section of salivary gland indicating mature virus penetrating the basal lamina to establish an infection in the acinar cells (A,C) and then penetrating the plasma membrane to fill the apical cavity and/or accumulate in the acinar cell junctions. Pink arrows trace arbovirus entry and movement in individual acinar cells leading to the virus’ eventual release through the apical cavity (yellow lines). Salivary gland escape (SGE)1 and SGE2 refer to potential barriers (entry into the apical cavity and acinar cell junctions, respectively) for virus escape from salivary glands (D) Immunofluorescence of salivary gland lobe infected with DENV-2 Merida (Cosmopolitan) at 7 dpi shows that mature viral particles clearly accumulated in acinar cell junctions (red arrow). (E) Uninfected salivary glands showed no immunofluorescence signal. N: nuclei. SD: salivary duct. AC: apical cavity for saliva storage. BL: basal lamina. Cyt: Cytoplasm. CM: Cell membrane. DLL: distal lateral lobe, ML: medial lobe. * Represents a virus cluster (1C).
Insects 12 00073 g001
Figure 2. Individual capillaries observed under the microscope for the presence of saliva and assigned saliva categories based on total volume. The volume of saliva was calculated using the cylinder volume formula (V = π (r2 × h) = 3.1416 × 0.01 × h = mm3 or µL) and measuring the height of the saliva with a digital fractional caliper (precision: ±0.02 mm), as described previously.
Figure 2. Individual capillaries observed under the microscope for the presence of saliva and assigned saliva categories based on total volume. The volume of saliva was calculated using the cylinder volume formula (V = π (r2 × h) = 3.1416 × 0.01 × h = mm3 or µL) and measuring the height of the saliva with a digital fractional caliper (precision: ±0.02 mm), as described previously.
Insects 12 00073 g002
Figure 3. Method for generating nested half-sib design used for calculating estimates of genetic and environmental variances and covariances.
Figure 3. Method for generating nested half-sib design used for calculating estimates of genetic and environmental variances and covariances.
Insects 12 00073 g003
Figure 4. Titer of Mapastepec Aedes aegypti mosquitoes 7 days after intrathoracically injecting with 500 plaque forming units (PFU) of virus. Titers were determined by plaque assay in the whole mosquito for each viral strain. Means and 95% confidence intervals are shown. * = p < 0.0001.
Figure 4. Titer of Mapastepec Aedes aegypti mosquitoes 7 days after intrathoracically injecting with 500 plaque forming units (PFU) of virus. Titers were determined by plaque assay in the whole mosquito for each viral strain. Means and 95% confidence intervals are shown. * = p < 0.0001.
Insects 12 00073 g004
Figure 5. (A) Average numbers of log10 (PFU/mL) per salivary gland for all seven viral genotypes. (B) Average numbers of log10 (PFU/mL) per saliva were calculated for all seven viral genotypes. Horizontal lines over pairs of genotypes in either graph indicate that average log10 (PFU/mL) per salivary gland did not differ significantly.
Figure 5. (A) Average numbers of log10 (PFU/mL) per salivary gland for all seven viral genotypes. (B) Average numbers of log10 (PFU/mL) per saliva were calculated for all seven viral genotypes. Horizontal lines over pairs of genotypes in either graph indicate that average log10 (PFU/mL) per salivary gland did not differ significantly.
Insects 12 00073 g005
Figure 6. Distribution of phenotypes log10 (PFU/mL) per salivary gland in half-sib families for the seven virus strains. Each column of points arising from the x-axis represents a sire family. In each column, the distribution of individual salivary glands titers arising from that sire appear as gray circles. Mean sire values and their 95% confidence intervals are displayed in red. The mean of dam families in each sire appear in green. (A) DENV2 American genotype, (B) DENV2 Asian genotype, (C) DENV2 Cosmopolitan genotype, (D) CHIKV ESCA-IOL genotype, (E) CHIKV Asian genotype, (F) ZIKV Asian genotype, and (G) ZIKV African genotype.
Figure 6. Distribution of phenotypes log10 (PFU/mL) per salivary gland in half-sib families for the seven virus strains. Each column of points arising from the x-axis represents a sire family. In each column, the distribution of individual salivary glands titers arising from that sire appear as gray circles. Mean sire values and their 95% confidence intervals are displayed in red. The mean of dam families in each sire appear in green. (A) DENV2 American genotype, (B) DENV2 Asian genotype, (C) DENV2 Cosmopolitan genotype, (D) CHIKV ESCA-IOL genotype, (E) CHIKV Asian genotype, (F) ZIKV Asian genotype, and (G) ZIKV African genotype.
Insects 12 00073 g006
Figure 7. Correlation analysis to test whether the average log10 (PFU/mL) virus in the salivary glands are negatively correlated with the magnitude of heritability.
Figure 7. Correlation analysis to test whether the average log10 (PFU/mL) virus in the salivary glands are negatively correlated with the magnitude of heritability.
Insects 12 00073 g007
Figure 8. Distribution of phenotypes log10 (PFU/mL) per saliva in half-sib families for the seven virus strains. Each column of points arising from the X-axis represents a sire family. In each column, the distribution of individual saliva titer arising from that sire appear as gray circles. Mean sire values and their 95% confidence intervals are displayed in red. The mean of dam families in each sire appear in green. (A) DENV2 American genotype, (B) DENV2 Asian genotype, (C) DENV2 Cosmopolitan genotype, (D) CHIKV ESCA-IOL genotype, (E) CHIKV Asian genotype, (F) ZIKV Asian genotype, and (G) ZIKV African genotype.
Figure 8. Distribution of phenotypes log10 (PFU/mL) per saliva in half-sib families for the seven virus strains. Each column of points arising from the X-axis represents a sire family. In each column, the distribution of individual saliva titer arising from that sire appear as gray circles. Mean sire values and their 95% confidence intervals are displayed in red. The mean of dam families in each sire appear in green. (A) DENV2 American genotype, (B) DENV2 Asian genotype, (C) DENV2 Cosmopolitan genotype, (D) CHIKV ESCA-IOL genotype, (E) CHIKV Asian genotype, (F) ZIKV Asian genotype, and (G) ZIKV African genotype.
Insects 12 00073 g008
Figure 9. Correlation between virus titer log10 (PFU/mL) per saliva and volume of saliva for the seven virus strains. Each column of points arising from the X-axis represents a volume of saliva assigned. In each column, the distribution of individual saliva titers appears as black circles, also the number of positive saliva/total number of saliva per volume and prevalence is in parenthesis. Mean titer values and their 95% confidence intervals are displayed in red. Red arrow indicates non-salivated samples (n = 5). Total number of saliva samples tested and p value per virus strain are indicated. (A) DENV2 American genotype, (B) DENV2 Asian genotype, (C) DENV2 Cosmopolitan genotype, (D) CHIKV ESCA-IOL genotype, (E) CHIKV Asian genotype, (F) ZIKV Asian genotype, and (G) ZIKV African genotype.
Figure 9. Correlation between virus titer log10 (PFU/mL) per saliva and volume of saliva for the seven virus strains. Each column of points arising from the X-axis represents a volume of saliva assigned. In each column, the distribution of individual saliva titers appears as black circles, also the number of positive saliva/total number of saliva per volume and prevalence is in parenthesis. Mean titer values and their 95% confidence intervals are displayed in red. Red arrow indicates non-salivated samples (n = 5). Total number of saliva samples tested and p value per virus strain are indicated. (A) DENV2 American genotype, (B) DENV2 Asian genotype, (C) DENV2 Cosmopolitan genotype, (D) CHIKV ESCA-IOL genotype, (E) CHIKV Asian genotype, (F) ZIKV Asian genotype, and (G) ZIKV African genotype.
Insects 12 00073 g009
Table 1. Infection rates for salivary gland infection barrier (SGIB) and salivary gland escape barrier (SGEB) for each of the seven viral strains.
Table 1. Infection rates for salivary gland infection barrier (SGIB) and salivary gland escape barrier (SGEB) for each of the seven viral strains.
Viral SpeciesDENV2
Qro94
DENV2
Guerrero
DENV2
Merida
CHIKV LR 2006-OPY1CHIKV R99659ZIKV PRVABC-59ZIKV
Dakar 41525
GenotypeAmericanAsianCosmopolitanECSA-IOLAsianAsianAfrican
Number of sire families24212122222121
Number paired samples691667681686735783623
Uninfected
Glands
38540000
SGIB%5.50%0.75%0.59%0.00%0.00%0.00%0.00%
Uninfected
Saliva
578480404327452408200
Infected
Saliva
113187277358282375423
SGEB%88.5%72.5%59.7%47.7%61.5%52.1%32.1%
Transmission11.5%27.5%40.3%52.3%38.5%47.9%67.9%
Table 2. Quantitative genetic analysis of salivary gland infection barrier in Aedes aegypti. Each analysis corresponds to one of the seven virus strains tested for log10(PFU)/salivary gland. We used PROC MIXED in SAS to obtain restricted maximum likelihood (REML) to estimate quantitative genetic parameters. The greater the variance in sire means, the larger the heritability.
Table 2. Quantitative genetic analysis of salivary gland infection barrier in Aedes aegypti. Each analysis corresponds to one of the seven virus strains tested for log10(PFU)/salivary gland. We used PROC MIXED in SAS to obtain restricted maximum likelihood (REML) to estimate quantitative genetic parameters. The greater the variance in sire means, the larger the heritability.
SGIB
VirusStrainGenotypeh2σ2s
DENV-2Qro94American0.8103 ***0.2700 **
Guerrero C932Asian0.5799 **0.1029 *
Merida BC17Cosmopolitan0.06850.0082
CHIKVLR2006-OPY-1ECSA-IOL0.5174 ***0.0337 **
R99659Asian0.03850.0012
ZIKVPRVABC59Asian0.1526 *0.0077
Dakar 42525African0.9808 ***0.0762 **
σ2s = average variance among offspring of a single sire; h2 = narrow-sense heritability (h2 = Va/Vp). * p = 0.03–0.05, ** p = 0.001, *** p < 0.00001.
Table 3. Quantitative genetic analysis of SGEB in Aedes aegypti. Each analysis corresponds to one of the seven virus strains tested for log10(PFU)/saliva. We used PROC MIXED in SAS to obtain restricted maximum likelihood (REML) to estimate quantitative genetic parameters. The greater the variance in sire means, the larger the heritability.
Table 3. Quantitative genetic analysis of SGEB in Aedes aegypti. Each analysis corresponds to one of the seven virus strains tested for log10(PFU)/saliva. We used PROC MIXED in SAS to obtain restricted maximum likelihood (REML) to estimate quantitative genetic parameters. The greater the variance in sire means, the larger the heritability.
SGEB
VirusStrainGenotypeh2σ2s
DENV-2Qro94American0.31860.0287
Guerrero C932Asian0.160.0110
Merida BC17Cosmopolitan0.28880.0258
CHIKVLR2006-OPY-1ECSA-IOL0.4611 ***0.0914 *
R99659Asian0.4372 ***0.1025 *
ZIKVPRVABC59Asian0.3286 ***0.0839 *
Dakar 42525African0.40500.1569
σ2s = average variance among offspring of a single sire; h2 = narrow-sense heritability (h2 = Va/Vp). * p = 0.01, *** p < 0.0001.
Table 4. Correlation analysis between the log10 (PFU/mL)/salivary gland and the log10 (PFU/mL)/saliva across all individual mosquitoes for each of the seven genotypes. Mosquitoes without an infected salivary gland were not included. p values that are bolded are significant at a threshold of p < 0.05.
Table 4. Correlation analysis between the log10 (PFU/mL)/salivary gland and the log10 (PFU/mL)/saliva across all individual mosquitoes for each of the seven genotypes. Mosquitoes without an infected salivary gland were not included. p values that are bolded are significant at a threshold of p < 0.05.
Virus StrainGenotypeCorr. Coeff. (R)Prob.
(p)
DENV2 Qro94American−0.05030.6003
DENV2 GuerreroC932Asian0.14440.0487
DENV2 Merida BC17Cosmopolitan0.02550.6726
CHIKV R99659Asian0.06490.2774
CHIKV LR2006OPY-1ECSA-IOL0.15910.0026
ZIKV PRVABC59Asian−0.04530.2118
ZIKV DakarAfrican0.1947<0.0001
Table 5. Correlation analysis between the log10 (PFU/mL)/saliva and saliva volume for each virus. p values that are bolded are significant at a threshold of p < 0.05.
Table 5. Correlation analysis between the log10 (PFU/mL)/saliva and saliva volume for each virus. p values that are bolded are significant at a threshold of p < 0.05.
Virus StrainGenotypeCorr. Coeff.
(R)
Prob.
(p)
DENV2 Qro94American0.1258 0.5001
DENV2 GuerreroC932Asian0.15280.3596
DENV2 Merida BC17Cosmopolitan0.15700.3466
CHIKV R99659Asian0.18560.1876
CHIKV LR2006OPY-1ECSA-IOL0.09240.4981
ZIKV PRVABC59Asian−0.03840.7728
ZIKV DakarAfrican0.15060.0018
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Sanchez-Vargas, I.; Olson, K.E.; Black, W.C. The Genetic Basis for Salivary Gland Barriers to Arboviral Transmission. Insects 2021, 12, 73. https://0-doi-org.brum.beds.ac.uk/10.3390/insects12010073

AMA Style

Sanchez-Vargas I, Olson KE, Black WC. The Genetic Basis for Salivary Gland Barriers to Arboviral Transmission. Insects. 2021; 12(1):73. https://0-doi-org.brum.beds.ac.uk/10.3390/insects12010073

Chicago/Turabian Style

Sanchez-Vargas, Irma, Ken E. Olson, and William C. Black. 2021. "The Genetic Basis for Salivary Gland Barriers to Arboviral Transmission" Insects 12, no. 1: 73. https://0-doi-org.brum.beds.ac.uk/10.3390/insects12010073

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop