Altitude and hillside orientation shapes the population structure of the Leishmania infantum vector Phlebotomus ariasi

Despite their role in Leishmania transmission, little is known about the organization of sand fly populations in their environment. Here, we used 11 previously described microsatellite markers to investigate the population genetic structure of Phlebotomus ariasi, the main vector of Leishmania infantum in the region of Montpellier (South of France). From May to October 2011, we captured 1,253 Ph. ariasi specimens using sticky traps in 17 sites in the North of Montpellier along a 14-km transect, and recorded the relevant environmental data (e.g., altitude and hillside). Among the selected microsatellite markers, we removed five loci because of stutter artifacts, absence of polymorphism, or non-neutral evolution. Multiple regression analyses showed the influence of altitude and hillside (51% and 15%, respectively), and the absence of influence of geographic distance on the genetic data. The observed significant isolation by elevation suggested a population structure of Ph. ariasi organized in altitudinal ecotypes with substantial rates of migration and positive assortative mating. This organization has implications on sand fly ecology and pathogen transmission. Indeed, this structure might favor the global temporal and spatial stability of sand fly populations and the spread and increase of L. infantum cases in France. Our results highlight the necessity to consider sand fly populations at small scales to study their ecology and their impact on pathogens they transmit.