Thermal niche evolution and geographical range expansion in a species complex of western Mediterranean diving beetles

Background: Species thermal requirements are one of the principal determinants of their ecology and biogeography, although our understanding of the interplay between these factors is limited by the paucity of integrative empirical studies. Here we use empirically collected thermal tolerance data in combination with molecular phylogenetics/phylogeography and ecological niche modelling to study the evolution of a clade of three western Mediterranean diving beetles, the Agabus brunneus complex. Results: The preferred mitochondrial DNA topology recovered A. ramblae (North Africa, east Iberia and Balearic islands) as paraphyletic, with A. brunneus (widespread in the southwestern Mediterranean) and A. rufulus (Corsica and Sardinia) nested within it, with an estimated origin between 0.60-0.25 Ma. All three species were, however, recovered as monophyletic using nuclear DNA markers. A Bayesian skyline plot suggested demographic expansion in the clade at the onset of the last glacial cycle. The species thermal tolerances differ significantly, with A. brunneus able to tolerate lower temperatures than the other taxa. The climatic niche of the three species also differs, with A. ramblae occupying more arid and seasonal areas, with a higher minimum temperature in the coldest month. The estimated potential distribution for both A. brunneus and A. ramblae was most restricted in the last interglacial, becoming increasingly wider through the last glacial and the Holocene. Conclusions: The A. brunneus complex diversified in the late Pleistocene, most likely in south Iberia after colonization from Morocco. Insular forms did not differentiate substantially in morphology or ecology, but A. brunneus evolved a wider tolerance to cold, which appeared to have facilitated its geographic expansion. Both A. brunneus and A. ramblae expanded their ranges during the last glacial, although they have not occupied areas beyond their LGM potential distribution except for isolated populations of A. brunneus in France and England. On the islands and possibly Tunisia secondary contact between A. brunneus and A. ramblae or A. rufulus has resulted in introgression. Our work highlights the complex dynamics of speciation and range expansions within southern areas during the last glacial cycle, and points to the often neglected role of North Africa as a source of European biodiversity.