Next Article in Journal
Closing the Gaps to Understand the Tick Transmission of Anaplasma marginale among Giant Anteaters (Myrmecophaga tridactyla) in Argentina
Next Article in Special Issue
Antitrypanosomal and Antileishmanial Activity of Chalcones and Flavanones from Polygonum salicifolium
Previous Article in Journal
Genomic Analysis of Aeromonas veronii C198, a Novel Mcr-3.41-Harboring Isolate from a Patient with Septicemia in Thailand
Previous Article in Special Issue
“Begging the Question”—Does Toxocara Infection/Exposure Associate with Multiple Sclerosis-Risk?
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Pathogens
Volume 9
Issue 12
10.3390/pathogens9121032
pathogens-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Integrative Approach to Phlebotomus mascittii Grassi, 1908: First Record in Vienna with New Morphological and Molecular Insights

by
Edwin Kniha
1,
Vít Dvořák
2,
Petr Halada
3,
Markus Milchram
4,
Adelheid G. Obwaller
5,
Katrin Kuhls
6,7,†,
Susanne Schlegel
6,8,†,
Martina Köhsler
1,
Wolfgang Poeppl
9,
Karin Bakran-Lebl
10,
Hans-Peter Fuehrer
10,
Věra Volfová
2,
Gerhard Mooseder
9,
Vladimir Ivovic
11,
Petr Volf
2 and
Julia Walochnik
1,*
1
Institute of Specific Prophylaxis and Tropical Medicine, Center for Pathophysiology, Infectiology and Immunology, Medical University of Vienna, 1090 Vienna, Austria
2
Department of Parasitology, Faculty of Science, Charles University Prague, 128 43 Prague, Czech Republic
3
BioCeV, Institute of Microbiology of the Czech Academy of Sciences, 252 50 Vestec, Czech Republic
4
Department of Integrative Biology and Biodiversity Research, Institute of Zoology, University of Natural Resources and Life Sciences Vienna, 1180 Vienna, Austria
5
Federal Ministry of Defence, Division of Science, Research and Development, 1090 Vienna, Austria
6
Division Molecular Biotechnology and Functional Genomics, Technical University of Applied Sciences Wildau, 15745 Wildau, Germany
7
Research Platform “Models & Simulation”, Leibniz Centre for Agricultural Landscape Research (ZALF), 15374 Müncheberg, Germany
8
Division Microsystems Engineering, Technical University of Applied Sciences Wildau, 15745 Wildau, Germany
9
Department of Dermatology and Tropical Medicine, Military Medical Cluster East, Austrian Armed Forces, 1210 Vienna, Austria
10
Department of Pathobiology, Institute of Parasitology, University of Veterinary Medicine Vienna, 1210 Vienna, Austria
11
Department of Biodiversity, FAMNIT, University of Primorska, 6000 Koper-Capodistria, Slovenia
 
add Show full affiliation list
 
remove Hide full affiliation list
*
Author to whom correspondence should be addressed.
†
Current affiliation.
Pathogens 2020, 9(12), 1032; https://doi.org/10.3390/pathogens9121032
Submission received: 17 November 2020 / Revised: 1 December 2020 / Accepted: 4 December 2020 / Published: 9 December 2020
(This article belongs to the Special Issue Translational Research for Zoonotic Parasites: New Findings toward Improved Diagnostics, Therapy and Prevention)
Browse Figures
Versions Notes

Abstract

:
Sand flies (Diptera: Psychodidae: Phlebotominae) are blood-feeding insects that transmit the protozoan parasites Leishmania spp. and various arthropod-borne (arbo) viruses. While in Mediterranean parts of Europe the sand fly fauna is diverse, in Central European countries including Austria mainly Phlebotomus mascittii is found, an assumed but unproven vector of Leishmania infantum. To update the currently understudied sand fly distribution in Austria, a sand fly survey was performed and other entomological catches were screened for sand flies. Seven new trapping locations of Ph. mascittii are reported including the first record in Vienna, representing also one of the first findings of this species in a city. Morphological identification, supported by fluorescence microscopy, was confirmed by two molecular approaches, including sequencing and matrix-assisted laser desorption/ionization-time of flight mass spectrometry (MALDI-TOF MS) protein profiling. Sand fly occurrence and activity were evaluated based on surveyed locations, habitat requirements and climatic parameters. Moreover, a first comparison of European Ph. mascittii populations was made by two marker genes, cytochrome c oxidase subunit 1 (COI), and cytochrome b (cytb), as well as MALDI-TOF mass spectra. Our study provides new important records of Ph. mascittii in Austria and valuable data for prospective entomological surveys. MALDI-TOF MS protein profiling was shown to be a reliable tool for differentiation between sand fly species. Rising temperatures and globalization demand for regular entomological surveys to monitor changes in species distribution and composition. This is also important with respect to the possible vector competence of Ph. mascittii.

Graphical Abstract

1. Introduction

Phlebotomine sand flies (Diptera: Psychodidae: Phlebotominae) occur in tropical, subtropical, as well as temperate regions. They are of significant medical importance as vectors of Leishmania spp., bacteria and several arthropod-borne (arbo) viruses in various regions of both Old and New World. Leishmaniasis is among the most prominent and yet neglected infectious diseases [1,2].
In Europe, sand flies are endemic throughout most of the Mediterranean countries while their occurrence north of the Alps and in Central Europe was overlooked for a long time. Phlebotomus mascittii Grassi, 1908, and Phlebotomus perniciosus Newstead, 1911, were recorded in Germany for the first time in 1999 and 2001, respectively [3,4]. A decade later, Ph. mascittii was recorded in Austria and its presence further confirmed by several entomological surveys [5,6,7], while a single Ph. mascittii specimen was also collected in Slovakia close to the Austrian border [8]. Very recently, Phlebotomus simici Nitzulescu, 1931, was recorded for the first time in Austria [9]. Apart from these findings, knowledge on sand fly distribution, species diversity and ecological factors determining their occurrence in Central Europe is scarce.
Ph. mascittii is the most widely distributed species in Europe and known to occur in Spain, France including Corsica, Italy, Switzerland, Germany, Belgium, Austria, Slovakia, Slovenia, Croatia, Hungary, and Serbia [3,6,8,10,11,12,13,14,15,16,17,18,19], as well as Algeria in North Africa [20]. It is also the sand fly species with the northernmost distribution in the Palearctic, occurring in Germany as far as 50° North [11,15,21].
Ph. mascittii is a suspected vector for Leishmania infantum [22,23], however, its vector competence has not yet been experimentally proven, albeit repeated reports of autochthonous leishmaniasis cases in Germany and Austria may indicate possible involvement of Ph. mascittii in transmission of Leishmania infantum [24,25,26,27].
To update the hitherto underreported sand fly distribution in Austria, dedicated entomological field studies were performed and bycatches from other entomological surveys were investigated. Species identification of obtained specimens was achieved by both, morphology and two different molecular approaches, including DNA sequencing of the two mitochondrial genes COI and cytb, and matrix-assisted laser desorption/ionization-time of flight mass spectrometry (MALDI-TOF MS) protein profiling. Ecological and climatic parameters of Ph. mascittii records were evaluated to promote prospective trapping success.

2. Results

2.1. Sand Fly Trapping and Identification

Altogether, 65 sampling sites in six federal states of Austria, namely Vienna, Lower Austria, Burgenland, Styria, Upper Austria, and Vorarlberg were surveyed (Table S1). At seven (10.8%) of these locations altogether 28 sand flies were trapped (Figure 1). Two (7.1%) were male and 26 (92.9%) were female, of which one specimen was engorged and three were gravid. The earliest capture was observed on 5 July and the latest capture was observed on 24 August.
In survey 1, two of the 46 (4.3%) sampled locations were positive. A single specimen was trapped in the 11th district of Vienna, the capital of Austria, on 17 July 2019 at a horse farm outside a barn used for hay storage. A second location was found positive in Laafeld, Styria, where one and two specimens were trapped in front of a chicken shed on 6 August 2018 and 5 August 2019, respectively. (Table 1). All locations surveyed in Upper Austria and Vorarlberg were negative (Figure 1).
In survey 2, only one location, namely Kaisersteinbruch in Lower Austria was sampled, where two specimens were caught on 11 and 26 July 2013 outside of a dog kennel (Table 1).
In survey 3, three of nine (33.3%) sampled locations were positive, which were all located in the federal district of Styria. In Hummersdorf, six and eight specimens were caught on 5 and 6 July 2015, respectively, inside and outside a chicken barn. In Bad Radkersburg, two specimens were caught on 6 July 2015 inside an old chicken shed and in Unterpurkla, four specimens were caught on 6 July 2015 inside a barn next to a chicken shed (Table 1).
In survey 4, screening bycatches of a mosquito monitoring, one of nine (11.1%) sampled locations were positive. Two sand fly specimens were trapped on 29 July and 24 August 2019 with a BG sentinel trap baited with CO2 in Neuhaus in the federal state of Burgenland (Table 1). The trap was set in the garden of a private property without animal barns or sheds. No animals were reported to live on the property and only rodents were observed to be present under a pile of wood close to the trap by the owners.
The majority of trapping sites, including those without sand fly catches, were animal farms with various domestic and farm animals as well as rodents present, being typical and suitable sand fly trapping sites.
Of all 28 caught specimens, 22 (78.6%) were caught with standard CDC light traps, four (14.3%) with standard CDC light traps with additional dry ice and two (7.1%) were caught as bycatch with BG sentinel traps using CO2 as bait. However, due to the low number of caught specimens and the fact that not all sampling techniques were applied at all locations, statistical analysis on trapping success by type of trap was not possible.
All specimens were morphologically identified as Phlebotomus mascittii. Further confirmation was achieved by fluorescence microscopy, all three applied light spectra providing detailed images of the spermathecae (Figure 2). Correct morphological identification as Ph. mascittii was confirmed by DNA sequencing of all specimens. Obtained sequences were compared to reference sequences stored in GenBank. Sequence similarities with sequences in GenBank ranged from 99.84% to 100% (KY848831.1, KX981913.1) for cytochrome c oxidase subunit 1 (COI) and from 98.44% to 100% (MG800324.1, KR336656.1) for cytochrome b (cytb).
MALDI-TOF MS protein profiling provided reproducible mass spectra with a high number of intense signals for all analyzed specimens (Figure 3), males and females producing similar protein profiles. All 5 spectra of specimens from Austria confirmed the species identification of Ph. mascittii. When these protein profiles were compared with those of Ph. mascittii originating from three European countries (France, Slovenia, and Serbia), the dendrogram generated by cluster analysis provided clear clustering according to the geographical origin with distinct clusters for all four countries. The protein profile of Ph. killicki, a closely related species of the subgenus Transphlebotomus, formed a sister cluster to all analyzed Ph. mascittii specimens (Figure 4).
Leishmania spp. DNA was not detected by PCR in any of the female specimens.

2.2. Pairwise Distances of Transphlebotomus Species

Overall, 27 COI sequences of all five Transphlebotomus species with a final length of 566 bp were included in the comparative analyses (Table S3). Mean intraspecific distances ranged from 0.09% to 0.6%, the lowest being calculated for Ph. mascittii (Table 2). In general, mean interspecific distances ranged from 6.6% to 15.4% and from 10.7% to 15.4% between Ph. mascittii and other Transphlebotomus species (Table 2).
Altogether, 31 cytb sequences of all five Transphlebotomus species with a final length of 475 bp were included in the comparative analysis (Table S4). Mean intraspecific distances ranged from 0.2% to 1.9% (Table 2). Overall, mean interspecific distances ranged from 8.7% to 13.9% and from 12.9% to 14.7% between Ph. mascittii and other Transphlebotomus species (Table 2).

2.3. Haplotype Analysis of Ph. Mascittii from Austria and Other Countries

Overall, 12 COI sequences with a final length of 613 bp were included in the analysis (Table S2) with an overall haplotype diversity (Hd) of 0.53 and an overall nucleotide diversity (π) of 0.53 (Table 3). The sequences of Ph. mascittii, originating from four different countries, namely Austria, Slovakia, Slovenia and Serbia, revealed two haplotypes (COI_1, COI_2) defined by one polymorphic site (Pos: 106, A/G) (Figure 5a). Interestingly, both haplotypes were found among the Austrian specimens. The specimens from the Austrian federal states Vienna and Styria shared one haplotype (COI_1) with specimens from Pernek, Slovakia, and Vojdovina, Serbia. The other Austrian specimens from Lower Austria, Burgenland, and Styria shared the other haplotype (COI_2) with specimens from Velike Zabjle and Cetore, Slovenia, as well as Krasava, Serbia (Table 3).
Altogether, 16 cytb sequences with a final length of 475 bp were used for the analysis (Table S3) with an overall haplotype diversity (Hd) of 0.45 and an overall nucleotide diversity (π) of 0.19 (Table 3). Sequences of Ph. mascittii originating from six different countries, namely France including Corsica, Belgium, Germany, Austria, Slovenia, and Serbia, revealed five haplotypes (Cytb_1–Cytb_5) defined by seven variable sites, of which three were parsimony informative (Figure 5b). In the cytb analysis, all Austrian specimens belonged to a single haplotype (Cytb_1), shared with specimens from Belgium, France including Corsica, Germany, and Serbia (Table 3).

2.4. Climatic Parameters

No significant differences were observed between annual mean temperature (Kruskal–Wallis test, p = 0.941), mean June temperature (Kruskal–Wallis test, p = 0.548), mean July temperature (Kruskal–Wallis test, p = 0.889) and mean August temperature (Kruskal–Wallis test, p = 0.883) at successful and unsuccessful trapping sites (Figure 6a–d).
Further, no significant differences were observed between mean day temperatures (Kruskal–Wallis test, p = 0.344) and mean night temperatures (Kruskal–Wallis test, p = 0.435) at successful and unsuccessful trapping sites (Figure 6e,f). Sand flies were found to be active between 14.1 and 22.9 °C (mean = 19.8, SD = 2.9) mean night temperature (Figure 6f). Minimum night temperatures during successful trap nights ranged from 13.5 to 21.4 °C.
However, significant differences based on successful and unsuccessful trap nights were observed between mean day relative humidity (RH) (Kruskal–Wallis test, p < 0.001) and mean night RH (Kruskal–Wallis test, p < 0.001) (Figure 6g,h). Sand flies were found to be active between 78.6% and 95.7% (mean = 88.0, SD = 7.8) mean night RH (Figure 6h).

3. Discussion

This study reports new findings of Phlebotomus mascittii in Austria, including the first record of this species in the capital Vienna, representing the northernmost record in Austria. In total, Ph. mascittii was found at seven new locations within four surveys, suggesting that its distribution in Austria may be wider than previously thought. We provide an integrative approach for its identification, using morphological as well as molecular methods.
Ph. mascittii is the type species of the subgenus Transphlebotomus Artemiev, 1984, which comprises four other species, namely Phlebotomus canaaniticus Adler and Theodor, 1931, Phlebotomus economidesi Léger, Depaquit, and Ferté, 2000, Phlebotomus killicki Dvořák, Votýpka and Volf, 2015, and Phlebotomus anatolicus Kasap, Depaquit, and Alten, 2015. A recent description of the latter two species from localities in the eastern Mediterranean region where composition of sand fly fauna had been repeatedly studied in the past underlines the challenges of species identification based on morphology within this subgenus [28]. Therefore, molecular approaches are useful to confirm species identification and the natrium dehydrogenase subunit 4 (NADH4) gene was presented as an informative marker for interspecific and intraspecific analysis within the Transphlebotomus subgenus [11].
Here, we applied three different techniques, including morphology, sequencing of the two commonly used mitochondrial marker genes COI and cytb and MALDI-TOF MS protein profiling. All specimens were morphologically identified to the species level, however, some individuals had hardly visible spermathecae, making morphological discrimination challenging. We further evaluated the usage of fluorescence microscopy to detect the autofluorescent spermathecae of Ph. mascittii, which we had already used for the identification of Phlebotomus simici [9]. By using different wavelengths of light the typical microstructures of the spermatheca of Ph. mascittii were visualized in high detail. Verification of morphological results was achieved with both mitochondrial marker genes used, COI and cytb. In general, interspecific distances exceeded 10% in both marker genes, with one exception, interspecific distances between Ph. canaaniticus and Ph. anatolicus were 6.6 and 8.7% for COI and cytb, respectively. This is in concordance with Kasap et al. [28], who observed similar interspecific distances with both marker genes and identified these two species as sister species.
MALDI-TOF MS protein profiling also confirmed species identification of Austrian specimens, providing distinct, species-specific protein profiles. This mass spectrometry method represents a currently emerging tool for species identification of various organisms including medically important arthropods that becomes favorable as a cost- and labor-effective alternative approach especially in field studies. However, in sand flies, it has been mostly used in studies focusing on the species diversity in rather small areas, often in foci of human leishmaniasis [29,30]. We attempted for the first time to compare protein profiles of distant populations of one species on a much larger geographical scale, analyzing specimens of Ph. mascittii from Austria, the French island of Corsica, Slovenia and Serbia. Our results clearly show that this method differentiates between these geographically distinct populations while successfully identifying all as Ph. mascittii when compared to the closely related species Ph. killicki. Specimens from a single locality in Corsica form a very compact cluster while specimens from Austria and Slovenia that originate from several localities show minor differences and form longer branches within their clusters. Specimens from Serbia were all collected at one locality, yet they also form longer branches, presumably due to the lower quality of obtained spectra because of compromised storage conditions.
Altogether, only seven out of the 64 surveyed locations were found positive for Ph. mascittii, even though most of the sampled locations were considered ideal for sand fly occurrence regarding annual mean temperature and potential breeding sites. Typically, only small trapping numbers were observed, as has been commonly reported from previous entomological surveys on sand flies in Austria [5,6] and bordering countries [8,31]. The attractiveness of light traps to Ph. mascittii is still a matter of discussion. However, eight individuals were trapped by a standard CDC light trap during a single night in Hummersdorf, indicating that light trapping is an appropriate technique for Ph. mascittii. In addition to this trapping method, trappings were also successful using CDC light traps with additional dry ice and BG sentinel traps baited with CO2 as an attractant. Successful trappings have also been reported with other methods, e.g., sticky traps [12,18,32], malaise traps [23] and emergence traps in a cellar in Switzerland [15], all nevertheless resulting in rather low trapping numbers. High numbers were reported from trappings with light traps in Sessa, southern Switzerland, where over 280 Ph. mascittii specimens were caught during a survey from 1987 to 1989 [15]. Comparative field studies are necessary to optimize trapping methods and attractants for Ph. mascittii. Based on our findings of Ph. mascittii in the frame of mosquito monitoring and the variety of suitable trapping methods for this species, the inspection of bycatch of entomological studies presents an opportunity to detect small overlooked Ph. mascittii populations.
The generally low population densities found might not only depend on the trapping method, but might be due to Ph. mascittii having specific ecological niches within typical habitats or might even depend on other yet unknown factors. As also observed in this study, Ph. mascittii in Central Europe is frequently found in old barns and sheds close to human dwellings and animals [3,5,6,8], even though Ph. mascittii has been proposed to be cavernicolous [11]. Ph. mascittii has been found to be active during winter time in a cave-like tunnel on the island of Corsica, where constant climatic conditions might increase population density and the time of activity [13]. A small number of individuals of this species have recently also been found in a small cave in Germany [31]. Despite these findings, records of such trapping sites are rather scarce. On the other hand, Ph. mascittii has been found to be present at other trapping sites such as wall crevices in Italy [12,32], dry stone walls at illegal waste sites in Slovenia [33] and in the basement of a house in Switzerland [15].
Despite the fact that global warming might promote the expansion of sand fly populations in Central Europe, the constant renovation and demolition of typical trapping sites such as old barns and sheds in Austria might be a potential limitation to the further dispersal, if other niches are not inhabited. This was observed while trapping in Upper Austria and Vorarlberg, where locating traditional old barns remained challenging, possibly one of the main reasons for not detecting sand flies in these areas so far, rather than climatic conditions.
In Vienna Ph. mascittii was found on a horse farm, of which there are many on the outskirts of Vienna but only very few towards the center. Even though climatic conditions are optimal in Vienna, typical breeding sites are scarce. Thus, alternative breeding sites should be evaluated and sampled within Vienna, to assess the potential spread of Ph. mascittii and the potential presence of other sand fly species within the city. The presence of sand flies in Vienna, with currently 1.9 million inhabitants, might be of public health relevance [22].
Temperature is one of the driving factors for sand fly occurrence and the 10 °C-annual-isotherm was proposed to be the boundary for northward dispersal of European sand fly species [34]. However, sand fly larvae overwinter in the ground and therefore we expect that temperatures during the summer season are more relevant for sand fly distribution. In this study, mean June, July, and August temperatures were slightly higher at successful trapping sites, but no significant differences compared to unsuccessful trapping sites were observed. Of all climatic parameters included in the analyses, only the relative humidity was significantly higher at successful sampling sites than at unsuccessful ones. Nonetheless, not all negative trapping sites can be explained by climatic conditions. Phlebotomus mascitti was observed to be active at night temperatures as low as 14.2 °C mean, which probably marks the lower boundary of sand fly activity and might affect population densities. The observed large activity spectrum between 14.2 °C and 22.9 °C is not surprising, as Ph. mascittii occurs in a variety of geographic regions including Central Europe, Southern Europe, and Northern Africa, where climatic conditions differ considerably. In addition, other factors such as relative humidity, precipitation and air pressure have an impact on sand fly activity during their active season and should be considered for field work [35,36].
Even though we did not detect Leishmania DNA in any of the sand fly females, such a monitoring is important. An anthropophilic behavior of Ph. mascittii was shown in Southern Switzerland [15] and the detection of Leishmania infantum DNA in an unfed female specimen in Austria [22] and on the Italian island of Montecristo [23] urge for further clarification of its suspected vector competence. Moreover, the frequent import of often asymptomatic Leishmania-infected dogs [37] adds another significant risk factor for the further spread of leishmaniasis.
The failure to trap Ph. mascittii in the federal states Vorarlberg and Upper Austria, of course, does not exclude that populations were overlooked. Further trappings over longer periods of time are necessary for confirmation. The absence of Ph. mascittii and other sand fly species in certain parts of Central Europe is most certainly a result of biogeographic events during and after the last glacial period [38]. During the last glacial period, Central Europe was most certainly free of sand flies and their return from Mediterranean refugial areas possibly started around 10,000 to 8000 years ago [39]. A maximum sand fly distribution in Central Europe was possibly reached during the Holocene optima around 6500 and 4500 years ago, respectively, and sand flies remained restricted to small favorable microclimatic habitats until today [34]. This hypothesis is supported by our haplotype networks which demonstrate close genetic relationships between European Ph. mascittii populations, indicating a rather recent dispersal. Although both marker genes, COI and cytb, are frequently used in phylogenetic studies [40], we could not detect substantial differences at the population level to further elucidate the origin of Ph. mascittii in Austria. Due to such very recent dispersal of Ph. mascittii within Europe, COI and cytb might not be ideal as marker genes and genes with higher mutation rates, such as the internal transcribed spacer 2, might be more suitable [41].
Overall, the Danube valley might constitute a corridor for sand fly dispersal, thereby connecting eastern and western populations. Whether eastern and western parts of Europe were post-glacially colonized from the same refugial area should be addressed in a detailed phylogeographic study, which would be crucial to understand the current sand fly distribution in Central Europe and to predict further sand fly dispersal in the future.

4. Material and Methods

4.1. Sand Fly Trapping and Available Material

The present study combines data of four different trapping surveys to update the current knowledge of sand fly distribution in Austria. Firstly, extensive field surveys to monitor sand flies were performed from 8 July to 21 August 2018 and 1 July to 25 August 2019, in six different federal states of Austria (survey 1). Trapping locations were chosen based on prevailing mean annual temperatures above 8 °C. In a second step, within these areas, suitable trapping sites were selected including farms, old barns, animal stables, private properties with and without animals present as well as two zoos. Every location was sampled for 2 to 3 consecutive nights based on weather conditions. Battery-operated CDC miniature light traps with fine gossamer collection bags (model #512, John W. Hock Company, Gainesville, FL, USA) were used. At outside locations, an additional source of CO2 was provided by a cooling bottle filled with ~700 g of dry ice/trap.
Secondly, material from prior sand fly trapping surveys in the years 2013 (survey 2) and 2015 (survey 3) with CDC light traps was included in this study.
Thirdly, bycatch of mosquito monitoring in the federal district Burgenland in July and August 2019 performed with BG sentinel traps baited with CO2 (Biogents AG, Regensburg, Germany) was screened for sand flies and included in this study (survey 4). Ethics Approval and Consent to Participate: Not applicable.

4.2. Climate Data

Mean annual, mean June, mean July and mean August temperatures of the trapping sites were accessed online (https://de.climate-data.org). Hourly temperature and relative humidity data of sampled locations were retrospectively obtained from the Central Institute for Meteorology and Geodynamics (ZAMG). Mean and minimum temperature and relative humidity were calculated for days (6 a.m. to 21 p.m.) and nights (22 p.m. to 5 a.m.).

4.3. Mapping of Sand Fly Distribution

Positive trapping sites as well as all previously published trapping locations in Austria were incorporated into a distribution map using Quantum GIS 3.4.11 [42].

4.4. Morphological Identification

Head and terminal segments of the abdomen of all caught sand fly specimens were dissected and slide-mounted in CMCP-10 mountant (Polysciences, Inc., Warrington, PA, USA). Identification was based on published morphological keys and descriptions of male genitalia, female spermatheca and pharyngeal armature [11,43]. Additionally, fluorescence microscopy was performed with a NIKON Eclipse E800 to identify hardly visible female spermatheca, as proposed in [9].

4.5. Molecular Identification by PCR and Sequencing

Molecular identification was performed by sequencing of the mitochondrial cytochrome c oxidase subunit 1 (COI) and cytochrome b (cytb) gene regions.
DNA was isolated from the remaining bodies with a QIAamp® DNA Mini Kit 250 (QIAGEN, Hilden, Germany). PCR amplifications of the COI and cytb gene regions were performed with the primer pairs LCO-1490/CoxUniEr and CytbEf1/CytbEr2, respectively, as published in [9].
All PCR amplifications were performed with an Eppendorf Mastercycler (Eppendorf AG, Hamburg, Germany). The PCR products were subjected to electrophoresis in 2% agarose gels stained with GelRed® Nucleic Acid Gel Stain (Biotium, Inc., Hayward, CA, USA). For further sequencing bands were analyzed with a Gel DocTM XR+ Imager (Bio-Rad Laboratories, Inc., Hercules, CA, USA), cut out from the gel and purified with the IllustraTM GFXTM PCR DNA and Gel Purification Kit (GE Healthcare, Buckinghamshire, UK). Sanger sequencing was performed with a Thermo Fisher Scientific SeqStudio (Thermo Fisher Scientific, Waltham, MA, USA). Sequences were obtained from both DNA strands and a consensus sequence was generated in GenDoc 2.7.0. Obtained sequences were compared to available sequences in the GenBank using the Basic Local Alignment Search Tool (BLAST) (https://blast.ncbi.nlm.nih.gov/Blast.cgi). One sequence of each location was submitted to GenBank: COI (MN812827.1–MN812830.1, MT332686.1–MT332688.1) and cytb (MN812832.1–MN812835.1, MT332689.1–MT332691.1).
To confirm the suitability of COI and cytb as marker genes to distinguish between Transphlebotomus species, mean interspecific and intraspecific Kimura-2-parameter (K2P) distances were calculated with MEGAX [44].

4.6. Identification by MALDI-TOF MS Protein Profiling

To confirm species identification of chosen specimens at surveyed localities (survey 3), MALDI-TOF MS protein profiling was applied according to previously optimized protocols for trapping and sample preparation [45,46]. Thoraxes of specimens trapped by CDC light traps were manually grinded by disposable pestles in 1.5 mL microtubes with 10 μL of 25% formic acid and briefly centrifuged at 10,000× g. Two μL of the homogenate were mixed with 2 μL of freshly prepared MALDI matrix, an aqueous 60% acetonitrile/0.3% TFA solution of sinapinic acid (30 mg/mL; Bruker Daltonics, Billerica, MA, USA). One μL of this mixture was applied on a steel MALDI plate in duplicates, air-dried and measured by Autoflex Speed MALDI-TOF spectrometer (Bruker Daltonics, Billerica, MA, USA) in a mass range of 4–25 kDa (20 × 300 laser shots from different positions of the sample spot). Obtained spectra were visualized by FlexAnalysis 3.4 software (Bruker Daltonics, Billerica, MA, USA), processed by MALDI Biotyper 3.1 and an in-house database that comprises reference spectra of 25 different sand fly species was searched to identify the species. Log score value (LSV) > 2.0 was decided as a threshold for an unambiguous identification. For MSP dendrogram creation, an individual main spectrum (MSP) was generated from each analyzed spectrum.
In total, 5 specimens from two Austrian localities (Hummersdorf and Unterpurkla) were subjected to MALDI-TOF MS protein profiling. For subsequent comparison of geographically distant European populations, we also included protein profiles of 4 specimens from another Austrian location (Rohrau) collected in 2018 [7] and three other European countries: 5 specimens (all males) collected in 2019 from a railway tunnel close to the road from Solenzara to Ste. Lucie de Porto Vecchio, Corsica, France [13], 3 specimens (1 male, 2 females) collected in 2016 from Krasava, Serbia [17] and 5 specimens (all females) collected in 2020 from Godovic, Ilirska Bistrica, Vrhnika, Zovnek and Velike Zablje, Slovenia. As an outgroup, a protein profile of Phlebotomus killicki, a closely related species of the subgenus Transphlebotomus, collected in 2019 in Xerokampos, Crete, Greece, was used [47] (Table S4).

4.7. Leishmania spp. Screening

Female specimens were screened by PCR for the presence of Leishmania parasites using the primers LITSR/L5.8S targeting the internal transcribed spacer 1 (ITS1) of the ribosomal DNA following the PCR protocol of protocol of El Tai et al. [48]. Five μL of extracted DNA from female sand flies was used in all PCR reactions. Microbial DNA free water (QIAGEN, Hilden, Germany) and 5 μL DNA of a male Ph. mascittii specimen were used as negative controls, 2 μL DNA of Leishmania infantum MHOM/TR/2000/OG-VL was used as a positive control.

4.8. Sequence Analyses

Obtained sequences were aligned with ClustalX 2.1 and edited with GeneDoc 2.7.0. for further analysis. DnaSP v.5 [49] was used to identify unique haplotypes. Median Joining Networks [50] were calculated and visualized with Popart v.1.7 [51].

4.9. Statistical Analyses

Data was analyzed with the R environment for Mac [52]. Descriptive statistics were performed with a one sample proportions test. A Shapiro–Wilk test was used to check for normality of the data. A nonparametric Kruskal–Wallis test was used to compare climate data. A two-sided p-value < 0.05 was considered statistically significant.

Supplementary Materials

The following are available online at https://www.mdpi.com/2076-0817/9/12/1032/s1, Table S1. Surveyed locations in Austria. Table S2. Included COI sequences of Transphlebotomus for inter- and intraspecific distance calculations. Table S3. Included cyt b sequences of Transphlebotomus for inter- and intraspecific distance calculations. Table S4. Included sand fly specimens to MALDI-TOF mass spectrometry.

Author Contributions

E.K., A.G.O., W.P., G.M., and J.W. designed the study. E.K., A.G.O., K.K., S.S., K.B.-L., H.-P.F., V.V., V.I., and P.V. conducted field work. E.K., V.D., and P.H. performed laboratory work. E.K., V.D., P.H., M.M., and M.K. analyzed the data. E.K., V.D., P.V., and J.W. wrote the manuscript. All authors have read and agreed to the published version of the manuscript.

Funding

The study was supported by the Austrian Federal Ministry of Defence and the Medical University of Vienna, Austria. The additional support from the European Regional Development Fund (projects BIOCEV CZ.1.05/1.1.00/02.0109 and CePaViP CZ.02.1.01/0.0/0.0/16_019/0000759) and from the Czech Infrastructure for Integrative Structural Biology (LM2015043 funded by MEYS CR) is gratefully acknowledged. Cooperative field work was funded by an Austrian–Czech Republic joint project (grant nr. CZ02/2020) granted by the Austrian Agency for International Cooperation in Education and Research (OEAD) and by the German Federal Ministry of Education and Research (BMBF) (grant nr. FKZ01DK14021). The funders had no role in the study design, data collection and analysis, decision to publish, or preparation of the manuscript. Edwin Kniha is a recipient of a DOC fellowship and funded by the Austrian Academy of Sciences (ÖAW). His stay at the Department of Parasitology, Charles University, was funded by the European Union’s Horizon 2020 research and innovation programme Infravec2, research infrastructures for the control of insect vector-borne diseases under grant agreement No 731060.

Acknowledgments

We would like to thank the Tiergarten Schönbrunn (Vienna), the Blumengarten Hirschstetten (Vienna), the Nationalpark Donauauen (Lower Austria/Vienna), and the Nationalpark Neusiedler See—Seewinkel (Burgenland) for their collaboration and their great support during fieldwork. We thank the owners of the properties were trappings have been performed. Special thanks to the ZAMG (Central Institute for Meteorology and Geodynamics) for collaboration and providing climate data. The authors are also grateful to Iveta Haefeli for technical support.

Conflicts of Interest

The authors declare that they have no competing interest.

Abbreviations

COI: cytochrome c oxidase subunit I, cytb: cytochrome b, Hd: haplotype diversity, MALDI-TOF MS: matrix-assisted laser desorption/ionization time of flight mass spectrometry, RH: relative humidity, SD: standard deviation.

References

  1. Akhoundi, M.; Kuhls, K.; Cannet, A.; Votýpka, 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, e0004349. [Google Scholar] [CrossRef] [PubMed]
  2. Ready, P.D. Biology of Phlebotomine Sand Flies as Vectors of Disease Agents. Annu. Rev. Entomol. 2013, 58, 227–250. [Google Scholar] [CrossRef] [PubMed]
  3. Naucke, T.J.; Pesson, B. Presence of Phlebotomus (Transphlebotomus) mascittii Grassi, 1908 (Diptera: Psychodidae) in Germany. Parasitol. Res. 2000, 86, 335–336. [Google Scholar] [CrossRef] [PubMed]
  4. Naucke, T.J.; Schmitt, C. Is leishmaniasis becoming endemic in Germany? Int. J. Med. Microbiol. 2004, 293, 179–181. [Google Scholar] [CrossRef]
  5. Naucke, T.J.; Lorentz, S.; Rauchenwald, F.; Aspöck, H. Phlebotomus (Transphlebotomus) mascittii Grassi, 1908, in Carinthia: First record of the occurrence of sandflies in Austria (Diptera: Psychodidae: Phlebotominae). Parasitol. Res. 2011, 109, 1161–1164. [Google Scholar] [CrossRef] [PubMed]
  6. Poeppl, W.; Obwaller, A.G.; Weiler, M.; Burgmann, H.; Mooseder, G.; Lorentz, S.; Rauchenwald, F.; Aspöck, H.; Walochnik, J.; Naucke, T.J. Emergence of sandflies (Phlebotominae) in Austria, a Central European country. Parasitol. Res. 2013, 112, 4231–4237. [Google Scholar] [CrossRef] [Green Version]
  7. Obwaller, A.G.; Poeppl, W.; Naucke, T.J.; Luksch, U.; Mooseder, G.; Aspöck, H.; Walochnik, J. Stable populations of sandflies ( Phlebotominae ) in Eastern Austria: A comparison of the trapping seasons 2012 and 2013. Trends Entomol. 2014, 2, 1–5. [Google Scholar]
  8. Dvořák, V.; Hlavackova, K.; Kocisova, A.; Volf, P. First record of Phlebotomus (Transphlebotomus) mascittii in Slovakia. Parasite 2016, 23, 48. [Google Scholar] [CrossRef] [Green Version]
  9. Kniha, E.; Dvořák, V.; Milchram, M.; Obwaller, A.G.; Koehsler, M.; Poeppl, W.; Antoniou, M.; Chaskopoulou, A.; Paronyan, L.; Stefanovska, J.; et al. Phlebotomus (Adlerius) simici Nitzulescu, 1931: First record in Austria and phylogenetic relationship with other Adlerius species. Parasit. Vectors 2020, (in press). [Google Scholar]
  10. Aransay, A.M.; Testa, J.M.; Morillas-Márquez, F.; Lucientes, J.; Ready, P.D. Distribution of sandfly species in relation to canine leishmaniasis from the Ebro Valley to Valencia, northeastern Spain. Parasitol. Res. 2004, 94, 416–420. [Google Scholar] [CrossRef]
  11. Depaquit, J.; Naucke, T.J.; Schmitt, C.; Ferté, H.; Léger, N. A molecular analysis of the subgenus Transphlebotomus Artemiev, 1984 (Phlebotomus, Diptera, Psychodidae) inferred from ND4 mtDNA with new northern records of Phlebotomus mascittii Grassi, 1908. Parasitol. Res. 2005, 95, 113–116. [Google Scholar] [CrossRef]
  12. Veronesi, E.; Pilani, R.; Carrieri, M.; Bellini, R. Trapping sand flies (Diptera: Psychodidae) in the Emilia-Romagna region of northern Italy. J. Vector Ecol. 2007, 32, 313–318. [Google Scholar] [CrossRef]
  13. Naucke, T.J.; Menn, B.; Massberg, D.; Lorentz, S. Winter activity of Phlebotomus (Transphlebotomus) mascittii, Grassi 1908 (Diptera: Psychodidae) on the island of Corsica. Parasitol. Res. 2008, 103, 477–479. [Google Scholar] [CrossRef]
  14. Bosnić, S.; Gradoni, L.; Khoury, C.; Maroli, M. A review of leishmaniasis in Dalmatia (Croatia) and results from recent surveys on phlebotomine sandflies in three southern counties. Acta Trop. 2006, 99, 42–49. [Google Scholar] [CrossRef] [PubMed]
  15. Grimm, F.; Gessler, M.; Jenni, L. Aspects of sandfly biology in southern Switzerland. Med. Vet. Entomol. 1993, 7, 170–176. [Google Scholar] [CrossRef] [PubMed]
  16. Praprotnik, E.; Zupan, S.; Ivović, V. Morphological and Molecular Identification of Phlebotomus mascittii Grassi, 1908 Populations From Slovenia. J. Med. Entomol. 2019, 56, 565–568. [Google Scholar] [CrossRef] [PubMed]
  17. Vaselek, S.; Dvořák, V.; Hlavackova, K.; Ayhan, N.; Halada, P.; Oguz, G.; Ivović, V.; Ozbel, Y.; Charrel, R.N.; Alten, B.; et al. A survey of sand flies (Diptera, Phlebotominae) along recurrent transit routes in Serbia. Acta Trop. 2019, 197, 105063. [Google Scholar] [CrossRef] [PubMed]
  18. Farkas, R.; Tánczos, B.; Bongiorno, G.; Maroli, M.; Dereure, J.; Ready, P.D. First Surveys to Investigate the Presence of Canine Leishmaniasis and Its Phlebotomine Vectors in Hungary. Vector-Borne Zoonotic Dis. 2011, 11, 823–834. [Google Scholar] [CrossRef]
  19. Dantas-Torres, F.; Tarallo, V.D.; Latrofa, M.S.; Falchi, A.; Lia, R.P.; Otranto, D. Ecology of phlebotomine sand flies and Leishmania infantum infection in a rural area of southern Italy. Acta Trop. 2014, 137, 67–73. [Google Scholar] [CrossRef]
  20. Berdjane-Brouk, Z.; Charrel, R.N.; Bitam, I.; Hamrioui, B.; Izri, A. Record of Phlebotomus (Transphlebotomus) mascittii Grassi, 1908 and Phlebotomus (Laroussius) chadli Rioux, Juminer & Gibily, 1966 female in Algeria. Parasite 2011, 18, 337–339. [Google Scholar]
  21. Melaun, C.; Krüger, A.; Werblow, A.; Klimpel, S. New record of the suspected leishmaniasis vector Phlebotomus (Transphlebotomus) mascittii Grassi, 1908 (Diptera: Psychodidae: Phlebotominae)—The northernmost phlebotomine sandfly occurrence in the Palearctic region. Parasitol. Res. 2014, 113, 2295–2301. [Google Scholar] [CrossRef]
  22. Obwaller, A.G.; Karakus, M.; Poeppl, W.; Töz, S.; Özbel, Y.; Aspöck, H.; Walochnik, J. Could Phlebotomus mascittii play a role as a natural vector for Leishmania infantum? New data. Parasit. Vectors 2016, 9, 458. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  23. Zanet, S.; Sposimo, P.; Trisciuoglio, A.; Giannini, F.; Strumia, F.; Ferroglio, E. Epidemiology of Leishmania infantum, Toxoplasma gondii, and Neospora caninum in Rattus rattus in absence of domestic reservoir and definitive hosts. Vet. Parasitol. 2014, 199, 247–249. [Google Scholar] [CrossRef] [PubMed]
  24. Bogdan, C.; Schönian, G.; Bañuls, A.L.; Hide, M.; Pratlong, F.; Lorenz, E.; Röllinghoff, M.; Mertens, R. Visceral leishmaniasis in a German child who had never entered a known endemic area: Case report and review of the literature. Clin. Infect. Dis. 2001, 32, 302–306. [Google Scholar] [CrossRef] [PubMed]
  25. Koehler, K.; Stechele, M.; Hetzel, U.; Domingo, M.; Schönian, G.; Zahner, H.; Burkhardt, E. Cutaneous leishmaniosis in a horse in southern Germany caused by Leishmania infantum. Vet. Parasitol. 2002, 109, 9–17. [Google Scholar] [CrossRef]
  26. Kollaritsch, H.; Emminger, W.; Zaunschirm, A.; Aspöck, H. Suspected Autochthonous Kala-azar in Austria. Lancet 1989, 1, 901–902. [Google Scholar] [CrossRef]
  27. Beyreder, J. Ein Fall von Leishmaniose in Niederösterreich. Wien. Med. Wochenschr. 1962, 115, 900–901. [Google Scholar]
  28. Kasap, O.E.; Dvořák, V.; Depaquit, J.; Alten, B.; Votypka, J.; Volf, P. Phylogeography of the subgenus Transphlebotomus Artemiev with description of two new species, Phlebotomus anatolicus n. sp. and Phlebotomus killicki n. sp. Infect. Genet. Evol. 2015, 34, 467–479. [Google Scholar] [CrossRef]
  29. Lafri, I.; Almeras, L.; Bitam, I.; Caputo, A.; Yssouf, A.; Forestier, C.L.; Izri, A.; Raoult, D.; Parola, P. Identification of Algerian Field-Caught Phlebotomine Sand Fly Vectors by MALDI-TOF MS. PLoS Negl. Trop. Dis. 2016, 10, e0004351. [Google Scholar] [CrossRef]
  30. Pareyn, M.; Dvořák, V.; Halada, P.; Van Houtte, N.; Medhin, G.; de Kesel, W.; Merdekios, B.; Massebo, F.; Leirs, H.; Volf, P. An integrative approach to identify sand fly vectors of leishmaniases in Ethiopia by morphological and molecular techniques. Parasit. Vectors 2020, 13, 580. [Google Scholar] [CrossRef]
  31. Oerther, S.; Jöst, H.; Heitmann, A.; Lühken, R.; Krüger, A.; Steinhausen, I.; Brinker, C.; Lorentz, S.; Marx, M.; Schmidt-Chanasit, J.; et al. Phlebotomine sand flies in Southwest Germany: An update with records in new locations. Parasit. Vectors 2020, 13, 173. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  32. Rossi, E.; Bongiorno, G.; Ciolli, E.; Di Muccio, T.; Scalone, A.; Gramiccia, M.; Gradoni, L.; Maroli, M. Seasonal phenology, host-blood feeding preferences and natural Leishmania infection of Phlebotomus perniciosus (Diptera, Psychodidae) in a high-endemic focus of canine leishmaniasis in Rome province, Italy. Acta Trop. 2008, 105, 158–165. [Google Scholar] [CrossRef]
  33. Ivović, V.; Kalan, K.; Zupan, S.; Bužan, E. Illegal waste sites as a potential micro foci of Mediterranean Leishmaniasis: First records of phlebotomine sand flies (Diptera: Psychodidae) from Slovenia. Acta Vet. Brno 2015, 65, 348–357. [Google Scholar]
  34. Trájer, A.J.; Sebestyén, V. The changing distribution of Leishmania infantum Nicolle, 1908 and its Mediterranean sandfly vectors in the last 140 kys. Sci. Rep. 2019, 9, 11820. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  35. Simsek, F.M.; Alten, B.; Caglar, S.S.; Ozbep, Y.; Aytekin, A.M.; Kaynas, S.; Belen, A.; Kasap, O.E.; Yaman, M.; Rastgeldi, S. Distribution and altitudinal structuring of phlebotomine sand flies (Diptera: Psychodidae) in southern Anatolia, Turkey: Their relation to human cutaneous leishmaniasis. J. Vector Ecol. 2007, 32, 269–279. [Google Scholar] [CrossRef]
  36. Tichy, H.; Kallina, W. Insect hygroreceptor responses to continuous changes in humidity and air pressure. J. Neurophysiol. 2010, 103, 3274–3286. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  37. Leschnik, M.; Löwenstein, M.; Edelhofer, R.; Kirtz, G. Imported non-endemic, arthropod-borne and parasitic infectious diseases in Austrian dogs. Wien. Klin. Wochenschr. 2008, 120, 59–62. [Google Scholar] [CrossRef] [PubMed]
  38. Aspöck, H.; Walochnik, J. When sandflies move north. Public Health 2009, 20, 24–31. [Google Scholar]
  39. Aspöck, H. Postglacial formation and fluctuations of the biodiversity of Central Europe in the light of climate change. Parasitol. Res. 2008, 103, 10–13. [Google Scholar] [CrossRef] [PubMed]
  40. Depaquit, J. Molecular systematics applied to Phlebotomine sandflies: Review and perspectives. Infect. Genet. Evol. 2014, 28, 744–756. [Google Scholar] [CrossRef] [PubMed]
  41. Depaquit, J.; Ferté, H.; Léger, N.; Lefranc, F.; Alves-Pires, C.; Hanafi, H.; Maroli, M.; Morillas-Márquez, F.; Rioux, J.A.; Svobodova, M.; et al. ITS 2 sequences heterogeneity in Phlebotomus sergenti and Phlebotomus similis (Diptera, Psychodidae): Possible consequences in their ability to transmit Leishmania tropica. Int. J. Parasitol. 2002, 32, 1123–1131. [Google Scholar] [CrossRef]
  42. QGIS Development Team QGIS Geographic Information System. Open Source Geospatial Foundation Project. Available online: http://qgis.osgeo.org (accessed on 29 October 2019).
  43. Lewis, D.J. A taxonomic review of the genus Phlebotomus (Diptera: Psychodidae). Bull. Br. Mus. Nat. Hist. 1982, 45, 121–209. [Google Scholar]
  44. Kumar, S.; Stecher, G.; Li, M.; Knyaz, C.; Tamura, K. MEGA X: Molecular evolutionary genetics analysis across computing platforms. Mol. Biol. Evol. 2018, 35, 1547–1549. [Google Scholar] [CrossRef] [PubMed]
  45. Dvořák, V.; Halada, P.; Hlacvackiva, K.; Dokianakis, E.; Antoniou, M.; Volf, P. Identification of phlebotomine sand flies (Diptera: Psychodidae) by matrix-assisted laser desorption/ionization time of flight mass spectrometry. Parasit. Vectors 2014, 7, 21. [Google Scholar] [CrossRef] [Green Version]
  46. Halada, P.; Hlavackova, K.; Risueño, J.; Berriatua, E.; Volf, P.; Dvořák, V. Effect of trapping method on species identification of phlebotomine sandflies by MALDI-TOF MS protein profiling. Med. Vet. Entomol. 2018, 32, 388–392. [Google Scholar] [CrossRef]
  47. Dvořák, V.; Tsirigotakis, N.; Pavlou, C.; Dokianakis, E.; Akhoundi, M.; Halada, P.; Volf, P.; Depaquit, J.; Antoniou, M. Sand fly fauna of Crete and the description of Phlebotomus (Adlerius) creticus n. sp. Parasit. Vectors 2020, 13, 547. [Google Scholar] [CrossRef]
  48. El Tai, N.O.; Osman, F.O.; El Far, M.; Presber, W.; Schönian, G. Genetic heterogeneity of ribosomal internal transcribed spacer in clinical samples of Leishmania donovani spotted on filter paper as revealed by single-strand conformation polymorphisms and sequencing. Trans. R. Soc. Trop. Med. Hyg. 2000, 94, 575–579. [Google Scholar] [CrossRef]
  49. Librado, P.; Rozas, J. DnaSP v5: A software for comprehensive analysis of DNA polymorphism data. Bioinformatics 2009, 25, 1451–1452. [Google Scholar] [CrossRef] [Green Version]
  50. Bandelt, H.-J.; Forster, P.; Röhl, A. Median-joining networks for inferring intraspecific phylogenies. Mol. Biol. Evol. 1999, 16, 37–48. [Google Scholar] [CrossRef]
  51. Leigh, J.W.; Bryant, D. Popart: Full-feature software for haplotype network construction. Methods Ecol. Evol. 2015, 6, 1110–1116. [Google Scholar] [CrossRef]
  52. R Core Team R: A Language and Environment for Statistical Computing. R Foundation for Statistical Computing: Vienna, Austria. Available online: https://www.R-project.org/ (accessed on 5 May 2020).
Pathogens 09 01032 g001 550
Figure 1. Distribution map of all published Ph. mascittii and Ph. simici records in Austria including all surveyed locations presented in this study. The climate map of Austria shows mean annual temperatures (1971–2000). Locations sampled within this study are marked with digits (Table S1), prior published records are marked with letters, (a) Naucke et al. 2011 [5], (b) Poeppl et al. 2013 [6], (c) Kniha et al. 2020 [9].
Figure 1. Distribution map of all published Ph. mascittii and Ph. simici records in Austria including all surveyed locations presented in this study. The climate map of Austria shows mean annual temperatures (1971–2000). Locations sampled within this study are marked with digits (Table S1), prior published records are marked with letters, (a) Naucke et al. 2011 [5], (b) Poeppl et al. 2013 [6], (c) Kniha et al. 2020 [9].
Pathogens 09 01032 g001
Pathogens 09 01032 g002 550
Figure 2. Phlebotomus mascittii characters for morphological identification. Pharynx with large irregular teeth (a), aedeagus with distal cup-like expansion (b), spermatheca not ringed, without neck indicated by black arrow (c), spermatheca under 330–380 nm light (d), spermatheca under 420–490 nm light (e), spermatheca under 585 nm light (f), tip of spermatheca with missing head is indicated by white arrow (d–f).
Figure 2. Phlebotomus mascittii characters for morphological identification. Pharynx with large irregular teeth (a), aedeagus with distal cup-like expansion (b), spermatheca not ringed, without neck indicated by black arrow (c), spermatheca under 330–380 nm light (d), spermatheca under 420–490 nm light (e), spermatheca under 585 nm light (f), tip of spermatheca with missing head is indicated by white arrow (d–f).
Pathogens 09 01032 g002
Pathogens 09 01032 g003 550
Figure 3. Matrix-assisted laser desorption/ionization-time of flight (MALDI-TOF) mass spectra of Ph. mascittii specimens originating from three different Austrian localities (Ro = Rohrau, Up = Unterpurkla, Hu = Hummersdorf).
Figure 3. Matrix-assisted laser desorption/ionization-time of flight (MALDI-TOF) mass spectra of Ph. mascittii specimens originating from three different Austrian localities (Ro = Rohrau, Up = Unterpurkla, Hu = Hummersdorf).
Pathogens 09 01032 g003
Pathogens 09 01032 g004 550
Figure 4. Dendrogram generated by cluster analysis of MALDI-TOF MS protein profiles of four geographically distant European populations of Ph. mascittii. As an outgroup, a mass spectrum of Phlebotomus killicki collected in Crete was used.
Figure 4. Dendrogram generated by cluster analysis of MALDI-TOF MS protein profiles of four geographically distant European populations of Ph. mascittii. As an outgroup, a mass spectrum of Phlebotomus killicki collected in Crete was used.
Pathogens 09 01032 g004
Pathogens 09 01032 g005 550
Figure 5. Haplotype networks of Ph. mascittii based on 12 COI sequences originating from 4 countries (a) and based on 16 cytb sequences originating from 6 countries (b).
Figure 5. Haplotype networks of Ph. mascittii based on 12 COI sequences originating from 4 countries (a) and based on 16 cytb sequences originating from 6 countries (b).
Pathogens 09 01032 g005
Pathogens 09 01032 g006 550
Figure 6. Comparison of climatic parameters at positive and negative surveyed locations analyzed by Kruskal–Wallis test. Mean annual temperature (a), mean June temperature (b), mean July temperature (c) and mean August temperature (d), parameters on temperature (e,f) and on relative humidity (RH) (g,h). Statistical differences between absence and presence are shown as not significant (n.s.), or as significant by asterisks (*** p < 0.001).
Figure 6. Comparison of climatic parameters at positive and negative surveyed locations analyzed by Kruskal–Wallis test. Mean annual temperature (a), mean June temperature (b), mean July temperature (c) and mean August temperature (d), parameters on temperature (e,f) and on relative humidity (RH) (g,h). Statistical differences between absence and presence are shown as not significant (n.s.), or as significant by asterisks (*** p < 0.001).
Pathogens 09 01032 g006
Table
Table 1. Sand fly positive locations surveyed in Austria.
Table 1. Sand fly positive locations surveyed in Austria.
IdLocation
(Federal State)		Latitude,
Longitude		Altitude
(m.a.s.l.)		Trapping Site
(Potential Host) a		Specimen
(Male/Female)		SpeciesTrap Date bTraps Set
(Type)
47Kaiser-steinbruch
(Burgenland)		47.9911,
16.7103		179 moutside, dog pound
(dog)		1
(0m/1f)		Ph. mascittii11.07.131
(CDC light)
1
(0m/1f)		Ph. mascittii26.07.131
(CDC light)
64Neuhaus/Klausenbach
(Burgenland)		46.8683,
16.0232		274 moutside, human dwelling
(rodents possibly)		1
(0m/1f)		Ph. mascittii30.07.191
(BG sentinel & CO2)
1
(0m/1f)		Ph. mascittii25.08.191
(BG sentinel & CO2)
16Neu Albern
(Vienna)		48.1676,
16.4801		157 moutside, roofed shed with hay
(horse, dog)		1
(1m/0f)		Ph. mascittii17.07.191
(CDC light & dry ice)
25Laafeld
(Styria)		46.6867,
16.0069		207 moutside, chicken shed
(chicken, geese, deer)		1
(0m/1f)		Ph. mascittii06.08.181
(CDC light & dry ice)
2
(0m/2f)		Ph. mascittii05.08.191
(CDC light & dry ice)
48Hummers-dorf
(Styria)		46.7076,
15.9812		209 minside and outside, barn
(chicken, dog)		6
(0m/6f)		Ph. mascittii05.07.151
(CDC light)
8
(1m/7f)		Ph. mascittii06.07.152
(CDC light)
49Bad Radkersburg
(Styria)		46.7012,
15.9756		209 minside, old chicken shed
(no)		2
(0m/2f)		Ph. mascittii06.07.151
(CDC light)
54Unterpurkla
(Styria)		46.7319,
15.9062		223 minside, barn
(chicken, cat)		4(0m/4f)Ph. mascittii06.07.151
(CDC light)
a potential host present within a 50 m radius. b date of the evening when trap was set.
Table
Table 2. Interspecific mean genetic distances (%) of the cytochrome c oxidase subunit 1 (COI) and cytochrome b (cytb) between the five Transphlebotomus species based on the Kimura-2-parameter model (COI/cytb). Diagonal bold values indicate intraspecific mean distances.
Table 2. Interspecific mean genetic distances (%) of the cytochrome c oxidase subunit 1 (COI) and cytochrome b (cytb) between the five Transphlebotomus species based on the Kimura-2-parameter model (COI/cytb). Diagonal bold values indicate intraspecific mean distances.
Species12345
1Phlebotomus mascittii0.09/0.3
2Phlebotomus canaaniticus10.7/14.7−/− a
3Phlebotomus anatolicus10.9/13.96.6/8.70.2/0.2
4Phlebotomus killicki12.2/12.911.4/13.712.1/12.70.6/1.4
5Phlebotomus economidesi15.4/13.813.2/13.314.5/12.813.4/9.9− a/1.9
a only one sequence available, included sequences are shown in Tables S2 and S3.
Table
Table 3. Included Ph. mascittii specimens to the haplotype network analysis based on cytochrome c oxidase subunit 1 (COI) and cytochrome b (cytb) DNA sequences.
Table 3. Included Ph. mascittii specimens to the haplotype network analysis based on cytochrome c oxidase subunit 1 (COI) and cytochrome b (cytb) DNA sequences.
COIcytb
Country, LocationGenBankHaplotypeGenBankHaplotypeReference
France, Cévennes--KR336654.1Cytb_2Kasap et al. (2015)
Corsica, Porto Vecchio 1--KR336655.1Cytb_3Kasap et al. (2015)
Corsica, Porto Vecchio 2--KR336656.1Cytb_1Kasap et al. (2015)
Belgium, Saint-Cécile--KR336656.1Cytb_1Kasap et al. (2015)
Germany, Neuenburg--KR336656.1Cytb_1Kasap et al. (2015)
France, Haute-Pyrénées--HQ023281.1Cytb_1Mahamdallie et al. (2011)
Serbia, VojdovinaKY848831.1COI_1--Vaselek et al. (2017)
Serbia, KrasavaMN003381.1COI_2MK991774.1Cytb_1Vaselek et al. (2019)
Slovakia, PernekKX963380.1COI_1--Dvořák et al. (2016)
Slovenia, Truske--MG800323.1Cytb_4Praprotnik et al. (2019)
Slovenia, CetoreKX981916.1COI_2MG800324.1Cytb_5Hlavackova (GenBank)/Praprotnik et al. (2019)
Slovenia, Velike ZabljeKX869078.1COI_2--Hlavackova (GenBank)
Austria, Burgenland, KaisersteinbruchMN812827.1COI_2MN812832.1Cytb_1present study
Austria, Burgenland, Neuhaus/KlausenMN812828.1COI_2MN812833.1Cytb_1present study
Austria, ViennaMN812829.1COI_1MN812834.1Cytb_1present study
Austria, Styria, LaafeldMN812830.1COI_2MN812835.1Cytb_1present study
Austria, Styria, HummersdorfMT332686.1COI_1MT332689.1Cytb_1present study
Austria, Styria, Bad RadkersburgMT332687.1COI_1MT332690.1Cytb_1present study
Austria, Styria, UnterpurklaMT332688.1COI_2MT332691.1Cytb_1present study
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Kniha, E.; Dvořák, V.; Halada, P.; Milchram, M.; Obwaller, A.G.; Kuhls, K.; Schlegel, S.; Köhsler, M.; Poeppl, W.; Bakran-Lebl, K.; et al. Integrative Approach to Phlebotomus mascittii Grassi, 1908: First Record in Vienna with New Morphological and Molecular Insights. Pathogens 2020, 9, 1032. https://doi.org/10.3390/pathogens9121032

AMA Style

Kniha E, Dvořák V, Halada P, Milchram M, Obwaller AG, Kuhls K, Schlegel S, Köhsler M, Poeppl W, Bakran-Lebl K, et al. Integrative Approach to Phlebotomus mascittii Grassi, 1908: First Record in Vienna with New Morphological and Molecular Insights. Pathogens. 2020; 9(12):1032. https://doi.org/10.3390/pathogens9121032

Chicago/Turabian Style

Kniha, Edwin, Vít Dvořák, Petr Halada, Markus Milchram, Adelheid G. Obwaller, Katrin Kuhls, Susanne Schlegel, Martina Köhsler, Wolfgang Poeppl, Karin Bakran-Lebl, and et al. 2020. "Integrative Approach to Phlebotomus mascittii Grassi, 1908: First Record in Vienna with New Morphological and Molecular Insights" Pathogens 9, no. 12: 1032. https://doi.org/10.3390/pathogens9121032

APA Style

Kniha, E., Dvořák, V., Halada, P., Milchram, M., Obwaller, A. G., Kuhls, K., Schlegel, S., Köhsler, M., Poeppl, W., Bakran-Lebl, K., Fuehrer, H. -P., Volfová, V., Mooseder, G., Ivovic, V., Volf, P., & Walochnik, J. (2020). Integrative Approach to Phlebotomus mascittii Grassi, 1908: First Record in Vienna with New Morphological and Molecular Insights. Pathogens, 9(12), 1032. https://doi.org/10.3390/pathogens9121032

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