Pronounced mito-nuclear discordance and various Wolbachia infections in the water ringlet Erebia pronoe have resulted in a complex phylogeographic structure
Erebia pronoe exhibits highly structured and strongly differentiated mitochondrial lineages, which are consistent with the distribution of previously described morphotaxa and analyses of Dincă et al.10 These genetic lineages are also reflected to varying degrees in the nuclear markers. The observed mito-nuclear discordances can be explained by different evolutionary rates of genetic markers, the effects of Wolbachia infections, and introgression. These aspects are discussed in more detail in the following sections on the phylogeographic history of this species complex.Mito-nuclear discordance and the systematic status of Erebia melas
Based on genital morphology and nuclear markers, E. melas represents a distinct group to E. pronoe. The common area of origin of both species was probably located in the eastern Alps, which is supported by a RASP analysis based on the nuclear markers. However, E. melas acts as an ingroup of E. pronoe based on the mitochondrial markers, and a RASP analysis indicates a common origin for both taxa in the Carpathian region. Since most Erebia species in Europe have at least parts of their distribution in the Alps21 and are adapted to Alpine environments and habitats22,23, we consider an eastern Alpine origin of the ancestor of E. pronoe and E. melas more likely. This hypothesis subsumes the assumption that the genetic proximity on the mitochondrial level was probably caused by hybridisation and introgression events, which could have occurred as a result of several eastward advances of E. pronoe to the Balkan Peninsula (see below). This seems plausible, because the ability and tendency of E. pronoe to hybridise with other Erebia species have been demonstrated repeatedly12,24,25.The existence of Wolbachia strain 2 in both species, and its distribution from the Pyrenees (in E. pronoe) to the Balkan Peninsula (in E. melas) also speaks for a common origin of both species. Thus, Wolbachia strain 2 might represent the ancient strain present in the common ancestor of this species group, surviving today at the geographic margins (i.e. Pyrenees, western Alps, Balkan Peninsula), but which at some time was replaced in the centre of the butterfly’s range (i.e. the eastern and central Alps) by strain 1. The link between co-occurrence in a common area and prevalence of one Wolbachia strain was also recently demonstrated in other Erebia species26 and might facilitate mitochondrial introgression27.Intraspecific differentiation and glacial refugia of Erebia pronoe
The Pyrenean region is inhabited by one of the oldest and most differentiated intraspecific lineages of E. pronoe. The high genetic diversity in the Pyrenees speaks for large effective population sizes throughout time, enabled by mostly altitudinal shifts in response to climatic cycles, and a lack of major genetic bottlenecks. Compared to the Pyrenean group, the genetic diversity of the western Alpine populations, also well differentiated from all other groups, is lower. This lower diversity was probably the result of repeated cold stage retreat to a geographically more restricted refugium at the foot of the south-western Alps, a well-known refugial area for numerous species28.We cannot say conclusively whether the populations in the Pyrenean region or in the western Alps differentiated first, due to the contradictory genetical markers. The higher evolutionary rate of the mitochondrial markers, the allopatric distribution, and the hybridisation with diverse Erebia species may have led to a greater differentiation of the Pyrenees and/or a loss of the genetic link between the western Alps and the Pyrenees. Since a link between the western Alps and the Pyrenees is still well reflected in the nuclear data set and by the shared Wolbachia strain 2, we consider the most likely scenario to be an early Pleistocene or even Pliocene expansion from the western Alps to the Pyrenees, with subsequent isolation and differentiation. Thus, the Pyrenees-western Alps populations might first have separated as one group from an eastern Alps group s.l., as suggested by nuclear information, and not in two independent events, as suggested by mitochondrial genes.Simultaneously to the split between western Alps and Pyrenees, a separation of the eastern Alpine group s.l. into a southern Alpine subgroup and an eastern Alpine subgroup should have occurred. The southern Alpine subgroup displays a high genetic diversity in their nuclear markers, but a significantly lower diversity in the mtDNA. This might be explained by the existence of a cold-stage refugial area in the southern Alps or their margin, supporting the constant survival of large populations, but also a reshaping of the mtDNA patterns through introgression from the eastern Alpine subgroup during secondary contact when both subgroups expanded into formerly glaciated east-central Alpine areas. The isolated occurrence of Wolbachia strain 1 and mitochondrial haplotypes H29 and H30 (shared with the eastern Alps subgroup) in the southern Alps further support the hypothesis of gene flow from the eastern Alpine region into the southern Alpine populations and vice versa.The eastern Alpine subgroup probably survived glacial periods in a large, cohesive refugium at the eastern edge of the Alps, as has been demonstrated for numerous other species28. This area is also seen as a potential centre of origin of the entire taxon. From there, a recent (most likely postglacial) dispersal must have taken place, which should be responsible at least partly for the star-like pattern of this group in both mitochondrial and nuclear haplotype networks. However, further dispersal events out of the eastern Alps during previous interglacials and maybe even going back to the Pliocene have to be postulated to explain the entire range dynamics in E. pronoe.Apparently, multiple advances out of the eastern Alps into the Balkan mountain systems have taken place from several independent glacial refugia in the region, as indicated by the different mtDNA lineages in Slovenia, western Balkan mountains, and eastern Balkan mountains. A separation between the eastern and western Balkans, and hence also separate glacial refugia in both areas, was frequently observed for mountain taxa28,31. This pattern may have resulted from a succession of independent dispersal events from the eastern Alps throughout the younger Pleistocene, with subsequent regional extinction events and/or independent dispersal events across the Carpathians, as has been demonstrated for numerous other species29.A similar pattern of two independent colonisation events also applies to the Carpathians. Thus, the highly isolated populations in the south-eastern Carpathians must go back to an older expansion out of the eastern Alps. This probably took place during one of the last interglacial phases. The route most likely followed the Carpathian arc, but only a few populations survived at their south-eastern edge. This underlines the phylogeographic independence of this part of the Romanian Carpathians, which is also supported by studies on numerous other mountain species30,31,32. On the other hand, the Tatra mountains, as the northernmost part of the Carpathians, were colonised very recently, most likely postglacially, out of the eastern Alpine area. The strong and rather recent link between these two areas is also supported by phylogeographic studies on many taxa30,33,34.Because of the slower evolutionary rate of nuclear DNA and the resulting incomplete lineage sorting, nuclear markers can contribute little to the reconstruction of these presumably recent events. In line with that, the Valais lineage also has little nuclear differentiation but is clearly distinguished from the western and eastern Alpine lineages by the exclusive mtDNA haplotype H17 and Wolbachia strain 3. The presence of a single, highly differentiated mtDNA haplotype and an exclusive Wolbachia strain indicates a selective sweep. This lineage most likely represents a chronological relict of an interglacial expansion of the eastern Alpine subgroup to the western-central Alps surviving since then in this area, finding glacial refugia in nearby unglaciated areas and becoming infested by a Wolbachia strain not present in any other E. pronoe lineage, hence accelerating its differentiation.Another selective sweep was probably the cause of the mito-nuclear unconformity in the southern Alps lineage. The occurrence of the mtDNA haplotypes H29 and H30 and the Wolbachia strain 1 indicate mitochondrial hybridisation between the eastern and southern Alpine lineages during an expansive interglacial phase. As a result, Wolbachia infection probably occurred, which might have impoverished the mitochondrial diversity of the southern Alps lineage.Consequences for subspecific differentiation in Erebia pronoe
In general, the support given by our data for the so-far described subspecies decreases from west to east. Erebia pronoe glottis Fruhstorfer, 1920, distributed in the Pyrenees, represents the best-supported subspecies. Fixed mitochondrial amino acid changes emphasize the distinctness of this taxon, which might be well advanced in the process of speciation; we cannot even exclude the possibility that it has already reached full species rank. The genetic separation of the western Alps from the Valais, geographically separated along the main Alpine ridge, justifies the recognition of the taxa E. pronoe vergy (Ochsenheimer, 1807) and E. pronoe psathura Fruhstorfer, 1920, respectively, and is supported by both marker sets as well as by the existence of two different Wolbachia strains. The eastern Alpine subgroup resembles the nominotypical E. pronoe pronoe. The existence of at least one lineage in the southern Alpine area is supported by both marker sets. A finer separation based on the mitochondrial markers is not possible, because of recent introgression events affecting east Alpine haplotypes, as also indicated by the existence of Wolbachia strain 1. This population group could be assigned to the taxon E. pronoe gardeina Schawerda, 1924, or to E. pronoe tarcenta Fruhstorfer, 1920, considering their ranges. Nevertheless, a final decision requires further regional studies. Erebia pronoe fruhstorferi Warren, 1933 was accepted to be widely distributed in the Balkan mountain systems. However, our data suggest independent lineages in the western and eastern Balkan mountain systems of which only the eastern populations can be assigned to this taxon. The lineage of the Slovenian Alps is primarily based on mitochondrial markers and morphological characteristics7. The existence of an independent lineage for the highly isolated populations in the southern Carpathians, justifies the subspecies status of E. pronoe regalis Hormuzachi, 1937. Both marker sets display a differentiation, which was more pronounced in the nuclear than in the mitochondrial DNA. More
