Abstract Wolbachia , a widespread bacterium that can reduce pathogen transmission in mosquitoes, has been detected within populations of Anopheles (An.) malaria vectors. In the An. gambiae complex, the primary vectors in Sub-Saharan Africa, Wolbachia strains are at low density and infection frequencies in wild populations. PCR-independent evidence is required to determine whether Wolbachia strains are true endosymbionts in Anopheles given most studies to date have used nested-PCR to identify strains. Here we report high-density strains found in geographically diverse populations of An. moucheti and An . demeilloni . Fluorescent in situ hybridization localized a heavy infection in the ovaries of An. moucheti and maternal transmission was observed. Genome sequencing of both strains obtained genome depths and coverages comparable to other known infections. Notably, homologs of cytoplasmic incompatibility factor ( cif ) genes were present indicating these strains possess the capacity to induce the phenotype cytoplasmic incompatibility which allows Wolbachia to spread through populations. The characteristics of these two strains suggest they are ideal candidates for Wolbachia biocontrol strategies in Anopheles .

Abstract Wolbachia is a genus of obligate bacterial endosymbionts that infect a diverse range of arthropod species as well as filarial nematodes, with its single described species, Wolbachia pipientis, divided into several ‘supergroups’ based on multilocus sequence typing. Wolbachia strains in mosquitoes have been shown to inhibit the transmission of human pathogens including Plasmodium malaria parasites and arboviruses. Despite their large host range, Wolbachia strains within the major malaria vectors of the Anopheles (A.) gambiae and A. funestus complexes appear at low density based solely on PCR-based methods. Questions have been raised as to whether this represents a true endosymbiotic relationship. However, recent definitive evidence for two distinct, high-density strains of supergroup B Wolbachia within A. demeilloni and A. moucheti has opened exciting possibilities to explore naturally occurring Wolbachia endosymbionts in Anopheles for biocontrol strategies to block Plasmodium transmission. Here we utilise genomic analyses to demonstrate that both Wolbachia strains have retained all key metabolic and transport pathways despite their smaller genome size. We further confirm the presence of cytoplasmic incompatibility factor genes, despite noticeably few prophage regions. Additionally, phylogenetic analysis indicates that these Wolbachia strains may have been introduced into these two Anopheles species via horizontal transmission events, and unlikely to be by ancestral acquisition and subsequent loss events in the Anopheles gambiae species complex. These are the first Wolbachia genomes that enable us to study the relationship between natural strains Plasmodium malaria parasites and their Anopheline hosts. Impact statement Wolbachia naturally infects a wide range of arthropod species, including insect vectors of human pathogens, where they may play a role in inhibiting their replication. These bacteria have been commonly found within Aedes (Ae.) albopictus and Culex pipiens mosquitoes but have been noticeably absent in the Anopheles mosquito genera, which includes all species responsible for malaria transmission. Recent PCR-based methods have suggested the potential for natural Wolbachia strains within the A. gambiae species complex, which includes major malaria vector species including A. gambiae s.s., A. coluzzii and A. arabiensis . We recently reported the presence of stable Wolbachia strains naturally occurring within two different Anopheles species ( A. demeilloni and A. moucheti) . In this study, we perform comparative genomic analysis of these two Wolbachia genomes against each other and published Wolbachia strains. The current assemblies are some of the smallest sequenced Wolbachia strains of insects, although their metabolic pathway repertoire is comparable to other strains. Interestingly, prophage fragments were identified within only one of the two strains. The findings of this study will be of significant interest to researchers investigating Wolbachia as a potential malaria biocontrol strategy, giving greater insight into the evolution and diversity of this obligate intracellular endosymbiont. Data summary Sequence data generated and used for this analysis are available in the National Centre for Biotechnology Information Sequence Read Archive (NCBI SRA bioproject number PRJNA642000). The two assembled Wolbachia genomes are available with genome accession numbers GCA_018491735.2 and GCA_018491625.2. Additional Wolbachia genomes used for comparative analysis are described in the supplementary material. The authors confirm all supporting data, code and protocols have been provided within the article or through supplementary data files. Additional supplementary data files used to generate several figures can be found at: https://figshare.com/projects/Wolbachia_endosymbionts_in_two_Anopheles_species_indicates_independent_acquisitions_and_lack_of_prophage_elements/126533

Abstract Background Resistance to major public health insecticides in Côte d’Ivoire has intensified and now threatens the long-term effectiveness of malaria vector control interventions. Methods This study evaluated the bioefficacy of conventional and next-generation long-lasting insecticidal nets (LLINs), determined resistance profiles, and characterized molecular and metabolic mechanisms in wild Anopheles coluzzii from South-East Côte d’Ivoire in 2019. Results Phenotypic resistance was intense: more than 25% of mosquitoes survived exposure to ten times the doses of pyrethroids required to kill susceptible populations. Similarly, 24-hour mortality to deltamethrin-only LLINs was very low and not significantly different to an untreated net. Sub-lethal pyrethroid exposure did not induce significant delayed vector mortality 72 hours later. In contrast, LLINs containing the synergist piperonyl butoxide (PBO), or new insecticides, clothianidin and chlorfenapyr, were highly toxic to An. coluzzii . Pyrethroid-susceptible An. coluzzii were significantly more likely to be infected with malaria, compared to those that survived insecticidal exposure. Pyrethroid resistance was associated with significant over-expression of CYP6P4, CPY6Z1 and CYP6P3 . Conclusions Study findings raise concerns regarding the operational failure of standard LLINs and support the urgent deployment of vector control interventions incorporating PBO, chlorfenapyr or clothianidin in areas of high resistance intensity in Côte d’Ivoire.

Abstract Background West Nile virus (WNV) is a zoonotic arbovirus of great medical and veterinary importance, threatening the health of humans and equines worldwide. Mosquitoes belonging to the Culex ( Cx .) pipiens complex are major vectors but numerous other mosquito species have also been implicated as vectors of WNV. Due to variations in blood-feeding behaviour, the different biotypes and hybrids of Cx. pipiens influence the transmission of WNV, from enzootic cycles (between mosquitoes and birds), to spill-over transmission to humans and equines. Methods In this study, mosquitoes were collected and analysed from two regional units (RUs) of Greece with reported cases of WNV within the past 4 years; Palaio Flairo and Argolida (in Attica and Peloponnese regions, respectively). Collections using different types of mosquito surveillance traps were undertaken in May-June 2018 during the early period of the WNV transmission season. Results A total of 1062 mosquitoes were collected, with Biogents Sentinel traps (BG traps) collecting both a greater number of mosquitoes across all species and Cx. pipiens complex individuals than Centres for Disease Control miniature light traps (CDC traps) or Heavy Duty Encephalitis Vector Survey traps (EVS traps). Identification of collected mosquitoes (using both morphological keys and molecular barcoding) confirmed the presence of additional species including Aedes (Ae.) albopictus, Ae. caspius and Culiseta (Cs.) longiareolata . The prevalence of Cx. pipiens biotypes in the RU of Palaio Faliro was 54.5% pipiens type, 20.0% molestus type and 25.5% hybrids. In the RU of Argolida, the collection comprised 68.1% pipiens type, 8.3% molestus type and 23.6% hybrids. Screening individual unfed female mosquitoes for WNV (molecular xenomonitoring) resulted in detection in three females of the pipiens type and in one hybrid; all collected from the RU of Argolida. Conclusions As hybrids play an important role in spill-over transmission of WNV to humans and equines, these findings highlight the importance of undertaking entomological surveillance programs incorporating molecular xenomonitoring at the onset of the transmission season to provide an early warning system for health authorities aiming to prevent WNV outbreaks in Greece.

Abstract Background Insecticide resistance among mosquito species is now a pervasive phenomenon, which threatens to jeopardise global malaria vector control efforts. Evidence of links between the mosquito microbiota and insecticide resistance is emerging, with significant enrichment of insecticide degrading bacteria and enzymes in resistant populations. Using 16S rRNA amplicon sequencing, we characterised and compared the microbiota of Anopheles ( An. ) coluzzii in relation to their deltamethrin resistance and exposure profiles. Results Comparisons between 2-3 day old deltamethrin resistant and susceptible mosquitoes, demonstrated significant differences in microbiota diversity (PERMANOVA, pseudo-F = 19.44, p=0.0015). Ochrobactrum, Lysinibacillus and Stenotrophomonas genera, each of which comprised insecticide degrading species, were significantly enriched in resistant mosquitoes. Susceptible mosquitoes had a significant reduction in alpha diversity compared to resistant individuals (Shannon index: H=13.91, q=0.0003, Faith’s phylogenetic diversity: H=6.68, q=0.01), with Asaia and Serratia dominating microbial profiles. There was no significant difference in deltamethrin exposed and unexposed 5-6 day old individuals, suggesting that insecticide exposure had minimal impact on microbial composition. Serratia and Asaia were also dominant in 5-6 day old mosquitoes, regardless of exposure or phenotype, and had reduced microbial diversity compared with 2-3 day old mosquitoes. Conclusions Our findings revealed significant alterations of An. coluzzii microbiota associated with deltamethrin resistance, highlighting the potential for identification of novel microbial markers for insecticide resistance surveillance. qPCR detection of Serratia and Asaia was consistent with 16S rRNA sequencing, suggesting that population level field screening of the bacterial microbiota may be feasibly integrated into wider resistance monitoring if reliable and reproducible markers associated with phenotype can be identified.

Abstract Culex (Cx.) tritaeniorhynchus is a mosquito species with an extensive and expanding inter-continental geographic distribution, currently reported in over 50 countries, across Asia, Africa, the Middle East, Europe and now Australia. It is an important vector of medical and veterinary concern, capable of transmitting multiple arboviruses which cause significant morbidity and mortality in human and animal populations. In regions endemic for Japanese encephalitis virus (JEV) in Asia, Cx. tritaeniorhynchus is considered the major vector and this species has also been shown to contribute to the transmission of several other significant zoonotic arboviruses, including Rift Valley fever virus and West Nile virus. Significant variation in vectorial capacity can occur between different vector populations. Obtaining knowledge of a species from across its geographic range is crucial to understanding its significance for pathogen transmission across diverse environments and localities. Vectorial capacity can be influenced by factors including the mosquito genetic background, composition of the microbiota associated with the mosquito and the co-infection of human or animal pathogens. In addition to enhancing information on vector surveillance and potential risks for pathogen transmission, determining the genetic and microbial diversity of distinct populations of a vector species is also critical for the development and application of effective control strategies. In this study, multiple geographically dispersed populations of Cx. tritaeniorhynchus from countries within Europe, Africa, Eurasia and Asia were sampled. Molecular analysis demonstrated a high level of genetic and microbial diversity within and between populations, including genetic divergence in the mosquito CO1 gene, as well as diverse microbiomes identified by 16S rRNA gene amplicon sequencing. Evidence for the detection of the endosymbiotic bacteria Wolbachia in some populations was confirmed using Wolbachia -specific PCR detection and sequencing of Wolbachia MLST genes; in addition to PCR-based detection of insect-specific viruses. Laboratory vector competence showed Cx. tritaeniorhynchus from a Greek population are likely to be competent vectors of JEV. This study expands understanding of the diversity of Cx. tritaeniorhynchus across its inter-continental range, highlights the need for a greater focus on this invasive vector species and helps to inform potential future directions for development of vector control strategies.

Wolbachia endosymbiotic bacteria are widespread throughout insect species and Wolbachia transinfected in Aedes mosquito species has formed the basis for biocontrol programs as Wolbachia strains inhibit arboviral replication and can spread through populations. Resident strains in wild Culicine mosquito populations (the vectors of most arboviruses) requires further investigation given resident strains can also affect arboviral transmission. As Madagascar has a large diversity of both Culicine species and has had recent arboviral outbreaks, an entomology survey was undertaken, in five ecologically diverse sites, to determine the Wolbachia prevalence. We detected diverse novel resident Wolbachia strains within the Aedeomyia, Culex, Ficalbia, Mansonia and Uranotaenia genera. Wolbachia prevalence rates and strain characterisation through Sanger sequencing with multilocus sequence typing (MLST) and phylogenetic analysis revealed significant diversity and we detected co-infections with the environmentally acquired bacterial endosymbiont Asaia. Mosquitoes were screened for major arboviruses to investigate if any evidence could be provided for their potential role in transmission and we report the presence of Rift Valley fever virus in three Culex species: Culex tritaeniorhynchus, Culex antennatus and Culex decens. The implications of the presence of resident Wolbachia strains are discussed and how the discovery of novel strains can be utilized for applications in the development of biocontrol strategies.

Outbreaks of mosquito-borne arboviral diseases including dengue virus (DENV), Zika virus (ZIKV), yellow fever virus (YFV) and chikungunya virus (CHIKV) have recently occurred in the Caribbean. The geographical range of the principle vectors responsible for transmission, Aedes (Ae.) aegypti and Ae. albopictus is increasing and greater mosquito surveillance is needed in the Caribbean given international tourism is so prominent. The island of Saint Lucia has seen outbreaks of DENV and CHIKV in the past five years but vector surveillance has been limited with the last studies dating back to the late 1970s. Natural disasters have changed the landscape of Saint Lucia and the island has gone through significant urbanisation. In this study, we conducted an entomological survey of Ae. aegypti and Ae. albopictus distribution across the island and analysed environmental parameters associated with the presence of these species. Although we collected Ae. aegypti across a range of sites across the island, no Ae. albopictus were collected despite traps being placed in diverse ecological settings. The number of Ae. aegypti collected was significantly associated with higher elevation and semi-urban settings yielded female mosquito counts per trap-day that were 5-fold lower than urban settings. Screening for arboviruses revealed a high prevalence of a novel insect-specific flavivirus closely related to cell fusing agent virus (CFAV). We discuss the implications that natural disasters, water storage and lack of mosquito surveillance have on arboviral outbreaks in Saint Lucia and implications for insect only flaviviruses on surveillance and detection of pathogenic flaviviruses.

Background: Wolbachia, a widespread bacterium that can influence mosquito-borne pathogen transmission, has recently been shown to infect Anopheles (An.) species that are malaria vectors in Sub-Saharan Africa. Although there are studies reporting strains in the An. gambiae complex, the apparent low density and low prevalence rates requires further investigation. In contrast, strains in other species appear higher density allowing a greater understanding of phylogenetics and more accurate determination of prevalence rates in wild mosquito populations. Methods: Anopheles mosquitoes were collected in the Faranah and Maferinyah regions of Guinea in June-July 2018. RNA was extracted from 542 females and reverse transcribed to determine Wolbachia prevalence rates, demonstrate gene expression and estimate relative strain densities using quantitative PCR. Molecular confirmation of mosquito species and Wolbachia multilocus sequence typing (MLST) was carried out to analyse phylogenetic relationships of newly discovered strains. Results: Low prevalence rates were detected in An. gambiae s.s. (0.0: 2.9%) and novel Wolbachia strains were discovered in two species: the wAnMe strain in An. melas (18/168: prevalence rate of 10.7%) and the wAnsX strain in a previously unidentified Anopheles species we have termed An. species X (1/1). Novel strains were phylogenetically diverse, with wAnMe clustering with Wolbachia Supergroup A strains, and wAnsX closest to Supergroup B strains. Significantly higher density strains were present in An. species X and An. melas compared to An. gambiae s.s. and were comparable in density to a novel Wolbachia strain discovered in Culex watti termed wWat. Conclusions: The discovery of novel Wolbachia strains provides further insight into the strain diversity within Anopheles species. Confirmation of gene expression for all MLST and wsp genes for wAnsX and wAnMe strains and qPCR analysis indicates these are candidate strains for transinfection to create stable infections in other Anopheles species. They have the potential for use in Anopheles biocontrol strategies either through population replacement or suppression control strategies.

Background: In recent years, the scale-up of long-lasting insecticidal nets (LLINs) and indoor residual spraying (IRS) has greatly reduced malaria transmission. However, malaria remains a global public health concern with the majority of disease burden in sub-Saharan Africa. Insecticide resistance is a growing problem among Anopheles vector populations, with potential implications for the continued effectiveness of available control interventions. Improved understanding of current resistance levels and underlying mechanisms is essential to design appropriate management strategies and to mitigate future selection for resistance. Methods: Anopheles gambiae s.l. mosquitoes were collected from three villages in Faranah Prefecture, Guinea and their levels of susceptibility to seven insecticides were measured using CDC resistance intensity bioassays. Synergist assays with piperonyl butoxide (PBO) were also undertaken to assess the role of elevated mixed-function oxidases in resistance. RNA was extracted from 563 individuals and PCR was performed on cDNA to determine vector species, presence of target site mutations (L1014F kdr, N1575Y and G119S Ace-1), Plasmodium falciparum infection, and relative expression of three metabolic genes (CYP6M2, CYP6P3 and GSTD3). Results: In Faranah, resistance to permethrin and deltamethrin was observed, as well as possible resistance to bendiocarb. All assayed vector populations were fully susceptible to alpha-cypermethrin, pirimiphos-methyl, clothianidin and chlorfenapyr. Plasmodium falciparum infection was detected in 7.3% (37/508) mosquitoes tested. The L1014F kdr mutation was found in 100% of a sub-sample of 60 mosquitoes, supporting its fixation in the region. The N1575Y mutation was identified in 20% (113/561) of individuals, with ongoing selection evidenced by significant deviations from Hardy-Weinberg equilibrium. The G119S Ace-1 mutation was detected in 62.1% (18/29) of mosquitoes tested and was highly predictive of bendiocarb bioassay survival. The metabolic resistance genes, CYP6M2, CYP6P3 and GSTD3, were found to be overexpressed in wild resistant and susceptible An. gambiae s.s. populations, compared to a susceptible G3 colony. Furthermore, CYP6P3 was significantly overexpressed in bendiocarb survivors, implicating its potential role in carbamate resistance in Faranah. Conclusions: Identification of intense resistance to permethrin and deltamethrin in Faranah, is of concern, as the Guinea National Malaria Control Program (NMCP) relies exclusively on the distribution of pyrethroid-treated LLINs for vector control. Study findings will be used to guide current and future control strategies in the region.