Next Article in Journal
Bovine Papillomavirus 1 Gets Out of the Flock: Detection in an Ovine Wart in Sicily
Previous Article in Journal
Serodynamic Analysis of the Piglets Born from Sows Vaccinated with Modified Live Vaccine or E2 Subunit Vaccine for Classical Swine Fever
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Host Species Determines the Composition of the Prokaryotic Microbiota in Phlebotomus Sandflies

by
Christos Papadopoulos
1,
Panagiotis A. Karas
1,
Sotirios Vasileiadis
1,
Panagiota Ligda
2,
Anastasios Saratsis
2,
Smaragda Sotiraki
2,* and
Dimitrios G. Karpouzas
1,*
1
Laboratory of Plant and Environmental Biotechnology, Department of Biochemistry and Biotechnology, University of Thessaly, 41500 Larissa, Greece
2
Group of Parasitology, Veterinary Research Institute, Hellenic Agricultural Organization-Demeter, 57100 Thessaloniki, Greece
*
Authors to whom correspondence should be addressed.
Pathogens 2020, 9(6), 428; https://doi.org/10.3390/pathogens9060428
Submission received: 9 May 2020 / Revised: 27 May 2020 / Accepted: 28 May 2020 / Published: 29 May 2020

Abstract

:
Phlebotomine sandflies are vectors of the humans’ and mammals’ parasite Leishmania spp. Although the role of gut microbiome in the biological cycle of insects is acknowledged, we still know little about the factors modulating the composition of the gut microbiota of sandflies. We tested whether host species impose a strong structural effect on the gut microbiota of Phlebotomus spp. Sandflies were collected from the island of Leros, Greece, and classified to P. papatasi, P. neglectus, P. tobbi, and P. similis, all being negative to Leishmania spp. The prokaryotic gut microbiota was determined via 16S rRNA gene amplicon sequencing. Phlebotomus species supported distinct microbial communities (p < 0.001). P. papatasi microbiota was the most distinct over-dominated by three Spiroplasma, Wolbachia and Paenibacillus operational taxonomic units (OTUs), while another Wolbachia OTU prevailed in P. neglectus. Conversely, the microbiota of P. tobbi and P. similis was composed of several less dominant OTUs. Archaea showed low presence with the dominant OTUs belonging to methanogenic Euryarcheota, ammonia-oxidizing Thaumarcheota, and Nanoarchaeota. We provide first insights into the composition of the bacterial and archaeal community of Phlebotomus sandflies and showed that, in the absence of Leishmania, host genotype is the major modulator of Phlebotomus sandfly gut microbiota.

1. Introduction

Phlebotomine sandflies (Diptera: Psychodidae, Phlebotominae) are insects of global health importance. This is based on their versatility as vectors of human pathogens including viruses (Phleboviruses) [1,2], bacteria (Bartonella bacilliformis) [3] and protozoa like Leishmania spp. [4]. In Europe and mainly in the Mediterranean region, sandflies of the genus Phlebotomus have been incriminated as vectors of Leishmania [5] transmitted by the bite of an infected female sandfly [6].
Collective field studies in the Mediterranean region [7] and neighboring countries like Albania [8] and Serbia [9,10], verified the presence of several Phlebotomus species like P. neglectus, P. tobbi, P. balkanicus, P. ariasi, P. perfiliewi, P. perniciosus, P. kandelakii, P. mascitii, P. alexandri, P. sergenti. P. similis, P. simici, P. patasi, all being proven or potential vectors of Leishmania. In Greece, monitoring surveys identified P. neglectus (32.8% of the population captured in the study), P. similis (30.3%), P. tobbi (16.7%), and P. perfiliewi (15.9%) in Ionian islands, whereas P. simici (50%), P. neglectus (24.5%), and P. tobbi (9.6%) predominated in the North Greece mainland (Xanthopoulou et al.) [11]. Recent studies in 11 islands in the Aegean Sea detected a rich sandfly fauna comprised of P. neglectus, P. tobbi, P. similis, P. simici, P. perfiliewi, P. alexandri and P. papatasi [12].
Insects live in association and interact with a diverse microbiota. These interactions have serious implications for the fitness, ecology and evolution of insects. Gut microbiota provides beneficial services to their insect host. They facilitate the nutrition of insects which feed on nutrient-deficient diets enabling their survival in oligotrophic environments [13]. Endosymbiotic vertically inherited bacteria like Wolbachia and Spiroplasma manipulate the reproduction system of their hosts to prevail in their population over endosymbiont-free individuals [14], alter the sex-ratio of the population via production of male-killing toxins [15] or protect their hosts (i.e., Drosophila melanogaster, Aeges aegypti) from parasitic wasps [16] and other insect parasites [17]. Insects also utilize microbiome as a rapid evolution mechanism to adapt to rapidly changing environmental conditions [18]. In this context, a symbiotic Citrobacter carrying an organophosphate hydrolase conferred resistance in Bactocera dorsalis to organophosphate insecticides [19]. Conversely, insect symbiotic bacteria could be detrimental to their hosts by conferring pathogenesis (i.e., Spiroplasma and honeybees) [20] or by interacting with plant defense systems exacerbating their effect on insect pests [21].
Although we have just started to appreciate the role of microorganisms, the association between the composition of insect gut microbiome and phylogenetic or ecologic traits of the host is not well defined yet. Several studies have investigated the structural role of diet, the local ecosystem and climatic conditions, the phylogeny and the life stage on insect gut microbiota. Kolasa et al. [22] showed that the composition of the gut microbial community in beetles was better explained by trophic guild. Muturi et al. [23] showed that sampling location exhibited the strongest structural effects on the gut microbial community of mosquitos, while Sanders et al. [24] showed a strong influence of host species in the composition of the bacterial community of turtle ants. In a survey including insects from various trophic guilds and taxa, Yun et al. [25] suggested that host habitat, diet, developmental stage and phylogeny, all contributed to structuring insect gut microbiota.
Special attention has been given to the microbiota of insects which are vectors of important parasites. Sandfly gut microbiota can be modulated (i) by their feeding habits: larvae feeding on soil dead organic matter, while female adults get blood and sugar meals (plant-fed) and (ii) through interactions with transmitting parasites [26,27]. To date the gut microbiota of Phlebotomus sandflies has been studied using culture-dependent approaches [28,29], while recent studies on Lutzomyia longipalpis have used culture-independent approaches [30,31]. The gut bacterial community of sandflies is dominated by Proteobacteria and Firmicutes. Initial studies showed that certain members of the sandfly gut microbiota exerted a negative effect on the development of Leishmania [32,33]. This was challenged recently by Louradour et al. [34] who showed that the gut microbiota enhances the survival of Leishmania by modulating optimal osmotic conditions for promastigotes in the midgut. Further studies by Dey et al. [35] suggested that gut microbes of L. longipalpis are egested into the host skin alongside Leishmania triggering neutrophil infiltration and facilitating parasite establishment. To date most studies have looked into the gut microbiota of sandflies through the prism of its interaction with the transmitting parasite or as a pool of symbiotic organisms with the potential to be used in paratransgenesis control strategies. However, we still lack data about the role of phylogeny on structuring the gut microbiota of sandflies. This is especially true for sandflies of the genus Phlebotomus, whose associated microbiota has been studied using less informative culture-dependent methods.
We posit that host phylogeny is a main determinant of the composition of the gut microbiota of Phlebotomus sandflies. To test this hypothesis, we collected wild sandfly specimens (over a two-year period) from a rather isolated population in a sandfly endemic area. Morphologically identified Phlebotomus species were used to perform amplicon sequencing analysis of the 16S rRNA gene of bacteria and archaea, particularly for the later whose presence in insect symbiome is underscored [36]. Specimens analyzed were all Leishmania-free, as verified by quantitative polymerase chain reaction (q-PCR), allowing us to study the structural role of sandfly phylogeny on the gut microbiota in the absence of any other interacting factors (i.e., Leishmania presence).

2. Results

2.1. Morphological Identification of Phlebotomus Species

In total 1235 sandflies belonging to 10 different species were collected from the island of Leros, Greece during the two-year sampling period (2017–2018), 566 (45.8%) males and 669 (54.2%) females. Collected sandflies were assigned to the following species in descending prevalence: P. papatasi (51.5%), P. neglectus (19.8%), P. tobbi (17.1%), Sergentomyia minuta (3.3%) P. similis (2.7%), P. simici (2.2%), P. perfiliewi (1.8%), S. dentata (1.4%), P. mascitti (0.2%) and P. alexandri (0.1%). Upon homogenization 37, 15, 10 and 10 samples of P. papatasi, P. neglectus, P. tobbi and P. similis, constituting the most prevalent Phlebotomus species, were processed for DNA extraction. From the sandflies collected none was found positive to L. infantum according to our q-PCR test (Supplementary Materials Figure S1). All these samples we further processed for determination of their prokaryotic gut microbiome.

2.2. The Composition of the Bacterial and Archaeal Microbiome in Sand Flies

From the 16S rRNA amplicon sequencing of all samples we obtained 1,503,942 quality sequences (19,281 per sample, range 2200–82,059) which were assigned to 3762 and 22 operational taxonomic units (OTUs) for bacteria and archaea, respectively. Our sequencing effort provided adequate coverage of the microbial diversity on the gut of sandflies as suggested by the rarefaction curves, which reached a plateau in all samples tested (Supplementary Materials Figure S2).
We observed significant differences in the α-diversity of bacteria among species. Hence, significantly lower values (p < 0.05) of Shannon and inverse Simpson diversity indices were observed in P. papatasi and P. neglectus compared to the other two species (Figure 1). In contrast, we did not observe significant differences between the different Phlebotomus species in the other α-diversity indices obtained.
The gut bacterial community of Phlebotomus was dominated by Proteobacteria (55.8%, 39.8%–77.3%), mostly of the classes of α-proteobacteria (39.3%) and γ-proteobacteria (15.9%), Firmicutes (14%, 7.5%–18.9%), Tenericutes (12.9%, 0.01%–50.0%), Actinobacteria (10.9%, 1.5%–19.0%) and Bacteroidetes (3.0%, 0.5%–6.5%) (Figure 2). The large representation of α-proteobacteria and Tenericutes was associated, almost entirely, with the presence of endosymbiotic Wolbachia and Spiroplasma, respectively. Non endosymbiotic proteobacteria like Pseudomonas (3.2%), Acinetobacter (2.3%), Methylobacterium (1.1%), Firmicutes like Staphylococcus (3.5%) and Actinobacteria like Cutibacterium (2.2%) were also present at lower relative abundance (RA) (<5%). The core microbiota of Phlebotomus species (all OTUs that participated with at least 0.1% RA in at least the 50% of the samples analysed) was composed of 14 OTUs belonging to α-proteobacteria (Wolbachia, Sphingomonadaceae, Rhizobiaceae), γ-protebacteria (Pseudomonas, Acinetobacter), Firmicutes (Staphylococcus), Actinobacteria (Corynebacteriaceae, Cutibacterium, Streptococcus) and Tenericutes (Spiroplasma) (Supplementary Materials Figure S3).
At the β-diversity level, canonical correspondence analysis (CCA) showed that P. papatasi supported a largely different bacterial community which was clearly separated along CCA1 from the bacterial communities (p < 0.001) in the other three Phlebotomus species which formed a second group (Figure 3). Within this second group, samples from P. similis and P. neglectus were also significantly separated (p < 0.01) along CCA2 (explaining 26% of the variance), while a weaker but still significant separation (p < 0.05) was evident between P. tobbi and the other two genera.
We further looked for bacterial OTUs that showed a significant association with specific Phlebotomus species and, hence, were responsible for the observed distinct clustering of the bacterial communities (Figure 4). OTUs 1, 2 and 159, belonging to Spiroplasma, Wolbachia and Spiroplasma respectively, showed significantly higher RA (p < 0.001) in P. papatasi (Figure 5), particularly the first two which were the most abundant OTUs in our study. OTUs 3, 263, 284, 391 and 720, belonging to Wolbachia, Pseudomonadaceae, Xanthomonadaceae (284 and 391) and Brachybacterium, respectively, showed significantly higher RA (p < 0.05) in P. neglectus (Figure 5). OTUs 18 and 295, assigned to γ-proteobacteria and Rhodobacteriaceae, respectively, showed significantly higher RA (p < 0.05) in P. similis compared to the other Phlebotomus species (Figure 6). Finally, OTUs 12, 154, 204 και 255, assigned to Pseudomonas, Flavobacterium, Saprospiraceae and Acinetobacter, respectively, showed significantly higher RA (p < 0.05) in P. tobbi (Figure 5). We further observed several OTUs that were either disfavored in a certain Phlebotomus species or they were associated with more than one species. In the former case, OTUs 5, 37 and 45, assigned to Staphylococcus, Micrococcus and Corynebacteriaceae, respectively, showed significantly lower RA (p < 0.05) in P. papatasi compared to the other three species (Figure 6). On the other hand, OTUs 7, 11 and 78, belonging to Cutlbacterium, Rhizobiaceae and Pseudomonas, respectively, showed significantly higher RA (p < 0.05) in P. similis and P. neglectus, while OTUs 220 and 471, belonging to the Actinobacterial genera Brevibacterium and Corynebacterium, respectively, were associated (p < 0.05) with P. neglectus and P. similis (Figure 5).
We also explored the presence of archaea in Phlebotomus species. Archaeal OTUs showed low overall RA compared to bacteria. The five most abundant OTUs belonged to Nanoarchaeota (Woesearchaeia), Thaumarchaeota (Nitrososphaeraceae), and Eyrarchaeota (Methanobacterium) with RA ranging from 0.002% to 0.006% (Figure 6). Correlation testing showed that OTU_503, assigned to Nitrosphaeraceae, being the second most abundant archaeal OTU, showed an exclusive presence in P. similis. Similarly, OTU_1309 assigned to Methanobrevibacter, showed significantly higher RA in P. neglectus compared to P. papatasi and P. tobbi.

3. Discussion

In this study we determined the composition of the gut prokaryotic microbiota of Phlebotomus species and explored the role of host phylogeny in structuring the gut bacterial and archaeal community. The sandfly population collected from the island of Leros was composed of 10 species, with P. papatasi, P. neglectus and P. tobbi being the most dominant, followed by S. minuta and P. simils, representing the sandfly fauna on the island of Leros, Greece. Previous studies in Serbia, Albania and other countries in the East Mediterranean basin have shown similar sandfly species distribution with P. neglectus most often identified as the dominant species [7,8,9,10]. A recent study in the same study region in Greece reported P. neglectus as the predominant species, followed by S. minuta, P. tobbi, P. simici and P. similis [12], which were also detected in our study. Additionally, in the current study we report the presence of P. papatasi, P. perfiliewi, P. mascitti and P. alexandri on the island. Discrepancies between monitoring surveys of sandfly populations performed in the same region at different time periods is a rather common observation and could be attributed to: (i) different sampling periods (i.e., in [12] samples were collected only during 3 sampling days in one sandfly season vs. several occasions in two consecutive sandfly seasons) and (ii) different sampling microhabitats (P. neglectus and P. papatasi have been reported to have different optima of relative humidity (RH) and temperature (RH 50%–60%, 27–29 ℃ vs. RH 30%–40%, 25–27 ℃) [12]. In line with this, Alten et al., [7] showed that the population density of P. neglectus in a given region could vary substantially from year to year as a function of varying environmental conditions.
The gut microbiota of the Phlebotomus species was dominated by Proteobacteria, mostly α- and γ-proteobacteria, Firmicutes, Tenericutes and Actinobacteria. This is in line with previous studies in wild caught sandflies [28,29,37] and other insects [25,38].
Most previous microbiota studies on Phlebotomus have used culture-dependent approaches [28,29,39,40], with their well-documented limitations [41], and did not explore the factors shaping the gut microbial community of wild populations of Phlebotomus sandflies. Previous studies have highlighted, among other factors, the role of host phylogeny on the composition of insect gut microbial community [24,42,43]. In this context, we explored the role of host species on the diversity of the gut microbiota of a wild population of Phlebotomus sandflies. In order to focus on host genotype effects, we only considered in our study Leishmania-free specimens. We noted a clear differentiation of the gut microbiota at the α-diversity level between P. similis and P. tobbi, dominated by several low-abundance members, and P. neglectus and P. papatasi whose bacteriome was dominated by a few highly prevalent OTUs. Multivariate analysis at the β-diversity level showed that all Phlebotomus species carried distinct microbial communities, with P. papatasi showing a more distinct bacterial community compared to the other three species, which also differed to each other to a lower extent. These findings verify our initial hypothesis that, in the absence of Leishmania, host phylogeny has a significant structural role in gut bacterial community assembly of Phlebotomus species.
The distinct composition of the bacterial community of P. papatasi was driven by three OTUs belonging to Spiroplasma and Wolbachia. These are common facultative endosymbionts of Phlebotomus [28,44], Lutzomyia [30,31] and other insects; it is now estimated that Wolbachia and Spiroplasma endosymbionts are present in up to 30% of all insects [45]. They are maternally inherited, and they affect host ecology, physiology and fitness [15,46,47]. Their extensive co-presence in P. papatasi contrasts previous studies in sandflies [44] and ants [48] which have suggested that their co-detection in the same host is a rather infrequent event. Similarly, the microbiota of P. neglectus was overwhelmed by another Wolbachia, shared with P. tobbi, and a Staphylococcus shared with P. tobbi and P. similis. Staphylococcus have been reported as common dwellers of the gut of L. longipalpis [49,50], L. evansi [37] and P. papatasi [28], although here they were specifically absent from this species. Staphylococcus have been incriminated as pathogens of various organisms including insects [51], although their role in the biology of insect hosts remains unknown. Overall, the high abundance of Wolbachia and Spiroplasma in the female Phlebotomus fauna in the study area might have serious, positive or negative, implications for the sandfly population depending on the type of effect these endosymbionts impose on their host. Further investigations will focus on the role of these symbionts on the biology of the Phlebotomus species considering the contrasting evidence for the role of endosymbiotic bacteria on the infectivity and survival of Leishmania [27,52].
The gut of P. tobbi and P. similis were co-colonized by bacteria showing low relative abundance (1%–4%) and belonged to Cutibacterium, Rhizobiaceae and Pseudomonas. Cutibacterium are common members of the insect gut microbiota of Bactocera oleae at early developmental stages [53]. They are mostly known as early colonizers of infants where they function as lactate-consumer and propionate-producer [54]. Their presence in the gut of sandflies might be associated with feeding and digestion mechanisms. Rhizobiaceae encompass a range of bacterial genera which are prevalent in soil and plant tissues, hence their presence in insect gut has been associated with plant feeding lifestyles [55,56]. Pseudomonas are common dwellers of the gut of wild caught sandflies [4,30,57]. Their common presence in larvae and adults of sandflies was attributed to acquisition by soil during the soil feeding larvae stage [37]. However recent studies with laboratory-fed sandflies showed that Pseudomonadaceae constitute significant members of the gut microbiota of L. longipalpis with their abundance depending on the meal, plant sugar or blood [26,27].
With the applied methods used we were also able to determine the archaeal community in sandflies gut. The contribution of archaea in the prokaryotic microbiota of sandflies was low (<0.01%), compared to termites (3%), whose archaeal symbiotic community has been thoroughly studied [58]. This underrepresentation of archaea in most studies, including ours, could be attributed either to the use of universal bacterial primers which was recently shown to be biased towards bacteria in the human gut microbiota [59], the lack of well populated 16S rRNA reference databases or indeed the low abundance of archaea in the studied samples [60]. Considering the lack of previous studies on sandflies, we provide first evidence for the presence of archaea belonging to Nanoarchaeota, methanogenic Euryarchaeota and Thaumarchaeota in the gut of Phlebotomus. Methanogenic Euryarchaeota belonging to Methanobrevibacter, Methanobactarium and Methanosarcina, all being obligate anaerobes, are common dwellers of the gut of termites [32,58] and Coleoptera [61], producing methane from CO2/H2 or formate. Their presence might suggest the establishment of anaerobic zones in the gut of Phlebotomus. Thaumarchaeota of the family Nitrososphaeraceae have been also detected in the gut of beetles [61]. They are ubiquitous in soil, plants and water ecosystems responsible for the oxidation of ammonium to nitrite [62]. However, their exact role in the gut of sandflies warrants further investigation. However, the most abundant archaeal OTU belonged to Nanoarchaeota, previously detected in the human lung [59], but reported for the first time in an insect gut. Nanoarchaeota are obligate symbionts of Crenarchaeota [63,64]. They are characterized by a reduced genome, hence relying on their host for central cellular biosyntheses [65]. Their detection in diverse aquatic and terrestrial ecosystems [66,67,68,69] suggests a rapidly expanding range of hosts and their dispersal in new environments. Nanoarchaeota could exert deleterious effects on their host like N. equitans, which prevents host replication [70], while terrestrial members do not seem to be deleterious to their host and might be beneficial under certain conditions [64].

4. Materials and Methods

4.1. Samples Collection

Sandflies were collected from the island of Leros in the south-eastern area of the Aegean Sea, Greece (Latitude: 37°08′60.00″ N, Longitude: 26°50′59.99″ E), a previously reported habitat of several Phlebotomus species [12]. The total area of the island is 54 km² and the natural environment is typical Mediterranean (Supplementary Materials Figure S4).
Sand fly collection was performed by using Centre for Disease Control (CDC) miniature light traps and BG Sentinel traps at regular intervals for two sandfly activity seasons (April to October of 2017 and 2018). Traps were placed in domestic and peridomestic environments nearby animal sheds with cattle, goats, sheep, or poultry to maximize trapping capacity [71]. The traps were operated all night and the sandflies were collected early in the morning, sorted, and placed in Eppendorf tubes in 80% ethanol before being transported to the laboratory for downstream processing and analyses.

4.2. Species-Level Identification of Phlebotomus Sandflies

The head and the rear part of the abdomen of all collected sandflies were mounted on permanent microscope slides and species identification was carried out based on the morphology of the pharynx, male genitalia or female spermathecae [72,73]. The stomach and the gut of each dissected specimen were then stored together at −80 °C for DNA extraction.

4.3. DNA Extraction

Upon species identification female sandflies of the same species, collected at the same time point from the same trap were homogenized in groups of 5. DNA was extracted from homogenized samples using the DNeasy Blood and Tissue Kit (Qiagen GmbH, Hilden, Germany) following the manufacturers’ protocol. The size and the integrity of the DNA extracted was verified by agarose gel electrophoresis (1%) and DNA was quantified using a Qubit fluorometer with a Quant-iT HS double-stranded DNA (dsDNA) assay kit (Invitrogen, Carlsbad, CA, USA).

4.4. Quantitative Polymerase Chain Reaction (q-PCR) Detection of Leishmania spp. in Sandfly Samples

All samples were tested for the presence of L. infantum with a TaqMan real time qPCR assay, targeting a 120 bp fragment of the kinetoplast minicircle DNA as reported before [74].

4.5. PCR Amplification of the 16S rRNA Gene and Amplicon Sequencing Analysis of the Microbiome

The composition of the bacterial and archaeal community (co-amplified with the primer sets and protocol used) was determined with amplicon sequencing of the 16S rRNA gene via HiSeq Illumina Rapid Mode 2 × 250 bp paired-end reads (Illumina Inc., San Diego, CA, USA) in the DNA Sequencing Center, Department of Biology, Brigham Young University (GSC-BYU, Provo, UT, US). Bacterial and archaeal 16S rRNA genes were amplified with the primer set 515f-806r [75,76], targeting the V4 region of the 16S rRNA gene, following the protocol of the Earth Microbiome Project [77]. For all PCR amplifications, the Q5® High-Fidelity DNA Polymerase (NEB, Massachusetts, USA) was used. All samples were initially amplified (28 amplification cycles) using the domain-specific primers mentioned above, followed by a PCR (7 amplification cycles) using the same primers but this time the forward primer 515r (5’-TTXXXXXXXXXGTGTGYCAGCMGCCGCGGTAA-3’) carried 5’ extensions comprising of linkers (italics) and indexes (underlined) used for samples barcoding for multiplex sequencing. PCR conditions are listed in Supplementary Materials Table S1.

4.6. Bioinformatic Analysis of Amplicon Sequencing Data

The amplicon sequencing data were analyzed as described by Katsoula et al., [78]. Briefly, the raw sequence data were demultiplexed with Flexbar v3.0 [79]. The reads were quality trimmed with Trimmomatic v0.32 [80] for a length cutoff of 150 bp and were assembled with FLASH v1.2.8 [81] using the default parameters unless otherwise stated. OTU calling at 97% identities was performed with the UPARSE v10.0.240 software [82]. Chimeric sequences were identified with the UCHIME v4.2 software [83] using the RDP Gold database vMicrobiomeutil-r20110519. Sequence classification was performed with Lambda v0.9.1 [84] against the Silva v128 small ribosomal subunit database [85]. Misclassified sequences were removed from downstream analysis.

4.7. Statistical Analysis

The OTU matrix of bacteria was used to assess the effect of host genotype on the α- and β-diversity. The impact on the α-diversity was determined via calculation of the diversity indices richness (S), Fisher alpha, inverse Simpson, Shannon [86] and Pielou’s evenness [87] and the data were subjected to one-way analysis of variance (ANOVA). Differences in the β-diversity of the bacterial community between different Phlebotomus species were modelled with CCA. We further used pairwise permutational multivariate analysis of variance (PERMANOVA) to identify significant differences in the structure of the bacterial communities between different sandfly species [88]. All statistical analyses were performed with the R v3.5.2 software [89]. The data were submitted to the Sequence Read Archive of the National Center for Biotechnological Information (NCBI) with bioproject accession number PRJNA630369.

5. Conclusions

The gut prokaryotic microbiota of a wild, Leishmania-free, population of Phlebotomus sandflies composed of four main species was structured according to host genotype. Endosymbiotic Wolbachia and Spiroplasma seemed to dominate the bacterial community of P. papatasi and P. neglectus, while P. tobbi and P. similis supported a more diverse bacterial community composed of several less-abundant bacteria. Archaea were minor members of the prokaryotic microbiota of sandflies with methanogenic Euryarchaeota, nitrifying Thaumarchaeota and obligate symbiotic Nanoarchaeota being most prevalent. Our study provides the first high-resolution analysis of the bacterial gut microbiota in Phlebotomus, which is structured (in the absence of Leishmania) according to host phylogeny, and reports pioneering evidence for the presence of archaea in the gut of sandflies. Further studies will explore the role of Leishmania in combination with biogeography and host genotype on the assembly of the gut microbial community of wild populations of Phlebotomus sandfies.

Supplementary Materials

The following are available online at https://www.mdpi.com/2076-0817/9/6/428/s1: Figure S1: The standard curve (a), amplification curve (b) and dissociation curves (c) of the quantitative polymerase chain reaction (q-PCR) measurements employed for the detection and quantification of Leishmania infantum in sandfly DNA samples. Only standard solutions used for the preparation of the standard curve were positive in amplification as indicated in the amplification and dissociation curves, Figure S2: Rarefaction curves indicating the depth of diversity coverage of the bacterial community in the sandfly samples analysed, Figure S3: Heatmap representing the OTUs composing the core microbiota of Phlebotomus species analysed in the current study. OTUs were characterized as core microbiota members if they participated with at least 0.1% in at least the 50% of the samples derived from the insect of interest, Figure S4: The geographic location of the study area where sandfly samples were collected, Table S1: PCR reagents and thermocycling conditions used for amplicon sequencing analysis.

Author Contributions

S.S. and D.G.K.: Conceptualization and supervision; D.G.K.: Writing—original draft; S.S., S.V., P.L., A.S.: Writing—review and editing; C.P., P.A.K., P.L., S.V. and A.S.: Investigation and Methodology; S.V.: Validation and Visualization. All authors have read and agreed to the published version of the manuscript.

Funding

No funding for this work is reported.

Acknowledgments

We would like to specifically acknowledge Michalis Kontrafouris, DVM in Leros island who supported our sandfly collection.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Pappa, A.; Konstantinou, G.; Pavlidou, V.; Antoniadis, A. Sandfly fever virus outbreak in Cyprus. Clin. Microbiol. Infec. 2006, 12, 192–194. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  2. Alwassouf, S.; Christodoulou, V.; Bichaud, L.; Ntais, P.; Mazeris, A.; Antoniou, M.; Charel, R.N. Seroprevalence of sand fly-borne phleboviruses belonging to three serocomplexes (sand fly fever Naples, sand fly fever Sicilian and Salehabad) in dogs from Greece and Cyprus using neutralization test. PLoS Negl. Trop Dis. 2016, 10, e0005063. [Google Scholar] [CrossRef] [PubMed]
  3. Billeter, S.A.; Levy, M.G.; Chomel, B.B.; Breitschwerdt, E.B. Vector transmission of Bartonella species with emphasis on the potential for tick transmission. Med. Vet. Entomol. 2008, 22, 1–15. [Google Scholar] [CrossRef] [PubMed]
  4. Akhoundi, M.; Kuhls, K.; Cannet, A.; Votypka, J.; Marty, P.; Delaunay, P.; Sereno, D. A historical overview of the classification, evolution, and dispersion of Leishmania parasites and sandflies. PLOS Negl. Trop Dis 2016, 10, e0004770. [Google Scholar] [CrossRef] [PubMed]
  5. Antoniou, M.; Gramiccia, M.; Molina, R.; Dvorak, V.; Volf, P. The role of indigenous phlebotomine sandflies and mammals in the spreading of leishmaniasis agents in the Mediterranean region. Euro Surveil. 2013, 18, e20540. [Google Scholar] [CrossRef] [Green Version]
  6. Bates, P.A.; Depquit, J.; Galati, E.A.B.; Kamhawi, S.; Maroli, M.; McDowell, M.A.; Picado, A.; Ready, P.D.; Salomon, O.D.; Shaw, J.J.; et al. Recent advances in phlebotomine sandfly research related to leishmaniasis control. Paras Vect 2015, 8, 131. [Google Scholar] [CrossRef] [Green Version]
  7. Alten, B.; Maia, C.; Afonso, M.O.; Campino, L.; Jiménez, M.; González, E.; Molina, R.; Banuls, A.-L.; Prudhomme, J.; Vergnes, B.; et al. Seasonal dynamics of Phlebotomine sand fly species proven vectors of Mediterranean Leishmaniasis caused by Leishmania infantum. PLoS Negl. Trop Dis. 2016, 10, e0004458. [Google Scholar] [CrossRef]
  8. Velo, E.; Bongiorno, G.; Kadriaj, P.; Myrseli, T.; Crilly, J.; Lika, A.; Mersini, K.; Di Muccio, T.; Bino, S.; Gramiccia, M.; et al. The current status of phlebotomine sand flies in Albania and incrimination of Phlebotomus neglectus (Diptera, Psychodidae) as the main vector of Leishmania infantum. PLoS ONE 2017, 12, e0179118. [Google Scholar] [CrossRef]
  9. Vaselek, S.; Ayhan, N.; Oguz, G.; Erisoz Kasap, O.; Savić, S.; Di Muccio, T.; Gradoni, L.; Ozbel, Y.; Alten, B.; Petric, D. Sand fly and Leishmania spp. survey in Vojvodina (Serbia): First detection of Leishmania infantum DNA in sand flies and the first record of Phlebotomus (Transphlebotomus) mascittii Grassi, 1908. Paras. Vect 2017, 10, 444. [Google Scholar] [CrossRef] [Green Version]
  10. Vaselek, S.; Dvorak, V.; Hlavackova, K.; Ayhan, N.; Halada, P.; Oguz, G.; Ivovic, V.; Ozbel, Y.; Charrel, R.N.; Alten, B.; et al. A survey of sand flies (Diptera, Phlebotominae) along recurrent transit routes in Serbia. Acta Tropica. 2019, 197, 105063. [Google Scholar] [CrossRef]
  11. Xanthopoulou, K.; Anagnostou, V.; Ivovic, V.; Djurkovic-Djiakovic, O.; Rogozi, E.; Sotiraki, S.; Papa, A. Distribution of sandflies (Diptera, Psychodidae) in two Ionian islands and northern Greece. Vector-Borne Zoon Dis. 2011, 11, 1591–1594. [Google Scholar] [CrossRef] [PubMed]
  12. Tsirigotakis, N.; Pavlou, C.; Christodoulou, V.; Dokianakis, E.; Kourouniotis, C.; Alten, B.; Antoniou, M. Phlebotomine sand flies (Diptera: Psychodidae) in the Greek Aegean Islands: Ecological approaches. Paras Vect. 2018, 11, 97. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  13. Douglas, A.E. Nutritional interactions in insect-microbial symbioses: Aphids and their symbiotic bacteria Buchnera. Ann. Rev. Entomol. 1998, 43, 17–37. [Google Scholar] [CrossRef] [Green Version]
  14. Dedeine, F.; Vavre, F.; Fleury, F.; Loppin, B.; Hochberg, M.; Bouletreau, M. Removing symbiotic Wolbachia bacteria specifically inhibits oogenesis in a parasitic wasp. Proc. Natl. Acad. Sci. USA 2001, 98, 6247–6252. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  15. Harumoto, T.; Lemaitre, B. Male-killing toxin in a Drosophila bacterial symbiont. Nature 2018, 557, 252–255. [Google Scholar] [CrossRef] [PubMed]
  16. Paredes, J.C.; Herren, J.K.; Schüpfer, F.; Lemaitre, B. The role of lipid competition for endosymbiont-mediated protection against parasitoid wasps in Drosophila. mBio 2016, 7, e01006-16. [Google Scholar] [CrossRef] [Green Version]
  17. Moreira, L.A.; Iturbe-Ormaetxe, I.; Jeffery, J.A.; Lu, G.; Pyke, A.T.; Hedges, L.M.; Rocha, B.C.; Hall-Mendelin, S.; Day, A.; Riegler, M.; et al. A Wolbachia symbiont in Aedes aegypti limits infection with dengue, Chikungunya, and Plasmodium. Cell 2009, 139, 1268–1278. [Google Scholar] [CrossRef] [Green Version]
  18. Guerrero, R.; Margullis, L.; Berlanga, M. Symbiogenesis: The holobiont as a unit of evolution. Int. Microbiol 2013, 16, 133–143. [Google Scholar]
  19. Cheng, D.; Guo, Z.; Riegler, M.; Xi, Z.; Liang, G.; Xu, Y. Gut symbiont enhances insecticide resistance in a significant pest, the oriental fruit fly Bactrocera dorsalis (Hendel). Microbiome 2017, 5, 13. [Google Scholar] [CrossRef] [Green Version]
  20. Schwarz, R.S.; Teixeira, E.W.; Tauber, J.P.; Birke, J.M.; Martins, M.F.; Fonseca, I.; Evans, J.D. Honey bee colonies act as reservoirs for two Spiroplasma facultative symbionts and incur complex, multiyear infection dynamics. Microbiol. Open 2014, 3, 341–355. [Google Scholar] [CrossRef]
  21. Mason, C.J.; Ray, S.; Shikano, I.; Peiffer, M.; Jones, A.G.; Luthe, D.S.; Hoover, K.; Felton, G.W. Plant defenses interact with insect enteric bacteria by initiating a leaky gut syndrome. Proc. Natl. Acad. Sci. USA 2019, 116, 15991–15996. [Google Scholar] [CrossRef] [Green Version]
  22. Kolasa, M.; Scibior, R.; Mazur, M.A.; Kubisz, D.; Dudek, K.; Kajtoch, L. How hosts taxonomy, trophy, and endosymbionts shape microbiome diversity in beetles. Microb. Ecol. 2018, 78, 995–1013. [Google Scholar] [CrossRef] [Green Version]
  23. Muturi, E.J.; Lagos-Kutz, D.; Dunlap, C.; Ramirez, J.L.; Rooney, A.P. Mosquito microbiota cluster by host sampling location. Paras Vect 2018, 11, 468. [Google Scholar] [CrossRef] [Green Version]
  24. Sanders, J.G.; Powell, S.; Kronauer, D.J.; Vasconcelos, H.L.; Frederickson, M.E.; Pierce, N.E. Stability and phylogenetic correlation in gut microbiota: Lessons from ants and apes. Molec Ecol. 2014, 23, 1268–1283. [Google Scholar] [CrossRef]
  25. Yun, J.-H.; Roh, S.W.; Whon, T.W.; Jung, M.-J.; Kim, M.-S.; Park, D.-S.; Yoon, C.; Nam, Y.-D.; Kim, Y.-J.; Choi, J.-H.; et al. Insect gut bacterial diversity determined by environmental habitat, diet, developmental stage, and phylogeny of host. Appl. Environ. Microbiol. 2014, 80, 5254–5264. [Google Scholar] [CrossRef] [Green Version]
  26. Pires, A.C.A.M.; Villegas, L.E.M.; Campolina, D.B.; Orfano, A.S.; Pimenta, P.F.P.; Secundino, N.F.C. Bacterial diversity of wild-caught Lutzomyia longipalpis (a vector of zoonotic visceral leishmaniasis in Brazil) under distinct physiological conditions by metagenomics analysis. Paras Vect 2017, 10, 627. [Google Scholar] [CrossRef] [Green Version]
  27. Kelly, P.H.; Bahr, S.M.; Serafim, T.D.; Ajami, N.J.; Petrosino, J.F.; Meneses, C.; Kirby, J.R.; Valenzuela, J.G.; Kamhawi, S.; Wilson, M.E. The gut microbiome of the vector Lutzomyia longipalpis is essential for survival of Leishmania infantum. mBio 2017, 8, e01121-16. [Google Scholar] [CrossRef] [Green Version]
  28. Fraihi, W.; Fares, W.; Perrin, P.; Dorkeld, F.; Sereno, D.; Barhoumi, V.; Sbissi, I.; Cherni, S.; Chelbi, I.; Durvasula, R.; et al. An integrated overview of the midgut bacterial flora composition of Phlebotomus perniciosus, a vector of zoonotic visceral leishmaniasis in the Western Mediterranean Basin. PLoS Negl. Trop Dis. 2017, 11, e0005484. [Google Scholar] [CrossRef]
  29. Karimian, F.; Vatandoost, H.; Rassi, Y.; Maleki-Ravasan, N.; Mohebali, M.; Shirazi, M.H.; Koosha, M.; Choubdar, N.; Oshagi, M.A. Aerobic midgut microbiota of sand fly vectors of zoonotic visceral leishmaniasis from northern Iran, a step toward finding potential paratransgenic candidates. Paras Vect. 2019, 12, 10. [Google Scholar] [CrossRef]
  30. Li, K.; Chen, H.; Jiang, J.; Li, X.; Xu, J.; Ma, Y. Diversity of bacteriome associated with Phlebotomus chinensis (Diptera: Psychodidae) sand flies in two wild populations from China. Sci. Rep. 2016, 6, 36406. [Google Scholar] [CrossRef] [Green Version]
  31. Vivero, R.J.; Villegas-Plazas, M.; Cadavid-Restrepo, G.E.; Herrera, C.X.M.; Uribe, S.I.; Junca, H. Wild specimens of sand fly phlebotomine Lutzomyia evansi, vector of leishmaniasis, show high abundance of Methylobacterium and natural carriage of Wolbachia and Cardinium types in the midgut microbiome. Sci. Rep. 2019, 9, 17746. [Google Scholar] [CrossRef] [Green Version]
  32. Santana, R.H.; Pires Catao, E.C.; Cardoso Lopes, F.A.; Constantino, R.; Chaves Barreto, C.; Kruger, R.H. The gut microbiota of workers of the litter-feeding termite Syntermes wheeleri (Termitidae: Syntermitinae): Archaeal, bacterial, and fungal communities. Microb. Ecol. 2015, 70, 545–556. [Google Scholar] [CrossRef]
  33. Schlein, Y.; Polacheck, I.; Yuval, B. Mycoses, bacterial infections and antibacterial activity in sandflies (Psychodidae) and their possible role in the transmission of leishmaniasis. Parasitology 1985, 90, 57–66. [Google Scholar] [CrossRef]
  34. Louradour, I.; Monteiro, C.C.; Inbar, E.; Ghosh, K.; Merkhofer, R.; Lawyer, P.; Paun, A.; Smelkinson, M.; Secundino, N.; Lewis, M.; et al. The midgut microbiota plays an essential role in sand fly vector competence for Leishmania major. Cell Microbiol. 2017, 19, e12755. [Google Scholar] [CrossRef] [Green Version]
  35. Dey, R.; Joshi, A.B.; Oliveira, F.; Pereira, L.; Guimaraes-Costa, A.B.; Serafim, T.D.; de Castro, W.; Coutinho-Abreu, I.V.; Bhattacharya, P.; Townsend, S.; et al. Microbes egested during bites of infected sand flies augment severity of Leishmaniasis via inflammasome-derived IL-1b. Cell Host Microbe 2018, 23, 134–143. [Google Scholar] [CrossRef] [Green Version]
  36. Gurung, K.; Wertheim, B.; Salles, J.F. The microbiome of pest insects: It is not just bacteria. Entomol Experim. Appl. 2019, 167, 156–170. [Google Scholar] [CrossRef] [Green Version]
  37. Vivero, J.; Jaramillo, G.; Cadavid-Restrepo, G.; Soto, S.; Herrera, C. Structural differences in gut bacteria communities in developmental stages of natural populations of Lutzomyia evansi from Colombia’s Caribbean coast. Paras Vect. 2016, 13, 496. [Google Scholar] [CrossRef] [Green Version]
  38. Paniagua Voirol, L.R.; Frago, E.; Kaltenpoth, M.; Hilker, M.; Fatouros, N.E. Bacterial symbionts in Lepidoptera: Their diversity, transmission, and impact on the host. Front. Microbiol. 2018, 9, 556. [Google Scholar] [CrossRef]
  39. Maleki-Ravasan, N.; Oshaghi, M.A.; Afshar, D.; Arandian, M.H.; Hajikhani, S.; Akhavan, A.A.; Yakhchali, B.; Shirazi, M.H.; Rassi, Y.; Jafari, R.; et al. Aerobic bacterial flora of biotic and abiotic compartments of a hyperendemic Zoonotic Cutaneous Leishmaniasis (ZCL) focus. Paras Vect. 2015, 8, 63. [Google Scholar] [CrossRef] [Green Version]
  40. Mukhopadhyay, J.; Braig, H.R.; Rowton, E.D.; Ghosh, K. Naturally occurring culturable aerobic gut flora of adult Phlebotomus papatasi, vector of Leishmania major in the Old World. PLoS ONE 2012, 7, e35748. [Google Scholar] [CrossRef] [Green Version]
  41. Lloyd, K.G.; Steen, A.D.; Ladau, J.; Yin, J.; Crosby, L. Phylogenetically novel uncultured microbial cells dominate Earth microbiomes. mSystems 2018, 3, e00055-18. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  42. Colman, D.R.; Toolson, E.C.; Takacs-Vesbach, C.D. Do diet and taxonomy influence insect gut bacterial communities? Mol. Ecol. 2012, 21, 5124–5137. [Google Scholar] [CrossRef] [PubMed]
  43. McLean, A.H.C.; Godfrey, H.C.J.; Ellers, J.; Henry, L.M. Host relatedness influences the composition of aphid Microbiomes. Environ. Microbiol. Rep. 2019, 11, 808–816. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  44. Karatepe, B.; Aksoy, S.; Karatepe, M. Investigation of Wolbachia spp. and Spiroplasma spp. in Phlebotomus species by molecular methods. Sci. Rep. 2018, 8, 10616. [Google Scholar]
  45. Duron, O.; Bouchon, D.; Boutin, S.; Bellamy, L.; Zhou, L.; Engelstadter, J.; Hurst, G.D. The diversity of reproductive parasites among arthropods: Wolbachia do not walk alone. BMC Biol. 2008, 6, 27. [Google Scholar] [CrossRef] [Green Version]
  46. Hedges, L.M.; Brownlie, J.C.; O’Neill, S.L.; Johnson, K.N. Wolbachia and virus protection in insects. Science 2008, 322, 702. [Google Scholar] [CrossRef]
  47. Xie, J.; Vilchez, I.; Mateos, M. Spiroplasma bacteria enhance survival of Drosophila hydei attacked by the parasitic wasp Leptopilina heterotoma. PLoS ONE 2010, 5, e12149. [Google Scholar] [CrossRef]
  48. Ballinger, M.J.; Moore, L.D.; Perlman, S.J. Evolution and diversity of inherited Spiroplasma symbionts in Myrmica ants. Appl. Environ. Microbiol. 2018, 84, e02299-17. [Google Scholar] [CrossRef] [Green Version]
  49. Gupta, A.K.; Rastogi, G.; Nayduch, D.; Sawant, S.S.; Bhonde, R.R.; Shouche, Y.S. Molecular phylogenetic profiling of gut-associated bacteria in larvae and adults of flesh flies. Med. Vet. Entomol. 2014, 28, 345–354. [Google Scholar] [CrossRef]
  50. Ngo, C.T.; Romano-Bertrand, S.; Manguin, S.; Jumas-Bilak, E. Diversity of the bacterial microbiota of anopheles mosquitoes from Binh Phuoc Province, Vietnam. Front. Microbiol. 2016, 7, 2095. [Google Scholar] [CrossRef] [Green Version]
  51. Haloi, K.; Kankana, M.; Nath, R.; Devi, D. Characterization and pathogenicity assessment of gut-associated microbes of muga silkworm Antheraea assamensis Helfer (Lepidoptera: Saturniidae). J. Invertebr. Pathol. 2016, 138, 73–85. [Google Scholar] [CrossRef]
  52. Sant’Anna, M.R.; Diaz-Albiter, H.; Aguiar-Martins, K.; Al Salem, W.S.; Cavalcante, R.R.; Dilon, V.M.; Bates, P.A.; Genta, F.A.; Dilon, R.J. Colonisation resistance in the sand fly gut: Leishmania protects Lutzomyia longipalpis from bacterial infection. Paras Vect. 2014, 7, 329. [Google Scholar] [CrossRef] [Green Version]
  53. Koskinioti, P.; Ras, E.; Augustinos, A.; Tsiamis, G.; Beukeboom, L.W.; Caceres, C.; Bourtzis, K. The effects of geographic origin and antibiotic treatment on the gut symbiotic communities of Bactrocera oleae populations. Entomol. Experim. Appl. 2019, 167, 197–208. [Google Scholar] [CrossRef] [Green Version]
  54. Martin, V.N.R.; Schwab, C.; Krych, L.; Voney, E.; Geirnaert, A.; Braegger, C.; Lacroix, C. Colonization of Cutibacterium avidum during infant gut microbiota establishment. FEMS Microbiol. Ecol. 2019, 95, fiy215. [Google Scholar]
  55. Jones, A.G.; Mason, C.J.; Felton, G.W.; Hoover, K. Host plant and population source drive diversity of microbial gut communities in two polyphagous insects. Sci. Rep. 2019, 9, 2792. [Google Scholar] [CrossRef] [Green Version]
  56. Muratore, M.; Prather, C.; Sun, Y. The gut bacterial communities across six grasshopper species from a coastal tallgrass prairie. PLoS ONE 2020, 15, e0228406. [Google Scholar] [CrossRef] [Green Version]
  57. Telleria, E.L.; Martins-da-Silva, A.; Tempone, A.J.; Traub-Csekö, Y.M. Leishmania, microbiota and sand fly immunity. Parasitology 2018, 145, 1336–1353. [Google Scholar] [CrossRef] [Green Version]
  58. Rahman, N.A.; Parks, D.H.; Willner, D.L.; Engelbrektson, A.L.; Goffredi, S.K.; Warnecke, F.; Scheffrahn, R.H.; Hugenholtz, P. A molecular survey of Australian and North American termite genera indicates that vertical inheritance is the primary force shaping termite gut microbiomes. Microbiome 2015, 3, 5. [Google Scholar] [CrossRef] [Green Version]
  59. Koskinen, K.; Pausan, M.-R.; Perras, A.K.; Beck, M.; Bang, C.; Mora, M.; Schilhabel, A.; Schmitz, R.; Moissl-Eichinger, C. First insights into the diverse human archaeome: Specific detection of Archaea in the gastrointestinal tract, lung, and nose and on skin. mBio 2017, 8, e00824-17. [Google Scholar] [CrossRef] [Green Version]
  60. Pausan, M.R.; Csorba, C.; Singer, G.; Till, H.; Schöpf, V.; Santigli, E.; Klug, B.; Högenauer, C.; Blohs, M.; Moissl-Eichinger, C. Exploring the archaeome: Detection of archaeal signatures in the human body. Front. Microbiol. 2019, 10, 2796. [Google Scholar] [CrossRef] [Green Version]
  61. Ziganshina, E.; Mohammed, W.; Shagimardanova, E.I.; Vankov, P.Y.; Gogoleva, N.E.; Ziganshin, A.M. Fungal, bacterial, and archaeal diversity in the digestive tract of several beetle larvae (Coleoptera). Biomed. Res. Intern. 2018. [Google Scholar] [CrossRef]
  62. Tourna, M.; Stieglmeier, M.; Spang, A.; Koenneke, M.; Schintlmeister, A.; Urich, T.; Engel, M.; Schloter, M.; Wagner, M.; Richter, A.; et al. Nitrososphaera viennensis, an ammonia oxidizing archaeon from soil. Proc. Natl. Acad. Sci. USA 2011, 108, 8420–8425. [Google Scholar] [CrossRef] [Green Version]
  63. Huber, H.; Hohn, M.J.; Rachel, R.; Fuchs, T.; Wimmer, V.C.; Stetter, K.O. A new phylum of Archaea represented by a nanosized hyperthermophilic symbiont. Nature 2002, 417, 63–67. [Google Scholar] [CrossRef]
  64. Jarett, J.K.; Nayfach, S.; Podar, M.L.; Inskeep, W.P.; Ivanova, N.N.; Munson-McGee, J.; Schulz, F.; Young, M.; Jay, Z.J.; Beam, J.P.; et al. Single-cell genomics of co-sorted Nanoarchaeota suggests novel putative host associations and diversification of proteins involved in symbiosis. Microbiome 2018, 6, 161. [Google Scholar] [CrossRef] [Green Version]
  65. Waters, E.; Hohn, M.J.; Ahel, I.; Graham, D.E.; Adams, M.D.; Barnstead, M.; Beeson, K.Y.; Bibbs, L.; Bolanos, R.; Keller, M.; et al. The genome of Nanoarchaeum equitans: Insights into early archaeal evolution and derived parasitism. Proc. Natl. Acad. Sci. USA 2003, 100, 12984–12988. [Google Scholar] [CrossRef] [Green Version]
  66. McCliment, E.A.; Voglesonger, K.M.; O’Day, P.A.; Dunn, E.E.; Holloway, J.R.; Cary, C. Colonization of nascent, deep-sea hydrothermal vents by a novel Archaeal and Nanoarchaeal assemblage. Environ. Microbiol. 2006, 8, 114–125. [Google Scholar] [CrossRef]
  67. Casanueva, A.; Galada, N.; Baker, G.C.; Grant, W.D.; Heaphy, S.; Jones, B.; Yanhe, M.; Ventosa, A.; Blamey, J.; Cowan, D.A. Nanoarchaeal 16S rRNA gene sequences are widely dispersed in hyperthermophilic and mesophilic halophilic environments. Extremophiles 2008, 12, 651–656. [Google Scholar] [CrossRef]
  68. Clingenpeel, S.; Kan, J.; Macur, R.E.; Woyke, T.; Lovalvo, D.; Varley, J.; Inskeep, W.P.; Nealson, K.; McDermott, T.R. Yellowstone lake nanoarchaeota. Front. Microbiol. 2013, 4, 274. [Google Scholar] [CrossRef] [Green Version]
  69. Wurch, L.; Giannone, R.J.; Belisle, B.S.; Swift, C.; Utturkar, S.; Hettich, R.L.; Reysenbach, A.L.; Podar, M. Genomics informed isolation and characterization of a symbiotic Nanoarchaeota system from a terrestrial geothermal environment. Nat. Commun. 2016, 7, 12115. [Google Scholar] [CrossRef] [Green Version]
  70. Giannone, R.J.; Wurch, L.L.; Heimerl, T.; Martin, S.; Yang, Z.; Huber, H.; Rachel, R.; Hettich, R.L.; Podar, M. Life on the edge: Functional genomic response of Ignicoccus hospitalis to the presence of Nanoarchaeum equitans. ISME J. 2015, 9, 101–114. [Google Scholar] [CrossRef] [Green Version]
  71. Feliciangeli, M.D. Natural breeding places of phlebotomine sandflies. Med. Vet. Entomol. 2004, 18, 71–80. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  72. Theodor, O. Psychodidae-Phlebotominae; Lindner, E., Ed.; Schweizerbart’sche Verlagsbuchhandlung: Stuttgart, Germany, 1958. [Google Scholar]
  73. Killick-Kendrick, R.; Tang, Y.; Killick-kendrick, M.; Sang, D.K.; Sirdar, M.K.; Ke, L.; Ashford, R.W.; Schorscher, J.; Johnson, R.W. The identification of female sandflies, of the subgenus Larroussius by the morphology of the spermathecal ducts. Parasitologia 1991, 33, 335–347. [Google Scholar]
  74. Francino, O.; Altet, L.; Sánchez-Robert, E.; Rodriguez, A.; Solano-Gallego, L.; Alberola, J.; Ferrer, L.; Sanchez, A.; Roura, X. Advantages of real-time PCR assay for diagnosis and monitoring of canine leishmaniosis. Vet. Parasitol. 2006, 137, 214–221. [Google Scholar] [CrossRef] [PubMed]
  75. Caporaso, J.G.; Lauber, C.L.; Walters, W.; Berg-Lyons, D.; Huntley, J.; Fierer, N.; Owens, S.M.; Betley, J.; Fraser, L.; Bauer, M.; et al. Ultra-high-throughput microbial community analysis on the Illumina HiSeq and MiSeq platforms. ISME J. 2012, 6, 1621–1624. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  76. Walters, W.; Hyde, E.R.; Berg-Lyons, D.; Ackermann, G.; Humphrey, G.; Parada, A.; Gilbert, J.A.; Jansson, J.K.; Caporaso, J.G.; Fuhrman, J.A.; et al. Improved bacterial 16S rRNA gene (V4 and V4-5) and fungal internal transcribed spacer marker gene primers for microbial community surveys. mSystems 2015, 1, e00009-15. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  77. Caporaso, J.G.; Ackermann, G.; Apprill, A.; Bauer, M.; Berg-Lyons, D.; Betley, J.; Fierer, N.; Fraser, L.; Gilbert, J.A.; Gormley, M.; et al. Earth Microbiome Project: EMP 16S Illumina Amplicon Protocol. 2018. Available online: https://www.protocols.io/view/emp-16s-illumina-amplicon-protocol-nuudeww (accessed on 29 May 2020).
  78. Katsoula, A.; Vasileiadis, S.; Sapountzi, M.; Karpouzas, D.G. The response of soil and phyllosphere microbial communities to repeated application of the fungicide iprodione: Accelerated biodegradation or toxicity? FEMS Microbiol. Ecol. 2020, fiaa056. [Google Scholar] [CrossRef]
  79. Dodt, M.; Roehr, J.T.; Ahmed, R.; Dieterich, C. FLEXBAR—Flexible barcode and adapter processing for next-generation sequencing platforms. Biology 2012, 1, 895–905. [Google Scholar] [CrossRef] [Green Version]
  80. Bolger, A.M.; Lohse, M.; Usadel, B. Genome analysis Trimmomatic: A flexible trimmer for Illumina sequence data. Bioinformatics 2014, 30, 2114–2120. [Google Scholar] [CrossRef] [Green Version]
  81. Magoc, T.; Salzberg, S.L. FLASH: Fast length adjustment of short reads to improve genome assemblies. Bioinformatics 2011, 27, 2957–2963. [Google Scholar] [CrossRef]
  82. Edgar, R.C. UPARSE: Highly accurate OTU sequences from microbial amplicon reads. Nat. Methods 2013, 647, 1–5. [Google Scholar] [CrossRef]
  83. Edgar, R.C.; Haas, B.J.; Clemente, J.C.; Quince, C.; Knight, R. UCHIME improves sensitivity and speed of chimera detection. Bioinformatics 2011, 27, 2194–2200. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  84. Hauswedell, H.; Singer, J.; Reinert, K. Lambda: The local aligner for massive biological data. Bioinformatics 2014, 30, 349–355. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  85. Yilmaz, P.; Parfrey, L.W.; Yarza, P.; Pruesse, E.; Quast, C.; Schweer, T.; Peplies, G.; Ludwig, W.; Glockner, F.O. The SILVA and “All-species Living Tree Project (LTP)” taxonomic frameworks. Nucleic Acids Res. 2014, 42, 643–648. [Google Scholar] [CrossRef] [Green Version]
  86. Jost, L. Entropy and diversity. Opinion 2006, 113, 363–375. [Google Scholar] [CrossRef]
  87. Pielou, E.C. Ecological Diversity, 8th ed.; John Wiley & Sons: Hoboken, NJ, USA, 1975. [Google Scholar]
  88. Martinez Arbizu, P. PairwiseAdonis: Pairwise Multilevel Comparison Using Adonis. R Package Version 0.3. 2019. Available online: https://github.com/pmartinezarbizu/pairwiseAdonis (accessed on 29 May 2020).
  89. R Core Team. R: A Language and Environment for Statistical Computing, Reference Index Version 3.3.3. 2017. Available online: http://www.r-project.org/ (accessed on 29 May 2020).
Figure 1. Analysis of the α-diversity of bacteria in the gut of P. tobbi, P. similis, P. papatasi and P. neglectus based on the calculation of the indices Shannon, Inverse Simpson, Fisher’s alpha, Evenness (S), Chao 1 and Abundance-based Coverage Estimator (ACE). Error bars represent the standard error. Within each index, bars designated by the same letter are not significantly different at the 5% level.
Figure 1. Analysis of the α-diversity of bacteria in the gut of P. tobbi, P. similis, P. papatasi and P. neglectus based on the calculation of the indices Shannon, Inverse Simpson, Fisher’s alpha, Evenness (S), Chao 1 and Abundance-based Coverage Estimator (ACE). Error bars represent the standard error. Within each index, bars designated by the same letter are not significantly different at the 5% level.
Pathogens 09 00428 g001
Figure 2. Stacked bar plot showing the relative abundance of bacteria in the gut of the different Phlebotomus species (N: P. neglectus, P: P. papatasi, S: P. similis and T: P. tobbi) presented at the phylum (a), class (b), order (c) and genus (d) taxonomic level. Error bars represent the standard deviation of three biological replicates. Taxa that participated with less than 1% in 80% or more of the samples analysed were grouped as “Rare”.
Figure 2. Stacked bar plot showing the relative abundance of bacteria in the gut of the different Phlebotomus species (N: P. neglectus, P: P. papatasi, S: P. similis and T: P. tobbi) presented at the phylum (a), class (b), order (c) and genus (d) taxonomic level. Error bars represent the standard deviation of three biological replicates. Taxa that participated with less than 1% in 80% or more of the samples analysed were grouped as “Rare”.
Pathogens 09 00428 g002
Figure 3. Canonical correspondence analysis (CCA) of the bacterial community in the gut of P. papatasi, P. tobbi, P. similis and P. neglectus. Ellipses encompass all samples of the same Phlebotomus species. The model testing the host genotype effect on the bacterial community structure was significant (p < 0.001) and explained 45.1% of the total variance.
Figure 3. Canonical correspondence analysis (CCA) of the bacterial community in the gut of P. papatasi, P. tobbi, P. similis and P. neglectus. Ellipses encompass all samples of the same Phlebotomus species. The model testing the host genotype effect on the bacterial community structure was significant (p < 0.001) and explained 45.1% of the total variance.
Pathogens 09 00428 g003
Figure 4. Heatmap analysis showing operational taxonomic units (OTUs) whose relative abundance exhibits significant correlation with certain Phlebotomus species. The level of significance is denoted by asterisks as follows: * p < 0.05, ** p < 0.01, *** p < 0.001). N: P. neglectus; P: P. papatasi; T: P. tobbi; S: P. similis.
Figure 4. Heatmap analysis showing operational taxonomic units (OTUs) whose relative abundance exhibits significant correlation with certain Phlebotomus species. The level of significance is denoted by asterisks as follows: * p < 0.05, ** p < 0.01, *** p < 0.001). N: P. neglectus; P: P. papatasi; T: P. tobbi; S: P. similis.
Pathogens 09 00428 g004
Figure 5. The relative abundance (RA) of bacterial OTUs showing species-specific patterns on the different Phlebotomus species. Error bars represent standard errors. For each OTU, bars designated by the same letter are not significantly different. * p < 0.05, ** p < 0.01, *** p < 0.001.
Figure 5. The relative abundance (RA) of bacterial OTUs showing species-specific patterns on the different Phlebotomus species. Error bars represent standard errors. For each OTU, bars designated by the same letter are not significantly different. * p < 0.05, ** p < 0.01, *** p < 0.001.
Pathogens 09 00428 g005
Figure 6. (a) The relative abundance (RA) of the most dominant archaeal OTUs detected in the sandfly samples analysed and (b) The relative abundance of archaeal OTUs 503 and 1309, which showed species-specific patterns, on the different Phlebotomus species studied. Error bars represent standard errors while bars designated by the same letter are not significantly different. * p < 0.05, ** p < 0.01.
Figure 6. (a) The relative abundance (RA) of the most dominant archaeal OTUs detected in the sandfly samples analysed and (b) The relative abundance of archaeal OTUs 503 and 1309, which showed species-specific patterns, on the different Phlebotomus species studied. Error bars represent standard errors while bars designated by the same letter are not significantly different. * p < 0.05, ** p < 0.01.
Pathogens 09 00428 g006

Share and Cite

MDPI and ACS Style

Papadopoulos, C.; Karas, P.A.; Vasileiadis, S.; Ligda, P.; Saratsis, A.; Sotiraki, S.; Karpouzas, D.G. Host Species Determines the Composition of the Prokaryotic Microbiota in Phlebotomus Sandflies. Pathogens 2020, 9, 428. https://doi.org/10.3390/pathogens9060428

AMA Style

Papadopoulos C, Karas PA, Vasileiadis S, Ligda P, Saratsis A, Sotiraki S, Karpouzas DG. Host Species Determines the Composition of the Prokaryotic Microbiota in Phlebotomus Sandflies. Pathogens. 2020; 9(6):428. https://doi.org/10.3390/pathogens9060428

Chicago/Turabian Style

Papadopoulos, Christos, Panagiotis A. Karas, Sotirios Vasileiadis, Panagiota Ligda, Anastasios Saratsis, Smaragda Sotiraki, and Dimitrios G. Karpouzas. 2020. "Host Species Determines the Composition of the Prokaryotic Microbiota in Phlebotomus Sandflies" Pathogens 9, no. 6: 428. https://doi.org/10.3390/pathogens9060428

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop