Skip to main content
Advertisement
  • Loading metrics

Exome-wide association study reveals largely distinct gene sets underlying specific resistance to dengue virus types 1 and 3 in Aedes aegypti

  • Laura B. Dickson ,

    Contributed equally to this work with: Laura B. Dickson, Sarah H. Merkling

    Roles Conceptualization, Data curation, Formal analysis, Investigation, Visualization, Writing – original draft, Writing – review & editing

    Affiliation Insect-Virus Interactions Unit, Institut Pasteur, UMR2000, CNRS, Paris, France

  • Sarah H. Merkling ,

    Contributed equally to this work with: Laura B. Dickson, Sarah H. Merkling

    Roles Conceptualization, Data curation, Investigation, Visualization, Writing – review & editing

    Affiliation Insect-Virus Interactions Unit, Institut Pasteur, UMR2000, CNRS, Paris, France

  • Mathieu Gautier,

    Roles Data curation, Formal analysis, Methodology, Writing – review & editing

    Affiliation CBGP, INRAE, CIRAD, IRD, Montpellier SupAgro, Univ. Montpellier, Montpellier, France

  • Amine Ghozlane,

    Roles Data curation, Formal analysis, Visualization, Writing – review & editing

    Affiliation Hub de Bioinformatique et Biostatistique–Département Biologie Computationnelle, Institut Pasteur, USR 3756 CNRS, Paris, France

  • Davy Jiolle,

    Roles Investigation, Writing – review & editing

    Affiliations Insect-Virus Interactions Unit, Institut Pasteur, UMR2000, CNRS, Paris, France, MIVEGEC, Univ. Montpellier, IRD, CNRS, Montpellier, France, Centre Interdisciplinaire de Recherches Médicales de Franceville, Franceville, Gabon

  • Christophe Paupy,

    Roles Investigation, Writing – review & editing

    Affiliations MIVEGEC, Univ. Montpellier, IRD, CNRS, Montpellier, France, Centre Interdisciplinaire de Recherches Médicales de Franceville, Franceville, Gabon

  • Diego Ayala,

    Roles Investigation, Writing – review & editing

    Affiliations MIVEGEC, Univ. Montpellier, IRD, CNRS, Montpellier, France, Centre Interdisciplinaire de Recherches Médicales de Franceville, Franceville, Gabon

  • Isabelle Moltini-Conclois,

    Roles Investigation, Writing – review & editing

    Affiliations Insect-Virus Interactions Unit, Institut Pasteur, UMR2000, CNRS, Paris, France, MIVEGEC, Univ. Montpellier, IRD, CNRS, Montpellier, France

  • Albin Fontaine,

    Roles Conceptualization, Investigation, Writing – review & editing

    Current address: Unité de Parasitologie et Entomologie, Département des Maladies Infectieuses, Institut de Recherche Biomédicale des Armées (IRBA), Marseille, France; Aix Marseille Université, IRD, AP-HM, SSA, UMR Vecteurs-Infections Tropicales et Méditerranéennes (VITROME), IHU-Méditerranée Infection, Marseille, France

    Affiliation Insect-Virus Interactions Unit, Institut Pasteur, UMR2000, CNRS, Paris, France

  • Louis Lambrechts

    Roles Conceptualization, Formal analysis, Funding acquisition, Methodology, Supervision, Visualization, Writing – original draft, Writing – review & editing

    louis.lambrechts@pasteur.fr

    Affiliation Insect-Virus Interactions Unit, Institut Pasteur, UMR2000, CNRS, Paris, France

Abstract

Although specific interactions between host and pathogen genotypes have been well documented in invertebrates, the identification of host genes involved in discriminating pathogen genotypes remains a challenge. In the mosquito Aedes aegypti, the main dengue virus (DENV) vector worldwide, statistical associations between host genetic markers and DENV types or strains were previously detected, but the host genes underlying this genetic specificity have not been identified. In particular, it is unknown whether DENV type- or strain-specific resistance relies on allelic variants of the same genes or on distinct gene sets. Here, we investigated the genetic architecture of DENV resistance in a population of Ae. aegypti from Bakoumba, Gabon, which displays a stronger resistance phenotype to DENV type 1 (DENV-1) than to DENV type 3 (DENV-3) infection. Following experimental exposure to either DENV-1 or DENV-3, we sequenced the exomes of large phenotypic pools of mosquitoes that are either resistant or susceptible to each DENV type. Using variation in single-nucleotide polymorphism (SNP) frequencies among the pools, we computed empirical p values based on average gene scores adjusted for the differences in SNP counts, to identify genes associated with infection in a DENV type-specific manner. Among the top 5% most significant genes, 263 genes were significantly associated with resistance to both DENV-1 and DENV-3, 287 genes were only associated with DENV-1 resistance and 290 were only associated with DENV-3 resistance. The shared significant genes were enriched in genes with ATP binding activity and sulfur compound transmembrane transporter activity, whereas the genes uniquely associated with DENV-3 resistance were enriched in genes with zinc ion binding activity. Together, these results indicate that specific resistance to different DENV types relies on largely non-overlapping sets of genes in this Ae. aegypti population and pave the way for further mechanistic studies.

Author summary

Compatibility between hosts and pathogens is often genetically specific in invertebrates but host genes underlying this genetic specificity have not been elucidated. We investigated the genetic architecture of dengue virus type-specific resistance in the mosquito vector Aedes aegypti. We used a natural population of Ae. aegypti from Bakoumba, Gabon, which is differentially resistant to dengue virus type 1 and dengue virus type 3. We surveyed genetic variation in protein-coding regions of the mosquito genome and compared the frequency of genetic polymorphisms between groups of mosquitoes that are either resistant or susceptible to each dengue virus type. We found that the Ae. aegypti genes associated with resistance to dengue virus type 1 or dengue virus type 3 were largely non-overlapping. This finding indicates that different sets of host genes, rather than different variants of the same genes, confer pathogen-specific resistance in this population. This study is an important step towards identification of mechanisms underlying the genetic specificity of invertebrate host-pathogen interactions.

Introduction

In many invertebrate host-pathogen systems, infection success depends on the specific pairing of host and pathogen genotypes [1]. Such genotype-by-genotype (G x G) interactions have been observed, for example, between crustaceans and bacteria [2], bumblebees and intestinal trypanosomes [3], nematodes and bacteria [4], anopheline mosquitoes and malaria parasites [5, 6] and butterflies and protozoan parasites [7]. In some instances, G x G interactions can result in extreme levels of host-pathogen specificity [8].

Understanding this genetic specificity in invertebrate host-pathogen systems has generated great enthusiasm because it may uncover new facets of invertebrate immunity [9, 10]. Although invertebrates lack the adaptive immunity of vertebrates, their immune system can generate a considerable diversity of immune receptors, revealing an unsuspected molecular complexity [11, 12]. Alternatively, G x G interactions could be mediated by variation in the host microbiota [13] or any other interplay between host and pathogen genomes. G x G interactions can be detected at the level of gene expression [14] and mapped to physical locations in the host genome [1517]. Accumulating observations on G x G interactions have led to some controversy [1820] because they challenge the prevailing view that invertebrate defense against pathogens relies on broad-spectrum recognition and effector mechanisms. A central aspect of the controversy is that the molecular mechanisms underlying G x G interactions between invertebrate hosts and pathogens are yet to be elucidated.

We previously documented significant G x G interactions between dengue virus (DENV) and its main mosquito vector Aedes aegypti [21, 22]. Both DENV and Ae. aegypti are genetically diverse in nature. DENV exists worldwide as four genetic types (DENV-1, DENV-2, DENV-3, and DENV-4) that loosely cluster antigenically [23] and are often referred to as serotypes. The mosquito Ae. aegypti consists of two major genetic subspecies that contain substantial genetic diversity [24]. In the last few decades, DENV has become a major public health threat worldwide with an estimated 390 million human infections per year [25]. The lack of an efficient vaccine and the failure of insecticide-based methods to control mosquito populations on the long term have stimulated basic research to develop novel vector control strategies [26]. One of these strategies aims at rendering mosquitoes resistant to pathogen infection [27]. The premise of this strategy makes it critical to understand G x G interactions because engineered mosquito resistance should be effective against all possible pathogen genotypes.

Dissecting the genetic basis of G x G interactions between DENV and Ae. aegypti was previously hampered by the large size of the Ae. aegypti genome and the paucity of genetic markers available [28]. Although G x G interactions between DENV and Ae. aegypti were statistically assigned to individual mosquito genetic markers [15, 17], the specific mosquito genes underlying these quantitative trait loci (QTL) have not been identified. In particular, it is unknown whether DENV type- or strain-specific resistance relies on allelic variants of the same genes or on distinct gene sets.

Here, we combined a natural phenotype of DENV type-specific resistance in Ae. aegypti and the power of high-throughput sequencing to provide insights into the genetic architecture of G x G interactions in this system. We previously established a laboratory colony from a wild Ae. aegypti population in Bakoumba, Gabon that is differentially susceptible to DENV types. Using dose-response experiments, we found the Bakoumba population to be strongly resistant to DENV-1 and only moderately resistant to DENV-3 infection. We took advantage of this field-derived colony to investigate the genetic basis for discriminating between different DENV types. We used a modified genome-wide association study design, in which replicate pools of individuals with contrasted phenotypes are pooled by phenotype and analyzed for differences in allele frequencies [29, 30]. We compared the genetic basis of Bakoumba resistance to DENV-1 and DENV-3 from extremes of the phenotypic distribution (i.e., resistant = uninfected at a high virus dose vs. susceptible = infected at a low virus dose). We estimated and compared allele frequencies among the phenotypic pools using the contrast statistic recently implemented in the software BayPass [31, 32]. Because the genome of Ae. aegypti is large and repetitive, we performed exome sequencing as a cost-effective genotyping method targeting the exons of all protein-coding genes [33], as was previously implemented in Ae. aegypti [3436]. We combined SNP-specific contrasts of allele frequency to compute gene-wide scores reflecting the statistical significance of association with DENV-1 and/or DENV-3 infection. The results of our exome-wide association study (EWAS) show that largely distinct gene sets underlie specific resistance to DENV-1 and DENV-3 infection in this Ae. aegypti population.

Results

Dose-response experiments

We initially compared the infection dose response of two outbred Ae. aegypti colonies from Cairns (Australia) and Bakoumba (Gabon) following oral exposure to either DENV-1 or DENV-3. In two separate experiments, female mosquitoes from both colonies were offered artificial blood meals containing one of three increasing infectious doses ranging from 1.0 x 104 to 3.3 x 107 FFUs/mL for DENV-1 and from 1.0 x 104 to 3.9 x 107 FFUs/mL for DENV-3 (Fig 1A). Ten days after virus exposure, infection status was determined by RT-PCR for 24–32 (mean 30.8; median 32) blood-fed females for each combination of experiment, population, virus and dose. In total, 366 individual females were tested for DENV-1 and 374 individual females were tested for DENV-3. As expected, infectious dose was a significant predictor of infection prevalence (Table 1). In addition, we observed a strong pattern of G x G interactions, whereby the Bakoumba population was significantly more resistant to DENV-1 than the three other population-virus combinations (Fig 1A). This G x G interaction was statistically supported by a significant effect (p = 0.0013) of the population x virus interaction on the infection phenotype (Table 1). Based on the regression curves of the combined dose responses, we estimated that the 50% oral infectious doses (OID50) of the Bakoumba population were 7.08 log10 FFUs/mL (95% confidence interval: 6.84–7.44) for DENV-1 and 5.54 log10 FFUs/mL (95% CI: 5.23–5.80) for DENV-3. The OID50 estimates of the Cairns population were 5.80 log10 FFUs/mL (95% CI: 5.58–6.00) for DENV-1 and 5.01 log10 FFUs/mL (95% CI: 4.68–5.26) for DENV-3. The dose-response analysis showed that the Bakoumba population was intrinsically more resistant to DENV-1 than to DENV-3 infection, although the DENV-1 and DENV-3 isolates were equally infectious to the mosquito population from Cairns. In the rest of the study, we took advantage of this natural resistance phenotype to dissect the genetic architecture of DENV type-specific resistance in Ae. aegypti. A direct comparison between the Bakoumba and the Cairns populations would be impractical because their genetic differentiation would likely mask any loci underlying DENV resistance that are shared or divergent between them. Indeed, the Bakoumba population presumably belongs to the African subspecies, Ae. aegypti formosus, whereas the Cairns population presumably belongs to the globally invasive subspecies, Ae. aegypti aegypti, therefore the two populations are expected to be genetically divergent [24]. In contrast, the comparison between phenotypically resistant or susceptible individuals within the Bakoumba population allowed us to control for the genetic background.

thumbnail
Fig 1. Investigating the DENV type-specific resistance phenotype of Bakoumba mosquitoes.

(A) Dose-response curves of Ae. aegypti colonies from Bakoumba (Gabon) or Cairns (Australia) challenged with DENV-1 and DENV-3. The percentage of DENV-infected mosquitoes at 10 days post exposure is shown as a function of the blood meal titer in log10-transformed focus-forming units (FFUs)/mL. The data was obtained in two separate dose-response experiments with three doses for each virus. The three doses of experiment 1 are covering a lower range than the three doses of experiment 2. Curves are logistic regressions of the data with their 95% confidence intervals indicated by shaded bands. Note that one data point (Bakoumba DENV-3 experiment 2) is masked due to the overlap with another value at 100% (Cairns DENV-3 experiment 2). (B) Experimental design for the exome-wide association analysis (EWAS) of DENV type-specific resistance of Bakoumba mosquitoes. Female mosquitoes of the Bakoumba population were exposed to low and high infectious doses of DENV-1 and DENV-3, expected to result in about one third and two thirds of infected mosquitoes, respectively. Bodies were used to screen for infection status by RT-PCR, and heads were collected for DNA extraction and exome sequencing. The EWAS analysis compared pools of individuals that were uninfected at the high dose (i.e., resistant) versus individuals infected at the low dose (i.e., susceptible) for each DENV type.

https://doi.org/10.1371/journal.pgen.1008794.g001

thumbnail
Table 1. Test statistics of the DENV dose-response analysis.

Infection prevalence was analyzed by logistic regression as a function of replicate experiment, dose, population, virus, and all their interactions. The table shows the minimal adequate model after sequentially removing non-significant effects. Oral infectious dose (blood meal titer) was log10-transformed prior to the analysis.

https://doi.org/10.1371/journal.pgen.1008794.t001

Exome-wide association study

We reasoned that if Bakoumba mosquitoes were intrinsically more resistant to DENV-1 than to DENV-3 infection, the genetic factors underlying resistance to each DENV type would be distinct. Our study was designed to determine whether DENV type-specific resistance in the Bakoumba population relied on different alleles of the same genes (i.e., if a DENV-1 resistance allele also confers susceptibility to DENV-3) or on genetic variation at distinct gene sets (i.e., if DENV-1 resistance and susceptibility alleles do not influence DENV-3 infection). To identify allelic variants associated with both DENV-1 and DENV-3 resistance or susceptibility, we sequenced the exome of four phenotypic pools of Bakoumba mosquitoes (Fig 1B): DENV-1 resistant (i.e., uninfected after exposure to a high DENV-1 dose), DENV-1 susceptible (i.e., infected after exposure to a low DENV-1 dose), DENV-3 resistant (i.e., uninfected after exposure to a high DENV-3 dose), and DENV-3 susceptible (i.e., infected after exposure to a low DENV-3 dose). Low and high infectious doses were adjusted to result in approximately one third and two thirds of infected mosquitoes, respectively. Due to the difference in the dose responses (Fig 1A), low and high infectious doses were different for each DENV type (Table 2). DNA was extracted from the head tissues of each individual mosquito whereas the infection phenotype was determined by RT-PCR testing of the head-less mosquito body 10 days post infectious blood meal. The phenotypic screen was performed in three replicate experiments.

thumbnail
Table 2. Infectious doses used for the EWAS phenotypic screen.

The table shows the measured titer of DENV-1 and DENV-3 infectious blood meals for each of the three experimental replicates of the exome-wide association study in the Bakoumba population.

https://doi.org/10.1371/journal.pgen.1008794.t002

Association analyses at the single-nucleotide level

The bodies of 668, 690 and 680 females was individually tested for DENV infection by RT-PCR in replicate experiments 1, 2 and 3, respectively (Table 3). Based on their resistant or susceptible phenotype, 182, 174 and 176 females were selected for DNA extraction from experiments 1, 2 and 3, respectively. DNA extracted from the head of individual females was combined into 12 phenotypic pools (3 experiments x 2 DENV types x 2 phenotypes) of 30–48 individuals (mean 44.3) to prepare 12 libraries for exome sequencing (Table 3). Individual DNA concentrations were adjusted prior to pooling so that each individual contributed the same amount of DNA to the library. Across all three chromosomes, a total of approximately 1.4 million SNPs were identified with an average sequencing depth ranging from 100.6X to 176.4X (mean 132.7X) among the 12 libraries. We first performed association analyses at the level of individual SNPs based on allele frequency differences between DENV-1 resistant and susceptible pools (DENV-1 comparison) and between DENV-3 resistant and susceptible pools (DENV-3 comparison). For both the DENV-1 and the DENV-3 comparisons, the largest number of significant SNPs was detected on chromosome 1 but each DENV type displayed a distinct profile (S1A Fig and S1B Fig). Among the top 0.001% most significant SNPs for each comparison, 5 SNPs were shared between the DENV-1 and DENV-3 comparisons, 142 were uniquely associated with DENV-1 infection status, and 146 were uniquely associated with DENV-3 infection status (S1 Table; S1C Fig). For the 5 shared SNPs, allelic effects were in the same direction (i.e., the same allele was associated with both DENV-1 and DENV-3 resistance).

thumbnail
Table 3. Summary of mosquitoes analyzed for each experimental replicate of the EWAS phenotypic screen.

The table shows the number of Bakoumba mosquitoes that tested DENV-positive or DENV-negative by RT-PCR (S2 Fig) in each of the experimental conditions of the screen. Cells highlighted in yellow are experimental conditions used as opposite phenotypes for the exome-wide association study. The number shown in brackets represents the number of mosquitoes selected for DNA extraction and subsequent sequencing.

https://doi.org/10.1371/journal.pgen.1008794.t003

Association analyses at the gene level

To identify candidate genes associated with DENV-1 and/or DENV-3 infection status, we performed a second series of association analyses based on an average score calculated for each annotated gene. An empirical p value was computed from the average gene scores adjusted for differences in SNP count per gene. Note that although the exome capture probes were designed using the AaegL3 genome build [28] because it was the only one available at the time, our gene-based analysis was performed using the more recent AaegL5 genome build [17], which is significantly better annotated. Of the 19,763 genes annotated in the AaegL5 reference sequence, our exome capture probes covered 11,076 genes. Overall, the chromosomal distribution of empirical p values computed for each gene and each DENV type (Fig 2) was consistent with the pattern observed with the SNP-based analysis (S1 Fig). The most significant genes were detected on chromosome 1 around the centromere for both DENV types. Although there are regions of similarity between the DENV-1 and DENV-3 comparisons, the genetic architecture underlying DENV-1 and DENV-3 infection phenotypes is markedly distinct (Fig 2). Among the top 5% most significant candidate genes for each DENV type, 287 genes are uniquely associated with DENV-1 infection status, 290 are uniquely associated with DENV-3 infection status, and 263 genes are shared (S2 Table; Fig 3A).

thumbnail
Fig 2. Genes associated with DENV-1 and DENV-3 infection in the Bakoumba population.

(A) Manhattan plot of empirical p values derived from average gene scores of DENV-1 resistance, distributed along the three chromosomes. The gene scores reflect variation in SNP frequencies (adjusted for differences in SNP counts) between phenotypic pools of mosquitoes that are either resistant of susceptible to DENV-1. (B) Manhattan plot of empirical p values derived from average gene scores of DENV-3 resistance, distributed along the three chromosomes. The gene scores reflect variation in SNP frequencies (adjusted for differences in SNP counts) between phenotypic pools of mosquitoes that are either resistant of susceptible to DENV-3. Each dot represents a single gene and is colored according to the statistical significance of the genotype-phenotype association. Blue dots represent the lower 5% of p values and red dots represent the lower 2.5% p values.

https://doi.org/10.1371/journal.pgen.1008794.g002

thumbnail
Fig 3. Limited overlap between genes associated with DENV-3 or DENV-1 infection in the Bakoumba population.

(A) Venn diagram of the 5% most significant genes associated with resistance to DENV-1, resistance to DENV-3, or both. (B) TreeMap visualization of molecular function enrichment for the 5% most significant genes. Each first level category is represented by a large rectangle filled with nested rectangles corresponding to the GO term annotation. The size of the rectangle is proportional to the number of genes in the list out of the total number of genes with similar roles in the Ae. aegypti genome. Bold font represents a higher category level.

https://doi.org/10.1371/journal.pgen.1008794.g003

Finally, we sought to determine whether the 5% most significant genes in the gene-based analysis belonged to different functional categories. We focused on the molecular mechanisms potentially underlying DENV type-specific resistance and accordingly we extracted gene ontology (GO) terms that exhibited a significant enrichment in the molecular function category. The 290 genes uniquely associated with resistance to DENV-3 infection were enriched in genes with zinc ion binding activity (Fig 3B). The 287 genes uniquely associated with resistance to DENV-1 were not significantly enriched in any particular molecular function. The 263 genes shared between the two comparisons were enriched in genes with ATP binding activity and sulfur compound transmembrane transporter activity (Fig 3B).

Discussion

Our EWAS revealed that the DENV type-specific resistance phenotype displayed by an Ae. aegypti population from Bakoumba, Gabon reflects a distinct genetic architecture of resistance to DENV-1 and DENV-3 infection. The midgut infection phenotype we focused on is an important component of vector competence, which also depends on subsequent virus dissemination from the midgut and virus release in saliva [37]. We unveiled a set of SNPs and genes that are significantly associated with DENV-1 and/or DENV-3 infection in the Bakoumba population, which open the way to future studies to functionally characterize the genetic basis of DENV type-specific resistance. The overlap between SNPs and genes associated with resistance to either DENV-1 or DENV-3 was only partial and demonstrated that DENV type-specific resistance relies on both shared and unique genes in this Ae. aegypti population. Although the mechanisms of DENV type-specific resistance in Ae. aegypti remain to be elucidated, this study is a leap forward to understand G x G specificity between Ae. aegypti and DENV at the molecular level.

G x G specificity between pathogens and their invertebrate hosts has been observed in many systems [1], but to date the molecular nature of this specificity is largely unknown. We previously detected G x G interactions between DENV and Ae. aegypti [21, 22] that could be statistically associated with QTL in the mosquito genome [15, 17]. However, the specific mosquito genes underlying these QTL have not been identified. One important question to elucidate the genetic basis of host-pathogen specificity is whether pathogen genotype-specific resistance relies on different alleles or different genes of the host. In the “matching-allele” model of host-pathogen interactions, a resistant allele against a given pathogen strain becomes a susceptible allele against another strain of the same pathogen [38]. This was observed, for example, by experimental viral evolution to specific mouse MHC genotypes [39] and by genetic crosses between genotypes of the crustacean Daphnia magna with different levels of resistance to the parasitic bacterium Pasteuria ramosa [40, 41]. In other host-pathogen systems, host resistance to different strains of the same pathogen relies on allelic variation at different genes. For example, genome-to-genome analysis of associations between human genetic variation and HIV-1 sequence found SNPs in different HLA class I genes [42]. Likewise, Arabidopsis thaliana SNPs exhibited genetic effects only when paired with certain Xanthomonas arboricola variants [43]. In the mosquito Anopheles gambiae, discrimination between human and rodent malaria species was shown to rely on different paralogs of the APL1 gene family [44]. Resistant and susceptible alleles of the TEP1 gene of An. gambiae explained resistance to some but not all strains of the human malaria parasite Plasmodium falciparum [45, 46].

We took advantage of the striking DENV type-specific phenotype of our Bakoumba Ae. aegypti colony to use it as a simplified model of G x G interactions that could be subjected to a genetic association study. Note that in this study, we focused on the host factors involved in the G x G interaction but the virus factors remain to be elucidated. We found that in the Bakoumba population, about half of the genes associated with resistance to DENV-1 infection did not confer resistance to DENV-3 infection and vice versa. Whether the genetic architecture of DENV type-specific resistance in the Bakoumba population is representative of G x G interactions between Ae. aegypti and DENV, and in other invertebrate host-pathogen systems in general, is unknown. Additional work will be necessary to extend these results to other mosquito populations, other DENV types and strains, and other host-pathogen systems. Moreover, our association study focused on the discovery of SNPs in exons, which represent only about 2% of the entire Ae. aegypti genome [17, 28]. It is possible that our exome-sequencing approach may have missed important genetic determinants underlying DENV type-specific resistance because they are not readily detected by SNPs (e.g., structural variants) and/or they occur in introns or intergenic regions. As sequencing technologies become more affordable, whole-genome sequencing will allow more comprehensive surveys of genetic variation in the Ae. aegypti genome. Nevertheless, detection of different exome variants associated with resistance to DENV-1 and resistance to DENV-3 is sufficient to support our conclusion that their respective genetic architecture is distinct.

It is worth noting that regardless of the biases introduced by exome sequencing, these biases are consistent across experimental conditions. In other words, comparing the genetic basis of resistance to DENV-1 infection to the genetic basis of resistance to DENV-3 infection remains valid despite shortcomings in the genotyping method, because it was performed with the same initial population. Our finding that only about one third of the top 5% genes was jointly associated with both DENV-1 resistance and DENV-3 resistance cannot be explained by differences in genetic ascertainment because for both DENV types the genes were identified with the same set of SNPs in the same population.

The bulk of significant SNPs and genes were found on chromosome 1, particularly around the centromere. This observation could be due to at least two, non-mutually exclusive explanations. First, most of the genetic loci underlying DENV susceptibility in the Bakoumba population could be truly located on chromosome 1. QTL associated with DENV susceptibility in Ae. aegypti were previously detected on all three chromosomes [15, 17, 47, 48], and the genetic architecture of pathogen susceptibility is known to vary across host populations [49]. Second, the centromeric region of chromosome 1 contains the M locus responsible for sex determination in Ae. aegypti [17, 50]. We previously reported that the M locus is surrounded by a large region of reduced recombination between the sex chromosomes, which is associated with sex-specific genetic differentiation and elevated linkage disequilibrium (LD) over about 100 Mb in natural Ae. aegypti populations [51]. Because statistical power to detect genotype-phenotype associations increases with LD, it is possible that elevated LD may have inflated the proportion of significant genotype-phenotype associations in the centromeric region of chromosome 1. Nevertheless, as mentioned above, comparisons using the same population and the same markers remain valid because any genotyping bias should be consistent across experimental conditions.

The candidate genes that we identified deserve further investigation at the functional level. Although a genetic association does not necessarily imply differential expression, it would be interesting to determine whether candidate genes also display allelic differences in gene expression. Gene knockdown assays could be used to confirm the functional implication of the candidate genes [52, 53], although gene knockouts would provide a more definitive answer [54] because they are less prone to false negatives. Mosquito genes uniquely associated with resistance to DENV-3 infection were enriched in genes with zinc ion binding activity, whereas genes associated with resistance to both DENV-1 and DENV-3 infection were enriched in genes with ATP binding activity and sulfur compound transmembrane transporter activity. Zinc is an essential cofactor that ensures the proper folding and functioning of not only cellular proteins but also viral proteins [55]. To our knowledge, there is no prior evidence for a link between zinc ion binding activity and mosquito-virus interactions, however host cellular systems controlling zinc balance are known to interfere with virus replication [55]. Likewise, ATP binding activity and sulfur compound transmembrane transporter activity have not been specifically reported to participate in mosquito-virus interactions, however the high dependence of viruses on the cellular machinery makes any molecular function potentially relevant to host-virus interactions.

Although the forces driving the evolution of DENV resistance in Ae. aegypti are largely unknown, one evolutionary implication of the distinct gene set underlying resistance to DENV-1 and DENV-3 is that the evolution of resistance to one DENV type is not expected to lead to the correlated evolution of resistance to another DENV type. This finding is consistent with the absence of virus cross-resistance [56] and the lack of genetic trade-offs between the levels of resistance to different viral genotypes [57] in Drosophila.

The epidemiological relevance of our results is difficult to assess because dengue epidemiology is poorly documented in Gabon. A recent study reported DENV-3 circulation in 2016–2017 [58] whereas previous dengue outbreaks were mainly associated with DENV-2 in 2007 and in 2010 [59], but this information is insufficient to make a link between the DENV-1 resistance phenotype of the Bakoumba population and the relative lack of this DENV type in recently reported outbreaks in Gabon. In our experiments, we used a DENV-1 isolate from Thailand and a DENV-3 isolate from Gabon but the geographical origin of the virus is unlikely to have influenced the results due to the lack of evidence for DENV local adaptation to Ae. aegypti populations [60].

In conclusion, our results support a model where the genetic basis of resistance to DENV infection in Ae. aegypti has a “universal” component that acts across types and strains, and a type- or strain-specific component [1517]. In the Bakoumba population, the DENV type-specific component of resistance consists, at least in part, of a distinct set of genes. Our gene enrichment analysis did not identify a functional overlap between the genes uniquely associated with DENV-1 resistance, uniquely associated with DENV-3 resistance, or associated with both. We speculate that Ae. aegypti resistance to DENV infection may involve different pathways and different molecular mechanisms. Probing this hypothesis will require follow-up studies to functionally test the effect of candidate genes by reverse genetics.

Materials and methods

Ethics statement

The Institut Pasteur animal facility has received accreditation from the French Ministry of Agriculture to perform experiments on live animals in compliance with the French and European regulations on care and protection of laboratory animals. This study was approved by the Institutional Animal Care and Use Committee at Institut Pasteur under protocol number 2015–0032.

Cells and virus isolates

C6/36 cells (derived from Ae. albopictus) were cultured in Leibovitz’s L-15 medium (Life Technologies) supplemented with 10% fetal bovine serum (FBS, Life Technologies), 1% non-essential amino acids (Life Technologies) and 0.1% Penicillin-Streptomycin (Life Technologies) at 28°C. DENV-1 isolate KDH0026A was originally obtained in 2010 from the serum of a dengue patient in Kamphaeng Phet, Thailand [15]. DENV-3 isolate GA28-7 was originally obtained in 2010 from the serum of a dengue patient in Moanda, Gabon [59]. Informed consent of the patient was not necessary because viruses isolated in laboratory cell culture are no longer considered human samples. Virus stock was prepared using C6/36 cells, and viral infectious titers were measured on C6/36 cells using a standard focus-forming assay (FFA) as previously described [61].

Mosquito rearing and exposure to DENV

All experiments were carried out with Ae. aegypti mosquito colonies derived from wild populations sampled in Bakoumba, Gabon in 2014, and Cairns, Australia in 2013. Both colonies had been maintained in the laboratory for 8 generations at the time of the experiments. Mosquitoes were reared and exposed to DENV as described previously [61]. Briefly, mosquitoes were reared under standard insectary conditions (28°C ± 1°C, 75 ± 5% relative humidity, 12:12 hour light-dark cycle) with a diet of fish food for larvae (Tetramin), and a 10% sucrose solution for adults. Twenty-four hours before experimental infection with DENV, 5- to 7-day-old female mosquitoes were deprived of sucrose solution. Fresh rabbit erythrocytes were washed and supplemented with ATP at a final concentration of 10 mM to stimulate blood uptake by mosquitoes. The blood meal consisted of a 2:1 mixture of erythrocytes and virus suspension, and was provided for 15 minutes using an artificial membrane feeding system (Hemotek) with pig intestine as a membrane. To expose mosquitoes to different virus concentrations in the dose-response experiments, virus titer in the blood meal was adjusted prior to its preparation by diluting the virus stock in cell culture medium. Aliquots of the blood meal were collected prior to feeding, and viral titers subsequently determined by FFA. Blood-fed females were transferred to 1-pint carton boxes, and kept under controlled conditions (28°C ± 1°C, 70 ± 5% relative humidity, 12:12 hour light-dark cycle) with access to a 10% sucrose solution.

DENV RNA detection

Mosquito bodies (abdomen and thorax) were ground for two rounds of 30 seconds at 5,000 rpm in a mixer mill (Precellys 24, Bertin Technologies) in 300 μL of squash buffer (10 mM Tris-Hcl pH 8.0, 50 mM NaCl, 1.25 mM EDTA diluted in Ultrapure water) supplemented with Proteinase K (Eurobio-Ingen) at a final concentration of 0.35 g/L. Lysate was heated to 56°C for 5 minutes, and 98°C for 10 minutes. Complementary DNA was synthetized using random hexamers and the M-MLV reverse transcriptase (Life Technologies). Five μL of lysate were included in a 20-μL reaction following the manufacturer’s instructions. Finally, a PCR assay using DreamTaq polymerase (Life Technologies), and following manufacturer’s instructions was used to generate amplicons from the conserved NS5 region of DENV-1 (forward 5’-CGAAGATCACTGGTTCAGCA-3’; reverse 5’-ACATCCATCACGGTTCCATT-3’) and DENV-3 (forward 5’-AGAAGGAGAAGGACTGCACA-3’; reverse 5’-ATTCTTGTGTCCCAACCGGCT-3’). Representative pictures of the electrophoresis gels and a rationale for scoring RT-PCR readouts are presented in S2 Fig.

Virus titration

Virus titration was performed by FFA modified from a published protocol [62], as previously described [61]. In brief, a 96-well plate was seeded sub-confluently with C6/36 cells, inoculated with the viral suspension, and covered with an overlay medium containing 1.6% carboxyl methylcellulose solution. After 5 days of incubation at 28°C, cells were fixed using 3.6% formaldehyde. Virus staining was performed using a primary mouse anti-DENV complex monoclonal antibody (MAB8705, Merck Millipore), and a secondary Alexa Fluor 488-conjugated goat anti-mouse antibody (Life Technologies). Infectious foci were counted using a fluorescent microscope and focus-forming units/mL (FFU/mL) subsequently determined.

DNA extraction

DNA was extracted from individual mosquito heads using a protocol originating from Dr. Pat Roman’s laboratory, and adapted from a previous publication [63]. First, the mosquito head was ground during two rounds of 30 seconds at 5,000 rpm in a mixer mill (Precellys 24, Bertin Technologies) in 300 μL of Pat Roman’s buffer (0.1 M NaCl, 0.2 M sucrose, 0.1 M Tris-Hcl pH 8.0, 0.05 M EDTA pH 8.0, 0.5% SDS, adjusted at pH 9.2). Lysate was incubated at 65°C for 35 minutes. While still warm, 42.5 μL of 8 M potassium acetate was added to each tube, followed by gentle mixing. After 30 minutes of incubation on ice, lysate was transferred to a new tube and centrifugated for 15 minutes at 12,000 g. Supernatant was transferred to a fresh tube, and mixed with 500 μL of 96% ethanol. DNA was left to precipitate overnight at -20°C. Then, tubes were centrifugated for 30 minutes at 12,000 g. The pellet was washed once with ice-cold 75% ethanol, and then with ice-cold 96% ethanol, each time followed by 5-minute centrifugation at 12,000 g. Finally, the pellet was resuspended in 30 μL of TE buffer (10 mM Tris-HCl pH 8.0, 1 mM EDTA pH 8.0).

Library preparation and sequencing

The amount of DNA isolated from each head was quantified using the Picogreen (Thermo Fisher Scientific) reagent following the kit protocol. To prepare the 12 sequencing libraries (2 DENV types x 2 phenotypes x 3 replicate experiments), individual DNA samples were combined to obtain 1.0 ng of DNA per pool. Individual DNA concentrations were adjusted prior to pooling so that each individual contributed the same amount of DNA to the pool. Library preparation and exome capture were performed following the detailed Nimblegen SeqCap EZ Library Protocol version 5 provided by Nimblegen (Roche). Briefly, genomic DNA was sheared using a Covaris to around 350-bp fragments and libraries were prepared using the Kapa HTP/LTP (Roche) library preparation kit following kit instructions. At the end of library preparation, the DNA concentration was quantified by Qubit (Thermo Fisher Scientific) and equal amounts of DNA from each library were pooled and allowed to hybridize with the exome-capture probes for 20 hours. Probes targeting the coding DNA sequences (CDS) were designed and synthesized by Nimblegen/Roche following previous studies [3436]. Probe design used the AaegL3 genome build of Ae. aegypti [28] as this was the latest version of the genome at the time of the study. Genome coordinates for probes are provided as the BED file in S1 File. The exome-enriched libraries were sequenced in paired end on an Illumina NextSeq 500 instrument with a high-output v2 kit 300 cycles (Illumina). Raw sequencing data were deposited to the European Nucleotide Archive repository under accession number PRJEB36810.

SNP detection

Sequencing reads were trimmed and cleaned of adaptor sequences using Trimmomatic [64] before mapping onto the AaegL5 genome assembly using BWA MEM [65, 66] The resulting BAM files were further sorted and indexed and PCR duplicate were removed using the SAMtools v1.2 sort, index and rmdup programs, respectively, with default options [67]. The 12 pooled-sequencing BAM files were processed using the mpileup program in SAMtools [67]. Variant calling was performed on the resulting mpileup file using Varscan v2.3.6 [68] with default options except for the minimum coverage threshold (-min-coverage option) set to 10 reads, the variant minor allele frequency threshold (-min-var-freq option) set to 0.01, and the average base quality score threshold (-min-avg-qual option) set to 20. The resulting VCF file was subsequently processed with custom awk scripts to retain only bi-allelic SNPs covered by more than 20 reads and less than 400 reads in each pool. In total, 1,498,923 SNPs were available for association testing with a median coverage ranging from 85X to 149X per pool.

Association testing

Exome-wide association testing of the 12 phenotypic pools was performed using the C2 contrast statistic implemented in the program BayPass 2.2 [31, 32]. This statistic compares allele frequencies between two groups of pools (resistant vs. susceptible), while properly accounting for the overall correlation structure resulting from the shared origin of the pools. Of note, the model implemented in BayPass allows to properly account for the sampling of reads from the (unobserved) allele count data in pooled-sequencing data when estimating and comparing the (unobserved) allele frequencies. Here, the 12 pools originated from the same colony derived from individuals sampled in Bakoumba, Gabon. To improve the statistical power of association tests, a diagonal structure was assumed for the scaled covariance matrix of pool allele frequencies as in the model initially proposed by Nicholson et al. [69]. The elements of the diagonal matrix were estimated as described in Gautier et al. [70] using the option -nicholsonprior newly introduced in version 2.2 of the BayPass program [31]. They can be interpreted as the inverse of the actual effective sample size of the pools and ranged from 2.50% to 14.7% (with a median equal to 3.62%). Three independent runs of BayPass (using the option -seed) were performed to estimate the contrasts for each SNP. The estimated C2 values together with the model hyper-parameters were highly correlated across the different runs (e.g., Pearson’s correlation >0.985 for C2 values) and only the results from the first BayPass run were retained. Support for association between each SNP and the DENV resistance phenotype was evaluated based on the p values derived from the corresponding contrast statistic. The p-value distribution was close to uniform except for extremely small p values. For each DENV type, gene-specific scores were calculated by adding the estimated C2 values of their underlying SNPs. Gene positions were obtained using the base feature annotation file for the AaegL5 genome assembly (available at https://www.vectorbase.org/organisms/aedes-aegypti). Using custom awk scripts, a total of 964,398 SNPs (excluding SNPs with a minor allele frequency <0.05) were assigned to 10,797 different genes with a median number of 62 SNPs per gene (ranging from 1 to 2,377). When a SNP mapped to several genes, it was arbitrarily assigned to only the first ranking one. For each DENV type, support for association of the gene scores with the resistance phenotype was evaluated based on empirical p values computed as follows. Let k represent the number of SNPs mapping to a given gene i and nk the number of loci of k consecutive SNPs (including gene i) over the whole genome. Locus-specific scores were calculated as the sum of the k C2 values estimated for the SNPs underlying the nk loci. The empirical p value of gene i was computed as the proportion of the nk loci with a higher score than the gene i score.

Gene enrichment analysis

Genes with the highest scores were analyzed for functional enrichment by submitting the gene lists to the DAVID Knowledgebase v6.8 [71]. Enriched terms according to the Ae. aegypti background were clustered with REVIGO [72]. A treemap for gene ontology (GO) terms enriched in the molecular function category was generated using an in-house R script. In the treemap, each term is represented by a rectangle whose size is proportional to the number of genes in the list out of the total number of genes with similar roles in the Ae. aegypti genome. Terms with a semantic similarity, e.g. genes with a proximity in the GO graph below 0.7 of similarity, are represented in the same color and the centroid term is emphasized.

Supporting information

S1 Fig. SNPs associated with DENV-1 and DENV-3 infection in the Bakoumba population.

(A) Manhattan plot of p values representing the association between SNP frequency and DENV-1 resistance, distributed along the three chromosomes. SNP frequency was compared between phenotypic pools of mosquitoes that are either resistant of susceptible to DENV-1. (B) Manhattan plot of p values representing the association between SNP frequency and DENV-3 resistance, distributed along the three chromosomes. SNP frequency was compared between phenotypic pools of mosquitoes that are either resistant of susceptible to DENV-3. In (A) and (B), each dot represents a single SNP and is colored according to the statistical significance of the genotype-phenotype association. Blue dots represent the lower 5% of p values and red dots represent the lower 2.5% p values. (C) Venn diagram of the 0.001% most significant SNPs associated with DENV-1 infection, DENV-3 infection, or both.

https://doi.org/10.1371/journal.pgen.1008794.s001

(TIF)

S2 Fig. Scoring of the mosquito body infection phenotype.

Photos of representative electrophoresis gels of RT-PCR products for DENV-1 (A) and DENV-3 (B) detection in mosquito bodies. The readout was based on five scores (exemplified in red font) as follows: 1 = clear and bright band at the right height; 2 = clear and moderately bright band at the right height; 3 = weak band at the right height; 4 = one or several bands at an unexpected height (sometimes accompanied by the right band); 5 = no band. A sample was only considered DENV-positive when its score was 1 or 2. For DENV-1 all five scores were typically present on the gels, whereas for DENV-3 scores 3 and 4 were typically absent. The + and–symbols denote positive and negative controls, respectively.

https://doi.org/10.1371/journal.pgen.1008794.s002

(TIF)

S1 Table. List of the 0.001% most significant SNPs associated with DENV-1 infection, DENV-3 infection, or both.

https://doi.org/10.1371/journal.pgen.1008794.s003

(XLSX)

S2 Table. List of the 5% most significant genes associated with DENV-1 infection, DENV-3 infection, or both.

https://doi.org/10.1371/journal.pgen.1008794.s004

(XLSX)

Acknowledgments

We are grateful to Gordana Rašić and Ary Hoffmann for providing the mosquito colony from Cairns. We thank Catherine Lallemand for assistance with mosquito rearing and three anonymous reviewers for helpful comments on an earlier version of the manuscript.

References

  1. 1. Lambrechts L. Dissecting the genetic architecture of host-pathogen specificity. PLoS Pathog. 2010;6(8):e1001019. pmid:20700450.
  2. 2. Carius HJ, Little TJ, Ebert D. Genetic variation in a host-parasite association: potential for coevolution and frequency-dependent selection. Evolution. 2001;55(6):1136–45. pmid:11475049.
  3. 3. Schmid-Hempel P. On the evolutionary ecology of host-parasite interactions: addressing the question with regard to bumblebees and their parasites. Naturwissenschaften. 2001;88(4):147–58. Epub 2001/08/02. pmid:11480702.
  4. 4. Schulenburg H, Ewbank JJ. Diversity and specificity in the interaction between Caenorhabditis elegans and the pathogen Serratia marcescens. BMC Evol Biol. 2004;4(1):49. pmid:15555070.
  5. 5. Harris C, Lambrechts L, Rousset F, Abate L, Nsango SE, Fontenille D, et al. Polymorphisms in Anopheles gambiae Immune Genes Associated with Natural Resistance to Plasmodium falciparum. PLoS Pathog. 2010;6(9). pmid:20862317.
  6. 6. Lambrechts L, Halbert J, Durand P, Gouagna LC, Koella JC. Host genotype by parasite genotype interactions underlying the resistance of anopheline mosquitoes to Plasmodium falciparum. Malar J. 2005;4(1):3. pmid:15644136.
  7. 7. de Roode JC, Altizer S. Host-parasite genetic interactions and virulence-transmission relationships in natural populations of Monarch butterflies. Evolution. 2009;64(2):502–14. pmid:19796153.
  8. 8. Luijckx P, Ben-Ami F, Mouton L, Du Pasquier L, Ebert D. Cloning of the unculturable parasite Pasteuria ramosa and its Daphnia host reveals extreme genotype-genotype interactions. Ecol Lett. 2011;14(2):125–31. pmid:21091597.
  9. 9. Little TJ, Hultmark D, Read AF. Invertebrate immunity and the limits of mechanistic immunology. Nat Immunol. 2005;6(7):651–4. pmid:15970937.
  10. 10. Schulenburg H, Boehnisch C, Michiels NK. How do invertebrates generate a highly specific innate immune response? Mol Immunol. 2007;44(13):3338–44. Epub 2007/03/30. pmid:17391764.
  11. 11. Dong Y, Taylor HE, Dimopoulos G. AgDscam, a hypervariable immunoglobulin domain-containing receptor of the Anopheles gambiae innate immune system. PLoS Biol. 2006;4(7):e229. Epub 2006/06/16. pmid:16774454.
  12. 12. Watson F, Püttmann-Holgado R, Thomas F, Lamar D, Hughes M, Kondo M, et al. Extensive diversity of Ig-superfamily proteins in the immune system of insects. Science. 2005;309(5742):1874–8. pmid:16109846
  13. 13. Koch H, Schmid-Hempel P. Gut microbiota instead of host genotype drive the specificity in the interaction of a natural host-parasite system. Ecol Lett. 2012;15(10):1095–103. Epub 2012/07/07. pmid:22765311.
  14. 14. Barribeau SM, Sadd BM, du Plessis L, Schmid-Hempel P. Gene expression differences underlying genotype-by-genotype specificity in a host-parasite system. Proc Natl Acad Sci U S A. 2014;111(9):3496–501. Epub 2014/02/20. pmid:24550506.
  15. 15. Fansiri T, Fontaine A, Diancourt L, Caro V, Thaisomboonsuk B, Richardson JH, et al. Genetic Mapping of Specific Interactions between Aedes aegypti Mosquitoes and Dengue Viruses. PLoS Genet. 2013;9(8):e1003621. Epub 2013/08/13. pmid:23935524.
  16. 16. Lambrechts L, Quillery E, Noel V, Richardson JH, Jarman RG, Scott TW, et al. Specificity of resistance to dengue virus isolates is associated with genotypes of the mosquito antiviral gene Dicer-2. Proc Biol Sci. 2013;280(1751):20122437. Epub 2012/11/30. pmid:23193131.
  17. 17. Matthews BJ, Dudchenko O, Kingan SB, Koren S, Antoshechkin I, Crawford JE, et al. Improved reference genome of Aedes aegypti informs arbovirus vector control. Nature. 2018;563(7732):501–7. Epub 2018/11/16. pmid:30429615.
  18. 18. Hauton C, Smith VJ. Adaptive immunity in invertebrates: a straw house without a mechanistic foundation. Bioessays. 2007;29(11):1138–46. Epub 2007/10/16. pmid:17935208.
  19. 19. Little TJ, Colegrave N, Sadd BM, Schmid-Hempel P. Studying immunity at the whole organism level. Bioessays. 2008;30(4):404–5; Epub 2008/03/19. pmid:18348252.
  20. 20. Schmid-Hempel P. Natural insect host-parasite systems show immune priming and specificity: puzzles to be solved. Bioessays. 2005;27(10):1026–34. pmid:16163710.
  21. 21. Lambrechts L. Quantitative genetics of Aedes aegypti vector competence for dengue viruses: towards a new paradigm? Trends Parasitol. 2011;27(3):111–4. pmid:21215699.
  22. 22. Lambrechts L, Chevillon C, Albright RG, Thaisomboonsuk B, Richardson JH, Jarman RG, et al. Genetic specificity and potential for local adaptation between dengue viruses and mosquito vectors. BMC Evol Biol. 2009;9:160. pmid:19589156.
  23. 23. Katzelnick LC, Fonville JM, Gromowski GD, Bustos Arriaga J, Green A, James SL, et al. Dengue viruses cluster antigenically but not as discrete serotypes. Science. 2015;349(6254):1338–43. Epub 2015/09/19. pmid:26383952.
  24. 24. Gloria-Soria A, Ayala D, Bheecarry A, Calderon-Arguedas O, Chadee DD, Chiappero M, et al. Global genetic diversity of Aedes aegypti. Mol Ecol. 2016;25(21):5377–95. pmid:27671732.
  25. 25. Bhatt S, Gething PW, Brady OJ, Messina JP, Farlow AW, Moyes CL, et al. The global distribution and burden of dengue. Nature. 2013;496(7446):504–7. Epub 2013/04/09. pmid:23563266.
  26. 26. Shaw WR, Catteruccia F. Vector biology meets disease control: using basic research to fight vector-borne diseases. Nat Microbiol. 2019;4(1):20–34. Epub 2018/08/29. pmid:30150735.
  27. 27. Kean J, Rainey SM, McFarlane M, Donald CL, Schnettler E, Kohl A, et al. Fighting Arbovirus Transmission: Natural and Engineered Control of Vector Competence in Aedes Mosquitoes. Insects. 2015;6(1):236–78. Epub 2015/10/16. pmid:26463078.
  28. 28. Nene V, Wortman JR, Lawson D, Haas B, Kodira C, Tu ZJ, et al. Genome sequence of Aedes aegypti, a major arbovirus vector. Science. 2007;316(5832):1718–23. pmid:17510324.
  29. 29. Schlotterer C, Tobler R, Kofler R, Nolte V. Sequencing pools of individuals—mining genome-wide polymorphism data without big funding. Nat Rev Genet. 2014;15(11):749–63. Epub 2014/09/24. pmid:25246196.
  30. 30. Sham P, Bader JS, Craig I, O'Donovan M, Owen M. DNA Pooling: a tool for large-scale association studies. Nat Rev Genet. 2002;3(11):862–71. Epub 2002/11/05. pmid:12415316.
  31. 31. Gautier M. Genome-Wide Scan for Adaptive Divergence and Association with Population-Specific Covariates. Genetics. 2015;201(4):1555–79. Epub 2015/10/21. pmid:26482796.
  32. 32. Olazcuaga L, Loiseau A, Parrinello H, Paris M, Fraimout A, Guedot C, et al. A whole-genome scan for association with invasion success in the fruit fly Drosophila suzukii using contrasts of allele frequencies corrected for population structure. Mol Biol Evol. msaa098. In press. Epub 2020/04/17. https://doi.org/10.1093/molbev/msaa098
  33. 33. Warr A, Robert C, Hume D, Archibald A, Deeb N, Watson M. Exome Sequencing: Current and Future Perspectives. G3. 2015;5(8):1543–50. Epub 2015/07/04. pmid:26139844.
  34. 34. Crawford JE, Alves JM, Palmer WJ, Day JP, Sylla M, Ramasamy R, et al. Population genomics reveals that an anthropophilic population of Aedes aegypti mosquitoes in West Africa recently gave rise to American and Asian populations of this major disease vector. BMC Biol. 2017;15(1):16. pmid:28241828.
  35. 35. Dickson LB, Campbell CL, Juneja P, Jiggins FM, Sylla M, Black WC. Exon-Enriched Libraries Reveal Large Genic Differences Between Aedes aegypti from Senegal, West Africa, and Populations Outside Africa. G3. 2017;7(2):571–82. pmid:28007834.
  36. 36. Juneja P, Ariani CV, Ho YS, Akorli J, Palmer WJ, Pain A, et al. Exome and Transcriptome Sequencing of Aedes aegypti Identifies a Locus That Confers Resistance to Brugia malayi and Alters the Immune Response. PLoS Pathog. 2015;11(3):e1004765. Epub 2015/03/31. pmid:25815506.
  37. 37. Black WC, Bennett KE, Gorrochotegui-Escalante N, Barillas-Mury CV, Fernandez-Salas I, de Lourdes Munoz M, et al. Flavivirus susceptibility in Aedes aegypti. Arch Med Res. 2002;33(4):379–88. pmid:12234528.
  38. 38. Agrawal A, Lively CM. Infection genetics: gene-for-gene versus matching-alleles models and all points in between. Evol Ecol Res. 2002;4:79–90.
  39. 39. Kubinak JL, Ruff JS, Hyzer CW, Slev PR, Potts WK. Experimental viral evolution to specific host MHC genotypes reveals fitness and virulence trade-offs in alternative MHC types. Proc Natl Acad Sci U S A. 2012;109(9):3422–7. Epub 2012/02/11. pmid:22323587.
  40. 40. Bento G, Routtu J, Fields PD, Bourgeois Y, Du Pasquier L, Ebert D. The genetic basis of resistance and matching-allele interactions of a host-parasite system: The Daphnia magna-Pasteuria ramosa model. PLoS Genet. 2017;13(2):e1006596. Epub 2017/02/22. pmid:28222092.
  41. 41. Luijckx P, Fienberg H, Duneau D, Ebert D. A matching-allele model explains host resistance to parasites. Curr Biol. 2013;23(12):1085–8. Epub 2013/05/28. pmid:23707426.
  42. 42. Bartha I, Carlson JM, Brumme CJ, McLaren PJ, Brumme ZL, John M, et al. A genome-to-genome analysis of associations between human genetic variation, HIV-1 sequence diversity, and viral control. eLife. 2013;2:e01123. Epub 2013/10/31. pmid:24171102.
  43. 43. Wang M, Roux F, Bartoli C, Huard-Chauveau C, Meyer C, Lee H, et al. Two-way mixed-effects methods for joint association analysis using both host and pathogen genomes. Proc Natl Acad Sci U S A. 2018;115(24):E5440–E9. Epub 2018/06/01. pmid:29848634.
  44. 44. Mitri C, Jacques JC, Thiery I, Riehle MM, Xu J, Bischoff E, et al. Fine pathogen discrimination within the APL1 gene family protects Anopheles gambiae against human and rodent malaria species. PLoS Pathog. 2009;5(9):e1000576. Epub 2009/09/15. pmid:19750215.
  45. 45. Molina-Cruz A, DeJong RJ, Ortega C, Haile A, Abban E, Rodrigues J, et al. Some strains of Plasmodium falciparum, a human malaria parasite, evade the complement-like system of Anopheles gambiae mosquitoes. Proc Natl Acad Sci U S A. 2012;109(28):E1957–62. Epub 2012/05/25. pmid:22623529.
  46. 46. White BJ, Lawniczak MK, Cheng C, Coulibaly MB, Wilson MD, Sagnon N, et al. Adaptive divergence between incipient species of Anopheles gambiae increases resistance to Plasmodium. Proc Natl Acad Sci U S A. 2010;108(1):244–9. pmid:21173248.
  47. 47. Bosio CF, Fulton RE, Salasek ML, Beaty BJ, Black WC. Quantitative trait loci that control vector competence for dengue-2 virus in the mosquito Aedes aegypti. Genetics. 2000;156(2):687–98. pmid:11014816.
  48. 48. Gomez-Machorro C, Bennett KE, del Lourdes Munoz M, Black WC. Quantitative trait loci affecting dengue midgut infection barriers in an advanced intercross line of Aedes aegypti. Insect Mol Biol. 2004;13(6):637–48. pmid:15606812.
  49. 49. Wilfert L, Schmid-Hempel P. The genetic architecture of susceptibility to parasites. BMC Evol Biol. 2008;8:187. pmid:18590517.
  50. 50. Hall AB, Basu S, Jiang X, Qi Y, Timoshevskiy VA, Biedler JK, et al. A male-determining factor in the mosquito Aedes aegypti. Science. 2015;348(6240):1268–70. Epub 2015/05/23. pmid:25999371.
  51. 51. Fontaine A, Filipovic I, Fansiri T, Hoffmann AA, Cheng C, Kirkpatrick M, et al. Extensive Genetic Differentiation between Homomorphic Sex Chromosomes in the Mosquito Vector, Aedes aegypti. Genome biology and evolution. 2017;9(9):2322–35. pmid:28945882.
  52. 52. Merkling SH, Raquin V, Dabo S, Henrion-Lacritick A, Blanc H, Moltini-Conclois I, et al. Tudor-SN Promotes Early Replication of Dengue Virus in the Aedes aegypti Midgut. iScience. 2020;23(2):100870. Epub 2020/02/15. pmid:32059176.
  53. 53. Raquin V, Merkling SH, Gausson V, Moltini-Conclois I, Frangeul L, Varet H, et al. Individual co-variation between viral RNA load and gene expression reveals novel host factors during early dengue virus infection of the Aedes aegypti midgut. PLoS Negl Trop Dis. 2017;11(12):e0006152. Epub 2017/12/21. pmid:29261661.
  54. 54. Kistler KE, Vosshall LB, Matthews BJ. Genome engineering with CRISPR-Cas9 in the mosquito Aedes aegypti. Cell Rep. 2015;11(1):51–60. Epub 2015/03/31. pmid:25818303.
  55. 55. Lazarczyk M, Favre M. Role of Zn2+ ions in host-virus interactions. J Virol. 2008;82(23):11486–94. Epub 2008/09/13. pmid:18787005.
  56. 56. Cogni R, Cao C, Day JP, Bridson C, Jiggins FM. The genetic architecture of resistance to virus infection in Drosophila. Mol Ecol. 2016;25(20):5228–41. Epub 2016/07/28. pmid:27460507.
  57. 57. Carpenter JA, Hadfield JD, Bangham J, Jiggins FM. Specific interactions between host and parasite genotypes do not act as a constraint on the evolution of antiviral resistance in Drosophila. Evolution. 2012;66(4):1114–25. Epub 2012/04/11. pmid:22486692.
  58. 58. Abe H, Ushijima Y, Loembe MM, Bikangui R, Nguema-Ondo G, Mpingabo PI, et al. Re-emergence of dengue virus serotype 3 infections in Gabon in 2016–2017, and evidence for the risk of repeated dengue virus infections. Int J Infect Dis. 2020;91:129–36. Epub 2019/12/11. pmid:31821892.
  59. 59. Caron M, Grard G, Paupy C, Mombo IM, Bikie Bi Nso B, Kassa Kassa FR, et al. First evidence of simultaneous circulation of three different dengue virus serotypes in Africa. PLoS ONE. 2013;8(10):e78030. Epub 2013/11/10. pmid:24205075.
  60. 60. Fansiri T, Pongsiri A, Klungthong C, Ponlawat A, Thaisomboonsuk B, Jarman RG, et al. No evidence for local adaptation of dengue viruses to mosquito vector populations in Thailand. Evol Appl. 2016;9(4):608–18. Epub 2016/04/22. pmid:27099625.
  61. 61. Fontaine A, Jiolle D, Moltini-Conclois I, Lequime S, Lambrechts L. Excretion of dengue virus RNA by Aedes aegypti allows non-destructive monitoring of viral dissemination in individual mosquitoes. Sci Rep. 2016;6:24885. Epub 2016/04/28. pmid:27117953.
  62. 62. Payne AF, Binduga-Gajewska I, Kauffman EB, Kramer LD. Quantitation of flaviviruses by fluorescent focus assay. J Virol Methods. 2006;134(1–2):183–9. pmid:16510196.
  63. 63. Black WC, DuTeau NM. RAPD-PCR and SSCP analysis for insect population genetic studies. In: Crampton JM, Beard CB, Louis C, editors. The Molecular Biology of Insect Disease Vectors: A Methods Manual. Dordrecht: Springer Netherlands; 1997. p. 361–73.
  64. 64. Bolger AM, Lohse M, Usadel B. Trimmomatic: a flexible trimmer for Illumina sequence data. Bioinformatics. 2014;30(15):2114–20. pmid:24695404.
  65. 65. Li H. Aligning sequence reads, clone sequences and assembly contigs with BWA-MEM2013. arXiv:1303.3997 [q-bio.GN].
  66. 66. Li H, Durbin R. Fast and accurate short read alignment with Burrows-Wheeler transform. Bioinformatics. 2009;25(14):1754–60. Epub 2009/05/20. pmid:19451168.
  67. 67. Li H, Handsaker B, Wysoker A, Fennell T, Ruan J, Homer N, et al. The Sequence Alignment/Map format and SAMtools. Bioinformatics. 2009;25(16):2078–9. pmid:19505943.
  68. 68. Koboldt DC, Chen K, Wylie T, Larson DE, McLellan MD, Mardis ER, et al. VarScan: variant detection in massively parallel sequencing of individual and pooled samples. Bioinformatics. 2009;25(17):2283–5. Epub 2009/06/23. pmid:19542151.
  69. 69. Nicholson G, Smith AV, Jonsson F, Gustafsson O, Stefansson K, Donnelly P. Assessing population differentiation and isolation from single-nucleotide polymorphism data. Journal of the Royal Statistical Society B. 2002;64:695–715.
  70. 70. Gautier M, Hocking TD, Foulley JL. A Bayesian outlier criterion to detect SNPs under selection in large data sets. PLoS One. 2010;5(8):e11913. Epub 2010/08/07. pmid:20689851.
  71. 71. Sherman BT, Huang da W, Tan Q, Guo Y, Bour S, Liu D, et al. DAVID Knowledgebase: a gene-centered database integrating heterogeneous gene annotation resources to facilitate high-throughput gene functional analysis. BMC bioinformatics. 2007;8:426. Epub 2007/11/06. pmid:17980028.
  72. 72. Supek F, Bosnjak M, Skunca N, Smuc T. REVIGO summarizes and visualizes long lists of gene ontology terms. PLoS One. 2011;6(7):e21800. Epub 2011/07/27. pmid:21789182.