Genetic consequences of cladogenetic vs. anagenetic speciation in endemic plants of oceanic islands

This paper presents for the first time a comparison of the genetic consequences of two different types of speciation in plants of an oceanic island. Genetic data, using two different DNA methods, were obtained from more than 4,000 plants from the two major islands of the Juan Fernández Archipelago (Chile). Results show that some immigrant populations undergo major splitting events and harbor limited genetic diversity within each evolving line. In contrast, other immigrant populations establish and enlarge, but they never split, hence accumulating higher levels of genetic diversity.


Introduction
Oceanic islands have long stimulated biologists to investigate patterns and processes of evolution (e.g. Darwin 1842; Wallace 1881; Whittaker and Ferná ndez-Palacios 2007;Bramwell and Caujapé -Castells 2011). These isolated land masses, far from continental source areas, offer opportunities for determining origins of immigrants and their evolutionary history after establishment. The low probability of long-distance dispersal and successful colonization, the reduction of genetic variation in founding populations and the challenges of adaptation to new environments are all features that combine to affect processes of evolution in island archipelagos, particularly speciation.
One dimension of speciation in island plants that has received considerable attention is adaptive radiation (Carlquist 1974;Whittaker and Ferná ndez-Palacios 2007;Rundell and Price 2009). This is a process that begins with dispersal from the original immigrant population into different habitats on the same or neighbouring island. This isolation leads to divergence of the new segregate populations, each becoming rapidly adapted to divergent habitats (Schluter 2001), such that eventually new species are recognized taxonomically. This general process of speciation is usually diagrammed (Fig. 1) as splitting events or cladogenesis (Rensch 1959). A number of dramatic species complexes have developed in oceanic islands through adaptive radiation, such as illustrated by the lobelioids (Givnish et al. 2009) and silverswords (Carlquist et al. 2003) in Hawaii, Aeonium (Liu 1989;Jorgensen and Olesen 2001) and Echium (Böhle et al. 1996) in the Canary Islands and Scalesia (Eliasson 1974) in the Gálapagos archipelago.
In addition to speciation via adaptive radiation (involving cladogenesis), another process, anagenesis (Fig. 1), has recently been emphasized (Stuessy et al. , 2006Whittaker et al. 2008). Some immigrant populations, especially when arriving on an island with limited ecological opportunity, proliferate in size and accumulate genetic diversity mainly through mutation and recombination. After many generations (perhaps over a million or more years), genetic changes result in different morphology that may be treated as a distinct species. This process has been labelled anagenetic speciation (Stuessy et al. 2006), being one type of progenitor-derivative speciation (Crawford 2010). It has been estimated that at least onequarter of all endemic plant species of oceanic islands have originated via anagenesis (Stuessy et al. 2006).
Some studies have been published on the genetic consequences of cladogenesis in endemic plants of different archipelagos. Bö hle et al. (1996) examined chloroplast sequence variation among endemic species of Echium (Boraginaceae) of the Canary Islands, showing very little nucleotide divergence even though the morphological variation is striking. Likewise, Baldwin (2003) examined internal transcribed spacer regions of nuclear ribosomal DNA (ITS) variation among species of the Hawaiian silverswords (Asteraceae) and again, limited sequence variation was seen. The general result from these, and other studies, is that during cladogenesis, the immigrant population becomes fragmented, with each segment containing a limited range of genetic variation in comparison with the continental progenitor population (Baldwin et al. 1998). Maximum morphological divergence occurs but with low levels of observable genetic diversity (Frankham 1997). There is some evidence (Perugganan et al. 2003) that the genetic changes responsible for the morphological adaptations involve alterations in regulatory rather than structural genes.
Results so far with anagenesis show a strikingly different pattern. Most of the investigations have been done on endemic species of Ullung Island, in which at least 88 % of the endemic species have originated anagenetically (Stuessy et al. 2006). The island is young (1.8 Ma; Kim 1985), of low elevation (,1000 m) and relatively ecologically uniform (Yim et al. 1981). Pfosser et al. (2005), using amplified fragment length polymorphisms (AFLPs), examined island and Japanese populations of Dystaenia takesimana and D. ibukiensis, respectively, and the results showed high levels of genetic variation within D. takesimana in comparison with D. ibukiensis. Similar results have been obtained in assessing the origin of Acer takesimensis and A. okomotoanum (Takayama et al. 2012(Takayama et al. , 2013a. Because there is no partitioning of the immigrant population, it survives and proliferates, during which time it accumulates genetic variation through mutation and recombination. Eventually, the level of genetic diversity may even equal (or surpass) that observed in parental source populations (Stuessy 2007).
Because the above studies have been done on different genera in different island archipelagos, it would be useful to compare the genetic consequences of both types of speciation within groups of the same island system, preferably within the same island. In this fashion, more direct comparisons can be made because the general environment is the same. Important, obviously, is to locate plant groups that have originated via both anagenesis and cladogenesis within the same archipelago. A good choice for examining the genetic consequences of anagenesis and cladogenesis in endemic plants of oceanic islands is the Juan Ferná ndez Archipelago, Chile. Approximately 64 % of the species have originated by cladogenesis and 36 % by anagenesis (Stuessy et al. 2006). From another perspective, it is estimated that 70 % of the colonists to the islands have diverged anagenetically, in contrast to only 30 % that have diverged via adaptive radiation ).
The Juan Fernández Archipelago consists of two major islands (Fig. 2): Robinson Crusoe (¼ Masatierra), located 667 km west of continental Chile at 338S latitude, and Alejandro Selkirk (¼ Masafuera) situated 181 km further westward into the Pacific Ocean. The former is known to be 4 million years old and the latter 1 -2 million years old (Stuessy et al. 1984). At present, these two islands are approximately the same size of 50 km 2 (Stuessy 1995). The flora is small, containing 78 native and 135 endemic vascular plant species (Danton et al. 2006). From a biogeographic standpoint, this setting is particularly favourable for generating initial hypotheses, because the near island (Robinson Crusoe) is also the older one, making it highly probable as the initial site for colonization of most groups. Furthermore, the older island is hypothesized to have been much larger when formed , making it a bigger target for dispersal from the mainland.
Numerous molecular markers now exist for assessing genetic variation within and among populations (Lowe et al. 2004). Amplified fragment length polymorphisms (Vos et al. 1995) have been used effectively to provide an overall evaluation of population genetic diversity (Tremetsberger et al. 2003;López-Sepú lveda et al. 2013a). These are treated as dominant markers and hence cannot be employed to determine allelic frequencies. An appropriate co-dominant and polymorphic marker that does allow allelic calculations are nuclear microsatellites or simple sequence repeats (SSRs). The challenge with this marker is to develop primers for locating sequences within the genome for comparison. Next-generation sequencing (NGS) methods are now available that allow this to be done much more easily and at reasonable cost (Takayama et al. 2011. Numerous successful applications of SSRs have shown their efficacy to reveal genetic variation at the population level (Gleiser et al. 2008;Kikuchi et al. 2009;López-Sepú lveda et al. 2013b).
Studies using AFLPs and SSRs have already been published on a number of endemic taxa of the Juan Fernández Archipelago, representing groups that have undergone speciation via cladogenesis and anagenesis. The largest (and endemic) genus that has been investigated is Robinsonia   Takayama et al. 2015), which has seven species on Robinson Crusoe Island that have originated cladogenetically and one on Alejandro Selkirk Island that has evolved anagenetically. The genus Erigeron (Asteraceae; Ló pez-Sepú lveda et al. 2015) has six species that evolved cladogenetically on the younger island, Alejandro Selkirk. These two genera were selected because Robinsonia has speciated primarily via cladogenesis on the older island, and Erigeron has done so on the younger island. Regarding anagenesis, studies have been The purposes of this article are to (i) summarize published data from AFLP and SSR investigations on endemic species of the genera Drimys, Myrceugenia, Rhaphithamnus, Robinsonia and Erigeron; (ii) compare and contrast differences in genetic diversity in groups that have undergone anagenetic or cladogenetic speciation and (iii) discuss the importance of considering modes of speciation for understanding levels of genetic diversity within endemic species of oceanic archipelagos.

Methods
The data summarized here (Table 1) provide the first comprehensive genetic comparisons (with AFLPs and SSRs) in the Juan Ferná ndez Archipelago of species that have evolved by anagenesis and cladogenesis, based on consistent samplings, laboratory methods and modes of analysis. A number of earlier studies utilizing isozymes and DNA sequences have examined genetic variation in endemic species of these islands (e.g. Crawford et al. 1998Crawford et al. , 2001a, but these investigations were not focussed on comparing modes of speciation. Genera in the present studies were selected for their representation of anagenesis and cladogenesis and for their occurrence on the two islands of different geological ages. The samples were collected during expeditions in February 2010 and 2011 from 1870 individuals in 163 populations in 15 endemic species, hence representing 14 % of the endemic angiosperms in the archipelago. The samples provide very good geographic coverage of populations over the landscape in both islands. The term population, as used here in the sense of sampling, refers to groups of individuals that were clearly delimited spatially in the field. The number of individuals analysed per population ranged from 1 to 31. The voucher data for these samples and details of data gathering and analysis are given in the respective publications. Briefly, the following approaches were used for AFLPs. Four or six selective primer combinations were chosen. Numerous (24 -85) primer trials were run with each genus to determine the best combination of primers for good resolution of individuals and populations. Data were obtained on an automated DNA sequencer (ABI 3130xl, Applied Biosystems, Waltham, MA, USA). Scoring was done using GeneMarker ver. 1.85 (SoftGenetics, State College, PA, USA). For analysis of AFLP data, the programs Arlequin 3.5.1.2 (Excoffier et al. 2005), FAMD ver. 1.25 (Schlü ter and Harris 2006), R-Script AFLPdat (Ehrich 2006) and SPSS ver. 15.0 (SPSS; IBM, Armonk, NY, USA) were employed to determine total number of fragments (TNB), percentage of polymorphic fragments (PPB), Shannon Diversity Index (SDI), average gene diversity over loci (AGDOL) and rarity index (RI).
For SSRs, NGS methods (Takayama et al. 2011) were used to generate 6-12 loci, selected for their repeatability and scoring convenience. Polymerase chain reaction-amplified fragments were also run on the same automated sequencer and scored with GeneMarker ver. 1.85. Data analysis involved using GENEPOP 4.0 (Raymond and Rousset 1995) The overall pattern of higher genetic diversities in anagenetically derived species in comparison with cladogenetically derived ones was examined by a Student's t-test (average TNB, PPB, SDI, AGDOL and RI in AFLPs, and H O , H E , N A and A R5 in SSRs) and shown in Table 2. To improve normality of H O and H E , a square-root transformation was applied. The overall patterns of higher genetic diversities in Robinson Crusoe Island (old) than Alejandro Selkirk Island (new) were also examined in the same way. The effects of two factors (speciation mode and island) and their interaction were analysed in a two-way ANOVA in R version 3.0.0 (R Core Team 2013) and shown in Table 3.
Data from both AFLPs and microsatellites were further analysed by assessing genetic distance (Nei et al. 1983) with the NeighborNet algorithm (Bryant and   For this article, to allow ease of visual comparisons of results among the species, emphasis has been placed on selected graphic presentations. SplitsTree Neighbor-Net was employed with the AFLP data, and the results are given in a series of graphs (Fig. 3). Neighbour-joining based on genetic distance was used for analysis of the SSRs, and simplified networks were used to show relationships among the populations (Fig. 4). For summary comparisons of genetic diversity among species, AGDOL was used with the AFLP data (Fig. 5). Not all calculated values for all original populations are presented or discussed in this review. The reader is referred to the original publications for additional methods and data.

Results
The results from the AFLP and SSR data analyses are given in Tables 1-4 and shown graphically in Figs 3-5. In general, the results from the two sources of genetic data are similar, with some exceptions, reinforcing confidence in the patterns seen. These data will be presented in context of the two modes of speciation, anagenesis and cladogenesis, but with attention also to the different ages of the islands. Robinson Crusoe Island is 4 million years old and Alejandro Selkirk 1 -2 million (Stuessy et al. 1984).

Anagenesis
The results from analysis of species that have evolved anagenetically include those from Myrceugenia fernandeziana, M. schulzei, Robinsonia masafuerae, D. confertifolia and R. venustus. The first species occurs only on the older island, the second and third species only on the younger island and the last two on both islands. A number of points seem evident. First, all anagenetically derived species show considerable levels of genetic diversity (Table 1,  . This is what might be expected from the predictions regarding anagenesis based on previous studies. Even more interesting, perhaps, is that the amount of genetic diversity differs in species on the two islands of different ages. In D. confertifolia, and R. venustus, which occur on both islands, one sees in both cases more genetic diversity (SDI) in populations on the older island than on the younger island except for estimates of SSRs in R. venustus (Table 1). The explanation of these data may relate to the time available for a genetic change to take place. Because Alejandro Selkirk Island is no more than 1 -2 million years old, this must be the maximum  time available for population divergence to take place. With anagenetically evolved species, all factors being equal, genetic variation increases through time, and this can be seen in the species investigated. One case of anagenesis in the archipelago also merits comment. Robinsonia masafuerae is a species that appears to have speciated from R. evenia, with which it has been closely associated in all studies so far (Crawford et al. 1993a;Sang et al. 1995;Takayama et al. 2015). Previous investigations on ITS 1 and 2 in Robinsonia (Sang et al. 1995) have shown sequence divergence between R. evenia and R. masafuerae as only 0.0063 (two base substitutions). Although one cannot place an absolute time on this divergence, it is the lowest level among any pair of species in the genus, which correlates well with the youthful geological age of Alejandro Selkirk Island. Genetic variation in R. masafuerae is much lower from AFLP data than in R. evenia from Robinson Crusoe (Table 1 and Fig. 5), but in SSRs, the pattern reverses with the anagenetically derived species, R. masafuerae, showing more variation than any single one of the cladogenetically originated species on Robinson Crusoe (Table 1).
It is also possible to make comparisons between populations of continental progenitors with endemic island derivatives. In the case of Myrceugenia schulzei, the closest continental congener is M. colchaguensis (Landrum 1981a, b;Ruiz et al. 2004). Although the sampling of populations on the continent is limited to two populations, the amount of genetic diversity is particularly low as shown by AFLP data, although somewhat higher with SSRs    submitted for publication). These results support the concept that over time, considerable genetic variation can accumulate in anagenetically derived populations, so much so that the degree of variation can approximate and even surpass that in the progenitor source populations.

Cladogenesis
Two of the largest genera of the archipelago are Robinsonia with eight endemic species and Erigeron with six. Both are in Asteraceae, although unrelated and placed in different tribes (Senecioneae vs. Astereae, respectively).
Robinsonia has adaptively radiated on Robinson Crusoe Island during the past 4 million years (maximum value) and Erigeron has done so on Alejandro Selkirk Island in the past 1 -2 million years.
Robinsonia is the second largest genus in the archipelago. The largest is Dendroseris, also of Asteraceae but from still another tribe (Cichorieae). This latter genus is of interest as it has derived cladogenetically on the older island with three independent dispersals to the younger island and three anagenetic speciations there (Sanders et al. 1987;Pacheco et al. 1991;Sang et al. 1994). Most of AoB PLANTS www.aobplants.oxfordjournals.org these species are quite rare now, however, which precluded our being able to obtain sufficient population data for genetic evaluation. Robinsonia has eight species, but two are presumed extinct (R. berteroi and R. megacephala; Danton et al. 2006). Our studies have focussed on five species having originated cladogenetically on the older island. Comments have already been made regarding the one anagenetically derived species (R. masafuerae) on Alejandro Selkirk Island. The results from AFLP data are shown in Fig. 3 and from SSRs in Fig. 4. Most notable from the SplitsTree graph in Fig. 3 is that the different species of Robinsonia are very distinct genetically. Divergence has obviously taken place during adaptive radiation and also during a maximum time available of 4 million years. The species R. gayana, R. thurifera and R. saxatilis form an evolutionary complex, which taxonomically has been regarded as sect. Robinsonia (Skottsberg 1922, as sect. Symphyolepis;Takayama et al. 2015). Robinsonia gracilis ties with R. evenia and its close anagenetic relative R. masafuerae in sect. Eleutherolepis (Skottsberg 1922). With SSR data (Fig. 4), the species are also very distinct and genetically more cohesive, with the anagenetic species R. masafuerae showing the greatest genetic diversity.
Another important point seen clearly in Figs 3 and 4 is that the range of genetic diversity within each of these cladogenetic species is limited in comparison with the anagenetically derived species discussed above.
Although Erigeron is not an endemic genus in the archipelago, six endemic species occur there having evolved via cladogenesis and adaptive radiation. The origin of this complex is unusual in that the colonist(s) presumably arrived directly to the younger island (Valdebenito et al. 1992). Amplified fragment length polymorphism and SSR data (Figs 3 and 4) reveal considerable genetic diversity within these endemic species, and each species is reasonably distinct. An exception is the Erigeron ingae complex consisting of E. ingae, E. luteoviridis and E. turricola. These species are sometimes difficult to distinguish morphologically. Solbrig (1962) and Marticorena et al. (1998), for example, placed E. turricola into synonymy with E. ingae, but Danton et al. (2006) kept them distinct. The molecular data parallel this morphological inconsistency. This may be a population complex in early stages of speciation, now undergoing divergence from within a pool of morphological and genetic variation. All of these species grow in the 'alpine zone' on the younger island (Skottsberg 1922), Figure 5. Summary of genetic diversities, AGDOL, within the endemic species of Drimys, Myrceugenia and Rhaphithamnus having originated by anagenesis, and Robinsonia and Erigeron having been derived through cladogenesis. Robinsonia masafuerae from the younger island is also an anagenetic derivative from the cladogenetic complex of Robinsonia on the older island. AS, Alejandro Selkirk Island; RC, Robinson Crusoe Island. White bar indicates an endemic species in RC, grey bar an endemic species in AS and black bar multiple species or islands combined. and we have not noticed any clear habitat differences among them. The species E. rupicola is confined to coastal rocks along the sea and also penetrates into the quebradas (ravines); its close relative, E. stuessyi, is also found on rocky ledges but residing inside the cool and deep ravines. Erigeron fernandezianus occurs in a broad altitudinal range (100-1200 m), and it inhabits mainly rocky areas in middle elevation plains, quebradas and ridges. This species also occurs on the older island, but it is found there in many plant communities and especially in disturbed sites. It appears, therefore, to be an example of back migration from the younger to the older island (Valdebenito et al. 1992;López-Sepú lveda et al. 2015).
Although most species of Erigeron on the younger island are distinct genetically, the degree of distinctness is much less than observed among species of Robinsonia on the older island (Figs 3 and 4). It may be that these species of Erigeron have had less time to diverge in comparison with those of Robinsonia. With the passage of time, therefore, the genetic profiles of species undergoing adaptive radiation may remain narrow due to strong directional selection in each different habitat. In both Erigeron and Robinsonia, however, the range of genetic variation seen is less than that in the anagenetically derived species.

Comparison of anagenesis and cladogenesis
Predictions from theory (Stuessy 2007) would suggest that higher levels of genetic diversity should be found within the anagenetically derived species. This is because the founding population increases in size over time, accumulating genetic diversity mainly through mutation and recombination. One would expect no (or very little) geographic partitioning over the landscape. Likewise, due to a lack of strong selection, one would not expect to find high levels of private alleles or bands, nor a high RI. With cladogenetic speciation, on the other hand, one would expect less overall genetic diversity within each species, but with more private alleles due to strong directional selection. As for impact from the age of the islands, one would predict less total genetic diversity within anagenetically derived species on the younger island because diversity increases through time. As for the cladogenetic species, one would predict less genetic divergence (distinctness) on the younger island in comparison with species on the older island, because directional selection continues over time and refines the genetic profile of each species as it adapts to the particular ecological zone.
Results from genetic analyses of 5 anagenetic species and 10 cladogenetic species allow comparisons between the two modes of speciation and the two islands of differing ages (Tables 1-3). A number of general points can be observed (Table 4 and Fig. 5). First, in anagenetic species, the level of genetic diversity tends to be higher per species than in the cladogenetic species, especially on Robinson Crusoe Island. This can be seen in percentage of polymorphic bands, SDI, AGDOL, observed proportion of heterozygotes, expected proportion of heterozygotes, number of alleles per locus and allelic richness. Second, in the . This is true on both islands of differing ages. This suggests that this pattern can develop easily within 1-2 million years and that it can persist for up to 4 million. This is consistent with the results reported for Ullung Island, Korea, which is known to be 1.8 million years old (Pfosser et al. 2005;Takayama et al. 2012Takayama et al. , 2013a. Third, the ability of an immigrant population to radiate adaptively has much to do with the properties of the colonists (and progenitors) and less with differences of habitat. Some colonists remain as a single larger population and are not responsive to adaptive change in different ecological zones, whereas others disperse well to micro-zones and quickly become modified morphologically and genetically. Fourth, perhaps most importantly, the total amount of genetic diversity within an anagenetically derived species in comparison with an entire adaptively radiating lineage is approximately the same (Fig. 5).

Genetics of speciation in endemic plants of oceanic islands
A number of previous studies have assessed levels of genetic variation within and among populations of endemic species of the Juan Ferná ndez Archipelago with other markers such as isozymes, random amplified polymorphic DNA (RAPDs) and inter simple sequence repeats (ISSRs). Isozymes have been analysed in Dendroseris (Crawford et al. 1987), Chenopodium sanctae-clarae (Crawford et al. 1988), Wahlenbergia , Robinsonia , Lactoris  and Myrceugenia (Jensen et al. 2002).
RAPDs have been investigated in Dendroseris (Esselman et al. 2000) and Lactoris (Brauner et al. 1992), and ISSRs also in Lactoris (Crawford et al. 2001b). Crawford et al. (2001a) summarized the results from isozyme studies on 29 endemic species of the Juan Ferná ndez Archipelago, and this represents the best set of observations to compare with the AFLP and SSR data summarized here. The most conspicuous result is that the mean genetic diversities at the species level are low (H es ¼ 0.065). Higher levels of diversity were seen in larger populations or in many small populations and also in outcrossing species in contrast to selfers. Of relevance for comparisons to the present study, isozymes have been analysed from four species of Robinsonia and in M. fernandeziana, E. fernandezianus and R. venustus. It is difficult to compare the results of the isozymes because they provide less detailed genetic information than from AFLPs and SSRs. Isozyme studies on the endemic Lactoris fernandezianus, for example Crawford et al. (1994), revealed virtually no variation, but ISSRs showed considerable variation within and among populations (Crawford et al. 2001b). Studies on isozymes (Crawford et al. 1987) and RAPDs (Esselman et al. 2000) from Dendroseris showed greater resolution of relationships from the latter. The isozyme data for the four cladogenetically derived species of Robinsonia show higher levels of genetic variation than in the anagenetic R. venustus (Crawford et al. 1993b) and Myrceugenia (Jensen et al. 2002), which would be in contrast to the trends documented here. It is important, therefore, that for questions involving population genetics in endemic plants of oceanic islands, rapidly evolving markers need to be used.
The employment of AFLPs and SSRs in the present study from 15 species of the Juan Ferná ndez Archipelago, therefore, does provide detailed genetic data at the population level for purposes of comparing consequences of different modes of speciation. A general review has recently been published on the general topic of interpretation of genetic variation within endemic species of oceanic islands , and the present data corroborate ideas summarized there. Clearly, the alternative modes of speciation, anagenesis and cladogenesis result in different genetic consequences. Interpretation of the evolutionary significance of levels of genetic diversity, therefore, must be done in context of type of speciation. As can be seen in the results of adaptive radiation in Erigeron and Robinsonia, on the young and older islands, respectively, the geological age of the island also matters, as this provides the time dimension in which the evolutionary processes unfold.
Another very significant impact on levels of genetic variation in populations of endemic plants of oceanic islands is that from human activity. Because oceanic islands often have agreeable climates and attractive beaches, people have come to live, play and build homes and apartments, all of which have caused pressures on the native vegetation. In the Juan Ferná ndez Archipelago, for example, people have been living continuously on Robinson Crusoe Island for .300 years (Woodward 1969;Wester 1991). It is not impossible that the species of Robinsonia on the older island have suffered some genetic loss due to human activity. Although these species occur either on high ridges or in deep forests, far removed from most persons who live at sea level in the village (San Juan Bautista), incursions into the native forest must have taken place and some plants destroyed. It is known that two species of Robinsonia, both on Robinson Crusoe Island, are now extinct (R. berteroi and R. megacephala; Danton and Perrier 2005;Danton et al. 2006). Assessing the level of human impact on the vegetation of an oceanic island, therefore, is challenging. At least in the Juan Ferná ndez Archipelago, there were no aboriginal peoples, and human activity could only have begun with discovery by Europeans (Juan Ferná ndez; Medina 1974) at the end of the 16th century. Since that time, however, considerable negative impact from human activity has been documented in the archipelago (Wester 1991;Matthei et al. 1993;Stuessy et al. 1997;Swenson et al. 1997;Cuevas and Leersum 2001;Greimler et al. 2002;Dirnbö ck et al. 2003;Cuevas et al. 2004;Ricci 2006;Vargas et al. 2011), especially from introduced animals, such as rats, rabbits and goats (e.g. Camus et al. 2008). These combined activities have surely had some impact on the levels of genetic variation within and among populations.

Sources of Funding
This work was supported by an FWF (Austrian Science Fund) grant (P21723-B16) to T.F.S. and a Japan Society for the Promotion of Science (JSPS) Postdoctoral Fellowship for Research Abroad (grant 526) to Ko.T.

Contributions by the Authors
Ko.T. conceived the idea behind the article; all authors participated in the field work except G.K. and Ka.T.; J.N., P.L.-S., G.K. and Ko.T. completed the laboratory work; J.N. coordinated the NGS data acquisition; T.F.S. and Ko.T wrote the initial draft and all authors contributed to subsequent drafts and offered comments for improvement.