
==== Front
Sci Rep
Sci Rep
Scientific Reports
2045-2322
Nature Publishing Group UK London

39232089
71507
10.1038/s41598-024-71507-y
Article
Characterization of the microbiome of Aedes albopictus populations in different habitats from Spain and São Tomé
Melo Tiago 1
Sousa Carla Alexandra 1
Delacour-Estrella Sarah 23
Bravo-Barriga Daniel 4
Seixas Gonçalo gseixas@ihmt.unl.pt

1
1 https://ror.org/02xankh89 grid.10772.33 0000 0001 2151 1713 Global Health and Tropical Medicine, GHTM, Associate Laboratory in Translation and Innovation Towards Global Health, LA-REAL, Instituto de Higiene e Medicina Tropical, Universidade Nova de Lisboa, UNL, Rua da Junqueira, 100, 1349-008 Lisboa, Portugal
2 https://ror.org/012a91z28 grid.11205.37 0000 0001 2152 8769 Animal Health Department, The AgriFood Institute of Aragon (IA2), School of Veterinary Medicine, University of Zaragoza, 50013 Zaragoza, Spain
3 Departamento de Investigación y Desarrollo (I+D), Quimera. B.S. Calle Olivo, 14, 50016 La Puebla de Alfindén, Spain
4 https://ror.org/05yc77b46 grid.411901.c 0000 0001 2183 9102 Departamento de Sanidad Animal, Grupo de Investigación en Salud Animal y Zoonosis (GISAZ), UIC Zoonosis y Enfermedades Emergentes (ENZOEM), Universidad de Córdoba, Córdoba, Spain
4 9 2024
4 9 2024
2024
14 205452 4 2024
28 8 2024
© The Author(s) 2024
2024
https://creativecommons.org/licenses/by-nc-nd/4.0/ Open Access This article is licensed under a Creative Commons Attribution-NonCommercial-NoDerivatives 4.0 International License, which permits any non-commercial use, sharing, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if you modified the licensed material. You do not have permission under this licence to share adapted material derived from this article or parts of it. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by-nc-nd/4.0/.
The mosquito microbiome significantly influences vector competence, including in Aedes albopictus, a globally invasive vector. Describing the microbiome and Wolbachia strains of Ae. albopictus from different regions can guide area-specific control strategies. Mosquito samples from Spain and São Tomé were analyzed using 16S rRNA gene sequencing and metagenomic sequencing. Wolbachia infection patterns were observed by sex and population. Female mosquitoes were blood-fed, a factor considered in analyzing their microbiota. Results revealed a dominance of dual Wolbachia infections, strains A and B, in the microbiome of both populations of Ae. albopictus, especially among females. Both populations shared a core microbiome, although 5 and 9 other genera were only present in Spain and São Tomé populations, respectively. Genera like Pelomonas and Nevskia were identified for the first time in Aedes mosquitoes. This study is the first to describe the Ae. albopictus bacteriome in Spain and São Tomé, offering insights for the development of targeted mosquito control strategies. Understanding the specific microbiome composition can help in designing more effective interventions, such as microbiome manipulation and Wolbachia-based approaches, to reduce vector competence and transmission potential of these mosquitoes.

Keywords

Microbiome
Wolbachia
Vector competence
Mosquito control
Metagenomics
Aedes albopictus
Subject terms

Biological techniques
Computational biology and bioinformatics
Ecology
Microbiology
Molecular biology
issue-copyright-statement© Springer Nature Limited 2024
==== Body
pmcIntroduction

Vector-borne diseases, which account for 20% of infectious diseases worldwide, pose a significant global health threat. Illnesses like dengue fever, yellow fever, chikungunya, and Zika put nearly 3.9 billion people1 at risk, and effective management and control of these disease vectors are crucial for reducing these risks of transmission.

Aedes albopictus, a mosquito species known for its ability to transmit diseases and for being a nuisance to humans2, has a worldwide distribution3 and is responsible for spreading dengue, chikungunya, and Zika viruses4–6. Its invasive nature7 and adaptability to different environments8 make it a significant threat to human health. Originally from Southeast Asia, Ae. albopictus has expanded its range to other continents, facilitated by global trade of used tires, ornamental plants, and road traffic4,9,10. Moreover, in recent decades, rapid urbanization has contributed to proliferation of Ae. albopictus populations, and this, along with inadequate or absent vector control measures, has further increased the risk of pathogen transmission by this species11. The rapid spread of Ae. albopictus mosquitoes in São Tomé and Príncipe12 and Spain13 underscores the need to understand and manage this species in these regions. São Tomé and Principe, an island nation in the Gulf of Guinea, is particularly vulnerable to vector-borne diseases due to its tropical climate and limited healthcare infrastructure but colonization events should be rare14. Aedes albopictus was first identified in São Tomé in 2016, although it is probable that the species had inhabited the island for a decade prior to this discovery12. In Spain, the first Ae. albopictus was detected in 2004, at Sant Cugat del Vallès15. Since then, the mosquito’s presence has proliferated across the country13,16–18 due to multiple introductions from abroad19. This expansion has escalated concerns among public health officials, emphasizing the need for increased surveillance and control measures20–22.

The host physiology of Aedes mosquitoes is heavily influenced by their microbiota, which has been shown to impact various aspects of reproduction, egg production, blood digestion23–25, regulation of immunity26, genetic diversity, vector competence, host–pathogen interactions and survival19,27,28. For instance, some bacteria can produce antimicrobial peptides that boost the mosquito’s immune system, while other might facilitate the digestion of blood meals. Additionally, the presence of certain microbial communities can influence the mosquito’s ability to transmit pathogens by either supporting or hindering the development of the pathogens within the mosquito. These factors collectively contribute to the vectorial capacity of mosquitoes and are key for pathogen transmission. Innovative vector control strategies have increasingly focused on manipulating the mosquito microbiome to reduce their capacity for disease transmission29 and control of populations. One promising approach involves introducing exogenous Wolbachia strains into Ae. albopictus populations, which has shown potential in disrupting pathogen transmission30. However, most of the Ae. albopictus populations are already naturally infected with this endosymbiont31 and different Wolbachia strains may exhibit various interactions when coexisting within a host30. Describing the native strains of our target-populations is the first step to assess possible interactions between the naturally occurring Wolbachia strains and the introduced type.

The genetic background of any given species plays a role in shaping its microbiome32. Different populations of mosquitoes may have genetic variations that affect their immune responses, susceptibility to certain microorganisms, and overall interactions with the microbiome. In the same way, different geographical locations may have unique microbial communities in soil, water, and vegetation, and the local climate and temperature can impact both the mosquito and the microorganisms it harbors33. Given the significance of understanding and managing vector-borne diseases, this study aimed to examine and to compare the microbiome composition and characterize the presence of Wolbachia strains in populations of Ae. albopictus from two countries with distinct climatic, habitat, epidemiological, and geographical conditions: São Tomé and Príncipe (Africa) and Spain (Europe). Investigation on the composition of the microbiome and Wolbachia strains in Ae. albopictus populations from these distinct geographical regions could provide valuable knowledge into potential control measures tailored to the specific problems of each area. By understanding the microbiome composition and the prevalence of different Wolbachia strains, researchers can develop targeted strategies, such as microbiome manipulation or Wolbachia-based interventions, to reduce the vector competence and transmission potential of these mosquitoes. This information could help develop new vector control strategies and contribute to reducing the global burden of vector-borne diseases.

Materials and methods

Mosquito collection and rearing

Eggs were collected using sixteen ovitraps installed in three distant regions of Spain: Valencia (Paterna), Huesca (Monzón), and Málaga (Rincón de La Victoria), as well as one region in São Tomé, Água Grande (São Tomé). This collection took place between late 2021 (September and October) and early 2022 (April), with the assistance of collaborators, as shown in Fig. 1. The ovitraps were designed as black plastic cups filled with 500 mL of dechlorinated water and lined with a strip of germination paper to provide a substrate for the mosquitoes to lay their eggs. The ovitraps were baited with a small amount of hay infusion (prepared by soaking 50 g of dried hay in 10 L of water for 3 days) to attract gravid female mosquitoes. The ovitraps were strategically placed in shaded and sheltered locations near human habitation and vegetation, ensuring optimal conditions for mosquito egg-laying. The germination papers were collected weekly, and new ones were placed in the ovitraps. The collected papers with eggs were transported to the laboratory for further processing.Fig. 1 Map indicating collection sites for Ae. albopictus. (A) Map of Spain showing the sampling locations in Valencia (Paterna), Huesca (Monzón), and Málaga (Rincón de La Victoria); (B) Regional map indicating the relative positions of Spain and São Tomé, per continent. (C) Map of São Tomé showing the sampling location in Água Grande. The map was generated with QGIS v3.34 software (https://qgis.org/en/site/).

Four Ae. albopictus colonies were subsequently established and maintained in the high-security insectary of the 'In Vivo Arthropod Security Facility' (VIASEF) available at IHMT, between March and May 2022. These colonies were initiated from a parental generation (F0) and maintained until the second generation (F2). Upon hatching in dechlorinated water, larvae were reared under controlled laboratory conditions (temperature: 26 ± 2 °C, relative humidity: 70 ± 5%, photoperiod: 12 h/12 h light/dark) and fed with Tetra Min flakes fish food (Tetra, Melle, Germany). The fish food flakes were processed into powder before being administered to the larvae. Adult mosquitoes were provided with a 10% glucose solution, and females were blood-fed on Mus musculus two to three times a week. The process of handling the animals used occurred under supervision and was carried out based on Community Council standards European Union of 24 November 1986 (86/609/EEC) and national legislation in force (Decree-Law 129/92 of June 2nd, Ordinance No. 100/92 of October 23rd). All animal experiments were based on protocols approved by the Direção-Geral de Veterinária, Ministério da Agricultura do Desenvolvimento Rural e das Pescas, Portugal (ID approvals: 023351 and 023355).

DNA extraction

DNA extractions from a total of 199 adult mosquitoes were carried out using two distinct methods, tailored to the specific requirements of our analyses. For the purpose of sequencing the 16S rRNA gene in a subset of 19 female samples, we employed the "NZY tissue gDNA isolation kit" provided by NZYtech, Portugal. Before DNA extraction, the 19 samples were treated with 10% bleach, followed by 70% alcohol, and finally distilled water to sterilize the mosquitoes. These selected samples originated from adult female mosquitoes belonging to the F0 generation, specifically targeted for our 16S rRNA gene sequencing study to provide a detailed microbiome analysis. The selection of these 19 samples was based on ensuring a representative subset that included mosquitoes from both geographic locations (Spain and São Tomé and Príncipe) and aimed to capture any potential differences in microbiota composition between these regions. The limited number of samples for 16S rRNA gene sequencing was due to resource constraints and the need to focus on a manageable subset for high-throughput sequencing and in-depth microbiome characterization.

For the remaining 180 samples, we utilized a ribosomal DNA extraction protocol that was adapted from the methodology described by Collins et al.34. This approach was applied to both male and female adult mosquitoes from the F0 and F1 generations to facilitate a comprehensive comparison of Wolbachia prevalence and strain distribution across the populations under study.

All DNA extractions were performed on individual mosquitoes, not pooled samples, to ensure accurate representation of each mosquito’s microbiome and Wolbachia infection status.

To ensure the integrity of our DNA extraction process and to eliminate the possibility of contamination, negative controls were systematically included in each batch of extractions.

16S metagenomic sequencing

The bacterial composition of Ae. albopictus microbiota was investigated by sequencing the V3 and V4 regions of the bacterial 16S rRNA gene in 19 samples extracted with the NZY tissue gDNA isolation kit. The sequencing of the amplified products was performed on the Illumina MiSeq platform using v3 chemistry and a 2 × 300 bp paired-end module. Library construction, sequencing, and bioinformatics analysis were performed by Eurofins Genomics Europe Sequencing GmbH (Konstanz, Germany). Briefly, bioinformatics pipeline began with the removal of sequences containing ambiguous bases, followed by the identification and exclusion of chimeric reads utilizing the UCHIME algorithm within the VSEARCH package35,36. The remaining high-quality sequences were then processed through Minimum Entropy Decomposition (MED), which efficiently partitions the marker gene datasets into Operational Taxonomic Units (OTUs)37,38. This method capitalizes on Shannon entropy, selectively considering only the informative nucleotide positions, thus enabling the delineation of sequences based on single nucleotide variances without the interference of random sequencing errors.

After OTU generation, taxonomic assignments were conducted through DC-MEGABLAST alignments against a curated reference database, ensuring the provision of the most specific taxonomic nomenclature to each OTU from the best-matching reference sequences. We applied the QIIME 1.9.1 software package for additional processing of OTUs and their taxonomic classifications, where the abundance data of bacterial taxa were normalized based on lineage-specific gene copy numbers to enhance accuracy39. A minimum sequence identity of 70% over at least 80% of the representative sequence was required to consider the sequence as a reference. Further processing of OTUs and taxonomic assignments was performed using QIIME software (version 1.9.1, http://qiime.org/). The abundances of bacterial taxonomic units were normalized using the 16S rRNA Gene Copy Number (GCN) correction. For this, the read numbers assigned to a species were divided by the known or assumed number of marker gene regions39. The analytical output comprised a comprehensive suite of files, elucidating the taxonomic structure and diversity within the Ae. albopictus microbiota. These outputs included summarized lists of identified taxonomic units per sample, representative sequences for each OTU, and matrices detailing the estimated abundances of OTUs and taxonomic units across samples. Alpha diversity was assessed using Shannon and Simpson diversity indices to evaluate the bacterial diversity within individual samples. Beta diversity was analyzed to compare the bacterial community composition between samples from Spain and São Tomé. These indices were calculated using the QIIME 1.91. software package36.

All sequences have been submitted to NCBI under the accession number for SRA data PRJNA1028981 (https://www.ncbi.nlm.nih.gov/bioproject/PRJNA1028981).

Wolbachia genotyping

PCR genotyping of Wolbachia infection in 180 samples used Wolbachia-specific primers 328Fw and 691Rv for wAlbA, and 183F and 691R for wAlbB40. PCR mix and thermal cycling conditions were standardized. Universal primers wsp 81Fw and 691Rv were used to confirm Wolbachia-negative samples. Primers and PCR conditions details are provided in Supplementary Table S1.

Statistical analysis

Mixed and single infection rates were calculated by sex and population41. Comparisons between mixed and single infections in females and males were performed using the chi-square test. A p-value less than 0.05 was considered statistically significant.

Results

Microbiome composition

To analyze the bacterial communities within Ae. albopictus from Spain and São Tomé, we conducted 16S rRNA gene sequencing of the V3–V4 hypervariable regions across 19 individual adult female samples. After pre-processing and quality control, five samples from São Tomé yielded 606.801 clean sequences, and 14 samples from Spain yielded 798.002 clean sequences (Table 1). All these clean sequences were classified into OTUs (Operational Taxonomic Units): 344 from São Tomé and 517 from Spain (Table 1). The OTUs represented different taxonomic levels, and a correction was made to the species found39 (Supplementary Tables S2 and S3). The Fig. 2 shows shared and non-shared bacterial genera among the populations.Table 1 Summarized statistics of Spain and São Tomé samples.

Sample origin	Total number of input sequences	Remaining sequences after processing	Total number of sequences assigned to OTUs/taxa	Number of OTUs assigned to taxa	
Spain	798,100	798,002	643,902	517	
São Tomé	606,877	606,801	525,776	344	

Fig. 2 Venn diagram showing shared and unshared genera between the Spain and São Tomé populations.

Compared to other studies on the microbiota of mosquitoes from the field, the number of taxa found per sample in our study is relatively low. To ensure that we captured the maximum diversity, rarefaction curves were generated for each sample. These curves, which plot the number of observed OTUs against the number of sequences sampled, indicated that the sequencing depth was sufficient to capture the majority of the bacterial diversity present in our samples (Supplementary Fig. S1). The rarefaction analysis confirmed that our sequencing effort was adequate, as the curves approached an asymptote, suggesting that additional sequencing would likely yield few additional OTUs.

Taxonomic units overview and microbiome diversity analysis

Tables 2 and 3 provide an overview of the taxonomic units identified in samples from São Tomé and Spain. Taxonomic units with readings below 0.1% were categorized as "other." OTU readings analysis showed a predominance of the genus Wolbachia in all samples, ranging between 92.4 and 98.8% in São Tomé samples and between 96.1 and 97.5% in Spanish samples.Table 2 Summary of the taxonomic composition of samples from São Tomé.

Province/district	Sample	Taxonomic unit	Taxonomic level	OTUs reads	Reads	Raw fraction	
Água Grande/São Tomé	STPF0.1	Wolbachia	Genus	191	101,778	98.1%	
Gammaproteobacteria	Class	1	157	0.2%	
Other	–			1.7%	
STPF0.2	Wolbachia	Genus	185	70,626	92.4%	
Asaia sp.	Species	1	1311	1.8%	
Gammaproteobacteria	Class	3	1035	1.4%	
Enterobacteriaceae	Family	1	876	1.2%	
Erwiniaceae	Family	1	242	0.3%	
Enterobacterales	Order	2	195	0.3%	
Acinetobacter	Genus	2	166	0,2%	
Dechloromonas	Genus	1	95	0.1%	
Other	–			2,3%	
STPF0.3	Wolbachia	Genus	146	126,113	98.8%	
Other	–			1.2%	
STPF0.4	Wolbachia	Genus	202	97,713	98.2%	
Other	–			1.8%	
STPF0.5	Wolbachia	Genus	206	73,196	96.7%	
Pelomonas	Genus	1	676	0.9%	
Brevundimonas	Genus	1	170	0.2%	
Pseudomonas	Genus	1	75	0,1%	
Other	–			2.1%	

Table 3 Summary of the taxonomic composition of samples from Spain.

Province/district	Population	Sample	Taxonomic unit	Taxonomic level	OTUs reads	Reads	Raw fraction	
Huesca	Monzón	MZF0.1	Wolbachia	Genus	196	39,959	97.0%	
Nevskia	Genus	1	84	0.2%	
Methylophilus	Genus	1	60	0.1%	
Other	–			2.7%	
MZF0.2	Wolbachia	Genus	192	42,677	97.8%	
Enterobacteriaceae	Family	1	64	0.1%	
Other	–			2.1%	
MZF0.3	Wolbachia	Genus	225	40,718	96.3%	
Acinetobacter	Genus	3	299	0.7%	
Nevskia	Genus	1	58	0.1%	
Comamonadaceae	Family	1	45	0.1%	
Other	–			2.8%	
MZF0.4	Wolbachia	Genus	200	45,034	96.9%	
Gammaproteobacteria	Class	1	47	0.1%	
Other	–			3.0%	
MZF0.5	Wolbachia	Genus	194	40,367	96.1%	
Nevskia	Genus	1	191	0.5%	
Methylophilus	Genus	1	137	0.3%	
Aquabacterium	Genus	1	109	0.3%	
Comamonadaceae	Family	1	61	0.1%	
Other	–			2.7%	
Valencia	Paterna	PAF0.1	Wolbachia	Genus	209	45,910	97.4%	
Sphingomonas	Genus	2	223	0.5%	
Other	–			2.1%	
PAF0.2	Wolbachia	Genus	202	49,022	97.4%	
Other	–			2.6%	
PAF0.3	Wolbachia	Genus	182	41,312	97.4%	
Sphingomonas	Genus	1	49	0.1%	
Pelomonas	Genus	1	46	0.1%	
Other	–			2.4%	
PAF0.4	Wolbachia	Genus	166	30,946	96.5%	
Pelomonas	Genus	1	120	0.4%	
Sphingomonas	Genus	1	97	0.3%	
Nevskia	Genus	1	44	0.1%	
Other	–			2.7%	
Málaga	Rincón de la Victoria	RCF0.1	Wolbachia	Genus	189	49,296	97.5%	
Other	–			2.5%	
RCF0.2	Wolbachia	Genus	179	31,177	96.9%	
Other	–			3.1%	
RCF0.3	Wolbachia	Genus	177	48,900	97.5%	
Other	–			2.5%	
RCF0.4	Wolbachia	Genus	197	44,330	97.3%	
Other	–			2.7%	
RCF0.5	Wolbachia	Genus	189	37,996	97.1%	
Other	–			2.9%	

Shannon and Simpson indices (Table 4) revealed similar bacterial diversity between Spanish and São Tomé populations, as well as among individual samples (Supplementary Table S4).Table 4 Average of the Shannon and Simpson indices of the Spain and São Tomé populations.

Population	Shannon (average)	Simpson (average)	
Spain	4.229	0.819	
São Tomé	4.071	0.773	

Wolbachia genotyping

Of the 180 adult mosquitoes analyzed (97 females and 83 males), 178 (98.89%) were infected with Wolbachia. Five (2.78%) were infected with the wAlbA strain and 34 (18.88%) with the wAlbB strain. Mixed infections with wAlbA and wAlbB strains occurred in 77.22% of the infected mosquitoes (Table 5). Single wAlbA infections were exclusive to females (p < 0.05), while single wAlbB infections occurred in both males and females but were more prevalent in males (p < 0.05). Mixed infections were more common in females than males (85.57% in females and 67.47% in males, of the total number of infected mosquitoes) (p < 0.05). A predominance of mixed infections was observed in all locations, with rates ranging from 62.50 to 82.76%. However, simple wAlbB infections were detected only in the three provinces of Spain (Valencia, Huesca, and Málaga) with rates ranging from 13.79 to 31.25% (Table 6) and simple wAlbA infections were found only in two provinces of Spain (Valencia and Huesca). In São Tomé, only mixed infections by wAlbAeB were found, all in females. Furthermore, it was the only place where negative samples (N = 2) were detected in males.Table 5 Rate of Wolbachia infection in Ae. albopictus by sex.

	Infected	wAlbA	wAlbB	wAlbAeB	Total	
Females	97 (100%)	5 (5.15%)	9 (9.28%)	83 (85.57%)	97	
Males	81 (97.59%)	0 (0.0%)	25 (30.12%)	56 (67.47%)	83	
Total infected mosquitoes	178 (98.88%)	5 (2.78%)	34 (18.88%)	139 (77.22%)	180	

Table 6 Rate of Wolbachia infection, wAlbA and wAlbB, by locality.

Country	Province/district	Total	Infected	wAlbA	wAlbB	wAlbAeB	
Spain	Valencia	16	16 (100%)	1 (6.25%)	5 (31.25%)	10 (62.50%)	
Huesca	116	116 (100%)	4 (3.45%)	16 (1379%)	96 (82.76%)	
Málaga	41	41 (100%)	0 (0.00%)	13 (31.71%)	28 (68.29%)	
São Tomé and Príncipe	Água Grande	7	5 (71.43%)	0 (0.00%)	0 (0.00%)	5 (71.43%)	

Discussion

Aedes albopictus microbiome

This study represents the inaugural attempt to provide a comprehensive insight into the composition of the Ae. albopictus bacterial microbiome in wild populations originating from Spain and São Tomé through the sequencing of the V3–V4 regions of the 16S rRNA gene. The main phyla identified in the Ae. albopictus bacterial microbiome were Proteobacteria and Firmicutes, consistent with previous studies42,43. Both populations exhibited unique bacterial genera (Spain: Aquabacterium, Methylophilus, Nevskia, Solimonas and Sphingomonas; São Tomé: Asaia, Dechloromonas, Aerococcus, Bacillus, Brevundimonas, Pseudomonas, Rhizobium and Stenotrophomonas), suggesting a geographic influence on the microbiome43. A shared bacterial microbiome was also observed, comprising Acinetobacter, Pelomonas, Ralstonia, and Wolbachia genera. In addition to Wolbachia, these results suggest the existence of a bacterial microbiome shared among populations of Ae. albopictus geographically very distant44. The presence of Pelomonas and Nevskia genera in the Aedes genus is a novel finding, and their role in the microbiome of hematophagous insects remains unknown, and further studies are needed to understand their interactions with hosts45–48. Pelomonas and Nevskia are often found in environmental samples, raising the question of whether these bacteria could be externally adhered to the mosquitoes. To address this, we ensured that the mosquitoes were thoroughly cleaned before DNA extraction. The cleaning process involved washing the mosquitoes with a solution of 0.1% bleach followed by 70% alcohol. This procedure is designed to remove external contaminants, thereby ensuring that the DNA extracted and analyzed represents the internal microbiome. Therefore, the presence of Pelomonas and Nevskia is likely indicative of their role within the mosquito microbiome rather than external contamination.

Consistent with prior research, the bacterial composition of Ae. albopictus was predominantly governed by the Wolbachia genus49. When Wolbachia infection is highly dominant, it can "mask" the DNA of other bacteria present42,50 leading to undetected and unsequenced sequences. Future research may consider individual tissue-based microbiome analysis51, to ensure that Wolbachia, which is naturally abundant in the ovaries, does not interfere with the sequencing data.

This highly dominant infection by Wolbachia in both populations may be related to interactions between microbiome components, such as interspecific competition for resources44,52 but also to laboratory rearing conditions49. Research into the microbiota of laboratory and wild Aedes mosquitoes has revealed that both are primarily characterized by a limited number of phyla44. Nonetheless, although the overall microbiota composition is comparable in laboratory-reared and wild mosquitoes, it was observed that the diversity of midgut bacterial communities was more extensive in mosquitoes collected from natural environments. In this study, female mosquitoes were blood-fed, which can substantially modify their microbiota. When compared to other studies on field-collected mosquitoes, the number of taxa found per sample in our study is relatively low. This discrepancy could be attributed to several factors. First, the laboratory rearing conditions, including the controlled diet and environment, may have influenced the microbial diversity, leading to a reduction in taxa compared to field-collected samples. Second, the dominance of Wolbachia in our samples might have overshadowed the detection of other bacterial taxa, potentially due to interspecific competition or the high abundance of Wolbachia masking other bacterial DNA. In fact, the diversity of female midgut bacteria in Ae. albopictus from the laboratory was reduced when compared to females from those from the field53. These disparities could be originated from variations in the origin of the blood meal, dietary factors, and habitat53,54. Furthermore, recent studies have indicated a more consistent and greater bacterial diversity in breeding water compared to larvae and adults, regardless of sample source, with a notable decrease in the microbial community diversity between the larvae and newly emerged adult mosquitoes53. The impact of blood meals and rearing within an insectary environment on the Ae. albopictus microbiome warrants thorough investigation and future studies should consider the analysis of larvae’s microbiome and the bacterial diversity of the breeding water.

Wolbachia genotyping

To better understand the Wolbachia populations naturally present in these Ae. albopictus, A and B strains31,55 were genotyped. Both strains were found to infect the Spanish and São Tomé populations in three different scenarios: mixed infection with both strains, simple infection with wAlbA, or simple infection with wAlbB. Mixed infection, which is predominant worldwide56–60 was also the predominant form in this study. Several factors could explain these results, including species susceptibility to infection, facilitation of secondary infections by an active Wolbachia infection, or the stable maintenance of dual infections in hosts61.

In line with previous studies31,62, the prevalence of mixed infections was significantly higher in females compared to males and this could be due to the maternal transmission31, ineffective transmission of one strain, or physiological mechanisms allowing Wolbachia infection in females but not in males63. The exact mechanisms of Wolbachia infection and dissemination in mosquito vectors remain unclear64. Maternal transmission is a key mechanism for the spread of Wolbachia, and high rates of maternal transmission have been documented for mixed infections31. However, some studies have identified simple infections in the progeny of females with mixed infections, indicating the possibility of inefficient transmission of one of the strains63. This phenomenon could also occur in the transmission from female progenitors to their male offspring, explaining the lower prevalence of mixed infections in males. Physiological differences between male and female mosquitoes also play a role. Wolbachia tends to infect the germline cells of the host, persisting in the ovaries of females but not infecting the sperm in males. This differential infection could be due to a tropism of Wolbachia for ovarian tissues, allowing stable transmission through the female line, while infection in males is not stable or necessary for the maintenance of the endosymbiont31,62.

The prevalence of wAlbB infection was higher than wAlbA, supporting prior findings for European populations55,57,62. This inequality is particularly pronounced in males and could be attributed to reduced wAlbA density over time57 or a lower reproduction rate63,65. Research has suggested that females may not efficiently transmit wAlbA to their progeny, resulting in a reduction in male fertility due to Cytoplasmic Incompatibility (CI)55.

In this study, males from São Tomé population were not infected by Wolbachia. In fact, low Wolbachia prevalence has been detected in natural Ae. albopictus infections in Cameroon, one of the possible origins of Ae. albopictus from São Tomé66. Therefore, is not possible to assume that males in this region do not host Wolbachia, as infected females were found. Thus, the absence of infected males in São Tomé may be related to the lower prevalence of endosymbiont infection in males and the small sample size (N = 2). This small sample size limits the ability to draw definitive conclusions and may not accurately represent the population. Cautious interpretation is required, and future studies with larger sample sizes are necessary to validate this finding and provide a more accurate representation of Wolbachia prevalence in the male mosquito population from São Tomé.

The number of Wolbachia-negative samples was minimal (1.11%), indicating that innate Wolbachia infection may provide selective advantages to hosts, such as higher fecundity or hatching rates67. Geographic variation in Wolbachia infection has been reported68, but this study did not find significant differences between samples from different locations in Spain, as well as other studies carried out in other locations of this country60. An uniform distribution pattern with a predominance of mixed infections has been suggested and supported by prior studies for Ae. albopictus populations across Europe62.

The rapid dissemination and uniformity of distribution in Europe could result from multiple introductions69,70 in different Mediterranean countries, leading to constant mixed infection of Ae. albopictus by Wolbachia. Maternal inheritance and cytoplasmic incompatibility may also play a role in promoting the spread of Wolbachia infection71.

Aedes albopictus vector control

To control Ae. albopictus populations, one approach involves creating a colony with a stable infection of specific exogenous Wolbachia strains that cause Cytoplasmic Incompatibility (CI) and non-viable offspring. The Incompatible Insect Technique (IIT) involves the release of males carrying these Wolbachia strains into the environment, leading to a reduction in the target population. This approach has been tested and implemented in Ae. albopictus, which naturally hosts two Wolbachia strains (wAlbA and wAlbB). Introducing a third strain, wPip (found in Culex pipiens molestus), resulted in a triple-infected colony that effectively decreased the number of Ae. albopictus females in the study area72. Recent strategies employed wPip transinfection in Ae. albopictus, accompanied by the removal of its native Wolbachia strains, yielding promising outcomes73. These findings underscore IIT's potential as a future tool for Ae. albopictus control.

To conclude, this study significantly enhances our understanding of the bacterial microbiome of Ae. albopictus populations from Spain and São Tomé. Our study provides valuable insights into the microbiota composition of Ae. albopictus populations from Spain and São Tomé, revealing the presence of Wolbachia and other bacterial endosymbionts, which establishes a foundational framework for future manipulations and the introduction of new strains for vector control. As the scientific community intensifies its search for innovative and effective methods to combat vector-borne diseases, the insights from this study become increasingly invaluable. Consistent with previous research, we found that the microbiota of Ae. albopictus is predominantly governed by the Wolbachia genus. The dominance of Wolbachia in both populations suggests a stable association, which could be considered a part of the core microbiota of Ae. albopictus. The concept of core microbiota refers to the set of microbial taxa that are consistently found across different populations of a host species, regardless of geographical location or environmental conditions. In our study, we identified several bacterial genera that were shared between the Spanish and São Tomé populations, such as Acinetobacter, Pelomonas, Ralstonia, and Wolbachia. This finding supports the idea that Ae. albopictus harbors a core microbiota, which may play crucial roles in the mosquito's physiology, including reproduction, immune regulation, and vector competence. Furthermore, the identification of novel genera like Pelomonas and Nevskia in Ae. albopictus highlights the dynamic nature of the mosquito microbiota and suggests potential new areas of research to explore their roles in mosquito biology and pathogen transmission. The presence of these core and unique bacterial taxa emphasizes the importance of understanding the microbiota composition for developing targeted vector control strategies.

Despite the valuable insights provided by this study, there are several limitations that should be acknowledged. The most significant limitation is the small sample size used for the microbiome analysis and the Wolbachia genotype study, particularly for mosquitoes from São Tomé. The limited number of samples may have restricted the detection of the full diversity of the microbiota and Wolbachia strains, potentially overlooking important taxa and variations. Additionally, the use of laboratory-reared mosquitoes, although necessary for controlled conditions, might not fully represent the natural microbiota found in field populations due to differences in environmental exposures and diet.

To ensure the robustness of our findings, future studies should aim to include larger sample sizes and consider a broader range of geographical locations. Furthermore, incorporating field-collected mosquitoes and longitudinal sampling could provide a more comprehensive understanding of the microbiota dynamics and their implications for vector competence and control strategies.

While the presence of Wolbachia in Ae. albopictus is well established, our study contributes to the understanding of Wolbachia strain distribution and prevalence in new geographical regions, specifically Spain and São Tomé. We observed a consistent dominance of Wolbachia across different populations, reinforcing the stability of this endosymbiont's association with its host. However, we also highlight the presence of dual infections with different Wolbachia strains, which could have implications for vector control strategies leveraging Wolbachia-induced pathogen interference.

Future research directions to expand upon our findings include examining the impact of the bacterial microbiome on Ae. albopictus vector competence, evaluating the potential of CRISPR/Cas9 gene editing in endosymbionts for vector control, investigating paratransgenic strategies, assessing the ecological consequences of microbiome manipulation, and using advanced sequencing technologies for a more detailed understanding of the mosquito microbiome.

Supplementary Information

Supplementary Information.

Supplementary Information

The online version contains supplementary material available at 10.1038/s41598-024-71507-y.

Acknowledgements

This work was initially conceived and designed during research stay of DBB at GHTM in 2019, as part of the Short-Term Scientific Missions (STSM) within the AIM-COST CA17108 Action. Subsequently, this study was funded by the Fundação para a Ciência e Tecnologia funding to GHTM-UID/04413/2020, LA-REAL-LA/P/0117/2020 and “healTh RIsk and social vulnerability to Arboviral Diseases in mainland Portugal” (TRIAD)—Ref. PTDC/GES-OUT/30210/2017.

Author contributions

Research design: G.S. and C.A.S.; Mosquito collections: C.A.S., D.B.B., S.D.E.; Molecular genotyping: T.M., G.S.; Data analysis: G.S., T.M., C.A.S.; Manuscript writing: all authors; all authors reviewed the manuscript.

Data availability

Sequence data that support the findings of this study have been deposited in the NCBI under the accession number for SRA data PRJNA1028981 (temporary submission ID: SUB13907152).SRA records will be accessible with the following link after the release date (2024-11-01): http://www.ncbi.nlm.nih.gov/bioproject/1028981.

Competing interests

The authors declare no competing interests.

Publisher's note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
==== Refs
References

1. WHO. Launch of the Global Arbovirus Initiative. https://www.who.int/news-room/events/detail/2022/03/31/default-calendar/global-arbovirus-initiative (2022).
2. Capinera, J. Encyclopedia of Entomology (Springer, 2008).
3. Manni M Genetic evidence for a worldwide chaotic dispersion pattern of the arbovirus vector, Aedes albopictus PLoS Negl. Trop. Dis. 2017 11 e0005332 10.1371/journal.pntd.0005332 28135274
Manni, M. et al. Genetic evidence for a worldwide chaotic dispersion pattern of the arbovirus vector, Aedes albopictus. PLoS Negl. Trop. Dis. 11, e0005332 (2017).28135274 10.1371/journal.pntd.0005332
4. Paupy C Delatte H Bagny L Corbel V Fontenille D Aedes albopictus, an arbovirus vector: From the darkness to the light Microbes Infect. 2009 11 1177 1185 10.1016/j.micinf.2009.05.005 19450706
Paupy, C., Delatte, H., Bagny, L., Corbel, V. & Fontenille, D. Aedes albopictus, an arbovirus vector: From the darkness to the light. Microbes Infect. 11, 1177–1185 (2009).19450706 10.1016/j.micinf.2009.05.005
5. Effler PV Dengue fever, Hawaii, 2001–2002 Emerg. Infect. Dis. 2005 10.3201/eid1105.041063 16318695
Effler, P. V. et al. Dengue fever, Hawaii, 2001–2002. Emerg. Infect. Dis.10.3201/eid1105.041063 (2005).16318695 10.3201/eid1105.041063
6. Ramchurn, S. K., Moheeput, K. & Goorah, S. S. An analysis of a short-lived outbreak of dengue fever in Mauritius. Euro Surveill. 14, (2009).
7. Medlock JM An entomological review of invasive mosquitoes in Europe Bull. Entomol. Res. 2015 10.1017/S0007485315000103 26271257
Medlock, J. M. et al. An entomological review of invasive mosquitoes in Europe. Bull. Entomol. Res.10.1017/S0007485315000103 (2015).26271257 10.1017/S0007485315000103
8. Swan, T. et al. A literature review of dispersal pathways of Aedes albopictus across different spatial scales: implications for vector surveillance. Parasites Vectors 15, (2022).
9. Eritja R Palmer JRB Roiz D Sanpera-Calbet I Bartumeus F Direct evidence of adult Aedes albopictus dispersal by Car Sci. Rep. 2017 7 14399 10.1038/s41598-017-12652-5 29070818
Eritja, R., Palmer, J. R. B., Roiz, D., Sanpera-Calbet, I. & Bartumeus, F. Direct evidence of adult Aedes albopictus dispersal by Car. Sci. Rep. 7, 14399 (2017).29070818 10.1038/s41598-017-12652-5
10. Gratz NG Critical review of the vector status of Aedes albopictus Med. Vet. Entomol. 2004 10.1111/j.0269-283X.2004.00513.x 15347388
Gratz, N. G. Critical review of the vector status of Aedes albopictus. Med. Vet. Entomol.10.1111/j.0269-283X.2004.00513.x (2004).15347388 10.1111/j.0269-283X.2004.00513.x
11. Hasan S Jamdar SF Alalowi M Beaiji AAA S. M. Dengue virus: A global human threat: Review of literature J. Int. Soc. Prev. Commun. Dent. 2016 10.4103/2231-0762.175416
Hasan, S., Jamdar, S. F., Alalowi, M., Beaiji, A. A. A. & S. M.,. Dengue virus: A global human threat: Review of literature. J. Int. Soc. Prev. Commun. Dent.10.4103/2231-0762.175416 (2016).10.4103/2231-0762.175416
12. Reis S Cornel AJ Melo M Pereira H Loiseau C First record of Aedes albopictus (Skuse 1894) on São tomé island Acta Trop. 2017 171 86 89 10.1016/j.actatropica.2017.03.035 28377219
Reis, S., Cornel, A. J., Melo, M., Pereira, H. & Loiseau, C. First record of Aedes albopictus (Skuse 1894) on São tomé island. Acta Trop. 171, 86–89 (2017).28377219 10.1016/j.actatropica.2017.03.035
13. Goiri F Progressive invasion of Aedes albopictus in Northern Spain in the period 2013–2018 and a possible association with the increase in insect bites Int. J. Environ. Res. Public Health 2020 17 1678 10.3390/ijerph17051678 32143518
Goiri, F. et al. Progressive invasion of Aedes albopictus in Northern Spain in the period 2013–2018 and a possible association with the increase in insect bites. Int. J. Environ. Res. Public Health 17, 1678 (2020).32143518 10.3390/ijerph17051678
14. Loiseau C High endemism of mosquitoes on São Tomé and Príncipe Islands: Evaluating the general dynamic model in a worldwide island comparison Insect. Conserv. Divers. 2019 12 69 79 10.1111/icad.12308
Loiseau, C. et al. High endemism of mosquitoes on São Tomé and Príncipe Islands: Evaluating the general dynamic model in a worldwide island comparison. Insect. Conserv. Divers. 12, 69–79 (2019).10.1111/icad.12308
15. Aranda C Eritja R Roiz D First record and establishment of the mosquito Aedes albopictus in Spain Med. Vet. Entomol. 2006 20 150 152 10.1111/j.1365-2915.2006.00605.x 16608499
Aranda, C., Eritja, R. & Roiz, D. First record and establishment of the mosquito Aedes albopictus in Spain. Med. Vet. Entomol. 20, 150–152 (2006).16608499 10.1111/j.1365-2915.2006.00605.x
16. European Centre for Disease Prevention and Control and European Food Safety Authority. Mosquito maps. European Centre for Disease Prevention and Control (ECDC). https://ecdc.europa.eu/en/disease-vectors/surveillance-and-disease-data/mosquito-maps (2024).
17. Garrido M Aedes albopictus in a recently invaded area in Spain: Effects of trap type, locality, and season on mosquito captures Sci. Rep. 2024 14 2131 10.1038/s41598-024-52040-4 38267495
Garrido, M. et al. Aedes albopictus in a recently invaded area in Spain: Effects of trap type, locality, and season on mosquito captures. Sci. Rep. 14, 2131 (2024).38267495 10.1038/s41598-024-52040-4
18. Collantes F Review of ten-years presence of Aedes albopictus in Spain 2004–2014: Known distribution and public health concerns Parasit Vectors 2015 8 1 11 10.1186/s13071-015-1262-y 25561160
Collantes, F. et al. Review of ten-years presence of Aedes albopictus in Spain 2004–2014: Known distribution and public health concerns. Parasit Vectors 8, 1–11 (2015).25561160 10.1186/s13071-015-1262-y
19. Lucati F Multiple invasions, Wolbachia and human-aided transport drive the genetic variability of Aedes albopictus in the Iberian Peninsula Sci. Rep. 2022 12 20682 10.1038/s41598-022-24963-3 36450768
Lucati, F. et al. Multiple invasions, Wolbachia and human-aided transport drive the genetic variability of Aedes albopictus in the Iberian Peninsula. Sci. Rep. 12, 20682 (2022).36450768 10.1038/s41598-022-24963-3
20. ECDC. Autochthonous Cases of Dengue in Spain and France (2019).
21. Monge S Characterization of the first autochthonous dengue outbreak in Spain (August–September 2018) Acta Trop. 2020 205 105402 10.1016/j.actatropica.2020.105402 32088276
Monge, S. et al. Characterization of the first autochthonous dengue outbreak in Spain (August–September 2018). Acta Trop. 205, 105402 (2020).32088276 10.1016/j.actatropica.2020.105402
22. Navero-Castillejos J Molecular characterization of imported and autochthonous dengue in northeastern spain Viruses 2021 13 1910 10.3390/v13101910 34696340
Navero-Castillejos, J. et al. Molecular characterization of imported and autochthonous dengue in northeastern spain. Viruses 13, 1910 (2021).34696340 10.3390/v13101910
23. Gaio ADO Contribution of midgut bacteria to blood digestion and egg production in Aedes aegypti (diptera: Culicidae) (L.) Parasit. Vectors 2011 4 1 10 10.1186/1756-3305-4-105 21205315
Gaio, A. D. O. et al. Contribution of midgut bacteria to blood digestion and egg production in Aedes aegypti (diptera: Culicidae) (L.). Parasit. Vectors 4, 1–10 (2011).21205315 10.1186/1756-3305-4-105
24. Coon KL Brown MR Strand MR Gut bacteria differentially affect egg production in the anautogenous mosquito Aedes aegypti and facultatively autogenous mosquito Aedes atropalpus (Diptera: Culicidae) Parasit. Vectors 2016 9 1 12 10.1186/s13071-016-1660-9 26728523
Coon, K. L., Brown, M. R. & Strand, M. R. Gut bacteria differentially affect egg production in the anautogenous mosquito Aedes aegypti and facultatively autogenous mosquito Aedes atropalpus (Diptera: Culicidae). Parasit. Vectors 9, 1–12 (2016).26728523 10.1186/s13071-016-1660-9
25. Sicard M Bonneau M Weill M Wolbachia prevalence, diversity, and ability to induce cytoplasmic incompatibility in mosquitoes Curr. Opin. Insect Sci. 2019 10.1016/j.cois.2019.02.005 31247412
Sicard, M., Bonneau, M. & Weill, M. Wolbachia prevalence, diversity, and ability to induce cytoplasmic incompatibility in mosquitoes. Curr. Opin. Insect Sci.10.1016/j.cois.2019.02.005 (2019).31247412 10.1016/j.cois.2019.02.005
26. Xi Z Ramirez JL Dimopoulos G The Aedes aegypti toll pathway controls dengue virus infection PLoS Pathog. 2008 4 e1000098 10.1371/journal.ppat.1000098 18604274
Xi, Z., Ramirez, J. L. & Dimopoulos, G. The Aedes aegypti toll pathway controls dengue virus infection. PLoS Pathog. 4, e1000098 (2008).18604274 10.1371/journal.ppat.1000098
27. Dickson LB Carryover effects of larval exposure to different environmental bacteria drive adult trait variation in a mosquito vector Sci. Adv. 2017 3 e1700585 10.1126/sciadv.1700585 28835919
Dickson, L. B. et al. Carryover effects of larval exposure to different environmental bacteria drive adult trait variation in a mosquito vector. Sci. Adv. 3, e1700585 (2017).28835919 10.1126/sciadv.1700585
28. Souza-Neto JA Powell JR Bonizzoni M Aedes aegypti vector competence studies: A review Infect. Genet. Evol. 2019 10.1016/j.meegid.2018.11.009 30465912
Souza-Neto, J. A., Powell, J. R. & Bonizzoni, M. Aedes aegypti vector competence studies: A review. Infect. Genet. Evol.10.1016/j.meegid.2018.11.009 (2019).30465912 10.1016/j.meegid.2018.11.009
29. Jupatanakul N Sim S Dimopoulos G The insect microbiome modulates vector competence for arboviruses Viruses. 2014 10.3390/v6114294 25415198
Jupatanakul, N., Sim, S. & Dimopoulos, G. The insect microbiome modulates vector competence for arboviruses. Viruses.10.3390/v6114294 (2014).25415198 10.3390/v6114294
30. Blagrove MSC Arias-Goeta C Di Genua C Failloux AB Sinkins SP A Wolbachia wMel transinfection in Aedes albopictus is not detrimental to host fitness and inhibits Chikungunya virus PLoS Negl. Trop. Dis. 2013 7 e2152 10.1371/journal.pntd.0002152 23556030
Blagrove, M. S. C., Arias-Goeta, C., Di Genua, C., Failloux, A. B. & Sinkins, S. P. A Wolbachia wMel transinfection in Aedes albopictus is not detrimental to host fitness and inhibits Chikungunya virus. PLoS Negl. Trop. Dis. 7, e2152 (2013).23556030 10.1371/journal.pntd.0002152
31. Ahmad NA Vythilingam I Lim YAL Zabari NZAM Lee HL Detection of Wolbachia in Aedes albopictus and their effects on chikungunya virus Am. J. Trop. Med. Hyg. 2017 96 148 156 10.4269/ajtmh.16-0516 27920393
Ahmad, N. A., Vythilingam, I., Lim, Y. A. L., Zabari, N. Z. A. M. & Lee, H. L. Detection of Wolbachia in Aedes albopictus and their effects on chikungunya virus. Am. J. Trop. Med. Hyg. 96, 148–156 (2017).27920393 10.4269/ajtmh.16-0516
32. Novakova, E. et al. Mosquito microbiome dynamics, a background for prevalence and seasonality of West Nile virus. Front. Microbiol. 8, (2017).
33. Dada N Considerations for mosquito microbiome research from the Mosquito Microbiome Consortium Microbiome. 2021 10.1186/s40168-020-00987-7 33522965
Dada, N. et al. Considerations for mosquito microbiome research from the Mosquito Microbiome Consortium. Microbiome.10.1186/s40168-020-00987-7 (2021).33522965 10.1186/s40168-020-00987-7
34. Collins FH Comparison of DNA-probe and isoenzyme methods for differentiating Anopheles gambiae and Anopheles arabiensis (Diptera: Culicidae) J. Med. Entomol. 1988 25 116 120 10.1093/jmedent/25.2.116 3280799
Collins, F. H. et al. Comparison of DNA-probe and isoenzyme methods for differentiating Anopheles gambiae and Anopheles arabiensis (Diptera: Culicidae). J. Med. Entomol. 25, 116–120 (1988).3280799 10.1093/jmedent/25.2.116
35. Edgar RC Haas BJ Clemente JC Quince C Knight R UCHIME improves sensitivity and speed of chimera detection Bioinformatics 2011 27 2194 2200 10.1093/bioinformatics/btr381 21700674
Edgar, R. C., Haas, B. J., Clemente, J. C., Quince, C. & Knight, R. UCHIME improves sensitivity and speed of chimera detection. Bioinformatics 27, 2194–2200 (2011).21700674 10.1093/bioinformatics/btr381
36. Rognes, T., Flouri, T., Nichols, B., Quince, C. & Mahé, F. VSEARCH: A versatile open source tool for metagenomics. PeerJ 2016, (2016).
37. Eren AM Oligotyping: di erentiating between closely related microbial taxa using 16s rRNA gene data Methods Ecol. Evol. 2013 4 1111 11119 10.1111/2041-210X.12114 24358444
Eren, A. M. et al. Oligotyping: di erentiating between closely related microbial taxa using 16s rRNA gene data. Methods Ecol. Evol. 4, 1111–11119 (2013).24358444 10.1111/2041-210X.12114
38. Eren AM Minimum entropy decomposition: Unsupervised oligotyping for sensitive partitioning of high-throughput marker gene sequences ISME J. 2015 9 968 979 10.1038/ismej.2014.195 25325381
Eren, A. M. et al. Minimum entropy decomposition: Unsupervised oligotyping for sensitive partitioning of high-throughput marker gene sequences. ISME J. 9, 968–979 (2015).25325381 10.1038/ismej.2014.195
39. Angly, F. E. et al. CopyRighter: A rapid tool for improving the accuracy of microbial community profiles through lineage-specific gene copy number correction. Microbiome 2, (2014).
40. Zhou W Rousset F O’Neill S Phylogeny and PCR-based classification of Wolbachia strains using wsp gene sequences Proc. R. Soc. B Biol. Sci. 1998 265 209 515 10.1098/rspb.1998.0324
Zhou, W., Rousset, F. & O’Neill, S. Phylogeny and PCR-based classification of Wolbachia strains using wsp gene sequences. Proc. R. Soc. B Biol. Sci. 265, 209–515 (1998).10.1098/rspb.1998.0324
41. Vassarstats. Vassarstats: Website for Statistical Computation). http://vassarstats.net/newcs.html.
42. Bennett KL Dynamics and diversity of bacteria associated with the disease vectors Aedes aegypti and Aedes albopictus Sci. Rep. 2019 9 12160 10.1038/s41598-019-48414-8 31434963
Bennett, K. L. et al. Dynamics and diversity of bacteria associated with the disease vectors Aedes aegypti and Aedes albopictus. Sci. Rep. 9, 12160 (2019).31434963 10.1038/s41598-019-48414-8
43. Zhao T Metagenome sequencing reveals the microbiome of Aedes albopictus and its possible relationship with dengue virus susceptibility Front. Microbiol. 2022 13 891151 10.3389/fmicb.2022.891151 35633691
Zhao, T. et al. Metagenome sequencing reveals the microbiome of Aedes albopictus and its possible relationship with dengue virus susceptibility. Front. Microbiol. 13, 891151 (2022).35633691 10.3389/fmicb.2022.891151
44. Scolari F Casiraghi M Bonizzoni M Aedes spp. and their microbiota: A review Front. Microbiol. 2019 10.3389/fmicb.2019.02036 31551973
Scolari, F., Casiraghi, M. & Bonizzoni, M. Aedes spp. and their microbiota: A review. Front. Microbiol.10.3389/fmicb.2019.02036 (2019).31551973 10.3389/fmicb.2019.02036
45. Vicente CSL Characterization of bacterial communities associated with the pine sawyer beetle Monochamus galloprovincialis, the insect vector of the pinewood nematode Bursaphelenchus xylophilus FEMS Microbiol. Lett. 2013 10.1111/1574-6968.12232 23927049
Vicente, C. S. L. et al. Characterization of bacterial communities associated with the pine sawyer beetle Monochamus galloprovincialis, the insect vector of the pinewood nematode Bursaphelenchus xylophilus. FEMS Microbiol. Lett.10.1111/1574-6968.12232 (2013).23927049 10.1111/1574-6968.12232
46. Zhukova M Sapountzis P Schiøtt M Boomsma JJ Diversity and transmission of gut bacteria in Atta and Acromyrmex leaf-cutting ants during development Front. Microbiol. 2017 8 1942 10.3389/fmicb.2017.01942 29067008
Zhukova, M., Sapountzis, P., Schiøtt, M. & Boomsma, J. J. Diversity and transmission of gut bacteria in Atta and Acromyrmex leaf-cutting ants during development. Front. Microbiol. 8, 1942 (2017).29067008 10.3389/fmicb.2017.01942
47. Santos-Garcia D Mestre-Rincon N Zchori-Fein E Morin S Inside out: Microbiota dynamics during host-plant adaptation of whiteflies ISME J. 2020 14 847 856 10.1038/s41396-019-0576-8 31896788
Santos-Garcia, D., Mestre-Rincon, N., Zchori-Fein, E. & Morin, S. Inside out: Microbiota dynamics during host-plant adaptation of whiteflies. ISME J. 14, 847–856 (2020).31896788 10.1038/s41396-019-0576-8
48. Tabbabi A Mizushima D Yamamoto DS Kato H Sand flies and their microbiota Parasitologia 2022 2 71 87 10.3390/parasitologia2020008
Tabbabi, A., Mizushima, D., Yamamoto, D. S. & Kato, H. Sand flies and their microbiota. Parasitologia 2, 71–87 (2022).10.3390/parasitologia2020008
49. Lin D Bacterial composition of midgut and entire body of laboratory colonies of Aedes aegypti and Aedes albopictus from Southern China Parasit. Vectors 2021 14 1 13 10.1186/s13071-021-05050-4 33388087
Lin, D. et al. Bacterial composition of midgut and entire body of laboratory colonies of Aedes aegypti and Aedes albopictus from Southern China. Parasit. Vectors 14, 1–13 (2021).33388087 10.1186/s13071-021-05050-4
50. Minard G Pyrosequencing 16S rRNA genes of bacteria associated with wild tiger mosquito Aedes albopictus: A pilot study Front. Cell Infect. Microbiol. 2014 4 59 10.3389/fcimb.2014.00059 24860790
Minard, G. et al. Pyrosequencing 16S rRNA genes of bacteria associated with wild tiger mosquito Aedes albopictus: A pilot study. Front. Cell Infect. Microbiol. 4, 59 (2014).24860790 10.3389/fcimb.2014.00059
51. Rosso F Reduced diversity of gut microbiota in two Aedes mosquitoes species in areas of recent invasion Sci. Rep. 2018 8 16091 10.1038/s41598-018-34640-z 30382151
Rosso, F. et al. Reduced diversity of gut microbiota in two Aedes mosquitoes species in areas of recent invasion. Sci. Rep. 8, 16091 (2018).30382151 10.1038/s41598-018-34640-z
52. Brinker P Fontaine MC Beukeboom LW Falcao Salles J Host, symbionts, and the microbiome: The missing tripartite interaction Trends Microbiol. 2019 10.1016/j.tim.2019.02.002 30857919
Brinker, P., Fontaine, M. C., Beukeboom, L. W. & Falcao Salles, J. Host, symbionts, and the microbiome: The missing tripartite interaction. Trends Microbiol.10.1016/j.tim.2019.02.002 (2019).30857919 10.1016/j.tim.2019.02.002
53. Scolari F Exploring changes in the microbiota of Aedes albopictus: Comparison among breeding site water, larvae, and adults Front. Microbiol. 2021 12 624170 10.3389/fmicb.2021.624170 33584626
Scolari, F. et al. Exploring changes in the microbiota of Aedes albopictus: Comparison among breeding site water, larvae, and adults. Front. Microbiol. 12, 624170 (2021).33584626 10.3389/fmicb.2021.624170
54. Tuanudom R Yurayart N Rodkhum C Tiawsirisup S Diversity of midgut microbiota in laboratory-colonized and field-collected Aedes albopictus (Diptera: Culicidae): A preliminary study Heliyon 2021 7 e08259 10.1016/j.heliyon.2021.e08259 34765765
Tuanudom, R., Yurayart, N., Rodkhum, C. & Tiawsirisup, S. Diversity of midgut microbiota in laboratory-colonized and field-collected Aedes albopictus (Diptera: Culicidae): A preliminary study. Heliyon 7, e08259 (2021).34765765 10.1016/j.heliyon.2021.e08259
55. Tortosa P Wolbachia age-sex-specific density in Aedes albopictus: A host evolutionary response to Cytoplasmic Incompatibility? PLoS One 2010 5 e9700 10.1371/journal.pone.0009700 20300514
Tortosa, P. et al. Wolbachia age-sex-specific density in Aedes albopictus: A host evolutionary response to Cytoplasmic Incompatibility?. PLoS One 5, e9700 (2010).20300514 10.1371/journal.pone.0009700
56. Lo N Casiraghi M Salati E Bazzocchi C Bandi C How many Wolbachia supergroups exist? Mol. Biol. Evol. 2002 10.1093/oxfordjournals.molbev.a004087 11861893
Lo, N., Casiraghi, M., Salati, E., Bazzocchi, C. & Bandi, C. How many Wolbachia supergroups exist?. Mol. Biol. Evol.10.1093/oxfordjournals.molbev.a004087 (2002).11861893 10.1093/oxfordjournals.molbev.a004087
57. Shaikevich, E., Bogacheva, A. & Ganushkina, L. Dirofilaria and Wolbachia in mosquitoes (Diptera: Culicidae) in central European Russia and on the Black Sea coast. Parasite. 26, (2019).
58. Hu, Y. et al. Identification and molecular characterization of Wolbachia strains in natural populations of Aedes albopictus in China. Parasit. Vectors. 13, (2020).
59. Puerta-Guardo, H. et al. Wolbachia in native populations of Aedes albopictus (Diptera: Culicidae) from Yucatan Peninsula, Mexico. J. Insect Sci. 20, (2020).
60. Bueno-Marí, R. et al. Infecciones por Wolbachia pipientis en poblaciones de Aedes albopictus en la ciudad de València (España): Implicaciones para el control de mosquitos. Rev. Esp. Salud. Publica. 97, (2023).
61. Werren JH Windsor D Guo L Distribution of Wolbachia among neotropical arthropods Proc. R. Soc. B Biol. Sci. 1995 262 197 204 10.1098/rspb.1995.0196
Werren, J. H., Windsor, D. & Guo, L. Distribution of Wolbachia among neotropical arthropods. Proc. R. Soc. B Biol. Sci. 262, 197–204 (1995).10.1098/rspb.1995.0196
62. Belo F Wolbachia Infection in European Populations of Aedes albopictus 2021 Universidade Nova de Lisboa Instituto de Higiene e Medicina Tropical
Belo, F. Wolbachia Infection in European Populations of Aedes albopictus (Universidade Nova de Lisboa Instituto de Higiene e Medicina Tropical, 2021).
63. Kittayapong P Baisley KJ Sharpe RG Baimai V O’Neill SL Maternal transmission efficiency of Wolbachia superinfections in Aedes albopictus populations in Thailand Am. J. Trop. Med. Hyg. 2002 66 103 107 10.4269/ajtmh.2002.66.103 12135258
Kittayapong, P., Baisley, K. J., Sharpe, R. G., Baimai, V. & O’Neill, S. L. Maternal transmission efficiency of Wolbachia superinfections in Aedes albopictus populations in Thailand. Am. J. Trop. Med. Hyg. 66, 103–107 (2002).12135258 10.4269/ajtmh.2002.66.103
64. Bi J Wang YF The effect of the endosymbiont Wolbachia on the behavior of insect hosts Insect Sci. 2020 10.1111/1744-7917.12731 31631529
Bi, J. & Wang, Y. F. The effect of the endosymbiont Wolbachia on the behavior of insect hosts. Insect Sci.10.1111/1744-7917.12731 (2020).31631529 10.1111/1744-7917.12731
65. Wiwatanaratanabutr I Kittayapong P Effects of crowding and temperature on Wolbachia infection density among life cycle stages of Aedes albopictus J. Invertebr. Pathol. 2009 102 220 224 10.1016/j.jip.2009.08.009 19686755
Wiwatanaratanabutr, I. & Kittayapong, P. Effects of crowding and temperature on Wolbachia infection density among life cycle stages of Aedes albopictus. J. Invertebr. Pathol. 102, 220–224 (2009).19686755 10.1016/j.jip.2009.08.009
66. Bamou R Wolbachia detection in field-collected mosquitoes from Cameroon Insects 2021 12 1133 10.3390/insects12121133 34940221
Bamou, R. et al. Wolbachia detection in field-collected mosquitoes from Cameroon. Insects 12, 1133 (2021).34940221 10.3390/insects12121133
67. Dobson SL Rattanadechakul W Marsland EJ Fitness advantage and cytoplasmic incompatibility in Wolbachia single- and superinfected Aedes albopictus Heredity (Edinb) 2004 93 135 142 10.1038/sj.hdy.6800458 15127087
Dobson, S. L., Rattanadechakul, W. & Marsland, E. J. Fitness advantage and cytoplasmic incompatibility in Wolbachia single- and superinfected Aedes albopictus. Heredity (Edinb) 93, 135–142 (2004).15127087 10.1038/sj.hdy.6800458
68. Joanne S Distribution and dynamics of Wolbachia infection in Malaysian Aedes albopictus Acta Trop. 2015 148 38 45 10.1016/j.actatropica.2015.04.003 25899523
Joanne, S. et al. Distribution and dynamics of Wolbachia infection in Malaysian Aedes albopictus. Acta Trop. 148, 38–45 (2015).25899523 10.1016/j.actatropica.2015.04.003
69. Pichler V Complex interplay of evolutionary forces shaping population genomic structure of invasive Aedes albopictus in southern Europe PLoS Negl. Trop. Dis. 2019 13 e0007554 10.1371/journal.pntd.0007554 31437154
Pichler, V. et al. Complex interplay of evolutionary forces shaping population genomic structure of invasive Aedes albopictus in southern Europe. PLoS Negl. Trop. Dis. 13, e0007554 (2019).31437154 10.1371/journal.pntd.0007554
70. Kotsakiozi P Population genomics of the Asian tiger mosquito, Aedes albopictus: Insights into the recent worldwide invasion Ecol. Evol. 2017 7 10143 10157 10.1002/ece3.3514 29238544
Kotsakiozi, P. et al. Population genomics of the Asian tiger mosquito, Aedes albopictus: Insights into the recent worldwide invasion. Ecol. Evol. 7, 10143–10157 (2017).29238544 10.1002/ece3.3514
71. Alphey L Genetic control of mosquitoes Annu. Rev. Entomol. 2014 59 205 224 10.1146/annurev-ento-011613-162002 24160434
Alphey, L. Genetic control of mosquitoes. Annu. Rev. Entomol. 59, 205–224 (2014).24160434 10.1146/annurev-ento-011613-162002
72. Mains JW Brelsfoard CL Rose RI Dobson SL Female adult Aedes albopictus suppression by Wolbachia-infected male mosquitoes Sci. Rep. 2016 6 33846 10.1038/srep33846 27659038
Mains, J. W., Brelsfoard, C. L., Rose, R. I. & Dobson, S. L. Female adult Aedes albopictus suppression by Wolbachia-infected male mosquitoes. Sci. Rep. 6, 33846 (2016).27659038 10.1038/srep33846
73. Caputo B A bacterium against the tiger: Preliminary evidence of fertility reduction after release of Aedes albopictus males with manipulated Wolbachia infection in an Italian urban area Pest Manag. Sci. 2020 76 1324 1332 10.1002/ps.5643 31603613
Caputo, B. et al. A bacterium against the tiger: Preliminary evidence of fertility reduction after release of Aedes albopictus males with manipulated Wolbachia infection in an Italian urban area. Pest Manag. Sci. 76, 1324–1332 (2020).31603613 10.1002/ps.5643
