Skip to main content

Genetic structuring of remnant forest patches in an endangered medicinal tree in North-western Ethiopia



Habitat loss and fragmentation may have detrimental impacts on genetic diversity, population structure and overall viability of tropical trees. The response of tropical trees to fragmentation processes may, however, be species, cohort or region-specific. Here we test the hypothesis that forest fragmentation is associated with lower genetic variability and higher genetic differentiation in adult and seedling populations of Prunus africana in North-western Ethiopia. This is a floristically impoverished region where all but a few remnant forest patches have been destroyed, mostly by anthropogenic means.


Genetic diversity (based on allelic richness) was significantly greater in large and less-isolated forest patches as well as in adults than seedlings. Nearly all pairwise FST comparisons showed evidence for significant population genetic differentiation. Mean FST values were significantly greater in seedlings than adults, even after correction for within population diversity, but varied little with patch size or isolation.


Analysis of long-lived adult trees suggests the formerly contiguous forest in North-western Ethiopia probably exhibited strong spatial patterns of genetic structure. This means that protecting a range of patches including small and isolated ones is needed to conserve the extant genetic resources of the valuable forests in this region. However, given the high livelihood dependence of the local community and the high impact of foreign investors on forest resources of this region, in situ conservation efforts alone may not be helpful. Therefore, these efforts should be supported with ex situ gene conservation actions.


Persistent dependence of human livelihoods on agriculture, tropical rainforests and products has remained the major cause of overharvesting, deforestation, forest degradation and land use alterations in most developing tropical regions [1, 2]. Such anthropogenic factors in turn have led to forest loss (the loss of a certain amount of forest habitat from a region) and fragmentation (the sub-division of a formerly large unbroken forest into several forest patches of smaller but variable sizes and spatial configurations) [3]. Thus habitat loss and fragmentation modify landscape composition and spatial configuration [3]. Also, following forest habitat loss and fragmentation, significant modifications in habitat qualities and biochemical conditions may be expected specifically around the edges of remnant forest patches [4]. The changes in landscape composition, spatial configuration, habitat qualities and biochemical conditions may cause negative effects on biotic interactions, reproductive processes, dispersal and species persistence in human modified landscapes [411]. Fragmentation may also disrupt gene flow among populations and ultimately contribute to the formation of small and isolated populations, which are more vulnerable to extinction due to demographic, genetic and environmental stochastic events [10, 12].

Conservation genetics studies are expected to provide information that enables a better understanding of the ways in which human disturbance and fragmentation processes affect population viability in fragmented landscapes and this may assist in devising more appropriate management options [11, 13]. Small and isolated populations in fragmented landscapes may experience reduced genetic connectivity, high rates of random genetic drift and elevated inbreeding [1315]. These may increase genetic differentiation among remnant populations and may erode the populations’ genetic variability. The erosion in genetic variation may in turn reduce individual fitness and a population’s ability to adapt to changing environmental conditions [14, 16]. These genetic effects of habitat loss and fragmentation could be unraveled by conducting species-specific empirical investigations in appropriate systems [17].

However, empirical research has documented contradictory genetic effects of habitat loss and fragmentation especially on long-lived tree species [18]. That is to say, some studies have shown adverse effects of habitat loss and fragmentation on the genetic diversity and structure of various tree species [1922] whereas others did not find genetic effects of forest fragmentation [2325]. The absence of genetic effects of fragmentation among long-lived tree species is most likely due to the recent history of human disturbance and fragmentation events whose genetic effects will only be reflected after multiple generations which might not have been accounted for by some researchers [26]. In spite of this, the historical and contemporary effects of fragmentation may be fairly assessed and compared by working on different life history stages of trees such as fruiting adults and early stage seedlings [27]. Genetic variation among adults may reflect pre or historical fragmentation effects whereas contemporary effects will be observable among seedlings. Accordingly, fruiting adults are expected to have higher genetic variability and lower genetic differentiation than seedlings [20, 27, 28]. Application of highly informative molecular tools such as microsatellites [29] may reveal these expected patterns.

Nevertheless, genetic effects of fragmentation may still remain undetected or absent depending on a number of other contexts. First, capacity of some tree species for extensive gene flow by pollen or seed can buffer a population against genetic effects of fragmentation [23, 30]. Second, recent debates [31] argued that enhanced gene flow can be observed in more fragmented and disturbed landscapes in at least two ways: 1) if dispersal vectors are not strictly dependent on intact forest habitats, a physical opening of canopies arising from disturbance and fragmentation may allow them to move freely and cover longer distances; 2) human disturbance and fragmentation cause reduction in density of potential mates at a site. As a result, there will be less competition between local and immigrant propaguales thus likely facilitating better gene flow.

The remnant forest patches of Awi Zone in North-western Ethiopia are mainly dominated by Prunus africana, Albizia schimperiana, Albizia gumifera and Celtis africana. According to the recent natural vegetation map of Eastern Africa [32], nearly all of these plant species potentially occur or are cited as indicator or characteristic species of Afromontane rainforests in Ethiopia. In recent years, the previously large and continuous rainforests of this region have been lost because of continued human exploitation and habitat fragmentation. North-western Ethiopia harbours several forest patches of this currently rare rainforest type. The surrounding matrix consists of open agricultural land with a few remnant tree species such as Prunus africana and Albizia schimperiana. The remnant forest patches are under heavy pressure due to uncontrolled grazing and illegal logging by the local inhabitants. The remnant forest patches are thus highly threatened due to human disturbance and progressive subdivision. Conducting empirical investigations on the genetic effects of human disturbance and habitat fragmentation on representative taxa may help to formulate valid recommendations for the conservation and restoration of these fragile ecosystems.

In this study, we chose Prunus africana (Hook.f.) Kalkman (Rosaceae) as a model species to assess the impacts of habitat loss and fragmentation in Awi Zone of North-western Ethiopia. P. africana is a long-lived monoecious tree species that has been listed as Endangered on Appendix II of CITES since 1995 and as Vulnerable in the IUCN Red List of Threatened Species since 1998 [33] but is locally abundant in the study region. The bark of P. africana is highly valued in the international medicinal plant trade as it is the major source of an extract used to treat urinary disorders, prostate gland hypertrophy and benign prostatic hyperplasia [34]. Because of this, either the bark or partially processed bark extract has been heavily and unsustainably exported from Africa (especially Cameroon, Madagascar, and Kenya) to Europe causing concerns for its conservation and sustainable utilization [35]. In Kenya, a fragmentation genetics study of P. africana at the Kagamega forest revealed significantly lower genetic variability among seedlings than adults [28]. However, data are lacking on genetic variability and the extent of population differentiation among the remnant forest patches of North-western Ethiopia. Conserving and rehabilitating these otherwise vanishing forest genetic resources require more region-specific and detailed studies. So, our findings, which are based on several new analyses than used in Farwig et al. [28], are worth to inform more appropriate management and conservation actions.

We used seven polymorphic microsatellite loci to test the hypothesis that recent and ongoing forest fragmentation explains patterns of genetic variation in remnant forest patches. Under this scenario, the relative magnitude of genetic differentiation among remnant patches would be associated with their degree of fragmentation. Accordingly, we predict a) less variation in smaller relative to larger forest patches; b) more evidence of bottleneck effects in smaller than larger patches; and c) more structure in seedlings than adults. Alternatively, if recent and ongoing fragmentation does not explain patterns of genetic variation in remnant patches, we predict either panmixia across the study area, or an isolation by distance (IBD) effect associated with historical structuring of the forest that predates its fragmentation.


Linkage and Hardy-Weinberg equilibrium

Analysis of data in Micro-checker revealed the potential occurrence of null alleles in one of the 8 loci. We analyzed data at both 8 and 7 loci and found some variation in statistical outputs between the two. As a result, we decided to omit the locus suspected to contain null alleles (EMPAS01) from further analyses. No significant linkage disequilibrium was detected for the various locus pair combinations consistently across all populations indicating independence of loci used. Significant Hardy-Weinberg disequilibria were detected for some locus-population pairs for locus EMPA001 and PS12A02 even after sequential Bonferroni corrections. However, these loci did not consistently deviate from Hardy-Weinberg equilibrium for all populations. We retained these loci for analyses because results were robust regardless of their exclusion or inclusion.

Genetic variability of remnant forest patches

All the 7 loci were highly polymorphic with the number of alleles sampled per locus ranging from 7 to 29 in adults and 5 to 24 in seedlings. The total number of alleles sampled at the 7 loci over all forest patches was 129, 123 and 96 for the combined data set, adults and seedlings, respectively. The mean number of alleles sampled per patch across 7 loci ranged from 5.6 to 11 (overall mean = 8.82) in adults and 5 to 9.1 (overall mean = 7.63) in seedlings. The difference between these overall means was significantly different from zero (Table 1). The overall mean of expected heterozygosity (HS) per patch across 7 loci was significantly greater in adults (HS = 0.77) than seedlings (HS = 0.71). Similarly, the overall mean of allelic richness per patch was significantly greater in adults (5.96) than seedlings (4.85) (Table 1). However, the overall mean of FIS per patch did not differ significantly between adult and seedling populations (Table 1).

Table 1 Comparison of genetic diversity estimates (averaged across 7 loci) between cohorts of Prunus africana populations (Note: Allelic richness per locus and population was based on a minimum sample size of 7 and 6 diploid individuals, respectively in adults and seedlings

For the overall data set, average allelic richness per forest patch was significantly (t-test: t = 3.14, DF = 6, P = 0.0201) greater for large (range: 7.91 to 8.67, mean = 8.38) than small (range: 6.78 to 7.89, mean = 7.44) forest patch groups (Table 2 and Additional file 1). Comparison of average allelic richness per patch also revealed significantly greater diversity for the large than small forest patches. Similarly, mean allelic richness was significantly (t-test: t = 2.57, DF = 6, P = 0.0425) greater for less-isolated (range: 7.76 to 8.67, mean = 8.34) than isolated (range: 6.78 to 7.91, mean = 7.48) forest patches (Table 2 and Additional file 1). Further, linear regression showed a positive relationship between allelic richness and patch size and a negative relationship with patch isolation (Additional file 2). These linear regression tests were, however, non-significant.

Table 2 Variation in mean allelic richness between large and small, and less-isolated and isolated populations of Prunus africana for the three data sets (Note: Allelic richness per locus and population was based on a minimum sample size of 7 and 6 diploid individuals, respectively in adults and seedlings)

The Wilcoxon’s signed rank test for mutation-drift equilibrium under TPM did not reveal significant population demographic bottlenecks or recent changes in effective population sizes of remnant forest patches (Additional file 3). In these bottleneck tests, distribution of allele frequencies followed a normal mode shift for all populations (data not shown).

Population differentiation and genetic structure

All pairwise FST estimates were significantly higher than random for the overall data set. Furthermore, most pairwise FST comparisons showed evidence for significant population genetic differentiation in both adults and seedlings. Average population-specific FST values were significantly (paired t-test, Adult-Seedling, p = 0.0004, t-value = -6.27, DF = 7) greater for seedlings (FST = 0.075) than adults (FST = 0.035) (Figure 1), but varied little (test results not shown) with patch size (large, FST = 0.038 vs small, FST = 0.062) or patch isolation (isolated, FST = 0.063 vs less-isolated, FST = 0.038). Even so, linear regression, though non-significant, revealed a positive relationship between FST and patch isolation and an inverse relationship with patch size (Additional file 2). The trends of average population-specific FST comparisons between the seedling and adult populations were similar to corrections using the standardized FST and Jost’s Dest (Figure 1). The Mantel test for correlation between genetic and spatial distances was significant (r = 0.569, p = 0.0025) for the combined data set demonstrating an isolation by distance (IBD) effect (Figure 2). A significant pattern of IBD was also detected in adults (r = 0.68, p = 0.0003) and seedlings (r = 0.43, p = 0.0214) (Table 3). Results from simple and partial Mantel tests showed patch size and nearest distance (between edges of patches) were not significantly correlated with genetic variance after accounting for geographic distance (Table 3).

Figure 1

Comparison of Jost’s D est , F ST and F ST between adult and seedling populations of Prunus africana (Note: error bars show the standard error of the mean; L = large; C = less-isolated).

Figure 2

Correlation of genetic vs geographical distances (significance testing at 100000 permutations for Mantel test) (Note: A = Seedlings, B = Adults, C = Overall).

Table 3 Simple and partial Mantel tests (100000 permutations) to assess effects of fragmentation variables on spatial patterns of genetic differentiation in Prunus africana populations

The GESTE analysis for pairs of environmental variables against population-specific FST values showed the highest posterior probabilities for geographical distance (from centers of patches) and nearest patch distance (edge to edge) consistently in the combined data set, adults and seedlings (Table 4). The next highest posterior probability was for the pair nearest patch distance and patch size and this was observed only in the seedling cohort (Table 4).

Table 4 Posterior probabilities of models supported by the pairwise GESTE analyses of fragmentation variables against population-specific F ST

Analysis of genetic data for first generation migrant detection revealed limitation in propagule dispersal. Only up to 10% of individual seedlings were migrants. Also, first generation migrant detection for the adult populations showed that less than 8% of adult individuals were migrants (Additional file 4).


Remnant patches of Prunus africana in NW Ethiopia exhibit a pattern of higher genetic diversity in larger and less isolated patches relative to smaller and more isolated patches. This trend is expected under conservation genetics theory [36]. Remarkably, analysis of population genetic structure among patches revealed possible limitation in both pollen and seed dispersal that is likely to reflect the pre-fragmentation state of formerly vast and contiguous forests of the study region. Reduced gene flow alone may not, however, explain the strong genetic structuring [37, 38]. Our results also showed that the extent of genetic differentiation was more than twice in seedlings than adults. This finding is, however, different from that of Farwig et al. [28] who reported only a very slightly greater genetic differentiation in seedlings than adults. Also, while Farwig et al. [28] did not find significant correlations between spatial and genetic distances in adults and seedlings, our results showed significant correlations in both cohorts. The differences between findings of the two studies might be attributed to several factors. First, our sampling design was slightly different, i.e., it involved sampling adults across transect lines and up to six seedlings per adult tree, whereas Farwig et al. [28] sampled two seedlings per adult tree. Second, there might be variation in the degree of polymorphism in the different sets of marker loci used; only 2 loci are shared between the two studies. Third, regional differences with respect to various ecological attributes such as edaphic and climatic factors, intensity of logging, and vegetation condition that can determine distribution and persistence of gene dispersal agents might explain the observed differences.

In addition to the findings of Farwig et al. [28], our results contribute to understanding of population structure in P. africana by 1) providing evidence for spatial patterns of genetic structuring and potential causes after controlling for correlating factors; 2) showing no evidence for recent reduction in effective population sizes; and 3) accounting for effects of genetic diversity when evaluating genetic differentiation.

Comparison of genetic diversity between seedling and adult cohorts

Despite the relatively larger sample size of seedlings, we found a significantly lower number of alleles, expected heterozygosity and allelic richness in the seedling than adult populations. The larger homozygosity among seedlings appears to be a result of mating amongst close relatives (spatial structure) as is common in trees, but this may decay via demographic thinning or selection against homozygotes. In addition, the sampling was performed on sets of seedlings gathered around mother trees. Therefore seedlings are likely related and share some alleles, thus reducing allelic richness and heterozygosity. Also, gene flow might have changed due to the changed spatial distribution of adult trees following human disturbance and fragmentation [3941]. Thus reduced gene flow among remnant patches might have contributed to the production of seedlings with lower genetic diversity. However, the observed significant difference in genetic diversity between the two cohorts might not be ecologically relevant if selection favours heterozygote survival. A more likely explanation is that the adult cohort represents overlapping generations while the juvenile cohort represents only a few reproductive events. Thus the adult cohort could be expected to possess higher genetic variability and this variation might not be adequately represented in the seedling cohort.

Although there are no previous data on the disturbance and fragmentation history of the study region, photographic evidence taken during field observations demonstrate the occurrence of extensive logging and felling of age-old P. africana trees in the remnant forest patches (see Additional file 5). As already described by Farwig et al. [28] and Hensen et al. [42], such destruction and removal of reproductive adult trees might result in fecundity variance [43], reduced effective population size and limited gene flow among the remnant populations and as a result might have caused reduced genetic diversity in the seedling cohort. These findings are in agreement with our prediction that the younger cohort exhibits genetic effects of fragmentation more strongly than the adult cohort. Results regarding the reduced genetic diversity in the juvenile cohort are also comparable to empirical research on populations of the same species in other regions [28] and other long-lived tree species such as Polylepis incana[42], Vateriopsis seychellarum[20], Swietenia humilis[7], Quercus ilex[22] and Elaeocarpus grandis[44].

Comparison of genetic diversity between small and large remnant forest patches

Regression analyses of genetic diversity with both patch size and isolation were non-significant. Despite these, the observed patterns conformed to our expectations that reduction in patch size and increase in isolation would cause lower genetic diversity in P. africana populations. This finding is in line with a simulation study [13], which showed that genetic diversity may decline following post-fragmentation reductions in population sizes after several generations. Additionally, based on the patch size and isolation category comparisons for the overall data set, large and less-isolated patches had significantly higher genetic diversity than small and isolated patches. Similar results have been reported for other species [7]. A possible explanation for our findings, however, could be that genetic variability between patches might reflect the legacy of pre-fragmentation genetic structuring (see next section below). In adults, however, only patch isolation was significant. Unexpectedly, both patch size and isolation (based on nearest neighbor distance) did not show significant differences in the seedling cohort but the trend still conformed to our expectation thus depicting very recent signatures of fragmentation genetic effects. As described in [6], such results might also arise from insufficient statistical power for comparisons of the remnant populations due to the relatively small sample sizes.

According to recent range-wide studies on the origin and historic migration routes of P. africana[45, 46], populations of the species in Ethiopia might be more genetically diverse than populations in other regions. In the current study, however, we found mixed results when comparing different genetic diversity parameters between populations in NW Ethiopia and Kakamega forest of Kenya [28]. For instance, the total number of alleles and the mean number of alleles per population was greater for the NW Ethiopian populations than those in Kakamega forest. In contrast, the mean expected heterozygosity per population was slightly lower in the NW Ethiopian than Kakamega populations. These differences can be attributed to differences in the quality of studied forests between the two regions, the size of populations and the levels of gene flow and connectivity between populations in each region. This comparison, however, should be viewed with caution because Farwig et al. [28] used a different set of loci (only two loci are shared).

Assessment of population bottlenecks among the remnant forest patches

There was no sign of severe reductions in effective population sizes (Ne) in the past, 2 Ne – 4 Ne generations. Given the ongoing severe deforestation and fragmentation processes, we might have expected to detect signs of bottlenecks in the remnant populations, especially the more isolated ones. This finding, however, may be additional evidence for the very recent human disturbance and that fragmentation ‘may not yet have left a population genetic signature’ [6] in the study region. Our informal discussion with some of the elderly local people also showed that fragmentation in NW Ethiopia in particular Guangua Wereda of Awi Zone was a very recent phenomenon (i.e., less than 100 years). According to the elderly local people, the most devastating forest destruction and fragmentation in this region occurred especially in the past 15-20 years. This reflects the absence of law enforcement exercises and also the lack of concerns and attention for conservation of forest resources of this region. To some extent, the failure to detect recent genetic bottlenecks might reflect the use of a relatively small number of loci (n = 7) as the recommended number of loci is at least 10, although the Wilcoxon test has been shown to be sufficiently powerful with as few as 4 loci [4749]. Also, the power to detect bottlenecks is often low when the decline is not instantaneous or when not many generations have passed. A similar decline occurred in butternut, a North American tree [50], but there was little power to detect the bottleneck due to the low number of markers and shape of decline [51].

Population genetic structure and spatial patterns of differentiation

Pairwise FST values between populations revealed significant genetic differentiation suggesting restricted gene flow or limited propagule (pollen and seed) dispersal of P. africana in NW Ethiopia. The limitation in propagule dispersal was supported by our first-generation migrant detection analysis in Geneclass2 [52], which demonstrated that only up to 10% of individual seedlings in each population were migrants to the respective sampled patches. As pollen and seed dispersal in P. africana are facilitated by animal dispersal vectors, possible disturbance and fragmentation effects on mobility, abundance and survival of these agents might have caused the limitation in pollen and seed dispersal [53]. In particular, the high genetic structuring might have resulted from limited spatial patterns of seed dispersal followed by seedling establishment and recruitment filters such as early seedling mortality [37] and selective logging of adult recruits [54]. In addition to the possibility of altered gene flow patterns, other demographic characteristics (e.g., differences in effective population sizes, divergence times and colonization history) as well as evolutionary processes such as mutation and genetic drift might have caused the significant genetic differentiation between populations [37, 38].

Average population-specific FST values did not significantly vary with patch size or isolation but observed patterns matched theoretical expectations perhaps implying slightly better gene flow rates during early or pre-fragmentation stages but reduced gene flow after severe fragmentation. These results could also be attributed to the lack of statistical power due to the relatively small sample sizes of patch isolation and patch size groups or the limited number of marker loci, or may indicate little or very recent impacts of fragmentation on population differentiation of P. africana in NW Ethiopia.

The mean population-specific FST values of the seedling cohort were significantly greater than the adult. One explanation could be that the seedling cohort was more genetically isolated than the adult as a result of recent human disturbance and habitat fragmentation. However, it may also reflect the differences in genetic diversity between the two cohorts. As FST measures within versus between population diversity, an increase in within population diversity will automatically reduce the FST value. The use of standardized FST (F’ST) and newer estimates of genetic distance, Jost’s Dest[55, 56], may be expected to account for this effect. The pattern of genetic differentiation between the seedling and adult cohort, however, remained similar even after using these corrections.

We would have expected a non isolation by distance effect, with only the fragmentation variables affecting the spatial patterns of genetic differentiation among the populations. However, both simple and partial Mantel tests revealed a strong effect of distance on genetic differentiation particularly in the adult populations. This pattern almost certainly predated forest fragmentation in this region, thus indicating limitation in both pollen and seed dispersal even in the previously vast and contiguous forests of the study region. Our first-generation migrant detection analysis for the adult populations also supports this prediction. The migrant analysis revealed that less than 8% of sampled adult individuals from a patch were migrants. As the main seed dispersal agents of P. africana are thought to be highly mobile animals such as frugivorous birds and mammals [57], better gene flow or seed dispersal among remnant populations would have been anticipated, of course, assuming a negligible impact of fragmentation on dispersers. Our work on bird community and functional trait analysis in the current study region (Yineger and Hughes, in review), however, demonstrated a significant impact of fragmentation on the diversity and abundance of avian frugivores. Such impacts of fragmentation on potential dispersers could limit propagule dispersal [58] and hence might be linked to the possibility of altered gene flow patterns. Estimating gene flow by pollen and seed might help to understand the species response to fragmentation. Also, mating system might be influenced by fragmentation (e.g. increased selfing with increased isolation) [26, 59] and increased selfing might increase FST[13]. We are currently working to understand fragmentation impacts on the mating system and probability of reproduction of P. africana and are hoping to communicate this in a separate paper.


Present day forest patches of Guangua Wereda in NW Ethiopia are remnants after loss of habitat from a previously large and contiguous forest. Analysis of long-lived adult trees suggests the former forest probably exhibited strong spatial patterns of genetic structure (IBD). This means that protecting a range of patches including small and isolated ones is needed to conserve genetic diversity in this degraded landscape. This should, however, be supplemented with conserving trees growing in the surrounding farmland and pasture because matrix trees have high reproductive output and potential to act as links between populations [60]. Nevertheless, given the high livelihood dependence of the local poor community and the high impact of foreign investors on forest resources of this region, in situ conservation efforts alone may not be helpful. Therefore, these efforts should be supported with ex situ gene conservation actions. Of course, it is known that seeds of P. africana are partially recalcitrant, which means that seeds may not be stored in gene banks for longer periods. Field gene banks appear better options when considering ex situ conservation. However, this by itself may not be easy to implement. Also it may not adequately represent the extant wild genetic variation of the species. Therefore, more emphasis should be given to the in situ conservation while increasing political commitment and conservation awareness among the local people.


Study species

Prunus africana is an endangered evergreen Afromontane medicinal tree species endemic to sub-Saharan Africa and nearby islands [61]. No reliable information is available regarding the longevity of Prunus africana[62] but we guess about 200 years based on our field observation of stumps and girth of old trees in NW Ethiopia. It grows to a height of 40 m and diameter of 1 m with elliptic-oblong leathery leaves, small flowers with creamy white petals, transversely ellipsoid dark reddish brown fruits, and dark brown to grey longitudinally fissured or scaly bark [63]. The flowering season is mainly between November and February with potential intermittent flowering that is observed all year round in some localities. Pollination of P. africana is assisted by insects and fruits start to develop within 4-6 months after pollination (ICRAF, undated). Fruit ripening usually occurs between September and November but in Kakamega forest of Kenya, Farwig et al. [57] observed individual trees with ripened fruits in October, December and March. The peak fruiting season in the current study area extends from February to April. Ripened fruits are dispersed by a number of bird species and monkeys [57]. Seeds of P. africana seem to be recalcitrant. Fresh seeds can readily germinate under shady conditions [34]. The species is, however, light demanding beyond the sapling stage [64].

Study sites and sampling design

The study sites were located in the fragmented landscapes of Awi Zone, Gojjam, NW Ethiopia, at about 10°45’-11°04’ N and 36°25’-36°48’ E. This is a floristically impoverished region where all but a few remnant forest patches have been destroyed, mostly by anthropogenic means. The major rainy season extends from June to October. Mean annual rainfall is about 1500 mm and mean annual temperature ranges between 19°C and 30°C. Elevation of remnant forest patches ranges between approximately 1750 and 2390 m. Matrix surrounding remnant patches is dominated by agricultural land cultivated in various shifts for millet, corn, tef (Eragrostis tef), wheat, barley and rarely beans. The matrix also has some range lands or pastures. Both the matrix of the pasture and agricultural land contained a few remnant tree species including Prunus africana.

Maps of the forest patches were digitized from high resolution satellite images of the study region. The digitized images were processed in Arc-Map10 (ESRI Inc.) to compute patch area and distance to the nearest patch (a measure of patch isolation). Following this, four large (650 -1389 ha) and four small (8 - 62 ha) remnant study patches were selected (see Figure 3). Nearest distances (edge to edge) for the large patches were between 4 m and 88 m except for one large patch, which was 291 m. Nearest distances for small patches were between 156 m and 1218 m. The unbalanced design for nearest distance was due to lack of the study species in some of the remnant forest patches. Pairwise geographical distances between centers of patches were estimated using the ‘distance measure’ facility in Arc-Map10 (ESRI Inc.). Estimated pair wise geographical distances among sampling sites ranged between 2.9 km and 15.5 km.

Figure 3

Location of study sites in NW Ethiopia.

In each patch, P. africana trees and seedlings were sampled across a transect line of variable shape. To minimize the probability of sampling clonally propagated individuals, fruiting P. africana trees growing at least 20 m apart were selected for leaf or bark sample collection. Selected fruiting trees were coded and tagged using Aluminium tags. We used an eight meter tele-pole pruner to collect healthy leaf samples. Bark samples were collected where access to leaf samples was difficult. Two hundred and ten adult trees were processed. The circumference (at breast height, i.e. 1.3 m from ground) of each fruiting tree was measured using a measuring tape. Height was visually estimated. Where available, up to 6 seedlings (estimated height 5 cm to 30 cm) per fruiting tree were sampled near the canopy cover of selected fruiting individuals. Two hundred and thirty two seedlings were processed. Readings of elevation and geographical position were taken for each sampled fruiting tree and seedling. Leaf and bark samples were dried in locked plastic bags containing crystals of silica gel. Samples were stored in a -80°C freezer until required for DNA extraction.

DNA extraction, amplification and genotyping

About 0.25 g of dried leaf and bark samples were machine-ground using a TissueLyser II and 3 mm tungsten beads (QIAGEN Inc.). The TissueLyser grinding was operated at a frequency of 30/s in two minute interval up to 7 times or until samples were finely powdered. Genomic DNA was extracted following a modified Cetyltrimethylammonium bromide (CTAB) protocol [65]. This protocol was slightly modified by applying a longer incubation period (for at least 8 hours) after adding and thoroughly mixing 750 μl of 2X CTAB buffer and 3μl of 2-mercaptoethanol to powdered samples in 1.5 ml Eppendrof tube.

Eight nuclear microsatellite DNA markers (EMPA001, EMPAS06, EMPA010, EMPAS10, EMPA016, PS12A02, EMPAS01 and UDP96-018) originally developed for other species of the same genus but later transferred to P. africana[66] were used to amplify DNA samples via polymerase chain reaction (PCR). A single PCR reaction (10 μl volume) consisted of: double distilled water, 10× buffer, 25 μM MgCl2, 10 μM dNTP, 10 μM forward primer, 10 μM reverse primer, 10 mg/ml BSA, 10 μM fluorescent-labelled M13(-21) primers (VIC, FAM, PET OR NED), 1 U/μl Red Taq DNA polymerase (PROMEGA) and about 10 ng of template DNA. The forward primer had the universal M13 tail where the labelled M13 primer can attach. Optimized PCR conditions consisted of a cycle of initial denaturing at 94°C for 5 minutes followed by 35 cycles of denaturing at 94°C for 30 s, annealing at 55°C for 45 s and extension at 72°C for 45 s. The final extension of amplicons was set at 72°C for 10 minutes. The PCR products were subjected to electrophoresis to isolate and size microsatellite fragments against an internal LIZ size standard (GeneScan500(-250)LIZ) on a 3130 Genetic Analayzer (APPLIED BIOSYSTEMS). Allele scoring was performed in GeneMapper 4.1 (APPLIED BIOSYSTEMS). Genotyped microsatellite data were tested (after 1000 randomizations) for occurrence of null alleles, stuttering or large allele dropout using MICRO-CHECKER version 2.2.3 [67]. These tests were crucial because errors introduced by null alleles, stuttering or large allele dropout may significantly bias population genetic estimates [68].

Statistical analyses

All statistical analyses were performed separately for adults and seedlings as well as for the overall genotype data set consisting of both adults and seedlings. To assess deviation of genotypes from Hardy-Weinberg equilibrium, exact tests [69] were performed by the Markov chain method using the GENEPOP software version 4.0.10 [70] at the default analysis parameters (i.e., after 1000 dememorizations, 100 batches and 1000 iterations per batch). GENEPOP was also used to test for genotypic linkage disequilibrium for each pair of loci in each population using the log likelihood ratio statistic at the default values of the Markov chain parameters.

Within population genetic diversity

We determined unbiased genetic diversity parameters (allelic richness and gene diversity) for each population at 7 polymorphic microsatellite markers as implemented in FSTAT version [71]. Allelic richness per locus and population was based on a minimum sample size of 16, 7 and 6 diploid individuals in the overall, adult and seedling populations, respectively. FSTAT was also used to test for differences in genetic diversity between populations from two forest patch size groups i.e. large (n = 4) vs small (n = 4). The p-values in these analyses (two-tailed tests) were obtained after 15000 permutations. Average allelic richness values per population were computed across the 7 loci and tested using linear regression analysis in SAS 9.2 (SAS Institute Inc.) against nearest distance (a measure of forest patch isolation) or forest patch size. Both environmental variables were log10 transformed prior to analysis to approximate normality. To meet the basic assumptions for the linear regression analyses we always performed residual diagnostics for normality, homoscedasticity, linearity and independence problems. Fortunately, we did not encounter any problem regarding these assumptions. The level of variance in genetic diversity based on allelic richness among the two patch size (i.e., large vs small) and patch isolation (less-isolated vs isolated) groups was also evaluated using the Student’s t-test in SAS 9.2 (SAS Institute Inc.). Patch size and isolation groupings were made using median values for size and isolation. The existence of inbreeding was evaluated by comparing both the inbreeding coefficient (FIS) and corrected genetic relatedness (Relc) among the two population groups as implemented in FSTAT. Mean allelic richness, gene diversity, number of alleles and inbreeding coefficient (FIS) were compared between cohorts using paired t-test in SAS 9.2 (SAS Institute Inc.).

Population differentiation and genetic structure analyses

Population differentiation and extent of genetic structuring was evaluated by computing pairwise population genetic distances (FST) and testing for significance after 10,100 permutations as implemented in ARLEQUIN version [72]. Average population-specific FST values were compared between cohorts using paired t-test and between patch size and isolation groups using the Student’s t-test in SAS 9.2 for the overall data set. The mean population-specific FST values (see below) were also linearly regressed against patch size or isolation in SAS 9.2 based on the overall data set. The spatial pattern of genetic differentiation among remnant forest patches was assessed by running a Mantel test [73] in ARLEQUIN. In this test 100,000 permutations were performed to evaluate significance of the correlation between matrices of the genetic (pair wise FST) and geographic distances (point distance measured in Km between centres of patches) between pairs of populations. Partial Mantel tests were performed in ARLEQUIN applying 100,000 permutations to check the correlation of matrices of the genetic and two environmental variables (patch size and geographical or patch nearest distances) while keeping constant one of these environmental variables (patch size or pairwise geographical distance between centres of sampled patches or patch nearest distance). Spatial genetic patterns were also assessed by applying a hierarchical Bayesian method to compute population-specific FST values and relate them to environmental variables through a generalized linear model in GESTE version 2.0 [74]. All parameters for this analysis were set to default values. Environmental variables tested in this analysis include patch nearest distance, patch size and mean pairwise geographical distances between centers of sampled patches.

Standardized F ST (F’ST) and Jost’s Dest[55, 56] were computed using RecodeData v.0.1 [56] and GeneAlex 6.5 [75, 76] to minimize the dependence of F ST on the genetic diversity of cohorts. Computation of these estimates was necessary because uncorrected FST comparison of seedlings and adults could be compromised by differences in their genetic diversity. Tests for recent reduction in effective population size were carried out in Bottleneck version 1.2.02 [49]. These analyses were performed by applying 10,000 simulations under the two-phase mutation model (TPM) with a 95% stepwise and 5% multi-step mutations, and 12% variance for the multiple steps as suggested for microsatellite data [49]. Because the numbers of loci employed in this study were fewer than recommended for Bottleneck analyses [47], statistical significance was based on the Wilcoxon’s signed-rank test for heterozygosity excess. Lastly, analysis for detection of first generation migrants was performed in Geneclass2 [52] using the ‘likelihood computation L = L_home/L_max’. The criterion for this computation was based on a Bayesian method [77]. Ten thousand simulations [78] were run setting a threshold p-value of 0.01.





Isolation by distance


africana Prunus africana:


Cetyltrimethylammonium bromide


Polymerase chain reaction




Effective population size


Two-phase mutation model


Institute of biodiversity conservation


Department of agriculture, fisheries and forestry


Gene diversity


Inbreeding coefficient


Mean number of alleles


Convention on international trade in endangered species of wild fauna and Flora


International union for conservation of nature


World agroforestry centre.


  1. 1.

    Kissinger GM, Herold Sy VD: Drivers of Deforestation and Forest Degradation: A Synthesis Report for REDD+ Policymakers. 2012, Vancouver Canada: Lexeme Consulting

    Google Scholar 

  2. 2.

    FAO: State of the world’s forests 2012. 2012, Rome, (accessed on 08/01/2013)

    Google Scholar 

  3. 3.

    Fahrig L: Effects of habitat fragmentation on biodiversity. Annu Rev Ecol Evol S. 2003, 34: 487-515. 10.1146/annurev.ecolsys.34.011802.132419.

    Article  Google Scholar 

  4. 4.

    Laurance WF, Camargo JLC, Luizao RCC, Laurance SG, Pimm SL, Bruna EM, Stouffer PC, Williamson GB, Benitez-Malvido J, Vasconcelos HL, Van Houtan KS, Zartman CE, Boyle SA, Didham RK, Andrade A, Lovejoy TE: The fate of Amazonian forest fragments: a 32-year investigation. Biol Conserv. 2011, 144 (1): 56-67. 10.1016/j.biocon.2010.09.021.

    Article  Google Scholar 

  5. 5.

    Leonardi S, Piovani P, Scalfi M, Piotti A, Giannini R, Menozzi P: Effect of Habitat Fragmentation on the Genetic Diversity and Structure of Peripheral Populations of Beech in Central Italy. J Hered. 2012, 103 (3): 408-417. 10.1093/jhered/ess004.

    Article  PubMed  Google Scholar 

  6. 6.

    Dixo M, Metzger JP, Morgante JS, Zamudio KR: Habitat fragmentation reduces genetic diversity and connectivity among toad populations in the Brazilian Atlantic Coastal Forest. Biol Conserv. 2009, 142 (8): 1560-1569. 10.1016/j.biocon.2008.11.016.

    Article  Google Scholar 

  7. 7.

    Rosas F, Quesada M, Lobo JA, Sork VL: Effects of habitat fragmentation on pollen flow and genetic diversity of the endangered tropical tree Swietenia humilis (Meliaceae). Biol Conserv. 2011, 144 (12): 3082-3088. 10.1016/j.biocon.2011.10.003.

    Article  Google Scholar 

  8. 8.

    Ramos FN, de Lima PF, Zucchi MI, Colombo CA, Solferini VN: Genetic Structure of Tree and Shrubby Species Among Anthropogenic Edges, Natural Edges, and Interior of an Atlantic Forest Fragment. Biochem Genet. 2010, 48 (3–4): 215-228.

    CAS  Article  PubMed  Google Scholar 

  9. 9.

    Aguilar R, Ashworth L, Galetto L, Aizen MA: Plant reproductive susceptibility to habitat fragmentation: review and synthesis through a meta-analysis. Ecol Lett. 2006, 9 (8): 968-980. 10.1111/j.1461-0248.2006.00927.x.

    Article  PubMed  Google Scholar 

  10. 10.

    Honnay O, Jacquemyn H, Bossuyt B, Hermy M: Forest fragmentation effects on patch occupancy and population viability of herbaceous plant species. New Phytol. 2005, 166 (3): 723-736. 10.1111/j.1469-8137.2005.01352.x.

    Article  PubMed  Google Scholar 

  11. 11.

    Apodaca JJ, Rissler LJ, Godwin JC: Population structure and gene flow in a heavily disturbed habitat: implications for the management of the imperilled Red Hills salamander (Phaeognathus hubrichti). Conserv Genet. 2012, 13 (4): 913-923. 10.1007/s10592-012-0340-3.

    Article  Google Scholar 

  12. 12.

    Ouborg NJ, Vergeer P, Mix C: The rough edges of the conservation genetics paradigm for plants. J Ecol. 2006, 94 (6): 1233-1248. 10.1111/j.1365-2745.2006.01167.x.

    Article  Google Scholar 

  13. 13.

    Lowe AJ, Boshier D, Ward M, Bacles CFE, Navarro C: Genetic resource impacts of habitat loss and degradation; reconciling empirical evidence and predicted theory for neotropical trees. Heredity. 2005, 95 (4): 255-273. 10.1038/sj.hdy.6800725.

    CAS  Article  PubMed  Google Scholar 

  14. 14.

    Young A, Boyle T, Brown T: The population genetic consequences of habitat fragmentation for plants. Trends Ecol Evol. 1996, 11 (10): 413-418. 10.1016/0169-5347(96)10045-8.

    CAS  Article  PubMed  Google Scholar 

  15. 15.

    Angeloni F, Ouborg NJ, Leimu R: Meta-analysis on the association of population size and life history with inbreeding depression in plants. Biol Conserv. 2011, 144 (1): 35-43. 10.1016/j.biocon.2010.08.016.

    Article  Google Scholar 

  16. 16.

    Gonzalez-Varo JP, Aparicio A, Lavergne S, Arroyo J, Albaladejo RG: Contrasting heterozygosity-fitness correlations between populations of a self-compatible shrub in a fragmented landscape. Genetica. 2012, 140 (1–3): 31-38.

    Article  PubMed  Google Scholar 

  17. 17.

    Storfer A, Murphy MA, Spear SF, Holderegger R, Waits LP: Landscape genetics: where are we now?. Mol Ecol. 2010, 19 (17): 3496-3514. 10.1111/j.1365-294X.2010.04691.x.

    Article  PubMed  Google Scholar 

  18. 18.

    Kramer AT, Ison JL, Ashley MV, Howe HF: The paradox of forest fragmentation genetics. Conserv Biol. 2008, 22 (4): 878-885. 10.1111/j.1523-1739.2008.00944.x.

    Article  PubMed  Google Scholar 

  19. 19.

    Jump AS, Penuelas J: Genetic effects of chronic habitat fragmentation in a wind-pollinated tree. Proc Natl Acad Sci U S A. 2006, 103 (21): 8096-8100. 10.1073/pnas.0510127103.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  20. 20.

    Finger A, Kettle CJ, Kaiser-Bunbury CN, Valentin T, Mougal J, Ghazoul J: Forest fragmentation genetics in a formerly widespread island endemic tree: Vateriopsis seychellarum (Dipterocarpaceae). Mol Ecol. 2012, 21 (10): 2369-2382. 10.1111/j.1365-294X.2012.05543.x.

    CAS  Article  PubMed  Google Scholar 

  21. 21.

    Dubreuil M, Riba M, Gonzalez-Martinez SC, Vendramin GG, Sebastiani F, Mayol M: Genetic Effects of Chronic Habitat Fragmentation Revisited: Strong Genetic Structure in a Temperate Tree, Taxus Baccata (Taxaceae), with Great Dispersal Capability. Am J Bot. 2010, 97 (2): 303-310. 10.3732/ajb.0900148.

    Article  PubMed  Google Scholar 

  22. 22.

    Ortego J, Bonal R, Munoz A: Genetic Consequences of Habitat Fragmentation in Long-Lived Tree Species: The Case of the Mediterranean Holm Oak (Quercus ilex, L.). J Hered. 2010, 101 (6): 717-726. 10.1093/jhered/esq081.

    CAS  Article  PubMed  Google Scholar 

  23. 23.

    White GM, Boshier DH, Powell W: Increased pollen flow counteracts fragmentation in a tropical dry forest: An example from Swietenia humilis Zuccarini. Proc Natl Acad Sci U S A. 2002, 99 (4): 2038-2042. 10.1073/pnas.042649999.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  24. 24.

    Moreira PA, Fernandes GW, Collevatti RG: Fragmentation and spatial genetic structure in Tabebuia ochracea (Bignoniaceae) a seasonally dry Neotropical tree. Forest Ecol Manag. 2009, 258 (12): 2690-2695. 10.1016/j.foreco.2009.09.037.

    Article  Google Scholar 

  25. 25.

    Oberprieler C, Meister J, Schneider C, Kilian N: Genetic structure of Anogeissus dhofarica (Combretaceae) populations endemic to the monsoonal fog oases of the southern Arabian Peninsula. Biol J Linn Soc. 2009, 97 (1): 40-51. 10.1111/j.1095-8312.2008.01173.x.

    Article  Google Scholar 

  26. 26.

    Aguilar R, Quesada M, Ashworth L, Herrerias-Diego Y, Lobo J: Genetic consequences of habitat fragmentation in plant populations: susceptible signals in plant traits and methodological approaches. Mol Ecol. 2008, 17 (24): 5177-5188. 10.1111/j.1365-294X.2008.03971.x.

    Article  PubMed  Google Scholar 

  27. 27.

    Vranckx G, Jacquemyn H, Muys B, Honnay O: Meta-Analysis of Susceptibility of Woody Plants to Loss of Genetic Diversity through Habitat Fragmentation. Conserv Biol. 2012, 26 (2): 228-237. 10.1111/j.1523-1739.2011.01778.x.

    Article  PubMed  Google Scholar 

  28. 28.

    Farwig N, Braun C, Boehning-Gaese K: Human disturbance reduces genetic diversity of an endangered tropical tree, Prunus africana (Rosaceae). Conserv Genet. 2008, 9 (2): 317-326. 10.1007/s10592-007-9343-x.

    Article  Google Scholar 

  29. 29.

    Selkoe KA, Toonen RJ: Microsatellites for ecologists: a practical guide to using and evaluating microsatellite markers. Ecol Lett. 2006, 9 (5): 615-629. 10.1111/j.1461-0248.2006.00889.x.

    Article  PubMed  Google Scholar 

  30. 30.

    Sork VL, Smouse PE: Genetic analysis of landscape connectivity in tree populations. Landscape Ecol. 2006, 21 (6): 821-836. 10.1007/s10980-005-5415-9.

    Article  Google Scholar 

  31. 31.

    Bacles CFE, Jump AS: Taking a tree's perspective on forest fragmentation genetics. Trends Plant Sci. 2011, 16 (1): 13-18. 10.1016/j.tplants.2010.10.002.

    CAS  Article  PubMed  Google Scholar 

  32. 32.

    Kindt R, van Breugel P, Lillesø J-PB, Bingham M, Sebsebe D, Dudley C, Friis I, Gachathi F, Kalema J, Mbago F, Minani V, Moshi H, Mulumba J, Namaganda M, Ndangalasi H, Ruffo C, Jamnadass R, Graudal LOV: Potential natural vegetation of eastern Africa. Volume 2: Description and tree species composition for forest potential natural vegetation types. Forest & Landscape Working Papers no. 62-2011. 2011

    Google Scholar 

  33. 33.

    World_Conservation_Monitoring_Centre: Prunus africana. IUCN 2012. IUCN Red List of Threatened Species. Version 2012.2. 1998, Downloaded on 26 March 2013

    Google Scholar 

  34. 34.

    Cunningham AB, Mbenkum FT: Sustainability of harvesting Prunus africana bark in Cameroon: a medicinal plant in international trade. People and Plants Working Paper 2. 1993, Paris: UNESCO

    Google Scholar 

  35. 35.

    Cunningham AB, Ayuk E, Franzel S, Duguma B, Asanga C: An economic evaluation of medicinal tree cultivation: Prunus africana in Cameroon. People and Plants working paper 10. 2002, Paris: UNESCO

    Google Scholar 

  36. 36.

    Allendorf FW, Luikart GH, Aitken SN: Conservation and the genetics of populations. 2012, UK: Wiley-Blackwell, 2

    Google Scholar 

  37. 37.

    Hampe A, El Masri L, Petit RJ: Origin of spatial genetic structure in an expanding oak population. Mol Ecol. 2010, 19 (3): 459-471. 10.1111/j.1365-294X.2009.04492.x.

    Article  PubMed  Google Scholar 

  38. 38.

    Marko PB, Hart MW: The complex analytical landscape of gene flow inference. Trends Ecol Evol. 2011, 26 (9): 448-456. 10.1016/j.tree.2011.05.007.

    Article  PubMed  Google Scholar 

  39. 39.

    Dick CW, Etchelecu G, Austerlitz F: Pollen dispersal of tropical trees (Dinizia excelsa: Fabaceae) by native insects and African honeybees in pristine and fragmented Amazonian rainforest. Mol Ecol. 2003, 12 (3): 753-764. 10.1046/j.1365-294X.2003.01760.x.

    Article  PubMed  Google Scholar 

  40. 40.

    Garcia C, Arroyo JM, Godoy JA, Jordano P: Mating patterns, pollen dispersal, and the ecological maternal neighbourhood in a Prunus mahaleb L. population. Mol Ecol. 2005, 14 (6): 1821-1830. 10.1111/j.1365-294X.2005.02542.x.

    CAS  Article  PubMed  Google Scholar 

  41. 41.

    Ismail SA, Ghazoul J, Ravikanth G, Shaanker RU, Kushalappa CG, Kettle CJ: Does long-distance pollen dispersal preclude inbreeding in tropical trees? Fragmentation genetics of Dysoxylum malabaricum in an agro-forest landscape. Mol Ecol. 2012, 21 (22): 5484-5496. 10.1111/mec.12054.

    CAS  Article  PubMed  Google Scholar 

  42. 42.

    Hensen I, Cierjacks A, Hirsch H, Kessler M, Romoleroux K, Renison D, Wesche K: Historic and recent fragmentation coupled with altitude affect the genetic population structure of one of the world's highest tropical tree line species. Global Ecol Biogeogr. 2012, 21 (4): 455-464. 10.1111/j.1466-8238.2011.00691.x.

    Article  Google Scholar 

  43. 43.

    Klein EK, Carpentier FH, Oddou-Muratorio S: Estimating the variance of male fecundity from genotypes of progeny arrays: evaluation of the Bayesian forward approach. Methods Ecol Evol. 2011, 2 (4): 349-361. 10.1111/j.2041-210X.2010.00085.x.

    Article  Google Scholar 

  44. 44.

    Rossetto M, Jones R, Hunter J: Genetic effects of rainforest fragmentation in an early successional tree (Elaeocarpus grandis). Heredity. 2004, 93 (6): 610-618. 10.1038/sj.hdy.6800585.

    CAS  Article  PubMed  Google Scholar 

  45. 45.

    Kadu CAC, Konrad H, Schueler S, Muluvi GM, Eyog-Matig O, Muchugi A, Williams VL, Ramamonjisoa L, Kapinga C, Foahom B, Katsvanga C, Hafashimana D, Obama C, Geburek T: Divergent pattern of nuclear genetic diversity across the range of the Afromontane Prunus africana mirrors variable climate of African highlands. Ann Bot-London. 2013, 111 (1): 47-60. 10.1093/aob/mcs235.

    Article  Google Scholar 

  46. 46.

    Kadu CAC, Schueler S, Konrad H, Muluvi GMM, Eyog-Matig O, Muchugi A, Williams VL, Ramamonjisoa L, Kapinga C, Foahom B, Katsvanga C, Hafashimana D, Obama C, Geburek T: Phylogeography of the Afromontane Prunus africana reveals a former migration corridor between East and West African highlands. Mol Ecol. 2011, 20 (1): 165-178. 10.1111/j.1365-294X.2010.04931.x.

    CAS  Article  PubMed  Google Scholar 

  47. 47.

    Cornuet JM, Luikart G: Description and power analysis of two tests for detecting recent population bottlenecks from allele frequency data. Genetics. 1996, 144 (4): 2001-2014.

    CAS  PubMed  PubMed Central  Google Scholar 

  48. 48.

    Luikart G, Cornuet JM: Empirical evaluation of a test for identifying recently bottlenecked populations from allele frequency data. Conserv Biol. 1998, 12 (1): 228-237. 10.1046/j.1523-1739.1998.96388.x.

    Article  Google Scholar 

  49. 49.

    Piry S, Luikart G, Cornuet JM: BOTTLENECK: A computer program for detecting recent reductions in the effective population size using allele frequency data. J Hered. 1999, 90 (4): 502-503. 10.1093/jhered/90.4.502.

    Article  Google Scholar 

  50. 50.

    Hoban SM, Borkowski DS, Brosi SL, McCleary TS, Thompson LM, McLachlan JS, Pereira MA, Schlarbaum SE, Romero-Severson J: Range-wide distribution of genetic diversity in the North American tree Juglans cinerea: a product of range shifts, not ecological marginality or recent population decline. Mol Ecol. 2010, 19 (22): 4876-4891. 10.1111/j.1365-294X.2010.04834.x.

    Article  PubMed  Google Scholar 

  51. 51.

    Hoban S, Gaggiotti O, Con GC, Bertorelle G: Sample Planning Optimization Tool for conservation and population Genetics (SPOTG): a software for choosing the appropriate number of markers and samples. Methods Ecol Evol. 2013, 4 (3): 299-303. 10.1111/2041-210x.12025.

    Article  Google Scholar 

  52. 52.

    Piry S, Alapetite A, Cornuet JM, Paetkau D, Baudouin L, Estoup A: GENECLASS2: A software for genetic assignment and first-generation migrant detection. J Hered. 2004, 95 (6): 536-539. 10.1093/jhered/esh074.

    CAS  Article  PubMed  Google Scholar 

  53. 53.

    Farwig N, Bleher B, Bohning-Gaese K: Consequences of forest fragmentation on frugivores, seed dispersal and genetic structure of Prunus africana populations in Kenya. J Ornithol. 2006, 147 (5): 114-114.

    Google Scholar 

  54. 54.

    Lacerda AEB, Kanashiro M, Sebbenn AM: Effects of Reduced Impact Logging on genetic diversity and spatial genetic structure of a Hymenaea courbaril population in the Brazilian Amazon Forest. Forest Ecol Manag. 2008, 255 (3–4): 1034-1043.

    Article  Google Scholar 

  55. 55.

    Meirmans PG, Hedrick PW: Assessing population structure: FST and related measures. Mol Ecol Resour. 2011, 11 (1): 5-11. 10.1111/j.1755-0998.2010.02927.x.

    Article  PubMed  Google Scholar 

  56. 56.

    Meirmans PG: Using the AMOVA framework to estimate a standardized genetic differentiation measure. Evolution. 2006, 60 (11): 2399-2402. 10.1554/05-631.1.

    Article  PubMed  Google Scholar 

  57. 57.

    Farwig N, Bohning-Gaese K, Bleher B: Enhanced seed dispersal of Prunus africana in fragmented and disturbed forests?. Oecologia. 2006, 147 (2): 238-252. 10.1007/s00442-005-0288-9.

    Article  PubMed  Google Scholar 

  58. 58.

    Herrera JM, Garcia D: Effects of Forest Fragmentation on Seed Dispersal and Seedling Establishment in Ornithochorous Trees. Conserv Biol. 2010, 24 (4): 1089-1098. 10.1111/j.1523-1739.2010.01459.x.

    Article  PubMed  Google Scholar 

  59. 59.

    Coates DJ, Williams MR, Madden S: Temporal and spatial mating-system variation in fragmented populations of Banksia cuneata, a rare bird-pollinated long-lived plant. Aust J Bot. 2013, 61 (4): 235-242. 10.1071/BT12244.

    Article  Google Scholar 

  60. 60.

    Aldrich PR, Hamrick JL: Reproductive dominance of pasture trees in a fragmented tropical forest mosaic. Science. 1998, 281 (5373): 103-105.

    CAS  Article  PubMed  Google Scholar 

  61. 61.

    Kalkman C: The Old world species of Prunus sub-genus Laurocerasus. Blumea. 1965, 13: 33-35.

    Google Scholar 

  62. 62.

    Hall JB, O'Brien EM, Sinclair FL: Prunus africana: a monograph. School of Agricultural and Forest Sciences Publication Number 18. 2000, Bangor: University of Wales, 104-

    Google Scholar 

  63. 63.

    Hedberg O: Rosaceae. Flora of Ethiopia and Eritrea. Volume 3. Pittosporaceae to Araliaceae. Edited by: Hedberg I, Edwards S. 1989, Addis Ababa, Ethiopia and Uppsala, Sweden: The National Herbarium, Department of Systematic Botany, 31-44.

    Google Scholar 

  64. 64.

    Kiama D, Kiyiapi J: Shade tolerance and regeneration of some tree species of a tropical rain forest in Western Kenya. Plant Ecol. 2001, 156 (2): 183-191. 10.1023/A:1012602930347.

    Article  Google Scholar 

  65. 65.

    Cota-Sanchez JH, Remarchuk K, Ubayasena K: Ready-to-use DNA extracted with a CTAB method adapted for herbarium specimens and mucilaginous plant tissue. Plant Mol Biol Rep. 2006, 24 (2): 161-167. 10.1007/BF02914055.

    CAS  Article  Google Scholar 

  66. 66.

    Cavers S, Munro RC, Kadu CAC, Konrad H: Transfer of Microsatellite Loci For The Tropical Tree Prunus africana (Hook. f.) Kalkman. Silvae Genet. 2009, 58 (5–6): 276-279.

    Google Scholar 

  67. 67.

    Van Oosterhout C, Hutchinson WF, Wills DPM, Shipley P: MICRO-CHECKER: software for identifying and correcting genotyping errors in microsatellite data. Mol Ecol Notes. 2004, 4 (3): 535-538. 10.1111/j.1471-8286.2004.00684.x.

    CAS  Article  Google Scholar 

  68. 68.

    DeWoody JA, Nason JD, Hipkins VD: Mitigating scoring errors in microsatellite data from wild populations. Mol Ecol Notes. 2006, 6 (4): 951-957. 10.1111/j.1471-8286.2006.01449.x.

    CAS  Article  Google Scholar 

  69. 69.

    Guo SW, Thompson EA: Performing the Exact Test of Hardy-Weinberg Proportion for Multiple Alleles. Biometrics. 1992, 48 (2): 361-372. 10.2307/2532296.

    CAS  Article  PubMed  Google Scholar 

  70. 70.

    Rousset F: Genepop'007: a complete reimplementation of the Genepop software for Windows and Linux. Mol Ecol Resour. 2008, 8: 103-106. 10.1111/j.1471-8286.2007.01931.x.

    Article  PubMed  Google Scholar 

  71. 71.

    Goudet J: FSTAT (Version 1.2): A computer program to calculate F-statistics. J Hered. 1995, 86 (6): 485-486.

    Google Scholar 

  72. 72.

    Excoffier L, Lischer HEL: Arlequin suite ver 3.5: a new series of programs to perform population genetics analyses under Linux and Windows. Mol Ecol Resour. 2010, 10 (3): 564-567. 10.1111/j.1755-0998.2010.02847.x.

    Article  PubMed  Google Scholar 

  73. 73.

    Mantel N: The detection of disease clustering and a generalized regression approach. Cancer Res. 1967, 27 (2): 209-220.

    CAS  PubMed  Google Scholar 

  74. 74.

    Foll M, Gaggiotti O: Identifying the environmental factors that determine the genetic structure of Populations. Genetics. 2006, 174 (2): 875-891. 10.1534/genetics.106.059451.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  75. 75.

    Peakall R, Smouse PE: GENALEX 6: genetic analysis in Excel, Population genetic software for teaching and research. Mol Ecol Notes. 2006, 6 (1): 288-295. 10.1111/j.1471-8286.2005.01155.x.

    Article  Google Scholar 

  76. 76.

    Peakall R, Smouse PE: GenAlEx 6.5: genetic analysis in Excel. Population genetic software for teaching and research-an update. Bioinformatics. 2012, 28 (19): 2537-2539. 10.1093/bioinformatics/bts460.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  77. 77.

    Rannala B, Mountain JL: Detecting immigration by using multilocus genotypes. Proc Natl Acad Sci U S A. 1997, 94 (17): 9197-9201. 10.1073/pnas.94.17.9197.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  78. 78.

    Paetkau D, Slade R, Burden M, Estoup A: Genetic assignment methods for the direct, real-time estimation of migration rate: a simulation-based exploration of accuracy and power. Mol Ecol. 2004, 13 (1): 55-65. 10.1046/j.1365-294X.2004.02008.x.

    CAS  Article  PubMed  Google Scholar 

Download references


This study was financially supported by the Rufford Small Grants Foundation (UK) and Griffith University (Australia). We thank the Guangua Wereda Department of Agriculture and Rural Development in particular Mr. Meseret Workie (the then head of the department) for officially granting us access to the remnant forest patches and permission to collect samples. We are grateful to Mr. Worku Fetene and Mr. Solomon Beyene who acted as the chief field assistants throughout the sampling stage of this work. We also thank Mr. Abaji Tessema, Mr. Alemnew Alamirew and Mr. Belew for their help as additional field assistants in some localities. Ms Kathryn Real and staff members of the Molecular Ecology Laboratory at Griffith University trained the first author in molecular techniques and offered kind technical assistance and monitoring during the DNA extraction, amplification, allele sizing and scoring stages. The first author was supported by a higher degree research scholarship from Griffith University. Finally we thank the Institute of Biodiversity Conservation (IBC) in Ethiopia for the export permit and the Department of Agriculture, Fisheries and Forestry (DAFF) Biosecurity in Australia for the import permit. We thank two anonymous reviewers for their constructive comments on earlier version of this paper.

Author information



Corresponding author

Correspondence to Haile Yineger.

Additional information

Competing interests

Authors declare that they have no competing interests.

Authors’ contributions

HY conceptualized, designed and conducted the field and laboratory works, analyzed data, and drafted the manuscript. DS conceptualized this work, guided in data analyses and critically improved the draft manuscript for intellectual content. JH conceptualized and designed the study, guided and commented all data analyses and critically improved the draft manuscript for intellectual content. All authors read and approved the final manuscript.

Electronic supplementary material

We show in this table environmental variables and genetic diversity estimates of

Additional file 1: Prunus africana in NW Ethiopia based on the overall data set.(PDF 67 KB)

In this file we illustrate linear regression of fragmentation variables against allelic richness and population-specific

Additional file 2: F ST of Prunus africana in NW Ethiopia (Note: Area and isolation were Log10 transformed; results of all the linear regression tests were non-significant).(PDF 127 KB)

Additional file 3: Here the results of the Wilcoxon’s test for population bottlenecks under TPM and SMM are provided.(PDF 108 KB)

This table shows sample size and proportion of first generation migrants in the seedling and adult populations of

Additional file 4: P. africana in NW Ethiopia.(PDF 53 KB)

This photographic file demonstrates human impact and the occurrence of extensive logging and felling of age-old

Additional file 5: P. africana trees in the remnant forest patches of NW Ethiopia.(TIFF 2 MB)

Authors’ original submitted files for images

Below are the links to the authors’ original submitted files for images.

Authors’ original file for figure 1

Authors’ original file for figure 2

Authors’ original file for figure 3

Rights and permissions

This article is published under license to BioMed Central Ltd. 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 work is properly credited. The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated.

Reprints and Permissions

About this article

Cite this article

Yineger, H., Schmidt, D.J. & Hughes, J.M. Genetic structuring of remnant forest patches in an endangered medicinal tree in North-western Ethiopia. BMC Genet 15, 31 (2014).

Download citation


  • Genetic diversity
  • Guangua Wereda
  • Microsatellite loci
  • Population genetic structure
  • Prunus africana
  • Spatial pattern of genetic differentiation