Phenotypic evidence of deltamethrin resistance and identification of selective sweeps in Anopheles mosquitoes on the Bijagós Archipelago, Guinea-Bissau.

Vector control in the Bijagós Archipelago of Guinea-Bissau currently relies on pyrethroid insecticide-treated nets. However, data on insecticide resistance in Guinea-Bissau is limited. This study identified deltamethrin resistance in the Anopheles gambiae sensu lato complex on Bubaque island using WHO tube tests in November 2022. Whole genome sequencing of An. gambiae sensu stricto mosquitoes identified six single nucleotide polymorphisms (SNPs) previously associated with, or putatively associated with, insecticide resistance: T791M, L995F, N1570Y, A1746S and P1874L in the vgsc gene, and L119V in the gste2 gene. Twenty additional non-synonymous SNPs were identified in insecticide-resistance associated genes. Four of these SNPs were present at frequencies over 5% in the population: T154S, I126F and G26S in the vgsc gene and A65S in ace1. Genome wide selection scans using Garud's H12 statistic identified two selective sweeps: one in chromosome X and one in chromosome 2R. Both selective sweeps overlap with metabolic genes previously associated with insecticide resistance, including cyp9k1 and the cyp6aa/cyp6p gene cluster. This study presents the first phenotypic testing for deltamethrin resistance and the first whole genome sequence data for Anophelesgambiae mosquitoes from the Bijagós, contributing data of significance for vector control policy in this region.

Genomic surveillance of Anopheles mosquitoes on the Bijagós Archipelago using custom targeted amplicon sequencing identifies mutations associated with insecticide resistance.

BACKGROUND: Insecticide resistance is reducing the efficacy of vector control interventions, consequently threatening efforts to control vector-borne diseases, including malaria. Investigating the prevalence of molecular markers of resistance is a useful tool for monitoring the spread of insecticide resistance in disease vectors. The Bijagós Archipelago (Bijagós) in Guinea-Bissau is a region of stable malaria transmission where insecticide-treated nets are the mainstay for malaria control. However, the prevalence of molecular markers of insecticide resistance in malaria vectors is not well understood. METHODS: A total of 214 Anopheles mosquitoes were analysed from 13 islands across the Bijagós. These mosquitoes were collected using CDC light traps in November 2019, during the peak malaria transmission season. High-throughput multiplex amplicon sequencing was used to investigate the prevalence of 17 different molecular markers associated with insecticide resistance in four genes: vgsc, rdl, ace1 and gste2. RESULTS: Of the 17 screened mutations, four were identified in mosquitoes from the Bijagós: vgsc L995F (12.2%), N1570Y (6.2%) and A1746S (0.7%) and rdl A269G (1.1%). This study is the first to report the L995F knock-down resistance (kdr)-west allele in Anopheles melas on the Archipelago. An additional eight non-synonymous single-nucleotide polymorphisms were identified across the four genes which have not been described previously. The prevalences of the vgsc L995F and N1570Y mutations were higher on Bubaque Island than on the other islands in this study; Bubaque is the most populous island in the archipelago, with the greatest population mobility and connection to continental Guinea-Bissau. CONCLUSIONS: This study provides the first surveillance data for genetic markers present in malaria vectors from islands across the Bijagós Archipelago. Overall prevalence of insecticide resistance mutations was found to be low. However, the identification of the vgsc L995F and N1570Y mutations associated with pyrethroid resistance warrants further monitoring. This is particularly important as the mainstay of malaria control on the islands is the use of pyrethroid insecticide-treated nets.

Faecal microbiota and cytokine profiles of rural Cambodian infants linked to diet and diarrhoeal episodes.

The gut microbiota of infants in low- to middle-income countries is underrepresented in microbiome research. This study explored the faecal microbiota composition and faecal cytokine profiles in a cohort of infants in a rural province of Cambodia and investigated the impact of sample storage conditions and infant environment on microbiota composition. Faecal samples collected at three time points from 32 infants were analysed for microbiota composition using 16S rRNA amplicon sequencing and concentrations of faecal cytokines. Faecal bacterial isolates were subjected to whole genome sequencing and genomic analysis. We compared the effects of two sample collection methods due to the challenges of faecal sample collection in a rural location. Storage of faecal samples in a DNA preservation solution preserved Bacteroides abundance. Microbiota analysis of preserved samples showed that Bifidobacterium was the most abundant genus with Bifidobacterium longum the most abundant species, with higher abundance in breast-fed infants. Most infants had detectable pathogenic taxa, with Shigella and Klebsiella more abundant in infants with recent diarrhoeal illness. Neither antibiotics nor infant growth were associated with gut microbiota composition. Genomic analysis of isolates showed gene clusters encoding the ability to digest human milk oligosaccharides in B. longum and B. breve isolates. Antibiotic-resistant genes were present in both potentially pathogenic species and in Bifidobacterium. Faecal concentrations of Interlukin-1alpha and vascular endothelial growth factor were higher in breast-fed infants. This study provides insights into an underrepresented population of rural Cambodian infants, showing pathogen exposure and breastfeeding impact gut microbiota composition and faecal immune profiles.

Whole genome sequence analysis of population structure and insecticide resistance markers in Anopheles melas from the Bijagós Archipelago, Guinea-Bissau.

BACKGROUND: Anopheles melas is an understudied malaria vector with a potential role in malaria transmission on the Bijagós Archipelago of Guinea-Bissau. This study presents the first whole-genome sequencing and population genetic analysis for this species from the Bijagós. To our knowledge, this also represents the largest population genetic analysis using WGS data from non-pooled An. melas mosquitoes. METHODS: WGS was conducted for 30 individual An. melas collected during the peak malaria transmission season in 2019 from six different islands on the Bijagós Archipelago. Bioinformatics tools were used to investigate the population structure and prevalence of insecticide resistance markers in this mosquito population. RESULTS: Insecticide resistance mutations associated with pyrethroid resistance in Anopheles gambiae s.s. from the Bijagós were absent in the An. melas population, and no signatures of selective sweeps were identified in insecticide resistance-associated genes. Analysis of structural variants identified a large duplication encompassing the cytochrome-P450 gene cyp9k1. Phylogenetic analysis using publicly available mitochondrial genomes indicated that An. melas from the Bijagós split into two phylogenetic groups because of differentiation on the mitochondrial genome attributed to the cytochrome C oxidase subunits COX I and COX II and the NADH dehydrogenase subunits 1, 4, 4L and 5. CONCLUSIONS: This study identified an absence of insecticide-resistant SNPs common to An. gambiae in the An. melas population, but did identify structural variation over insecticide resistance-associated genes. Furthermore, this study presents novel insights into the population structure of this malaria vector using WGS analysis. Additional studies are required to further understand the role of this vector in malaria transmission.

Genetic diversity of Plasmodium vivax isolates from pregnant women in the Western Brazilian Amazon: a prospective cohort study.

Background: Each year, 92 million pregnant women are at risk of contracting malaria during pregnancy, with the underestimation of the mortality and morbidity burden associated with Plasmodium vivax. During pregnancy, P. vivax infection is associated with low birth weight, maternal anaemia, premature delivery, and stillbirth. In the State of Acre (Brazil), high transmission leaves pregnant women at greater risk of contracting malaria and having a greater number of recurrences. The study of genetic diversity and the association of haplotypes with adverse pregnancy effects is of great importance for the control of the disease. Here we investigate the genetic diversity of P. vivax parasites infecting pregnant women across their pregnancies. Methods: P. vivax DNA was extracted from 330 samples from 177 women followed during pregnancy, collected in the State of Acre, Brazil. All samples were negative for Plasmodium falciparum DNA. Sequence data for the Pvmsp1 gene was analysed alongside data from six microsatellite (MS) markers. Allelic frequencies, haplotype frequencies, expected heterozygosity (HE) were calculated. Whole genome sequencing (WGS) was conducted on four samples from pregnant women and phylogenetic analysis performed with other samples from South American regions. Findings: Initially, the pregnant women were stratified into two groups-1 recurrence and 2 or more recurrences-in which no differences were observed in clinical gestational outcomes or in placental histological changes between the two groups. Then we evaluated the parasites genetically. An average of 18.5 distinct alleles were found at each of the MS loci, and the HE calculated for each marker indicates a high genetic diversity occurring within the population. There was a high percentage of polyclonal infections (61.7%, 108/175), and one haplotype (H1) occurred frequently (20%), with only 9 of the haplotypes appearing in more than one patient. Interpretation: Most pregnant women had polyclonal infections that could be the result of relapses and/or re-infections. The high percentage of H1 parasites, along with the low frequency of many other haplotypes are suggestive of a clonal expansion. Phylogenetic analysis shows that P. vivax population within pregnant women clustered with other Brazilian samples in the region. Funding: FAPESP and CNPq - Brazil.

Detection of insecticide resistance markers in Anopheles funestus from the Democratic Republic of the Congo using a targeted amplicon sequencing panel.

Vector control strategies have been successful in reducing the number of malaria cases and deaths globally, but the spread of insecticide resistance represents a significant threat to disease control. Insecticide resistance has been reported across Anopheles (An.) vector populations, including species within the An. funestus group. These mosquitoes are responsible for intense malaria transmission across sub-Saharan Africa, including in the Democratic Republic of the Congo (DRC), a country contributing > 12% of global malaria infections and mortality events. To support the continuous efficacy of vector control strategies, it is essential to monitor insecticide resistance using molecular surveillance tools. In this study, we developed an amplicon sequencing ("Amp-seq") approach targeting An. funestus, and using multiplex PCR, dual index barcoding, and next-generation sequencing for high throughput and low-cost applications. Using our Amp-seq approach, we screened 80 An. funestus field isolates from the DRC across a panel of nine genes with mutations linked to insecticide resistance (ace-1, CYP6P4, CYP6P9a, GSTe2, vgsc, and rdl) and mosquito speciation (cox-1, mtND5, and ITS2). Amongst the 18 non-synonymous mutations detected, was N485I, in the ace-1 gene associated with carbamate resistance. Overall, our panel represents an extendable and much-needed method for the molecular surveillance of insecticide resistance in An. funestus populations.

Towards environmental detection, quantification, and molecular characterization of Anopheles stephensi and Aedes aegypti from experimental larval breeding sites.

The invasion and establishment of An. stephensi mosquitoes in the Horn of Africa represents a significant regional threat, which may jeopardise malaria control, particularly in urban areas which were formally free from disease transmission. Novel vector surveillance methods are urgently needed, both agnostic to mosquito larval morphology, and simple to implement at the sampling stage. Using new multiplex TaqMan assays, specifically targeting An. stephensi and Ae. aegypti, we validated the use of environmental DNA (eDNA) for simultaneous vector detection in shared artificial breeding sites. Study findings demonstrated that An. stephensi and Ae. aegypti eDNA deposited by as few as one second instar larva in 1L of water was detectable. Characterization of molecular insecticide resistance mechanisms, using novel amplicon-sequencing panels for both vector species, was possible from eDNA shed by as few as 16-32 s instar larvae in 50 ml of water. An. stephensi eDNA, derived from emergent pupae for 24 h, was remarkably stable, and still detectable ~ 2 weeks later. eDNA surveillance has the potential to be implemented in local endemic communities and at points of country entry, to monitor the spread of invasive vector species. Further studies are required to validate the feasibility of this technique under field conditions.

Identification of two insecticide resistance markers in Ethiopian Anopheles stephensi mosquitoes using a multiplex amplicon sequencing assay.

Since its first detection in 2012 in Djibouti, Anopheles stephensi has invaded and established in the Horn of Africa, and more recently Nigeria. The expansion of this vector poses a significant threat to malaria control and elimination efforts. Integrated vector management is the primary strategy used to interrupt disease transmission; however, growing insecticide resistance is threatening to reverse gains in global malaria control. We present a next-generation amplicon-sequencing approach, for high-throughput monitoring of insecticide resistance genes (ace1, GSTe2, vgsc and rdl), species identification and characterization of genetic diversity (its2 and cox1) in An. stephensi. Ninety-five An. stephensi mosquitoes, collected in Ethiopia, were screened, identifying 104 SNPs, including the knock-down mutation L958F (L1014F in Musca domestica), and for the first time in this vector species, the A296S substitution (A301S in Drosophila melanogaster) in the rdl locus. Two other amino acid substitutions (ace1-N177D, GSTe2-V189L) were also identified but have not been previously implicated in insecticide resistance. Genetic diversity in the mitochondrial cox1 gene revealed shared haplotypes between Ethiopian An. stephensi with samples from Pakistan, Sudan, and Djibouti. Overall, we present a reliable, cost-effective strategy using amplicon-sequencing to monitor known insecticide resistance mutations, with the potential to identify new genetic variants, to assist in the high-throughput surveillance of insecticide resistance in An. stephensi populations.

A next generation targeted amplicon sequencing method to screen for insecticide resistance mutations in Aedes aegypti populations reveals a rdl mutation in mosquitoes from Cabo Verde.

Aedes mosquito vectors transmit many viruses of global health concern, including dengue, chikungunya and Zika. These vector-borne viral diseases have a limited number of treatment options, and vaccines vary in their effectiveness. Consequently, integrated vector management is a primary strategy for disease control. However, the increasing emergence and spread of insecticide resistance is threatening the efficacy of vector control methods. Identifying mutations associated with resistance in vector populations is important to monitor the occurrence and evolution of insecticide resistance and inform control strategies. Rapid and cost-effective genome sequencing approaches are urgently needed. Here we present an adaptable targeted amplicon approach for cost-effective implementation within next generation sequencing platforms. This approach can identify single nucleotide polymorphisms (SNPs) and small insertions and deletions (indels) in genes involved in insecticide resistance in Aedes aegypti mosquitoes. We designed and tested eleven amplicons, which included segments of the ace-1 (carbamate target), the Voltage-Gated Sodium Channel (vgsc; pyrethroids, DDT and organochlorines), and rdl (dieldrin) genes; thereby covering established knockdown resistance (kdr) mutations (e.g., S989P, I1011M/V, V1016G/I and F1534C), with the potential to identify novel ones. The amplicon assays were designed with internal barcodes, to facilitate multiplexing of large numbers of mosquitoes at low cost, and were sequenced using an Illumina platform. Our approach was evaluated on 152 Ae. aegypti mosquitoes collected in Cabo Verde, an archipelago with a history of arbovirus outbreaks. The amplicon sequence data revealed 146 SNPs, including four non-synonymous polymorphisms in the vgsc gene, one in ace-1 and the 296S rdl mutation previously associated with resistance to organochlorines. The 296S rdl mutation was identified in 98% of mosquitoes screened, consistent with the past use of an organochlorine compound (e.g., DDT). Overall, our work shows that targeted amplicon sequencing is a rapid, robust, and cost-effective tool that can be used to perform high throughput monitoring of insecticide resistance.

Rapid MinION profiling of preterm microbiota and antimicrobial-resistant pathogens.

The MinION sequencing platform offers near real-time analysis of DNA sequence; this makes the tool attractive for deployment in fieldwork or clinical settings. We used the MinION platform coupled to the NanoOK RT software package to perform shotgun metagenomic sequencing and profile mock communities and faecal samples from healthy and ill preterm infants. Using Nanopore data, we reliably classified a 20-species mock community and captured the diversity of the immature gut microbiota over time and in response to interventions such as probiotic supplementation, antibiotic treatment or episodes of suspected sepsis. We also performed rapid real-time runs to assess gut-associated microbial communities in critically ill and healthy infants, facilitated by NanoOK RT software package, which analysed sequences as they were generated. Our pipeline reliably identified pathogenic bacteria (that is, Klebsiella pneumoniae and Enterobacter cloacae) and their corresponding antimicrobial resistance gene profiles within as little as 1 h of sequencing. Results were confirmed using pathogen isolation, whole-genome sequencing and antibiotic susceptibility testing, as well as mock communities and clinical samples with known antimicrobial resistance genes. Our results demonstrate that MinION (including cost-effective Flongle flow cells) with NanoOK RT can process metagenomic samples to a rich dataset in < 5 h, which creates a platform for future studies aimed at developing these tools and approaches in clinical settings with a focus on providing tailored patient antimicrobial treatment options.