http://orcid.org/0000-0003-2212-3362
Figures
Abstract
The scattered eastern African high mountains harbor a renowned and highly endemic flora, but the taxonomy and phylogeographic history of many plant groups are still insufficiently known. The high-alpine populations of the Geranium arabicum/kilimandscharicum complex present intricate morphological variation and have recently been suggested to comprise two new endemic taxa. Here we aim to contribute to a clarification of the taxonomy of these populations by analyzing genetic (AFLP) variation in range-wide high-alpine samples, and we address whether hybridization has contributed to taxonomic problems. We identified only two genetic groups. One corresponded to G. kilimandscharicum, which has been reported as exclusively high-alpine and confined to the eastern Rift mountains in East Africa. The other corresponded to G. arabicum, reported from lower altitudes on the same mountains as well as from a wide altitudinal span in Ethiopia and on the western Rift mountains in East Africa. The four populations analyzed of a recently described species from the Bale Mts in Ethiopia were admixed, indicating that they result from recent long-distance dispersal of G. kilimandscharicum from East Africa followed by hybridization with local G. arabicum in naturally disturbed habitats. Some admixture between the two genetic groups was also inferred on other mountains, supporting earlier suggestions of introgression based on morphology. We did not find support for recognition of the recently suggested new subspecies of G. arabicum in Ethiopia. Interestingly, the high-alpine G. kilimandscharicum lacked clear geographic structuring, suggesting a recent history of colonization of the different mountains or extensive intermountain gene flow.
Citation: Wondimu T, Gizaw A, Tusiime FM, Masao CA, Abdi AA, Hou Y, et al. (2017) Genetic structuring, dispersal and taxonomy of the high-alpine populations of the Geranium arabicum/kilimandscharicum complex in tropical eastern Africa. PLoS ONE 12(5): e0178208. https://doi.org/10.1371/journal.pone.0178208
Editor: Tzen-Yuh Chiang, National Cheng Kung University, TAIWAN
Received: January 27, 2017; Accepted: May 9, 2017; Published: May 26, 2017
Copyright: © 2017 Wondimu et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.
Data Availability: All relevant data are within the paper and its Supporting Information files.
Funding: This work was funded by The Norwegian Programme for Development, Research and Higher Education (NUFU; project no 2007/1058) to S. Nemomissa and C. Brochmann.
Competing interests: The authors have declared that no competing interests exist.
Introduction
The flora of the scattered African high mountains is renowned in biogeography for its peculiar life forms, many endemics, and close relationships to temperate floras in other parts of the world. Although the tropical afro-alpine flora is quite poor in terms of species number, with only 521 species recognized in a recent enumeration [1], not only the bio- and phylogeographic history but also the taxonomy of many plant groups are still insufficiently known. There is a clear need for renewed efforts including more extensive field sampling to clarify the taxonomy of such groups, to test delimitation of species based on genetic data, and to address the relative importance of processes such as intermountain divergence, long-distance dispersal and hybridization in this fascinating system of fragmented ‘sky islands’ [2–4].
Recent phylogeographic studies have shown that afro-alpine plant species or species complexes can show more or less distinct genetic structuring corresponding to individual mountains or mountain groups, suggesting a long history of isolation after initial colonization(Erica trimera, [2]; Deschampsia cespitosa, [3]; Carex monostachya, [5]). Several studies have shown however that intermountain gene flow may be more common than previously thought(Lobelia gibberoa,[6]; Trifolium cryptopodium, [4]). In some species there is no or only little geographic structuring of the genetic diversity, suggesting recent colonization and/or extensive recent dispersal(Erica arborea,[2]; Koeleria capensis, [3]) and in others, there is clear evidence for intermountain dispersal followed by hybridization between divergent lineages (Carex spp., [5]; Carduus schimperi, [4]).Because the species and species complexes studied to date show a surprising variety of phylogeographic histories, more case studies are needed to assess to what degree there are general patterns in the history of the enigmatic afro-alpine flora, and to assess the influence of processes such as intermountain dispersal and hybridization on the morphology of afro-alpine plants.
A typical example showing intricate morphological variation and poorly resolved taxonomy is provided by the eastern African populations of the Geranium arabicum/kilimandscharicum complex. Typical G. arabicumhas pentagonal leaves and flowers in pairs, whereas G. kilimandscharicum has reniform leaves and solitary flowers. However, morphological intermediates do occur and are thought to result from introgression between the two species [7], and it has been suggested that it might be most appropriate to recognize them at the subspecies level [8].These small perennial herbs grow in montane and alpine grasslands and on rocky ground with open vegetation, and show the typical Geranium mode of short-distance dispersal by ejecting their mericarps after explosive curling of the awns [8].The complex morphological variation in this group resulted in description of several species by early taxonomists, but only two were accepted in the treatments of [9] and [8]:They reported G. kilimandscharicum Engl. as exclusively high-alpine and confined to the eastern Rift mountains in East Africa (Mt Kilimanjaro (type), Mt Meru, Mt Kenya, Aberdare Mts, Mt Elgon; [9]; [7]. The other species, G. arabicum Forssk., is typified from Yemen and reported as widespread in tropical African highlands between 1000 and 4000 m [7]; [8]. The two species co-occur on the eastern Rift mountains in East Africa, where G. arabicumusually occurs at lower altitudes (recorded up to 3950 m, but usually well below this altitude) and G. kilimandscharicum usually at higher altitudes (3200–4400 m; [7]).
Recently, a third species of the complex was tentatively described as endemic to the Bale mountains in Flora of Ethiopia and Eritrea, listed as ‘Geranium sp.’ [8]. This species also has reniform leaves and solitary flowers, but it appears conspicuously different from G. kilimandscharicum by its acaulescent stem and by its exceptionally overlapping leaf lobes, which make the leaves look ‘double’. It is only known from a few high-alpine (4150–4200 m) localities where the ground has been disturbed by activities of the giant mole rat or by frost heaving. In addition, a tentative new endemic subspecies of G. arabicum was listed in the Flora of Ethiopia and Eritrea, referred to as subspecies ‘Ash 1711’ based on a collection made in the Bale Mts [8]. This subspecies was reported from afro-alpine grasslands also on several other mountains in Ethiopia and Eritrea, typically at higher altitudes (2900–4000 m) than subspecies arabicum(1300–3650 m). Interestingly, [8] pointed out that some of the material referred to this subspecies was morphologically similar to the East African G. kilimandscharicum. A third subspecies, G. arabicum ssp.latistipulatum (A.Rich.) Kokwaro, has been recognized from low altitudes (1000–3150 m) in Erica forests both in East Africa and Ethiopia. The main leaf lobes of ssp. latistipulatum are deeply pinnatisect into narrowly oblong segments separated by wide sinuses. The other two subspecies have less dissected leaf lobes and narrow sinuses, but ssp. ‘Ash 1711’ differs from ssp. arabicum by having more reniform (vs pentagonal) leaf outline and wedge-shaped (vs rhombic) outline of the main leaf lobes [8].
Further studies are clearly needed to clarify the variation in the eastern African populations of the Geranium arabicum/kilimandscharicum complex. Here we aim to contribute to a clarification of the taxonomy of the high-alpine populations of this complex by analyzing genetic (AFLP) variation in nearly range-wide samples, and also to address their phylogeographic history. In particular, we ask whether hybridization has contributed to taxonomic problems in this group. For logistic reasons we were not able to cover the lower-altitude subspecies in the widespread G. arabicum, because our sampling was carried out as part of a larger project restricted to high-alpine areas.
Materials and methods
Materials
We carried out field work in 10 mountain systems in Ethiopia, Kenya, Tanzania and Uganda (Fig 1, Table 1).Permits were obtained from the responsible authorities in each country: The Ethiopian Wildlife Conservation Authority (EWCA) for the Simen and Bale Mountains National Parks; the Kenyan Wildlife Service (KWS) for Mount Kenya; the Tanzanian Wildlife Authority (TAWA) for Mt. Kilimanjaro and Mt. Arusha; and the Ugandan Wildlife Authority (UWA) for Mt. Ruwenzori and Mt. Muhavura. Leaf samples were collected from five individual plants within an area of100 m × 100 m, taken to represent a single population, and dried in silica gel. We aimed to cover the morphological variation observed in the alpine zone of each mountain and collected a total of 289 plants from 51populations. On the basis of their morphology, 28 populations were tentatively referred to G. kilimandscharicum, 10 to G. arabicumssp.arabicum, six to G. arabicum ssp. ‘Ash 1711’, and seven to Geranium sp.sensu [8]. We were not able to collect G. arabicum ssp. latistipulatum, which is recorded only from a few subalpine areas in East Africa(below 2800 m, [7]) and Ethiopia (below 3150 m; [8]). Our sampling represented all mountains recorded for the eastern East African endemic G. kilimandscharicum [9]; [7], the single mountain range (Bale Mts) recorded for the tentative Ethiopian endemic Geranium sp. [8], and two of the mountain ranges (Bale and Choke Mts) recorded for the tentative Ethiopian endemic G. arabicum ssp. ‘Ash 1711’ (Fig 1). For G. arabicum ssp. arabicum, our sampling included its high-alpine occurrences along the western branch of the Rift Valley in East Africa (Mt Ruwenzori and Mt Muhavura) and in Ethiopia (Bale and Simen), outside the recorded range of G. kilimandscharicum. In addition, we sampled one population referred to G. arabicum ssp. arabicum from the montane forest zone in Mt Kilimanjaro (< 3000 m; Fig 1). Three of the five plants from each population were pressed and deposited in the following herbaria: one in the Natural History Museum, University of Oslo (O), Norway; one in the National Herbarium, Addis Ababa University (ETH), Ethiopia; and the third voucher was deposited according to country of collection, i.e., in the East African Herbarium (EA), Kenya, at the Sokoine University of Agriculture (SUA), Tanzania, or Makerere University Herbarium (MHU), Uganda.
Colours represent the two main genetic groups inferred in the STRUCTURE analyses (blue: G. arabicum group, red: G. kilimandscharicum group).
P(%): Percentage of polymorphic loci; D: Nei’s gene diversity following [10]; DW: frequency-down-weighted marker value as a measure of genetic rarity following [11].The three populations tentatively identified as Geranium sp. sensu [8] with numbers in bold consisted of plants belonging to both genetic groups. Database number refers to the DNA Bank at the Natural History Museum, University of Oslo.
AFLP analysis
DNA was extracted from silica-dried leaf tissue using a GeneMole® robot and the Mole StripTM Plant DNA Kit (QIAGEN, Nordic, Oslo, Norway) or DNeasyTM Plant Mini Kit (QIAGEN, Valencia, CA). Leaf tissue was ground in 2.0 μL tubes with two tungsten carbide beads for 2 min at 15 Hz in a mixer mill (MM301, Retsch GmbH & Co., Haan, Germany), after which 250 μL of lysis buffer was added, vortexed, spun briefly, incubated on a heat block for 10 min at 65°C, and centrifuged at 14000 Hz for 2 min. 200 μL of the lysate was transferred to new tubes and loaded to the robot, which was set to produce a final elution volume of 100 μL.
About 10% of the samples were extracted twice to test for reproducibility of the markersfollowing [12]. AFLP data was generated following [13],except that the PCR reaction volumes were reduced by 50% and pre-selective PCR products were diluted ten times. Twelve primers pairs were initially tested on one plant from each of eight mountains, and the following three primer combinations, which resulted in many polymorphic and well separated bands, were selected for the final analysis:6FAM-EcoRI-ATG/MseI-CGA, VIC-EcoRI-ACA/MseI-CAC, and NED-EcoRI-AGC/MseI-CTG. For each sample, 2.0 μL 6-FAM, 2.0μL VIC and 3.0μL NED labeled selective PCR products were mixed and added to a master mix of 11.7 μL formamide and 0.3μL GENESCAN ROX 500 internal-lane size standard, denatured at 95°C for 5 min and cooled on ice before run on an ABI3100 sequencer (Applied Biosystems, Foster City, USA).
Data analyses
Markers in the size range 50–500 base pairs (bp)were scored as present (1) or absent (0) using GeneMapper® version 4.0 (Applied Biosystems, Foster City, USA).Peaks of low intensity were only scored when unambiguous. Error rate calculation and data cleaning [12] were done separately for each primer combination, and duplicates were removed and matrices combined prior to further analyses. Nei’s gene diversity (D; estimated as the average proportion of pairwise differences among genotypes; [10]), proportion of polymorphic markers (P%), and genetic distinctiveness or rarity (DW; estimated as frequency-down-weighted markers; [11]) were calculated using AFLPdat [14]. Pairwise genetic similarity between AFLP phenotypes was estimated using Dice’s coefficient of similarity in NTSYSpc 2.1 [15] and visualized using Principal Coordinate Analyses (PCoAs).
Bayesian clustering was carried out in STRUCTURE v 2.3.3 [16]. We used the recessive allele model to accommodate the dominant nature of the AFLP markers [17] and compared the no admixture model with uncorrelated allele frequencies vs the admixture model with correlated allele frequencies. Based on the result from the preliminary analysis, we selected the admixture model with correlated allele frequency for the final analysis. Analyses were performed at the Lifeportal, University of Oslo (http://www.lifeportal.uio.no). Analyses were run with K ranging from 1–10, and for each K, 10 replicate runs with a burn-in period of 200,000 and 1,000,000 iterations were used. We used the R-script STRUCTURE-SUM to summarize the results. Log probability of the data, L(K), as a function of K ranging from 1to 10 and the rate of change in the probability between successive runs, ΔK, were calculated according to [18]. In addition, similarity among different runs for the same K was estimated according to[19]. We used the program CLUMPP [20] to estimate the average individual admixture value among the replicated runs for the selected optimal K, and DISTRUCT [21] to graphically visualize the clustering.
A Neighbor-Net diagram [22] was constructed based on uncorrected p-distance using SplitsTree4 v. 4.12.6 [23], and support for branches was estimated from 1000 bootstrap replicates using TreeCon version 1.3b [24]. Analyses of Molecular Variance (AMOVAs) were performed to investigate partitioning of genetic variation at hierarchical and non-hierarchical levels using ARLEQUIN v. 3 [25]. Significance of genetic differentiation and pairwise population differentiation among mountains were estimated with 1000 permutations. We tested the correlation between geographic and genetic (FST) pairwise distances for the total dataset and for each of the two genetic groups separately, using Mantel test implemented in GeneAlex 6.5[26], using 1000 iterations.
Results
A total of 211 samples and 425 markers, of which 306 (72%) were polymorphic, were kept after cleaning the data. Reproducibility of the markers was 98.1%. In the STRUCTURE analyses, the change in the log probability of the data, L(K), showed highest increment from K = 1 to K = 2 (S1A Fig). The rate of change in the probability between successive Ks, DeltaK, indicated a single, highest peak at K = 2 (S1B Fig), and the similarity among the replicated runs also revealed highest convergence at this point (S1C Fig). We therefore inferredthe optimal partitioning of the data to be intotwo genetic groupsfrom the STRUCTURE analyses (K = 2; Figs 1–3).The G. kilimandscharicum group included all populations referred to the eastern East African endemic G. kilimandscharicum (all populations from Mt. Kenya, Mt. Aberdare, Mt. Elgon and Mt. Meru, and all but one from Mt. Kilimanjaro), as well as some of the plants from the Ethiopian Bale Mts referred to Geranium sp. The G. arabicum group comprised all populations referred to the widespread G. arabicum (all populations from the western East African mountains Muhavura and Ruwenzori and the Ethiopian Simen and Choke mountains, most populations from the Ethiopian Bale Mts, and the single montane forest populationfrom Mt Kilimanjaro), as well as the remaining Bale Mts plants referred to Geranium sp. The plants from the four populations (ET0677, ET0954, ET0974 and ET1030) referred to Geranium sp. were accordingly divided between the two genetic groups, with one complete population (ET0677) and four individual plants from the three other populations placed in the G. arabicumgroup, and with the remaining plants placed in the G. kilimandscharicum group (Table 1). All of the individual Geranium sp. plants showed however some degree of admixture (Fig 1). Some degree of admixture was also found between the two genetic groups in other mountains.
a) Total datasetincluding all 45 populations, with colours representing the two main genetic groups inferred in the STRUCTURE analyses (blue: G. arabicum group, red: G. kilimandscharicum group). b) Subset of the 30 populations belonging to the G. kilimandschariumgroup. c) Subset of the 15 populations belonging to the G. arabicumgroup.
Colours represent the two main genetic groups inferred in the STRUCTURE analyses (blue: G. arabicum group, red: G. kilimandscharicum group). No major branch obtained bootstrap support greater than 50%.
Similar genetic structuring was observed in the PCoA and Neighbor-Net analyses. The two genetic groups inferred from the STRUCTURE analyses were distinguished along the first axis of the PCoA of the total dataset, explaining 13.8% of the total variation (Fig 2A). Also in this analysis, the Bale Mts plants referred to Geranium sp. were divided between the two groups. Separate PCoAs of each genetic group showed large and continuous variation within each of them. In the G. kilimandscharicum group, there was no clear differentiation among mountains, and the Geranium sp. plants grouped closely with plants from Mt Aberdare (Fig 2B). In the G. arabicum group, the two subspecies could not beconsistently separated. The Bale plants referred to ssp. ‘Ash 1711’ were placed at one extreme of axis 1, whereas the Choke plants referred to the same subspecies grouped close to ssp. arabicumplants from the Simen Mts (Fig 2C). The mountains were however more clearly differentiated in this genetic group. Axis 1 corresponded to a division across the Rift Valley in Ethiopia (Bale Mts vs Simen/Choke Mts), and most of the East African plants were separated along axis 2. The Neighbor-Net analysis (Fig 3) was more or less star-shaped and without strongly supported major branches, but with the most distinct division corresponding to the two genetic groups inferred from the STRUCTURE analyses.
The G. arabicum genetic group had higher average within-population gene diversity (D = 0.072) and higher rarity (DW = 2.68) than the G. kilimandscharicumgenetic group (D = 0.051; DW = 1.75; Table 2). The highest diversity and rarity were recorded in the Mt Ruwenzori and Mt Muhavura populations of G. Arabicum(D = 0.061–0.066; DW = 3.20–3.62). Very little diversity was recorded in the Mt Kilimanjaro and Mt Kenya populations of G. Kilimandscharicum(D = 0.031–0.035; Table 2). In the non-hierarchical AMOVA analysis, most of the genetic variation was found within populations (62.32%; Table 3). In hierarchical AMOVAs, 21.53% of the variation was found between the two genetic groups, 24.42% among populations within groups, and 54.05% within populations. When dividing the data according to the three species inferred from morphology, 20.05% of the variation was found among the three groups, 24.65% among populations, and 55.30% within populations.Within each genetic group, much higher proportions of the genetic variation werefound within (66.87% and 71.14%) than among populations (33.13% and 28.86%, FST = 0.3313 and 0.2886). We found highly significant (P< 0.001) genetic differentiation among the ten mountains, with the strongest divergence between the Simen Mts and Mt Elgon (Fig 4). We also found a significant pattern of isolation by distance in the Mantel test for the total dataset (r = 0.4113, S2A Fig) and within each of the two genetic groups (r = 0.5710and r = 0.4719, for G. kilimandscharicum S2B Fig and G. arabicum S2C Fig, respectively, P< 0.001 for all values).
F-statistics are provided for genetic differentiation among groups (FCT), among populations within groups (FSC), and among all populations (FST). All values were significant at P<0.0001.
Discussion
Using genome-wide genetic markers and range-wide sampling, we have shown that the high-alpine populations of the Geranium arabicum/kilimandscharicum complex in the eastern African mountains present quite simple genetic structuring, in spite of their intricate morphological variation which recently led to the description of two new endemic taxa. We identified only two genetic groups, of which one is widespread (G. arabicum), and the other seem to have evolved as a specialized high-alpine ecotype in one group of mountains but not elsewhere (G. kilimandscharicum in eastern East Africa). It appears that the morphological complexity and taxonomic problems in the group have been caused by a combination of four processes: by differentiation into one low-alpine and one high-alpine ecotype in eastern East Africa, by local hybridization between these ecotypes, by at least one episode of long-distance dispersal of the high-alpine ecotype followed by hybridization and establishment of admixed plants in naturally disturbed habitats, and by phenotypic plasticity or parallel evolution in some leaf characters. In the following, we discuss the available evidence for the influence of each of these processes.
Morphological differentiation among the low-altitude and high-altitude populations growing in the same mountains in eastern East Africa has long been recognized e.g. [9]; [7], and the observation of occasional morphological intermediates led [7] to suggest that introgression between them occurs at intermediate altitudes. His hypothesis is supported by our data. Although our sampling was concentrated at high altitudes as part of a larger project restricted to the alpine zone, we sampled one montane forest population on Mt Kilimanjaro that clearly belonged to G. arabicum based on both morphology and AFLP data (Figs 1–3). It nevertheless showed some admixture with G. kilimandscharicum. Several of the G. kilimandscharicum plants sampled from each of the other eastern East African mountains also showed some degree of admixture with G. arabicum (Fig 1). The fact that only modest morphological and genetic divergence occurs between the two taxonomically recognized species, and due to their sympatric occurrence on the eastern East African mountains as well as the apparently frequent introgressive hybridization, we tend to agree with the suggestion of [8] to recognize the high-alpine G. kilimandscharicum as a subspecies of the more widespread and morphologically variable G. arabicum.
The conspicuous and puzzling morphology of the high-alpine populations referred to the tentative new species endemic to the Bale Mountains in Ethiopia seems to be caused by hybridization between G. arabicum and G. kilimandscharicum. These populations are restricted to habitats that are naturally disturbed, typically by the giant mole rat. While other mole rats feed underground, the giant mole rat that is endemic to the Bale Mountains mostly forages above ground, heavily disturbing the soil and clearing the vegetation around their tunnel openings [27]. This creates a type of habitat well-known to be suitable for establishment of plant hybrids, e.g. [28]. We found that all individual plants referred to Geranium sp. showed some degree of admixture between the local G. arabicum and the East African G. kilimandscharicum, suggesting that they originated by hybridization after long-distance dispersal of G. kilimandscharicum from East Africa to the Bale Mountains. Because G. kilimandscharicum is exclusively high-alpine in East Africa, we consider gradual migration of this ecotype across the vast lower-lying Kenyan gap to Ethiopiato be unlikely even under colder climatic periods, when the alpine belt may have been shifted downwards by 1000 m [29]. Notably, long-distance dispersal across the Kenyan gap has also been inferred based on genetic data for several other afro-alpine species ([4–5];). The plants referred to Geranium sp. in our sampling do not seem to represent F1 hybrids, but rather segregating later-generation hybrids or backcrosses, as inferred from their variable assignments to the two genetic groups (Fig 1). Long-distance dispersal of the high-alpine G. kilimandscharicumseems also to have occurred to other mountains, as inferred from the occurrence of some admixed plants in the western Rift mountains (Ruwenzori and Muhavura) and the Ethiopian Simen Mts (Fig 1).
The high-alpine Ethiopian populations referred to the tentative new subspecies of G. arabicum, ssp. ‘Ash 1711’, grouped with ssp. arabicum according to their geographic origin rather than morphology (Fig 2C). This pattern suggests that the leaf characters used to separate this subspecies from ssp. arabicum, leaf outline and shape of leaf lobes, are to some degree phenotypically plastic or subjected to parallel evolution, adding further taxonomic confusion to this plant complex.
The lack of geographic structuring we observed in the high-alpine G. kilimandscharicum across the widely separated eastern Rift mountains (Fig 2B) suggests a recent history of colonization of the different mountains or extensive intermountain gene flow. This finding, as well as the evidence for admixture with G. arabicum in other distant mountain groups, suggest that G. kilimandscharicum has a considerable capacity for long-distance dispersal even if it is obviously adapted to short-distance dispersal in terms of morphology. High capacity for long-distance dispersal in plants adapted to short-distance dispersal has also been inferred in some other afro-alpine plants (Wondimu et al., 2014) and in many arctic-alpine plant species [30–31]. In contrast, genetic diversity was more distinctly structured according to mountains in G. arabicum (Fig 2C), which has similar dispersal adaptations. The Rift Valley in Ethiopia appears as a barrier in G. arabicum, as also found for other species [32]; [29]; [33]; [14]. This result is surprising because G. arabicum is able to grow at lower altitudes and therefore could have been thought to cross lowland barriers via gradual migration during colder climates. It appears that the mode of dispersal and lower limit of the altitudinal range cannot always be used as predictors of the strength of intermountain barriers against gene flow in the afro-alpine flora.
Supporting information
S1 Fig. Summary results of STRUCTURE analyses.
a) Log likelihood of the data,L(K), as a function of K ranging from 1 to 10. b) Mean DeltaK for the rate of change in the probability between successive runs, DeltaK as a function of K, calculated according to [33]. c) Average similarity coefficients for the pairwise comparisons among 10 runs for a given K.
https://doi.org/10.1371/journal.pone.0178208.s001
(DOCX)
S2 Fig. Mantel test for correlation between pairwise genetic-distance (Fst) and geographic-distance (P< 0.001; 1000 permutation).
a) for the total dataset, b) G. kilimandscharicum and c) G. arabicum genetic groups.
https://doi.org/10.1371/journal.pone.0178208.s002
(DOCX)
Acknowledgments
This study formed part of the project ‘AFROALP-II—Afro-alpine ‘sky islands’: genetic versus taxonomic biodiversity, climate change, and conservation’ funded by The Norwegian Programme for Development, Research and Higher Education (NUFU; project no 2007/1058) to S. Nemomissa and C. Brochmann. We thank the other members of the AFROALP-II team for discussions and help during the field work, and Manuel Pimentel who determined/verified the identity of our collections.
Author Contributions
- Conceptualization: TW SN CB.
- Data curation: TW AG YH FT CM AA.
- Formal analysis: TW AG.
- Funding acquisition: CB SN.
- Investigation: TW.
- Methodology: TW AG.
- Project administration: CB SN.
- Resources: CB SN.
- Software: TW AG.
- Supervision: SN CB.
- Validation: CB.
- Visualization: TW.
- Writing – original draft: TW.
- Writing – review & editing: TW AG CB.
References
- 1. Gehrke B, Linder HP. Species richness, endemism and species composition in the tropical Afroalpine flora. Alp Bot 2014; 124: 165–177.
- 2. Gizaw A, Kebede M, Nemomissa S, Ehrich D, Bekele B, Mirré V, et al. Phylogeography of the heathers Erica arborea and E. trimera in the afro-alpine ‘sky islands’ inferred from AFLPs and plastid DNA sequences. Flora 2013; 208: 453–463.
- 3. Masao CA, Gizaw A, Piñeiro R, Tusiime FM, Wondimu T, Abdi AA, et al. Phylogeographic history and taxonomy of some afro-alpine grasses assessed based on AFLPs and morphometry: Deschampsia cespitosa, D. angusta and Koeleria capensis. Alp Bot 2013; 123: 107–122.
- 4. Wondimu T, Gizaw A, Tusiime FM, Masao CA, Abdi AA, Gussarova G, et al. Crossing barriers in an extremely fragmented system: two case studies in the afro-alpine sky island flora. Plant Syst Evol 2014; 300: 415–430.
- 5. Gizaw A, Wondimu T, Mugizi TF, Masao CA, Abdi AA, Popp M, et al. Vicariance, dispersal, and hybridization in a naturally fragmented system: the afro-alpine endemics Carex monostachya and C. runssoroensis (Cyperaceae). Alp Bot 2016; 126: 59–71.
- 6. Kebede M, Ehrich D, Taberlet P, Nemomissa S, Brochmann C. Phylogeography and conservation genetics of a giant lobelia (Lobelia giberroa) in Ethiopian and Tropical East African mountains. Mol Ecol 2007; 16: 1233–1243. pmid:17391409
- 7.
Kokwaro JO. Geraniaceae. In: Milne-Redhead E, Polhill RM eds. Flora of Tropical East Africa. London, Great Britain: Whitefriars 1971; 1: 1–10.
- 8.
Gilbert MG, Vorster P. Geraniaceae. In: Edwards S, Tadesse M, Demissew S, Hedberg I eds. Flora of Ethiopia & Eritrea, Vol. 2, Part 1. Addis Ababa: The National Herbarium, Addis Ababa University, Ethiopia. 2000; 346–369.
- 9. Hedberg O. Afroalpine vascular plants. A taxonomic revision. Symb Bot Upps 1957; 15: 1–411.
- 10. Kosman E. Nei’s gene diversity and the index of average differences are identical measures of diversity within populations. Plant Pathol 2003;52: 533–535.
- 11. Schönswetter P, Tribsch A. Vicariance and dispersal in the alpine perennial Bupleurum stellatum L. (Apiaceae). Taxon 2005; 54: 725–732.
- 12. Bonin A, Bellemain E, Eidesen PB, Pompanon F, Brochmann C, Taberlet P. How to track and assess genotyping errors in population genetics studies. Mol Ecol 2004; 13: 3261–3273. pmid:15487987
- 13. Gaudeul M, Taberlet P, Till-Bottraud I. Genetic diversity in an endangered alpine plant, Eryngium alpinum L. (Apiaceae), inferred from amplified fragment length polymorphism markers. Mol Ecol 2000; 9: 1625–1637. pmid:11050557
- 14. Ehrich D. AFLPdat: a collection of R functions for convenient handling of AFLP data. MolEcol Notes 2006; 6: 603–604.
- 15.
Rohlf F. NTSYSpc: Numerical Taxonomy and Multivariate Analysis System. Version 2.11a. 2000; Setauket, New York: Exeter Software.
- 16. Pritchard JK, Stephens M, Donnelly P. Inference of population structure using multilocus genotype data. Genetics 2000; 155: 945–959. pmid:10835412
- 17. Falush D, Stephens M, Pritchard J. Inference of population structure using multilocus genotype data: dominant markers and null allele. Mol Ecol Notes 2007; 7: 574–578. pmid:18784791
- 18. Evanno G, Regnaut S, Goudet J. Detecting the number of clusters of individuals using the software STRUCTURE: a simulation study. Mol Ecol 2005;14: 2611–2620. pmid:15969739
- 19. Nordborg M, Hu TT, Ishino Y, Jhaveri J, Toomajian C, Zheng H, et al. The pattern of polymorphism in Arabidopsis thaliana. PLoS Biol 2005; 3: e196. pmid:15907155
- 20. Jakobsson M, Rosenberg NA. CLUMPP: a cluster matching and permutation program for dealing with label switching and multimodality in analysis of population structure. Bioinformatics 2007; 23: 1801–1806. pmid:17485429
- 21. Rosenberg NA. DISTRUCT: a program for the graphical disply of population structure. Mol EcolNotes. 2004; 4: 137–138.
- 22. Bryant D, Moulton V. NeighborNet: an agglomerative algorithm for the construction of phylogenetic networks. Mol Biol Evol 2004; 21: 255–265. pmid:14660700
- 23. Huson DH, Bryant D. Application of phylogenetic networks in evolutionary studies. Mol Biol Evol 2006; 23: 254–267. pmid:16221896
- 24. Van de Peer Y, De Wachter R. TREECON for Windows: a software package for the construction and drawing of evolutionary trees for the Microsoft Windows environment. Comput Appl Biosci 1994; 10: 569–570. pmid:7828077
- 25. Excoffier L. & Lischer H.E.L. Arlequin suite ver. 3.5: A new series of programs to perform population genetics analyses under Linux and Windows. Mol Ecol Resour 2010; 10: 564–567. pmid:21565059
- 26. Peakall R, Smouse PE. GenAlEx 6: genetic analysis in Excel. Population genetic software for teaching and research. Mol. Ecol. Notes 2006; 6:288–295.
- 27. Yaba M, Mekonnen T, Bekele A, Malcolm J. Food selection and feeding behavior of giant mole rat (Tachyoryctes macrocephalus, Ruppell, 1842) from the Sanetti Plateau of bale Mountains National Park, Ethiopia. Asian J App Sci 2011; 4: 735–740.
- 28. Brochmann C, Borgen L, Stabbetorp OE. Multiple diploid hybrid speciation of the Canary Island endemic Argyranthemum sundingii (Asteraceae).Plant Syst Evol 2000; 220: 77–92.
- 29. Gottelli D, Marino J, Sillero-Zubiri C, Funk SM. The effect of the last glacial age on speciation and population genetic structure of the endangered Ethiopian wolf (Canis simensis). Mol Ecol 2004; 13: 2275–2286. pmid:15245401
- 30. Alsos IG, Ehrich D, Eidesen PB, Solstad H, Westergaard KB, Schönswetter P, et al. Long-distance plant dispersal to North Atlantic islands: colonization routes and founder effect. AoB Plants 2015; 7: plv036 pmid:25876627
- 31. Gussarova G, Allen GA, Mikhaylova Y, McCormick LJ, Mirré V, Marr KL, et al. Vicariance, long-distance dispersal, and regional extinction–recolonization dynamics explain the disjunct circumpolar distribution of the arctic-alpine plant Silene acaulis. Am J Bot 2015; 102: 1703–1720. pmid:26437887
- 32. Arctander P, Johansen C, Coutellec-Vreto MA. Phylogeography of three closely related African bovids (tribe Alcelaphini). Mol Biol Evol 1999; 16: 1724–1739. pmid:10605114
- 33. Evans BJ, Bliss SM, Mendel SA, Tinsley R. The Rift Valley is a major barrier to dispersal of African clawed frogs (Xenopus) in Ethiopia. Mol Ecol 2011;20: 4216–4230. pmid:21933293