Molecular phylogeny of heritable symbionts and microbiota diversity analysis in phlebotominae sand flies and Culex nigripalpus from Colombia

BACKGROUND: Secondary symbionts of insects include a range of bacteria and fungi that perform various functional roles on their hosts, such as fitness, tolerance to heat stress, susceptibility to insecticides and effects on reproduction. These endosymbionts could have the potential to shape microbial communites and high potential to develop strategies for mosquito-borne disease control. METHODOLOGY/PRINCIPAL FINDINGS: The relative frequency and molecular phylogeny of Wolbachia, Microsporidia and Cardinium were determined of phlebotomine sand flies and mosquitoes in two regions from Colombia. Illumina Miseq using the 16S rRNA gene as a biomarker was conducted to examine the microbiota. Different percentages of natural infection by Wolbachia, Cardinium, and Microsporidia in phlebotomines and mosquitoes were detected. Phylogenetic analysis of Wolbachia shows putative new strains of Lutzomyia gomezi (wLgom), Brumptomyia hamata (wBrham), and a putative new group associated with Culex nigripalpus (Cnig) from the Andean region, located in Supergroup A and Supergroup B, respectively. The sequences of Microsporidia were obtained of Pi. pia and Cx. nigripalpus, which are located on phylogeny in the IV clade (terrestrial origin). The Cardinium of Tr. triramula and Ps. shannoni were located in group C next to Culicoides sequences while Cardinium of Mi. cayennensis formed two putative new subgroups of Cardinium in group A. In total were obtained 550 bacterial amplicon sequence variants (ASVs) and 189 taxa to the genus level. The microbiota profiles of Sand flies and mosquitoes showed mainly at the phylum level to Proteobacteria (67.6%), Firmicutes (17.9%) and Actinobacteria (7.4%). High percentages of relative abundance for Wolbachia (30%-83%) in Lu. gomezi, Ev. dubitans, Mi. micropyga, Br. hamata, and Cx. nigripalpus were found. ASVs assigned as Microsporidia were found in greater abundance in Pi. pia (23%) and Cx. nigripalpus (11%). An important finding is the detection of Rickettsia in Pi. pia (58,8%) and Bartonella sp. in Cx. nigripalpus. CONCLUSIONS/SIGNIFICANCE: We found that Wolbachia infection significantly decreased the alpha diversity and negatively impacts the number of taxa on sand flies and Culex nigripalpus. The Principal Coordinate Analysis (PCoA) is consistent, which showed statistically significant differences (PERMANOVA, F = 2.4744; R2 = 0.18363; p-value = 0.007) between the microbiota of sand flies and mosquitoes depending on its origin, host and possibly for the abundance of some endosymbionts (Wolbachia, Rickettsia).