Range expansion and ancestral niche reconstruction in the Mediterranean diving beetle genus Meladema (Coleoptera, Dytiscidae)

Species of the genus Meladema (Dytiscidae, Colymbetinae) are some of the largest macroinvertebrates in the western Palearctic region, being top predators in fishless streams. Two of the three described species, Meladema imbricata (Wollaston, 1871) and Meladema lanio (Fabricius, 1775) are Macaronesian endemics from the Canary Islands and Madeira, respectively, while the third, Meladema coriacea Laporte, 1835, is widely distributed from Morocco and the Iberian Peninsula to Turkey, including the Canary Islands. Previous phylogenetic analysis using only mitochondrial markers revealed the existence of two cryptic lineages within M.coriacea, one restricted to Corsica and the other including the rest of sampled populations. We reconstruct here the evolutionary history of the species of Meladema using a more comprehensive sampling covering its whole geographical range, adding nuclear markers and Bayesian molecular dating. Using environmental niche modelling, we test for possible differences in climatic preferences among lineages and reconstruct their ancestral climatic niche. Our results strongly supported the existence of four monophyletic lineages represented by the three recognized species plus a fourth cryptic lineage with populations of M.coriacea from the Tyrrhenian islands (Corsica, Sardinia and Montecristo). This pattern is not likely to be the result of mitochondrial artefacts due to Wolbachia infection, as all 11 tested individuals were negative for this parasite. Dating analysis placed the origin of Meladema in the Middle Miocene although diversification among extant Meladema lineages started in the early Pleistocene and took place in a relatively short time period. Phylogeographic analysis inferred a continental origin of Meladema, with an independent colonization of the Macaronesian and Mediterranean islands. From the south-western Mediterranean region, the continental M.coriacea expanded its range up to Bulgaria and Turkey in the northern basin and to Tunisia in the southern. Results of niche modelling showed that seasonality is the critical factor in shaping the current distribution of Meladema. Island lineages (M.imbricata, M.lanio and the Tyrrhenian lineage of M.coriacea) occur in sites with low seasonality, within the range of the reconstructed ancestral climatic niche of the genus. On the contrary, continental M.coriacea expanded its range to localities outside the ancestral climatic range of the genus, with a higher seasonality and aridity.