Identifying insecticide resistance genes in mosquito by combining AFLP genome scans and 454 pyrosequencing

Neglected impacts of plant protection products on invertebrate aquatic biodiversity: a focus on eco-evolutionary processes

The application of plant protection products (PPPs) may have delayed and long-term non-intentional impacts on aquatic invertebrates inhabiting agricultural landscapes. Such effects may induce population responses based on developmental and transgenerational plasticity, selection of genetic resistance, as well as increased extirpation risks associated with random genetic drift. While the current knowledge on such effects of PPPs is still scarce in non-target aquatic invertebrate species, evidences are accumulating that support the need for consideration of evolutionary components of the population response to PPPs in standard procedures of risk assessment. This mini-review, as part of a contribution to the collective scientific assessment on PPP impacts on biodiversity and ecosystem services performed in the period 2020-2022, presents a brief survey of the current results published on the subject, mainly in freshwater crustaceans, and proposes some research avenues and strategies that we feel relevant to fill this gap.

Role of plasmid plasticity and mobile genetic elements in the entomopathogen Bacillus thuringiensis serovar israelensis

Bacillus thuringiensis is a well-known biopesticide that has been used for more than 80 years. This spore-forming bacterium belongs to the group of Bacillus cereus that also includes, among others, emetic and diarrheic pathotypes of B. cereus, the animal pathogen Bacillus anthracis and the psychrotolerant Bacillus weihenstephanensis. Bacillus thuringiensis is rather unique since it has adapted its lifestyle as an efficient pathogen of specific insect larvae. One of the peculiarities of B. thuringiensis strains is the extent of their extrachromosomal pool, with strains harbouring more than 10 distinct plasmid molecules. Among the numerous serovars of B. thuringiensis, 'israelensis' is certainly emblematic since its host spectrum is apparently restricted to dipteran insects like mosquitoes and black flies, vectors of human and animal diseases such as malaria, yellow fever, or river blindness. In this review, the putative role of the mobile gene pool of B. thuringiensis serovar israelensis in its pathogenicity and dedicated lifestyle is reviewed, with specific emphasis on the nature, diversity, and potential mobility of its constituents. Variations among the few related strains of B. thuringiensis serovar israelensis will also be reported and discussed in the scope of this specialised insect pathogen, whose lifestyle in the environment remains largely unknown.

Bacterial microbiota of Aedes aegypti mosquito larvae is altered by intoxication with Bacillus thuringiensis israelensis

Background

Insect microbiota is a dynamic microbial community that can actively participate in defense against pathogens. Bacillus thuringiensis (Bt) is a natural entomopathogen widely used as a bioinsecticide for pest control. Although Bt’s mode of action has been extensively studied, whether the presence of microbiota is mandatory for Bt to effectively kill the insect is still under debate. An association between a higher tolerance and a modified microbiota was already evidenced but a critical point remained to be solved: is the modified microbiota a cause or a consequence of a higher tolerance to Bt?

Methods

In this study we focused on the mosquito species Aedes aegypti, as no work has been performed on Diptera on this topic to date, and on B. thuringiensis israelensis (Bti), which is used worldwide for mosquito control. To avoid using antibiotics to cure bacterial microbiota, mosquito larvae were exposed to an hourly increasing dose of Bti during 25 hours to separate the most susceptible larvae dying quickly from more tolerant individuals, with longer survival.

Results

Denaturing gradient gel electrophoresis (DGGE) fingerprinting revealed that mosquito larval bacterial microbiota was strongly affected by Bti infection after only a few hours of exposure. Bacterial microbiota from the most tolerant larvae showed the lowest diversity but the highest inter-individual differences. The proportion of Bti in the host tissue was reduced in the most tolerant larvae as compared to the most susceptible ones, suggesting an active control of Bti infection by the host.

Conclusions

Here we show that a modified microbiota is associated with a higher tolerance of mosquitoes to Bti, but that it is rather a consequence of Bti infection than the cause of the higher tolerance. This study paves the way to future investigations aiming at unraveling the role of host immunity, inter-species bacterial competition and kinetics of host colonization by Bti that could be at the basis of the phenotype observed in this study.

Electronic supplementary material

The online version of this article (10.1186/s13071-018-2741-8) contains supplementary material, which is available to authorized users.

Application of cDNA-AFLP to biomarker exploration in a non-model species Grandidierella japonica

Biomarkers of exposure can be used to identify specific contaminants that are adversely affecting aquatic organisms. However, it remains prohibitively costly to investigate multiple novel biomarkers of exposure in a non-model species, despite the development of next-generation sequencing technology. In this study, we focused on the use of cDNA-amplified fragment length polymorphism (AFLP) as a cost-effective biomarker discovery tool to test whether it could identify biomarkers of exposure in the non-model amphipod species Grandidierella japonica. Loci were identified that were differentially expressed in amphipods exposed to reference chemicals (Cu, Zn, and nicotine) and to an environmental sample (road dust) at sublethal concentrations. Eight loci were shown to respond consistently to nicotine at different concentrations, but not to Cu or Zn. Some of the loci also responded to an environmental road dust sample containing nicotine. These findings suggest that loci identified using cDNA-AFLP could be used as biomarkers of nicotine exposure in environmental samples with complex matrices. Further studies with other organisms and toxicants are needed, but we have demonstrated that the use of cDNA-AFLP to identify biomarkers for ecotoxicological studies of non-model species is at least feasible.

Identification of protease m1 zinc metalloprotease conferring resistance to deltamethrin by characterization of an AFLP marker in Culex pipiens pallens

BACKGROUND

Continuous and excessive application of deltamethrin (DM) has resulted in the rapid development of insecticide resistance in Culex pipiens pallens. The quantitative trait loci (QTL) responsible for resistance to DM had previously been detected in Cx. pipiens pallens. But locating the QTLs on the chromosomes remained difficult. An available approach is to first characterize DNA molecular markers linked with the phenotype, and then identify candidate genes.

METHODS

In this study, the amplified fragment length polymorphism (AFLP) marker L3A8.177 associated with the QTL, was characterized. We searched for potential candidate genes in the flank region of L3A8.177 in the genome sequence of the closely related Cx. pipiens quinquefasciatus and conducted mRNA expression analysis of the candidate gene via quantitative real-time PCR. Then the relationship between DM resistance and the candidate gene was identified using RNAi and American CDC Bottle Bioassay in vivo. We also cloned the ORF sequences of the candidate gene from both susceptible and resistant mosquitoes.

RESULTS

The genes CYP6CP1 and protease m1 zinc metalloprotease were in the flank region of L3A8.177 and had significantly different expression levels between susceptible and resistant strains. Protease m1 zinc metalloprotease was significantly up-regulated in the susceptible strains compared with the resistant and remained over-expressed in the susceptible field-collected strains. For deduced amino acid sequences of protease m1 zinc metalloprotease, there was no difference between susceptible and resistant mosquitoes. Knockdown of protease m1 zinc metalloprotease not only decreased the sensitivity of mosquitoes to DM in the susceptible strain but also increased the expression of CYP6CP1, suggesting the role of protease m1 zinc metalloprotease in resistance may be involved in the regulation of the P450 gene expression.

CONCLUSION

Our study represents an example of candidate genes derived from the AFLP marker associated with the QTL and provides the first evidence that protease m1 zinc metalloprotease may play a role in the regulation of DM resistance in Cx. pipiens pallens.

Identification of QTLs Conferring Resistance to Deltamethrin in Culex pipiens pallens

Culex pipiens pallens is the most abundant Culex mosquito species in northern China and is an important vector of bancroftian filariasis and West Nile virus. Deltamethrin is an insecticide that is widely used for mosquito control, however resistance to this and other insecticides has become a major challenge in the control of vector-borne diseases that appear to be inherited quantitatively. Furthermore, the genetic basis of insecticide resistance remains poorly understood. In this study, quantitative trait loci (QTL) mapping of resistance to deltamethrin was conducted in F₂ intercross segregation populations using bulked segregation analysis (BSA) and amplified fragment length polymorphism markers (AFLP) in Culex pipiens pallens. A genetic linkage map covering 381 cM was constructed and a total of seven QTL responsible for resistance to deltamethrin were detected by composite interval mapping (CIM), which explained 95% of the phenotypic variance. The major QTL in linkage group 2 accounted for 62% of the variance and is worthy of further study. 12 AFLP markers in the map were cloned and the genomic locations of these marker sequences were determined by applying the Basic Local Alignment Search Tool (BLAST) tool to the genome sequence of the closely related Culex quinquefasciatus. Our results suggest that resistance to deltamethrin is a quantitative trait under the control of a major QTL in Culex pipiens pallens. Cloning of related AFLP markers confirm the potential utility for anchoring the genetic map to the physical map. The results provide insight into the genetic architecture of the trait.

The genetic architecture of a complex trait: Resistance to multiple toxins produced by Bacillus thuringiensis israelensis in the dengue and yellow fever vector, the mosquito Aedes aegypti

The bacterial insecticide Bacillus thuringiensis subsp. israelensis (Bti) is an increasingly popular alternative to chemical insecticides for controlling mosquito populations. Because Bti toxicity relies on the action of four main toxins, resistance to Bti is very likely a complex phenotype involving several genes simultaneously. Dissecting the underlying genetic basis thus requires associating a quantitative measure of resistance to genetic variation at many loci in a segregating population. Here, we undertake this task using the dengue and yellow fever vector, the mosquito Aedes aegypti, as a study model. We conducted QTL (Quantitative Trait Locus) and admixture mapping analyses on two controlled crosses and on an artificial admixed population, respectively, all obtained from resistant and susceptible lab strains. We detected 16 QTL regions, among which four QTLs were revealed by different analysis methods. These four robust QTLs explained altogether 29.2% and 62.2% of the total phenotypic variance in the two QTL crosses, respectively. They also all showed a dominant mode of action. In addition, we found six loci showing statistical association with Bti resistance in the admixed population. Five of the supercontigs highlighted in this study contained candidate genes as suggested by their function, or by prior evidence from expression and/or outlier analyses. These genomic regions are thus good starting points for fine mapping of resistance to Bti or functional analyses aiming at identifying the underlying genes and mutations. Moreover, for the purpose of this work, we built the first Ae. aegypti genetic map based on markers associated with genes expressed in larvae. This genetic map harbors 229 SNP markers mapped across the three chromosomes for a total length of 311.9cM. It brought to light several assembly discrepancies with the reference genome, suggesting a high level of genome plasticity in Ae. aegypti.

Towards the identification of the loci of adaptive evolution

1. Establishing the genetic and molecular basis underlying adaptive traits is one of the major goals of evolutionary geneticists in order to understand the connection between genotype and phenotype and elucidate the mechanisms of evolutionary change. Despite considerable effort to address this question, there remain relatively few systems in which the genes shaping adaptations have been identified. 2. Here, we review the experimental tools that have been applied to document the molecular basis underlying evolution in several natural systems, in order to highlight their benefits, limitations and suitability. In most cases, a combination of DNA, RNA and functional methodologies with field experiments will be needed to uncover the genes and mechanisms shaping adaptation in nature.

Neutral and adaptive drivers of microgeographic genetic divergence within continuous populations: the case of the neotropical tree Eperua falcata (Aubl.)

BACKGROUND

In wild plant populations, genetic divergence within continuous stands is common, sometimes at very short geographical scales. While restrictions to gene flow combined with local inbreeding and genetic drift may cause neutral differentiation among subpopulations, microgeographical variations in environmental conditions can drive adaptive divergence through natural selection at some targeted loci. Such phenomena have recurrently been observed in plant populations occurring across sharp environmental boundaries, but the interplay between selective processes and neutral genetic divergence has seldom been studied.

METHODS

We assessed the extent of within-stand neutral and environmentally-driven divergence in the Neotropical tree Eperua falcate Aubl. (Fabaceae) through a genome-scan approach. Populations of this species grow in dense stands that cross the boundaries between starkly contrasting habitats. Within-stand phenotypic and candidate-gene divergence have already been proven, making this species a suitable model for the study of genome-wide microgeographic divergence. Thirty trees from each of two habitats (seasonally flooded swamps and well-drained plateaus) in two separate populations were genotyped using thousands of AFLPs markers. To avoid genotyping errors and increase marker reliability, each sample was genotyped twice and submitted to a rigorous procedure for data cleaning, which resulted in 1196 reliable and reproducible markers.

RESULTS

Despite the short spatial distances, we detected within-populations genetic divergence, probably caused by neutral processes, such as restrictions in gene flow. Moreover, habitat-structured subpopulations belonging to otherwise continuous stands also diverge in relation to environmental variability and habitat patchiness: we detected convincing evidence of divergent selection at the genome-wide level and for a fraction of the analyzed loci (comprised between 0.25% and 1.6%). Simulations showed that the levels of differentiation for these outliers are compatible with scenarios of strong divergent selection.

Continental-scale footprint of balancing and positive selection in a small rodent (Microtus arvalis)

Genetic adaptation to different environmental conditions is expected to lead to large differences between populations at selected loci, thus providing a signature of positive selection. Whereas balancing selection can maintain polymorphisms over long evolutionary periods and even geographic scale, thus leads to low levels of divergence between populations at selected loci. However, little is known about the relative importance of these two selective forces in shaping genomic diversity, partly due to difficulties in recognizing balancing selection in species showing low levels of differentiation. Here we address this problem by studying genomic diversity in the European common vole (Microtus arvalis) presenting high levels of differentiation between populations (average F_ST = 0.31). We studied 3,839 Amplified Fragment Length Polymorphism (AFLP) markers genotyped in 444 individuals from 21 populations distributed across the European continent and hence over different environmental conditions. Our statistical approach to detect markers under selection is based on a Bayesian method specifically developed for AFLP markers, which treats AFLPs as a nearly codominant marker system, and therefore has increased power to detect selection. The high number of screened populations allowed us to detect the signature of balancing selection across a large geographic area. We detected 33 markers potentially under balancing selection, hence strong evidence of stabilizing selection in 21 populations across Europe. However, our analyses identified four-times more markers (138) being under positive selection, and geographical patterns suggest that some of these markers are probably associated with alpine regions, which seem to have environmental conditions that favour adaptation. We conclude that despite favourable conditions in this study for the detection of balancing selection, this evolutionary force seems to play a relatively minor role in shaping the genomic diversity of the common vole, which is more influenced by positive selection and neutral processes like drift and demographic history.

Transcriptome profiling of pyrethroid resistant and susceptible mosquitoes in the malaria vector, Anopheles sinensis

BACKGROUND

Anopheles sinensis is a major malaria vector in China and other Southeast Asian countries, and it is becoming increasingly resistant to the insecticides used for agriculture, net impregnation, and indoor residual spray. Very limited genomic information on this species is available, which has hindered the development of new tools for resistance surveillance and vector control. We used the 454 GS FLX system and generated expressed sequence tag (EST) databases of various life stages of An. sinensis, and we determined the transcriptional differences between deltamethrin resistant and susceptible mosquitoes.

RESULTS

The 454 GS FLX transcriptome sequencing yielded a total of 624,559 reads (average length of 290 bp) with the pooled An. sinensis mosquitoes across various development stages. The de novo assembly generated 33,411 contigs with average length of 493 bp. A total of 8,057 ESTs were generated with Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) annotation. A total of 2,131 ESTs were differentially expressed between deltamethrin resistant and susceptible mosquitoes collected from the same field site in Jiangsu, China. Among these differentially expressed ESTs, a total of 294 pathways were mapped to the KEGG database, with the predominant ESTs belonging to metabolic pathways. Furthermore, a total of 2,408 microsatellites and 15,496 single nucleotide polymorphisms (SNPs) were identified.

CONCLUSIONS

The annotated EST and transcriptome databases provide a valuable genomic resource for further genetic studies of this important malaria vector species. The differentially expressed ESTs associated with insecticide resistance identified in this study lay an important foundation for further functional analysis. The identified microsatellite and SNP markers will provide useful tools for future population genetic and comparative genomic analyses of malaria vectors.

Insecticide-driven patterns of genetic variation in the dengue vector Aedes aegypti in Martinique Island

Effective vector control is currently challenged worldwide by the evolution of resistance to all classes of chemical insecticides in mosquitoes. In Martinique, populations of the dengue vector Aedes aegypti have been intensively treated with temephos and deltamethrin insecticides over the last fifty years, resulting in heterogeneous levels of resistance across the island. Resistance spreading depends on standing genetic variation, selection intensity and gene flow among populations. To determine gene flow intensity, we first investigated neutral patterns of genetic variability in sixteen populations representative of the many environments found in Martinique and experiencing various levels of insecticide pressure, using 6 microsatellites. Allelic richness was lower in populations resistant to deltamethrin, and consanguinity was higher in populations resistant to temephos, consistent with a negative effect of insecticide pressure on neutral genetic diversity. The global genetic differentiation was low, suggesting high gene flow among populations, but significant structure was found, with a pattern of isolation-by-distance at the global scale. Then, we investigated adaptive patterns of divergence in six out of the 16 populations using 319 single nucleotide polymorphisms (SNPs). Five SNP outliers displaying levels of genetic differentiation out of neutral expectations were detected, including the kdr-V1016I mutation in the voltage-gated sodium channel gene. Association tests revealed a total of seven SNPs associated with deltamethrin resistance. Six other SNPs were associated with temephos resistance, including two non-synonymous substitutions in an alkaline phosphatase and in a sulfotransferase respectively. Altogether, both neutral and adaptive patterns of genetic variation in mosquito populations appear to be largely driven by insecticide pressure in Martinique.

Investigating the genetics of Bti resistance using mRNA tag sequencing: application on laboratory strains and natural populations of the dengue vector Aedes aegypti

Mosquito control is often the main method used to reduce mosquito-transmitted diseases. In order to investigate the genetic basis of resistance to the bio-insecticide Bacillus thuringiensis subsp. israelensis (Bti), we used information on polymorphism obtained from cDNA tag sequences from pooled larvae of laboratory Bti-resistant and susceptible Aedes aegypti mosquito strains to identify and analyse 1520 single nucleotide polymorphisms (SNPs). Of the 372 SNPs tested, 99.2% were validated using DNA Illumina GoldenGate® array, with a strong correlation between the allelic frequencies inferred from the pooled and individual data (r = 0.85). A total of 11 genomic regions and five candidate genes were detected using a genome scan approach. One of these candidate genes showed significant departures from neutrality in the resistant strain at sequence level. Six natural populations from Martinique Island were sequenced for the 372 tested SNPs with a high transferability (87%), and association mapping analyses detected 14 loci associated with Bti resistance, including one located in a putative receptor for Cry11 toxins. Three of these loci were also significantly differentiated between the laboratory strains, suggesting that most of the genes associated with resistance might differ between the two environments. It also suggests that common selected regions might harbour key genes for Bti resistance.

Genomic distribution of AFLP markers relative to gene locations for different eukaryotic species

BACKGROUND

Amplified fragment length polymorphism (AFLP) markers are frequently used for a wide range of studies, such as genome-wide mapping, population genetic diversity estimation, hybridization and introgression studies, phylogenetic analyses, and detection of signatures of selection. An important issue to be addressed for some of these fields is the distribution of the markers across the genome, particularly in relation to gene sequences.

RESULTS

Using in-silico restriction fragment analysis of the genomes of nine eukaryotic species we characterise the distribution of AFLP fragments across the genome and, particularly, in relation to gene locations. First, we identify the physical position of markers across the chromosomes of all species. An observed accumulation of fragments around (peri) centromeric regions in some species is produced by repeated sequences, and this accumulation disappears when AFLP bands rather than fragments are considered. Second, we calculate the percentage of AFLP markers positioned within gene sequences. For the typical EcoRI/MseI enzyme pair, this ranges between 28 and 87% and is usually larger than that expected by chance because of the higher GC content of gene sequences relative to intergenic ones. In agreement with this, the use of enzyme pairs with GC-rich restriction sites substantially increases the above percentages. For example, using the enzyme system SacI/HpaII, 86% of AFLP markers are located within gene sequences in A. thaliana, and 100% of markers in Plasmodium falciparun. We further find that for a typical trait controlled by 50 genes of average size, if 1000 AFLPs are used in a study, the number of those within 1 kb distance from any of the genes would be only about 1-2, and only about 50% of the genes would have markers within that distance.

CONCLUSIONS

The high coverage of AFLP markers across the genomes and the high proportion of markers within or close to gene sequences make them suitable for genome scans and detecting large islands of differentiation in the genome. However, for specific traits, the percentage of AFLP markers close to genes can be rather small. Therefore, genome scans directed towards the search of markers closely linked to selected loci can be a difficult task in many instances.

Natural selection and neutral evolution jointly drive population divergence between alpine and lowland ecotypes of the allopolyploid plant Anemone multifida (Ranunculaceae)

Population differentiation can be driven in large part by natural selection, but selectively neutral evolution can play a prominent role in shaping patters of population divergence. The decomposition of the evolutionary history of populations into the relative effects of natural selection and selectively neutral evolution enables an understanding of the causes of population divergence and adaptation. In this study, we examined heterogeneous genomic divergence between alpine and lowland ecotypes of the allopolyploid plant, Anemone multifida. Using peak height and dominant AFLP data, we quantified population differentiation at non-outlier (neutral) and outlier loci to determine the potential contribution of natural selection and selectively neutral evolution to population divergence. We found 13 candidate loci, corresponding to 2.7% of loci, with signatures of divergent natural selection between alpine and lowland populations and between alpine populations (Fst  = 0.074-0.445 at outlier loci), but neutral population differentiation was also evident between alpine populations (FST  = 0.041-0.095 at neutral loci). By examining population structure at both neutral and outlier loci, we determined that the combined effects of selection and neutral evolution are associated with the divergence of alpine populations, which may be linked to extreme abiotic conditions and isolation between alpine sites. The presence of outlier levels of genetic variation in structured populations underscores the importance of separately analyzing neutral and outlier loci to infer the relative role of divergent natural selection and neutral evolution in population divergence.

Challenges and pitfalls in the characterization of anonymous outlier AFLP markers in non-model species: lessons from an ocellated lizard genome scan

In the last few years, dozens of studies have documented the detection of loci influenced by selection from genome scans in a wide range of non-model species. Many of those studies used amplified fragment length polymorphism (AFLP) markers, which became popular for being easily applicable to any organism. However, because they are anonymous markers, AFLPs impose many challenges for their isolation and identification. Most recent AFLP genome scans used capillary electrophoresis (CE), which adds even more obstacles to the isolation of bands with a specific size for sequencing. These caveats might explain the extremely low number of studies that moved from the detection of outlier AFLP markers to their actual isolation and characterization. We document our efforts to characterize a set of outlier AFLP markers from a previous genome scan with CE in ocellated lizards (Lacerta lepida). Seven outliers were successfully isolated, cloned and sequenced. Their sequences are noncoding and show internal indels or polymorphic repetitive elements (microsatellites). Three outliers were converted into codominant markers by using specific internal primers to sequence and screen population variability from undigested DNA. Amplification in closely related lizard species was also achieved, revealing remarkable interspecific conservation in outlier loci sequences. We stress the importance of following up AFLP genome scans to validate selection signatures of outlier loci, but also report the main challenges and pitfalls that may be faced during the process.