Comprehensive characterization of internal and cuticle surface microbiota of laboratory-reared F1 Anopheles albimanus originating from different sites.

BACKGROUND: Research on mosquito-microbe interactions may lead to new tools for mosquito and mosquito-borne disease control. To date, such research has largely utilized laboratory-reared mosquitoes that typically lack the microbial diversity of wild populations. A logical progression in this area involves working under controlled settings using field-collected mosquitoes or, in most cases, their progeny. Thus, an understanding of how laboratory colonization affects the assemblage of mosquito microbiota would aid in advancing mosquito microbiome studies and their applications beyond laboratory settings. METHODS: Using high throughput 16S rRNA amplicon sequencing, the internal and cuticle surface microbiota of F1 progeny of wild-caught adult Anopheles albimanus from four locations in Guatemala were characterized. A total of 132 late instar larvae and 135 2-5 day-old, non-blood-fed virgin adult females that were reared under identical laboratory conditions, were pooled (3 individuals/pool) and analysed. RESULTS: Results showed location-associated heterogeneity in both F1 larval internal (p = 0.001; pseudo-F = 9.53) and cuticle surface (p = 0.001; pseudo-F = 8.51) microbiota, and only F1 adult cuticle surface (p = 0.001; pseudo-F = 4.5) microbiota, with a more homogenous adult internal microbiota (p = 0.12; pseudo-F = 1.6) across collection sites. Overall, ASVs assigned to Leucobacter, Thorsellia, Chryseobacterium and uncharacterized Enterobacteriaceae, dominated F1 larval internal microbiota, while Acidovorax, Paucibacter, and uncharacterized Comamonadaceae, dominated the larval cuticle surface. F1 adults comprised a less diverse microbiota compared to larvae, with ASVs assigned to the genus Asaia dominating both internal and cuticle surface microbiota, and constituting at least 70% of taxa in each microbial niche. CONCLUSIONS: These results suggest that location-specific heterogeneity in filed mosquito microbiota can be transferred to F1 progeny under normal laboratory conditions, but this may not last beyond the F1 larval stage without adjustments to maintain field-derived microbiota. These findings provide the first comprehensive characterization of laboratory-colonized F1 An. albimanus progeny from field-derived mothers. This provides a background for studying how parentage and environmental conditions differentially or concomitantly affect mosquito microbiome composition, and how this can be exploited in advancing mosquito microbiome studies and their applications beyond laboratory settings.

Co-occurrence of kdr Mutations V1016I and F1534C and Its Association With Phenotypic Resistance to Pyrethroids in Aedes aegypti (Diptera: Culicidae) Populations From Costa Rica.

Aedes aegypti (Linnaeus, 1762) is considered the most important mosquito vector species for several arboviruses (e.g., dengue, chikungunya, Zika) in Costa Rica. The primary strategy for the control and prevention of Aedes-borne diseases relies on insecticide-based vector control. However, the emergence of insecticide resistance in the mosquito populations presents a significant threat to these prevention actions. The characterization of the mechanisms driving the insecticide resistance in Ae. aegypti is vital for decision making in vector control programs. Therefore, we analyzed the voltage-gated sodium channel (VGSC) gene for the presence of the V1016I and F1534C kdr mutations in Ae. aegypti populations from Puntarenas and Limon provinces, Costa Rica. The CDC bottle bioassays showed that both Costa Rican Ae. aegypti populations were resistant to permethrin and deltamethrin. In the case of kdr genotyping, results revealed the co-occurrence of V1016I and F1534C mutations in permethrin and deltamethrin-resistant populations, as well as the fixation of the 1534C allele. A strong association between these mutations and permethrin and deltamethrin resistance was found in Puntarenas. Limon did not show this association; however, our results indicate that the Limon population analyzed is not under the same selective pressure as Puntarenas for the VGSC gene. Therefore, our findings make an urgent call to expand the knowledge about the insecticide resistance status and mechanisms in the Costa Rican populations of Ae. aegypti, which must be a priority to develop an effective resistance management plan.

Pyrethroid exposure alters internal and cuticle surface bacterial communities in Anopheles albimanus.

A deeper understanding of the mechanisms underlying insecticide resistance is needed to mitigate its threat to malaria vector control. Following previously identified associations between mosquito microbiota and insecticide resistance, we demonstrate for the first time, the effects of pyrethroid exposure on the microbiota of F1 progeny of field-collected Anopheles albimanus. Larval and adult mosquitoes were exposed to the pyrethroids alphacypermethrin (only adults), permethrin, and deltamethrin. While there were no significant differences in bacterial composition between insecticide-resistant and insecticide-susceptible mosquitoes, bacterial composition between insecticide-exposed and non-exposed mosquitoes was significantly different for alphacypermethrin and permethrin exposure. Along with other bacterial taxa not identified to species, Pantoea agglomerans (a known insecticide-degrading bacterial species) and Pseudomonas fragi were more abundant in insecticide-exposed compared to non-exposed adults, demonstrating that insecticide exposure can alter mosquito bacterial communities. We also show for the first time that the cuticle surfaces of both larval and adult An. albimanus harbor more diverse bacterial communities than their internal microbial niches. Together, these findings demonstrate how insecticide pressure could be selecting for certain bacteria within mosquitoes, especially insecticide-metabolizing bacteria, thus potentially contributing to insecticide resistance.

Development of molecular assays to detect target-site mechanisms associated with insecticide resistance in malaria vectors from Latin America.

BACKGROUND: Malaria remains an important public health problem in Latin America, and the development of insecticide resistance in malaria vectors poses a major threat to malaria elimination efforts. Monitoring of insecticide susceptibility and the determination of the mechanisms involved in insecticide resistance are needed to effectively guide the deployment of appropriate vector control measures. Here, molecular assays have been developed to screen for mutations associated with insecticide resistance on the voltage-gated sodium channel (VGSC) and acetylcholinesterase-1 (Ace-1) genes in four malaria vectors from Latin America. METHODS: Degenerate primers were designed to amplify a partial fragment on the VGSC and Ace-1 genes. Wild-caught individuals for Anopheles albimanus (also historical samples and individuals from a laboratory strain), Anopheles darlingi, Anopheles vestitipennis and Anopheles pseudopunctipennis were used to optimize the PCR assays. All samples were sequenced to validate the PCR results and DNA alignments were constructed for each gene using the unique haplotypes observed. RESULTS: Primers designed successfully amplified the VGSC gene in An. albimanus, An. darlingi, An. vestitipennis and An. pseudopunctipennis, and the Ace-1 gene in both An. albimanus and An. darlingi. DNA sequencing revealed that compared with Anopheles gambiae, there were a total of 29, 28, 21 and 24 single nucleotide polymorphisms (SNPs) on the VGSC gene for An. albimanus (308 bp), An. darlingi (311 bp), An. pseudopunctipennis (263 bp) and An. vestitipennis (254 bp), respectively. On the 459 bp fragment of the Ace-1 gene, a total of 70 SNPs were detected in An. darlingi and 59 SNPs were detected in An. albimanus compared with An. gambiae. The SNPs detected on the VGSC gene were all synonymous. On the Ace-1 gene, non-synonymous substitutions were identified on three different codons. All species showed the homozygous wild-type kdr allele (coding for leucine) at codon 995 (formerly reported as codon 1014) on the VGSC gene, but one sample was heterozygous at codon 280 (formerly reported as codon 119) on the Ace-1 gene, coding for both the resistant (serine) and susceptible (glycine) amino acids. CONCLUSIONS: New molecular assays to amplify and screen the regions of the VGSC and Ace-1 genes associated with insecticide resistance are reported for An. albimanus, An. darlingi, An. vestitipennis, and An. pseudopunctipennis. The development of these PCR assays presents an important advance in the analysis of target-site resistance in malaria vectors in the Americas, and will further facilitate the characterization of insecticide resistance mechanisms in these species.

Contrasting patterns of gene expression indicate differing pyrethroid resistance mechanisms across the range of the New World malaria vector Anopheles albimanus.

Decades of unmanaged insecticide use and routine exposure to agrochemicals have left many populations of malaria vectors in the Americas resistant to multiple classes of insecticides, including pyrethroids. The molecular basis of pyrethroid resistance is relatively uncharacterised in American malaria vectors, preventing the design of suitable resistance management strategies. Using whole transcriptome sequencing, we characterized the mechanisms of pyrethroid resistance in Anopheles albimanus from Peru and Guatemala. An. albimanus were phenotyped as either deltamethrin or alpha-cypermethrin resistant. RNA from 1) resistant, 2) unexposed, and 3) a susceptible laboratory strain of An. albimanus was sequenced and analyzed using RNA-Seq. Expression profiles of the three groups were compared based on the current annotation of the An. albimanus reference genome. Several candidate genes associated with pyrethroid resistance in other malaria vectors were found to be overexpressed in resistant An. albimanus. In addition, gene ontology terms related to serine-type endopeptidase activity, extracellular activity and chitin metabolic process were also commonly overexpressed in the field caught resistant and unexposed samples from both Peru and Guatemala when compared to the susceptible strain. The cytochrome P450 CYP9K1 was overexpressed 14x in deltamethrin and 8x in alpha-cypermethrin-resistant samples from Peru and 2x in deltamethrin-resistant samples from Guatemala, relative to the susceptible laboratory strain. CYP6P5 was overexpressed 68x in deltamethrin-resistant samples from Peru but not in deltamethrin-resistant samples from Guatemala. When comparing overexpressed genes between deltamethrin-resistant and alpha-cypermethrin-resistant samples from Peru, a single P450 gene, CYP4C26, was overexpressed 9.8x (p<0.05) in alpha-cypermethrin-resistant samples. In Peruvian deltamethrin-resistant samples, the knockdown resistance mutation (kdr) variant alleles at position 1014 were rare, with approximately 5% frequency, but in the alpha-cypermethrin-resistant samples, the frequency of these alleles was approximately 15-30%. Functional validation of the candidate genes and the kdr mutation as a resistance marker for alpha-cypermethrin will confirm the role of these mechanisms in conferring pyrethroid resistance.

Molecular evidence for historical presence of knock-down resistance in Anopheles albimanus, a key malaria vector in Latin America.

BACKGROUND: Anopheles albimanus is a key malaria vector in the northern neotropics. Current vector control measures in the region are based on mass distributions of long-lasting insecticidal nets (LLINs) and focal indoor residual spraying (IRS) with pyrethroids. Resistance to pyrethroid insecticides can be mediated by increased esterase and/or multi-function oxidase activity and/or mutations in the voltage-gated sodium channel gene. The aim of this work was to characterize the homologous kdr region of the voltage-gated sodium channel gene in An. albimanus and to conduct a preliminary retrospective analysis of field samples collected in the 1990's, coinciding with a time of intense pyrethroid application related to agricultural and public health insect control in the region. METHODS: Degenerate primers were designed to amplify the homologous kdr region in a pyrethroid-susceptible laboratory strain (Sanarate) of An. albimanus. Subsequently, a more specific primer pair was used to amplify and sequence the region that contains the 1014 codon associated with pyrethroid resistance in other Anopheles spp. (L1014F, L1014S or L1014C). RESULTS: Direct sequencing of the PCR products confirmed the presence of the susceptible kdr allele in the Sanarate strain (L1014) and the presence of homozygous-resistant kdr alleles in field-collected individuals from Mexico (L1014F), Nicaragua (L1014C) and Costa Rica (L1014C). CONCLUSIONS: For the first time, the kdr region in An. albimanus is described. Furthermore, molecular evidence suggests the presence of kdr-type resistance in field-collected An. albimanus in Mesoamerica in the 1990s. Further research is needed to conclusively determine an association between the genotypes and resistant phenotypes, and to what extent they may compromise current vector control efforts.