1,379
Views
2
CrossRef citations to date
0
Altmetric
Research Article

Population genetic analysis of two Polylepis microphylla (Wedd.) Bitter (Rosaceae) forests in Ecuador

ORCID Icon, & ORCID Icon
Pages 184-197 | Received 19 Oct 2020, Accepted 13 Apr 2021, Published online: 07 May 2021

ABSTRACT

Habitat fragmentation is a severe problem that negatively affects forests around the world. This phenomenon has consequences at a genetic level in populations such as genetic diversity and gene flow reduction. As a result, population fitness decreases, and the possibility of extinction increases. Polylepis (Rosaceae) forests are endemic to the Andes and their distribution reflects hundreds of years of human activities. Polylepis microphylla is the most endangered species of the genus in Ecuador and is represented by only two known forests in Chimborazo province, separated from each other by 1,6 km. The present study focused on i) exploring the genetic structure and connectivity, ii) analyzing the genetic diversity patterns in these two forests of P. microphylla. Parameters such as Fixation index, molecular variance (AMOVA) and genetic diversity estimators were calculated from the haplotypic sequences of one chloroplastic DNA region (trnG intron), one nuclear DNA region (ITS), and the concatenation of both. We found that DNA regions used were highly polymorphic with a total of nine, four and thirteen segregation sites, respectively, and nine, six and fifteen haplotypes, respectively. Also, we detected that both forest work as a single population unit due to low fixation index values (Fst = 0,03822, −0,03883 and −0,00469, respectively), which indicate high gene flow between these two forests. In addition, we detected high genetic diversity reflected by haplotype diversity in both forest (hd = 0,54204, 0,63664 and 0,81682, respectively). Likewise, haplotype relationships were visualized in networks, which demonstrated limited spatial structuring between haplotypes, corroborating the low genetic differentiation among these forests. In conclusion, our results suggest there is a single population unit of P. microphylla in Ecuador. Habitat fragmentation does not seem to have negative effects on the genetic diversity of this species. Nevertheless, the species may be vulnerable to rapid environmental fluctuations despite having relatively high genetic diversity. Its long generational time causes increased susceptibility of P. microphylla to environmental changes and may eventually lead to the extinction of the species in Ecuador.

Introduction

Habitat fragmentation caused by human activities is a major problem that negatively affects forests around the world, leading to isolated and small forest patches [Citation1–3]. On a large scale, habitat fragmentation can be evidenced by decline of biodiversity, changes in community composition and alteration of ecological processes [Citation3]. As a consequence of spatial isolation due to habitat fragmentation, genetic structure and connectivity can be altered, leading to genetic diversity and gene flow reduction, inbreeding depression, and genetic erosion [Citation4–7]. Thus, population fitness might decrease, and populations may become more vulnerable to strong environmental changes, increasing their probability of extinction [Citation4,Citation5].

Genetic diversity, structure and connectivity are closely related with occupied geographical space and can be affected by habitat fragmentation [Citation7,Citation8]. According to Young et al. [Citation9], genetic consequences of habitat fragmentation such as reduced variability, increased inbreeding and reduced gene flow have serious implications for species conservation. The reduction and isolation of forests can diminish the genetic pool of populations and reduce their fitness and capacity to quickly and effectively adapt to environmental changes [Citation4,Citation5,Citation10–12].

The tropical Andes are considered as one of the most important biodiversity hotspots in the world [Citation13,Citation14], which offers a wide variety of ecosystem services to the region [Citation15 p. 1–7, 13–19]. Due to high species richness and endemism, it has been categorized as one of the priority regions for conservation around the world [Citation13]. However, human activities such as agriculture, which has been carried out in the region for hundreds of years [Citation16], and fires [Citation17], have fragmented Andean ecosystems [Citation18–20].

Trees and shrubs of Polylepis Ruiz & Pav. (Rosaceae), endemic to the Andean region, exemplify habitat fragmentation related to anthropogenic impact. Polylepis forests have a patchy distribution that is the result of hundreds of years of constant human activities [Citation21]. Logging, livestock, agricultural activities and, especially, fire are determining factors that cause an accelerated reduction in the size of these forests and are a possible explanation of their currently fragmented distribution [Citation22–25].

Fjeldså [Citation26] considered Polylepis forests as small oases of life compared with the characteristic and dominant páramo grasslands in which they occur. They are also considered as biodiversity islands, providing shelter to many endemic animal and plants species [Citation27–29]. In addition, they are valuable habitats that offer several ecosystem services [Citation27]. Despite these characteristics, Polylepis forests are considered as one of the most endangered ecosystems around the world [Citation30].

The genus Polylepis comprises about 45 species [Citation31–34] that are distributed from northern Venezuela to northern Argentina and Chile [Citation35,Citation36] and can be found between 2700 and 5200 m a.s.l [Citation35–37]. In Ecuador, there are ten native species ranging from 2700 to 4300 m a.s.l [Citation36,Citation38]. All species are wind pollinated.

Polylepis microphylla (Wedd.) Bitter is the most endangered species of the genus in Ecuador, being known in the country from only one location [Citation39, p. 748–750]. In agreement with the forest definition of the United Nations Framework Conservation on Climate Change (UNFCCC) [Citation40], the species is represented by two forest remnants near Achupallas village in Chimborazo province, in an open grassland area dominated by agricultural activities between 3500 and 4000 m a.s.l. In 2004 IUCN Red List of Threatened Species, P. microphylla is recorded as vulnerable (VU) [Citation41], but the 2011 Red Book of the Endemic Plants of Ecuador [Citation42] categorizes it as critically endangered (CR) [39, p.748–750]. Based on our field observations, we consider that the biggest threats for this species are accelerated expansion of rural settlements, logging and, especially, fires ().

Figure 1. Photograph of burnt Polylepis microphylla shrubs in the study area

Figure 1. Photograph of burnt Polylepis microphylla shrubs in the study area

Available information about P. microphylla is limited and is mainly taxonomic. Simpson [Citation35] considered the species a synonym of P. weberbaueri Pilg. but Bitter [Citation43], Romoleroux [Citation36] and Kessler & Schmidt-Lebuhn [Citation31] recognized P. microphylla as a separate species. Schmidt‐Lebuhn et al. [Citation44] estimated the genome size of P. microphylla by flow cytometry and presumed that it is a diploid species.

It is generally assumed that long-distance pollen dispersion is operative in wind-pollinated species [Citation11,Citation45], implying high levels of genetic diversity and low genetic differentiation among populations [Citation11,Citation46]. However, Hensen et al. [Citation20] found high genetic differentiation among populations in wind pollinated Polylepis incana Kunth, indicating that habitat fragmentation has a negative effect on the genetic diversity in this species. Given the limited distance (1,6 km) between the two forest stands of P. microphylla, we hypothesize that these two stands share a common genetic pool and work as a single population. Thus, the present study will focus on: i) exploring the genetic structure in the two known forests of P. microphylla in Ecuador, and ii) analyzing the genetic diversity patterns of P. microphylla.

Materials and methods

Study species and study area. Since the first review of the genus [Citation43], various species concepts have been applied to Polylepis. Here we applied the general lineage concept by which species are considered a segment of an evolutionary lineage [Citation33,Citation47], because this definition is applied at the population level and takes ecological aspects into account.

Polylepis microphylla has been recorded from one location in Ecuador and three in Peru, ranging from 3200 to 4000 m a.s.l [Citation48,Citation49]. Following the taxonomical description of Romoleroux [Citation36], the species is a shrub or tree to 7 m high (); leaves are alternate, compound and imparipinnate formed by 3–6 pairs of leaflets, always accompanied by a pair of villous to glabrescent stipules that are connate and adnate to the petiole; leaflets are oblong or slightly ovate, 0,3–0,7 × 0,2–0,6 cm, the lower surface is densely villous with sparse small red and black hairs, while the upper surface is sparsely villous; inflorescences are simple racemes; flowers are apetalous with oval and green sepals; stamens 15–20, with elongate filaments; hypanthium with irregular crests densely covered by lanose pubescence; and fruits have small, flattened spines.

Figure 2. Tree habit of Polylepis microphylla of approximately 4 m high located in forest 1

Figure 2. Tree habit of Polylepis microphylla of approximately 4 m high located in forest 1

Figure 3. Detail of mature flower and leaves of Polylepis microphylla

Figure 3. Detail of mature flower and leaves of Polylepis microphylla

The study was carried out in the surroundings of Totoras indigenous community along the road Ozogoche-Achupallas, about 18 km SE of Alausí (Chimborazo province). The two forest of P. microphylla occur here on relatively dry slopes, in open areas dominated by grassland and crops, separated by a valley of ≈1,6 km () ( near here). According to UNFCCC [Citation40], forests are defined as having a minimum area of 0,05–1,0 hectares with vegetation coverage of more than 10–30%, made up of trees with a capacity to reach 2–5 m height at maturity. Forest 1 was located at S 02° 16ʹ 16,79ʹ’ W 78° 42ʹ 16,88ʹ’, between 3556 and 3800 m a.s.l., with a total area of ≈19,06 ha and an area of P. microphylla coverage of ≈10,56 ha (). Forest 2 was located at S 02° 15ʹ 04,53ʹ’ W 78° 41ʹ 30,77ʹ’, between 3964 and 3673 m a.s.l., and had a total area of ≈49,05 ha and an area of P. microphylla coverage of ≈30,14 ha ().

Figure 4. Map of two Polylepis microphylla forests in the study area in Chimborazo province, Ecuador

Figure 4. Map of two Polylepis microphylla forests in the study area in Chimborazo province, Ecuador

Figure 5. Photograph of forest 1. The open grassland area is dominated by agricultural activities where Polylepis microphylla grows

Figure 5. Photograph of forest 1. The open grassland area is dominated by agricultural activities where Polylepis microphylla grows

Figure 6. Photographs of forest 2. A) Polylepis microphylla represented by shrubs of maximum 1,5 m high. Area dominated by grassland. B) A slope of P. microphylla coverage

Figure 6. Photographs of forest 2. A) Polylepis microphylla represented by shrubs of maximum 1,5 m high. Area dominated by grassland. B) A slope of P. microphylla coverage

Sampling. In each forest, we randomly sampled 50 adult individuals of P. microphylla at minimally 20 m distance from each other to limit genetic closeness [Citation20]. Leaf samples were collected from each individual and stored in paper envelopes with silica gel. A GPS waypoint was noted for each individual, and vouchers were deposited in Herbarium QCA.

DNA extraction, PCR amplification and sequencing. Total DNA extraction was carried out with 50 mg of grinded leaf tissue using Doyle & Doyle [Citation50] protocol with some modifications. β-mercaptoethanol was placed independently of the CTAB solution, chloroform and isoamyl alcohol were added separately too. Also, DNA was precipitated for 24 hours. Finally, the extracted DNA was quantified in NanoDrop® 1000 Spectrophotometer and stored at concentrations of 25 ng/µL in 0.1 X TE solution on the freezer at −20 °C.

After testing a total of 3 chloroplast DNA (cpDNA) regions (trnL-trnF, rpoB-trnC, and trnS-trnG), we decided to select the intron trnG, which is part of the trnS-trnG intergenic spacers for the present study. We selected this cpDNA region based on the informational potential reflected as number of segregation sites. In the case of nuclear DNA, in order to make the best use of resources, we selected the ITS region, which has proven to be highly informative at the intraspecific level [Citation51–56]. The ITS region amplified included the spacer I and II and the 5.8 subunit of nuclear ribosomal DNA (nrDNA) [Citation51].

The final reaction volume was 25 µL consisting of 12,5 µL of GoTaq® Green Master Mix (Promega©), 4 µL of DNA (100 ng), 1 µL of each primer (10 pmol) and 6,5 µL of UltraPure™ DNase/RNase-Free Distilled Water (Thermo Fisher Scientific©).

Polymerase chain reaction (PCR) protocol for both regions was carried out under the following conditions of time and temperature: 95 °C x 3 min; 35 cycles (94 °C x 1 min; 54 °C x 1 min; 72 °C x 1 min); 72 °C x 7 min. The primers used of each region are presented in . ( near here)

Table 1. DNA regions used. The type of region, primers used, their sources and temperatures of annealing (TA) are indicated

The PCR products were verified by electrophoresis in 1% agarose gel in a solution of Tris-borate-EDTA (TBE) 1X with 4 µL of SYBR safe 10 000X in DMCO (Thermo Fisher Scientific©). Finally, we ran it at 100 v in the electrophoresis chamber for one hour. PCR positive samples were sequenced in Macrogen Korea commercial laboratory.

Data analysis.- In each region analyzed, the forward and reverse sequences received from Macrogen Korea were reviewed and aligned using the Geneious Alignment function [Citation57] in the Geneious Prime® 2019.2.1 program. Each consensus sequence obtained was reviewed and edited manually. For the alignment of consensus sequences between individuals of the same population we used the MUSCLE Alignment function [Citation58,Citation59] of the same program.

The regions were analyzed independently and concatenated. To concatenate them we used the SuperMatrix function from evobiR package in RStudio v.1.2.1335 program [Citation60].

The estimators of genetic diversity used at the intrapopulation and interpopulation levels were: Number of segregation sites (S), number of haplotypes (h), haplotypic diversity, (hd) and nuceotide diversity (π). Haplotypes are the genetic constitution of a single chromosome [Citation61], while haplotype diversity is defined as the probability that two randomly sampled alleles are different. On the other hand nucleotide diversity represents the average number of nucleotide differences per site in the comparison between different sequences of DNA [Citation62]. These genetic estimators were obtained using the DNA Sequence Polymorphism v.5 program.

To determine the percentage of variation between and within the forests, the analysis of the molecular variance (AMOVA) was performed. Excoffier et al. [Citation63] proposed AMOVA as a method to test the hierarchical population structure. They suggested that under the null hypothesis, random samples taken from the complete pool belong to a global population. This means that the analyzed populations are structured as one. While, to determine the genetic differentiation between populations, the value corresponding to fixation index (Fst) [Citation64] was used. This parameter is widely used to determine population structure type, and is defined as the probability of sampling randomly alleles from within a subpopulation in relation to the whole population [Citation65]. These values were calculated using the Arlequin v.3 program.

To visualize the genealogical relationships of DNA sequences at intraspecific level and population structure, haplotypic networks were generated by means of the Median Joining method [Citation66–68] in Population Analysis with Reticulate Trees (PopART) v.5. program. Finally, node names of haplotype networks were added using Inkscape v 0.92.4 program.

Results

DNA extraction, amplification and sequencing.- DNA extraction was successfully completed in 19 samples of forest 1 (F1) and 18 samples from forest 2 (F2), which represent the 37% of all the individuals sampled. The amplification of the trnG intron and ITS region was positive for all samples.

Reviewing, editing and concatenation of regions.- Alignments of the trnG intron and ITS region resulted in sequences of 737 and 618 bp, respectively. The trnG intron was the most polymorphic region with 9 segregation sites (S), while the ITS region showed 4 (). The concatenated region resulted in 37 sequences, 19 from F1 and 18 from F2 with a final length of 1355 bp and a total of 13 segregating sites (). AMOVA.- trnG intron. The analysis of the molecular variance showed that highest percent variation occurred within forests (96,18%), compared to among forests (3,82%) (P > 0,05). The Fst value revealed low population differentiation (Fst = 0,03822) ().

Table 2. Segregation sites and haplotype frequencies of trnG intron with their position in the final alignment and their frequencies in the two forests. (-) represents the same nucleotide in relation to the first sequence

Table 3. Segregation sites and haplotype frequencies of ITS region with their position in the final alignment and their frequencies in the two forests. (-) represents the same nucleotide in relation to the first sequence

Table 4. Segregation sites and haplotype frequencies of the concatenated regions with their position in the final alignment and their frequencies in the two forests. (-) represents the same nucleotide in relation to the first sequence

Table 5. Analysis of molecular variance (AMOVA) in trnG intron, ITS and concatenated regions in the two forests (P > 0,05). (d.f) degrees of freedom, (Va) Variance of components among forests, (Vb) Variance of components within forests, (Fst) Weir & Cockerham (1984) fixation index

ITS. AMOVA showed the same pattern for ITS: percent variation within forests (103,88%) was higher than among forests (−3,88%) (P > 0,05). The Fst value revealed low population differentiation (Fst = −0,03883) ().

Concatenated regions. Main percent of variation was located within forests (100,47%) compared to among forests (−0,47%) (P > 0,05). The Fst value revealed a low population differentiation (Fst = −0,00469) ().

Intrapopulation genetic diversity estimators.- trnG intron. Within F1, 6 segregation sites were identified, generating 6 haplotypes: H3, H5, H6, H7, H8 and H9. Haplotypes H3, H6, H7 and H8 were unique for this forest, while H9 was the most frequent. Haplotypic diversity was relatively high (>0,5) (hd = 0,53,801), nucleotide diversity was low (<0,5) (π = 0,00099) (). In the case of F2, five segregating sites were identified, generating 5 haplotypes: H1, H2, H4, H5 and H9. Haplotypes H1, H2 and H4 were unique for this forest, and H9 was the most frequent. Haplotypic diversity was relatively high (hd = 0,55556) while nucleotide diversity was low (π = 0,00139) ().

Table 6. Genetic diversity estimators, distribution and frequency of haplotypes within and both forests in trnG intron, ITS and concatenated regions. (n) sample size; (S) number of segregation sites; (h): number of haplotypes; (hd) haplotype diversity; (π) nucleotide diversity, * represents unique haplotypes

ITS. Within F1, 3 segregation sites were identified, generating a total of 5 haplotypes: H1, H2, H3, H5 and H6. Haplotypes H1, H3 and H5 turned out to be unique for this forest and H6 proved to be the most frequent. Haplotypic diversity was high (hd = 0,69591) while the nucleotide diversity was low (π = 0,00194) (). In the case of F2, 3 segregation sites were identified, generating a total of 3 haplotypes: H2, H4 and H6. Haplotype H4 was unique and H2 was the most frequent. Likewise, haplotypic diversity was high (hd = 0,58170) while the nucleotide diversity was low (π = 0.00191) ().

Concatenated regions. 9 segregation sites were identified in F1, generating a total of 9 haplotypes: H4, H7, H9, H10, H11, H12, H13, H14 and H15. Haplotypes H4, H7, H9, H10, H11 and H14 were unique for this forest, while H13 was the most frequent. Haplotypic diversity turned out to be high (hd = 0,84795) while the nucleotide diversity was low (π = 0,00142) (). In the case of F2, we found 8 segregation sites, and 8 haplotypes: H1, H2, H3, H5, H6, H8, H12 and H13. Haplotypes H1, H2, H3, H5, H6 and H8 were unique for this forest, while H12 and H13 were the most frequent ones. Haplotypic diversity was high (hd = 0,80392) and nucleotide diversity low (π = 0,00163) ().

Interpopulation genetic diversity estimators.- trnG intron. Nine segregation sites were detected in both forests, generating nine haplotypes (). Haplotypic diversity value was relatively high (hd = 0,54204) while nucleotide diversity was low (π = 0.00121) ().

ITS. A total of 4 segregation sites were detected in both forests, generating 6 haplotypes (). Haplotypic diversity value was high (hd = 0,63664) while nucleotide diversity was low (π = 0,00189) ().

Concatenated regions. Thirteen segregation sites were detected in the two forests, generating 15 haplotypes (). Haplotypic diversity value was high (hd = 0,80392) while nucleotide diversity was low (π = 0,00152) ().

Haplotype networks.- trnG intron. The network corroborates the low differentiation between forests. As shown in , the most common haplotype H9 is shared by the two forests with a frequency of 25. Other haplotypes are separated from the most common by one or two nucleotide substitutions. Haplotype 5, the second most common one, is also shared by both forests ().

Figure 7. Relationships of the 9 haplotypes generated from trnG intron represented by a median joining network. Each circle represents a different haplotype and the crossed lines represent nucleotide differences

Figure 7. Relationships of the 9 haplotypes generated from trnG intron represented by a median joining network. Each circle represents a different haplotype and the crossed lines represent nucleotide differences

ITS. Haplotypes H2 and H6, being the most frequent ones with frequencies of 15 and 17 m respectively, are present in both forests (). Unique haplotypes have low frequencies: H1(1), H3(1) and H5(2) in forest 1 and H4(1) in forest 2. These haplotypes are different from the most frequent for only one nucleotide substitution, showing limited spatial structuring between haplotypes ().

Figure 8. Relationships of the 6 haplotypes generated from ITS region represented by a median joining network. Each circle represents a different haplotype and the crossed lines represent nucleotide differences

Figure 8. Relationships of the 6 haplotypes generated from ITS region represented by a median joining network. Each circle represents a different haplotype and the crossed lines represent nucleotide differences

Concatenated regions. The patterns described above are repeated here, the most frequent haplotypes, H12 and H13, being shared by both forests with almost the same frequencies, respectively 11 and 12. All other haplotypes have frequencies of 1 with exception of H10 which has a frequency of 2 ().

Figure 9. Relationships of the 15 haplotypes generated from the concatenated regions represented by a median joining network. Each circle represents a different haplotype and the crossed lines represent nucleotide differences. Median vector could be interpreted as an unsampled or extinct haplotype

Figure 9. Relationships of the 15 haplotypes generated from the concatenated regions represented by a median joining network. Each circle represents a different haplotype and the crossed lines represent nucleotide differences. Median vector could be interpreted as an unsampled or extinct haplotype

Discussion

The results obtained in this study represent the first analysis of genetic diversity and structure of P. microphylla forests in Ecuador. They should be handled with caution because of the possibility of different ploidy levels within the Polylepis species [Citation44] and the relatively low number of samples used in the study.

According to Brown and Biggs [Citation69], theoretically 60–70% of alleles are captured by sampling 10% of the individuals in a population. Therefore, the number of samples to be analyzed in population genetics studies is not defined and may depend on several factors, such as the number of alleles captured by the type of molecular markers used [Citation70]. Zhao et al. [Citation71] found that 90% of population diversity as reflected by the expected heterozygosity or genetic diversity, can be captured with 25–30 samples using AFLP’s markers and 10–15 samples using SSRs. On the other hand, Hale et al. [Citation72] mentioned that minimally 25–30 individuals are needed using microsatellite markers.

The AMOVA results showed that the greatest source of variation is explained within forests, compared to among forests (). This points to a single population structure of P. microphylla in Ecuador. Due to the significance of the P value exceeding the tolerable range to take the alternative hypothesis [Citation73], the null hypothesis of the AMOVA, which implies that all samples belong to a global population, was selected [Citation63]. Our results are not in line with those reported by Hensen et al. [Citation20], who found that habitat fragmentation has a negative effect on the connectivity and indicate genetic differentiation among populations in P. incana using AFLP markers. Our results are also not in complete agreement with those of Vranckx et al. [Citation74], who indicated that severe habitat fragmentation can affect to wind-pollinated tree species and may lead to pollen limitation and limited gene flow.

As the studied P. microphylla forests were at a distance from each other of approximately 1,6 km, we suggest that the spatial separation was not a barrier for completely random reproductive processes. Since Polylepis trees are wind pollinated, their pollen grains can be dispersed over long distances [Citation75] and have reproductive success between spatially separated patches. For instance, Seltmann (unpubl.) found P. australis pollen grains separated approximately by 80 km from the nearest forest. Additionally, she mentioned that these pollen grains had a viability rate of up to 5 days after being released.

Values corresponding to Fst index can be interpreted as low genetic differentiation between the two forests, indicating high levels of gene flow and random mating [Citation76] (). Negative values in the ITS region and the concatenated analysis can occur when populations are not or are weakly differentiated [Citation77]. The results obtained in this study agree with studies on population genetics of Polylepis rugulosa Bitter and P. tarapacana Philippi [Citation78], P. australis Bitter [Citation12,Citation79], P. pauta Hieron [Citation80] and P. besseri Hieron [Citation81], all of which revealed high levels of gene flow and low population differentiation.

Because of the low genetic differentiation between the two forests of P. microphylla, we suggest that the two forests act as a single population unit and have a panmictic reproduction system, implying random mating between individuals [Citation82].

High diversity, reflected as haplotypic diversity, was detected in P. microphylla (>0.5) in each of the three genetic regions analyzed. This is consistent with results obtained in P. australis in Argentina [Citation79] and P. pauta in Ecuador [Citation80], where high genetic diversity was detected based on analysis of different molecular markers.

In both the individual and the concatenated analysis of the analyzed genetic regions, a well-defined pattern of genetic diversity estimators was found. The high haplotype diversity (>0.5) and low nucleotide diversity (<0,5) may be a consequence of relatively short sequences and few haplotypes in the population [Citation83] (). This pattern can indicate a low genetic differentiation between the observed haplotypes, separated from each other by few nucleotide substitutions () [Citation84,Citation85].

According to Hamrick [Citation11], populations at geographical and ecological border situations are generally more susceptible to environmental changes. Due to the major, human-induced environmental changes [Citation86], P. microphylla could be severely affected and even become locally extinct, despite having a high genetic variability. Its slow life cycle might not allow P. microphylla to adapt quickly to the changing habitat conditions [Citation87]. Because of its highly restricted geographical distribution in Ecuador, together with our results that points to a single population in the country, P. microphylla is vulnerable to environmental changes and susceptible to extinction [Citation11,Citation88].

Acknowledgments

This study would not have been possible without the support of the Pontificia Universidad Católica del Ecuador (PUCE), through the PUCE research project named “Dinámica y análisis poblacional de los bosques de Polylepis en el Ecuador”. Similarly, we extend our thanks to Ministerio de Ambiente y Agua de Ecuador and the project “Caracterización genética y análisis bioinformático de la biodiversidad de los Parques Nacionales Podocarpus, Sangay, Yacuri y Llanganates” for support. Thanks to the research projects “Diversidad de Plantas Vasculares”, 012-2017-IC-FLO-DNB/MA and “Diversidad Genética de la Flora Ecuatoriana”, MAE-DNB-CM-2018-0082 carried out in conjunction with PUCE, it was possible to collect the botanical and leaf tissue samples. Finally, we thank to Dr R. Gradstein and the two reviewers for their valuable comments.

Disclosure statement

No potential conflict of interest was reported by the author(s).

Additional information

Funding

This work was supported by Pontificia Universidad Católica del Ecuador through the research project: “Dinámica y análisis poblacional de los bosques de Polylepis en el Ecuador” Grant O13106/QINV0103;Pontificia Universidad Católica del Ecuador (EC) [O13106/QINV0103].

References

  • Skole D, Tucker C. Tropical deforestation and habitat fragmentation in the Amazon: satellite data from 1978 to 1988. Science. 1993;260(5116):1905–1910.
  • Sanderson EW, Jaiteh M, Levy MA, et al. The human footprint and the last of the wild: the human footprint is a global map of human influence on the land surface, which suggests that human beings are stewards of nature, whether we like it or not. BioScience. 2002;52(10):891–904.
  • Haddad N, Brudving L, Clobert J, et al. Habitat fragmentation and its lasting impact on earth’s ecosystems. Sci Adv. 2015;1(2):e1500052.
  • Amos W, Balmford A. When does conservation genetics matter? Heredity (Edinb). 2015;87(3):257–265. 2001.
  • Bijlsma R, Loeschcke V. Genetic erosion impedes adaptive responses to stressful environments. Evol Appl. 2012;5(2):117–129.
  • Wang R, Compton SG, Chen XY. Fragmentation can increase spatial genetic structure without decreasing pollen‐mediated gene flow in a wind-pollinated tree. Mol Ecol. 2011;20(21):4421–4432.
  • Teich I Análisis de la estructura genética espacial de especies arbóreas y su asociación con la variabilidad fenotípica y ambiental. Buenos Aires. Facultad de Ciencias Exactas y Naturales. Universidad de Buenos Aires; 2012.
  • Caicedo AL, Schaal BA. Population structure and phylogeography of Solanum pimpinellifolium inferred from a nuclear gene. Mol Ecol. 2004;13(7):1871–1882.
  • Young A, Boyle T, Brown T. The population genetic consequences of habitat fragmentation for plants. Trends Ecol Evol. 1996;11(10):413–418.
  • DeSalle R, Amato G. The expansion of conservation genetics. Nat Rev Genet. 2004;5(9):702.
  • Hamrick JL. Response of forest trees to global environmental changes. For Ecol Manage. 2004;197(1–3):323–335.
  • Seltmann P, Hensen I, Renison D, et al. Biparental inbreeding depression, genetic relatedness and progeny vigour in a wind-pollinated treeline species in Argentina. Plant Ecol. 2009b;205(1):155–164.
  • Myers N, Mittermeier RA, Mittermeier CG, et al. Biodiversity hotspots for conservation priorities. Nature. 2000;403(6772):853.
  • Hutter CR, Guayasamín JM, Wiens JJ. Explaining Andean megadiversity: the evolutionary and ecological causes of glassfrog elevational richness patterns. Ecol Lett. 2013;16(9):1135–1144.
  • Anderson EP, Marengo J, Villalva R, et al. Consequences of climate change for ecosystems and ecosystem services in the tropical Andes. In: Herzog SK, Martínez R, Jørgensen, editors. Climate change and biodiversity in the tropical Andes. São José dos Campos-Brazil: Inter-American Institute for Global Change Research (IAI) and Scientific Committee on Problems of the Environment (SCOPE);2011.1–14
  • Åkesson CM, Matthews-Bird F, Bitting M, Åkesson CM, Matthews-Bird F, Bitting M, et al. 2,100 years of human adaptation to climate change in the High Andes. Nat Ecol Evol. 2020;4(1):66–74.
  • Lippok D, Beck SG, Renison D, et al. Forest recovery of areas deforested by fire increases with elevation in the tropical Andes. For Ecol Manage. 2013;295:69–76.
  • Arango-Vélez N, Kattan GH. Effects of forest fragmentation on experimental nest predation in Andean cloud forest. Biol Conserv. 1997;81(1–2):137–143.
  • Armenteras D, Gast F, Villareal H. Andean forest fragmentation and the representativeness of protected natural areas in the eastern Andes, Colombia. Biol Conserv. 2003;113(2):245–256.
  • Hensen I, Cierjacks A, Hirsch H, et al. 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.
  • Kessler M. The “Polylepis problem”: where do we stand? Ecotropica. 2002;8(2):97–110.
  • Kessler M Present and potential distribution of Polylepis (Rosaceae) forests in Bolivia. Biodiversity and Conservation of Neotropical Montane Forests. 1995;281–294.
  • Renison D, Cingolani AM, Suarez R. Efectos del fuego sobre un bosquecillo de Polylepis australis (Rosaceae) en las montañas de Córdoba, Argentina. Rev Chil Hist Nat. 2002;75(4):719–727.
  • Teich I, Cingolani AM, Renison D, et al. Do domestic herbivores retard Polylepis australis Bitt. woodland recovery in the mountains of Córdoba Argentina? For Ecol Manage. 2005;219(2):229–241.
  • Cierjacks A, Salgado S, Wesche K, et al. Post‐fire population dynamics of two tree species in high-altitude Polylepis forests of central Ecuador. Biotropica. 2008;40(2):176–182.
  • Fjeldså J. Key areas for conserving the avifauna of Polylepis forests. Ecotropica. 2002;8(2):125–131.
  • Gareca EE, Hermy M, Fjeldså J, et al. Polylepis woodland remnants as biodiversity islands in the Bolivian high Andes. Biodivers Conservat. 2010;19(12):3327–3346.
  • Enrico L, Funes G, Cabido M. Regeneration of Polylepis australis Bitt. in the mountains of central Argentina. For Ecol Manage. 2004;190(2–3):301–309.
  • Romoleroux K, Cárate-Tandalla D, Erler R, et al. Plantas vasculares de los bosques de Polylepis en los páramos de Oyacachi. 1st. Quito-Ecuador: Centro de Publicaciones, Pontificia Universidad Católica del Ecuador; 2016.
  • Kessler M, Bosques de Polylepis. In: Moraes M, Øllgaard B, Kvist LP, editors. Botánica Económica de los Andes Centrales. La Paz-Bolivia: Universidad Mayor de San Andrés; 2006. p. 110–120.
  • Kessler M, Schmidt-Lebuhn AN. Taxonomical and distributional notes on Polylepis (Rosaceae). Organ Divers Evol. 2006;6(1):67–70.
  • Seltmann P, Cocucci A, Renison D, et al. Mating system, outcrossing distance effects and pollen availability in the wind-pollinated treeline species Polylepis australis Bitt. (Rosaceae) Basic Appl Ecol. 2009a;10(1):52–60.
  • Boza Espinoza TE, Quispe-Melgar HR, Kessler M. Taxonomic reevaluation of the Polylepis sericea complex (Rosaceae), with the description of a new species. Systemat Bot. 2019;44(2):324–334.
  • Boza Espinoza TE, Kessler M A monograph of the genus Polylepis (Rosaceae). Phytokeys. in prep.
  • Simpson BBA. Revision of the genus Polylepis (Rosaceae: sanguisorbeae). Smithson Contrib Bot. 1979;43:1–62.
  • Romoleroux K. Rosaceae 79. In: Harling G, Andersson L, editors. Flora of Ecuador. Göteborg: University of Gothenburg; 1996. p. 71–89.
  • Argollo J, Soliz C, Villalba R. Potencialidad dendrocronológica de Polylepis tarapacana en los Andes Centrales de Bolivia. Ecología en Bolivia. 2004;39(1):5–24.
  • Espinoza TEB, Romoleroux K, Kessler M. Taxonomic revaluation of the Polylepis pauta and P. sericea (Rosaceae) from Ecuador. Phytotaxa. 2020;454(2):111–126.
  • Romoleroux K. Rosaceae. In: León-Yánez S, Valencia R, Pitman N, editors. Libro Rojo de las Plantas Endémicas del Ecuador. Quito-Ecuador: Publicaciones del Herbario QCA, Pontificia Universidad Católica del Ecuador; 2011. p. 748–750.
  • United Nations Framework Convention on Climate Change (UNFCCC). Report of the conference of the parties on its seventh session, held at Marrakesh from 29 October to 10 November 2001, addendum part two: action taken by the conference of the parties, Volume I; 2001.
  • Polylepis microphylla. The IUCN red list of threatened species 2004 [internet]. Romoleroux K, Pitman N. [ cited 2021 Feb 2]. Available from: https://doi.org/10.2305/IUCN.UK.2004.RLTS.T36206A9987050
  • León-Yánez S, Valencia R, Pitman N, et al. Libro Rojo de las Plantas Endémicas del Ecuador. 2nd ed. Quito-Ecuador: Publicaciones del Herbario QCA, Pontificia Universidad Católica del Ecuador; 2011.
  • Bitter G. Revision der Gattung Polylepis. Botanische Jahrbücher für Systematik. 1911;45:564–656.
  • Schmidt‐Lebuhn AN, Fuchs J, Hertel D, et al. An Andean radiation: polyploidy in the tree genus Polylepis (Rosaceae, Sanguisorbeae). Plant Biol. 2010;12(6):917–926.
  • Kramer AT, Ison JL, Ashley MV, et al. The paradox of forest fragmentation genetics. Conserv Biol. 2008;22(4):878–885.
  • Llorens TM, Tapper SL, Coates DJ, et al. Does population distribution matter? Influence of a patchy versus continuous distribution on genetic patterns in a wind‐pollinated shrub. J Biogeograph. 2017;44(2):361–374.
  • De Queiroz K. The general lineage concept of species, species criteria, and the process of speciation: a conceptual unification and terminological recommendations. In: Howard DJ, Berlocher SH, editors. Endless forms: species and speciation. New York-USA: Oxford University Press; 1998. p. 57–75.
  • Jørgensen PM, León-Yánez S, editors. Catálogo de las Plantas Vasculares del Ecuador. St. Louis-USA: Missouri Botanical Garden Press; 1999. p. 633–668.
  • Mendoza W, Cano A. Diversidad del género Polylepis (Rosaceae, Sanguisorbeae) en los Andes peruanos. Revista Peruana de Biología. 2011;18(2):197–200.
  • Doyle J, Doyle J. Preservation of plant samples for DNA restriction endonuclease analysis. Taxon. 1987;36(4):715–722.
  • Harris DJ, Crandall KA. Intragenomic variation within ITS1 and ITS2 of freshwater crayfishes (Decapoda: cambaridae): implications for phylogenetic and microsatellite studies. Mol Biol Evol. 2000;17(2):284–291.
  • Gao Q, Zhang D, Duan Y, et al. Intraspecific divergences of Rhodiola alsia (Crassulaceae) based on plastid DNA and internal transcribed spacer fragments. Bot J Linnean Soc. 2012;168(2):204–215.
  • González-Tortuero E, Rusek J, Maayan I, et al. Genetic diversity of two Daphnia-infecting microsporidian parasites, based on sequence variation in the internal transcribed spacer region. Parasit Vectors. 2016;9(1):293.
  • Hapsari L, Azrianingsih R, Arumingtyas EL. Genetic variability and relationship of banana cultivars (Musa L.) from East Java, Indonesia based on the internal transcribed spacer region nrDNA sequences. J Trop Biol Conserv. 2018;15:101–120.
  • White TJ, Bruns T, Lee S, et al. Amplification and direct sequencing of fungal ribosomal RNA genes for phylogenetics. In: Innis M, Gaelfand D, Sninsky J, et al., editors. PCR protocols: a guide to methods and applications. San Diego-USA: Academic Press; 1990. 315–322.
  • Yang JY, Pak JH. Phylogeny of Korean Rubus (Rosaceae) based on ITS (nrDNA) and trnL/F intergenic region (cpDNA). J Plant Biol. 2006;49(1):44–54.
  • Smith TF, Waterman MS. Identification of common molecular subsequences. J Mol Biol. 1981;147(1):195–197.
  • Edgar RC. MUSCLE: multiple sequence alignment with high accuracy and high throughput. Nucleic Acids Res. 2004a;32(5):1792–1797.
  • Edgar RC. MUSCLE: a multiple sequence alignment method with reduced time and space complexity. BMC Bioinformatics. 2004b;5(1):113.
  • evobiR: comparative and population genetic analyses [Internet]. Arlington(TX): R-project; [cited 2020 Jul 15]. Available from: https://CRAN.R-project.org/package=evobiR. 2015.
  • Marsjan PA, Oldenbroek JK Molecular markers, a tool for exploring genetic diversity. In: food and agriculture organization (FAO). The State of the World’s Animal Genetic Resources for Food and Agriculture, first draft. Rome-Italy: FAO; 2006. p. 374.
  • Nei M. Molecular evolutionary genetics. New York (NY): Columbia University Press; 1987.
  • Excoffier L, Smouse PE, Quattro JM. Analysis of molecular variance inferred from metric distances among DNA haplotypes: application to human mitochondrial DNA restriction data. Genetics. 1992;131(2):479–491.
  • Weir BS, Cockerham CC. Estimating F-statistics for the analysis of population structure. Evolution. 1984;38(6):1358–1370.
  • Holsinger KE, Weir BS. Genetics in geographically structured populations: definition, estimation and interpretation of F ST. Nat Rev Genet. 2009;10(9):639.
  • Bandelt HJ, Forster P, Röhl A. Median-joining networks for inferring intraspecific phylogenies. Mol Biol Evol. 1999;16(1):37–48.
  • Mardulyn P. Trees and/or networks to display intraspecific DNA sequence variation? Mol Ecol. 2012;21(14):3385–3390.
  • Paradis E, Pearse W. Analysis of haplotype networks: the randomized minimum spanning tree method. Meth Ecol Evolut. 2018;9(5):1308–1317.
  • Brown AH, Briggs JD. Sampling strategies for genetic. In: Falk AD, Holsinger KE, editors. Genetics and conservation of rare plants. New York: Oxford University Press; 1991. p. 99–119.
  • Ward SM, Jasieniuk M. Review: sampling weedy and invasive plant populations for genetic diversity analysis. Weed Sci. 2009;57(6):593–602.
  • Zhao R, Cheng Z, Lu W, et al. Estimating genetic diversity and sampling strategy for a wild soybean (Glycine soja) population based on different molecular markers. Chinese Sci Bull. 2006;51(10):1219–1227.
  • Hale ML, Burg TM, Steeves TE. Sampling for microsatellite-based population genetic studies: 25 to 30 individuals per population is enough to accurately estimate allele frequencies. PloS One. 2012;7(9):e45170.
  • Sánchez OJ Introducción al diseño experimental. Quito-Ecuador: Giro Creativo; 2017
  • Vranckx GU, Jacquemyn H, Muys B, et al. Meta‐analysis of susceptibility of woody plants to loss of genetic diversity through habitat fragmentation. Conserv Biol. 2012;26(2):228–237.
  • Schmidt-Lebuhn AN, Kessler M, Kumar M. Promiscuity in the Andes: species relationships in Polylepis (Rosaceae, Sanguisorbeae) based on AFLP and morphology. Systemat Bot. 2006;31(3):547–559.
  • Wright S. The genetical structure of populations. Ann Eugenics. 1951;15(1):323–354.
  • Willing EM, Dreyer C, Van Oosterhout C. Estimates of genetic differentiation measured by FST do not necessarily require large sample sizes when using many SNP markers. PLoS One. 2012;7(8):e42649.
  • Schmidt-Lebuhn AN, Kumar M, Kessler M. An assessment of the genetic population structure of two species of Polylepis Ruiz & Pav. (Rosaceae) in the Chilean Andes. Flora-Morphol Distrib Funct Ecol Plants. 2006;201(4):317–325.
  • Julio N, Sobral A, Dueñas JR, et al. RAPD and ISSR markers indicate diminished gene flow due to recent fragmentation of Polylepis australis woodlands in central Argentina. Biochem Syst Ecol. 2008;36(5–6):329–335.
  • Aragundi S, Hamrick JL, Parker K. Genetic insights into the historical distribution of Polylepis pauta (Rosaceae) in the northeastern Cordillera Oriental of Ecuador. Conserv Gene. 2011;12(3):607–618.
  • Gareca EE, Breyne P, Vandepitte K, et al. Genetic diversity of Andean Polylepis (Rosaceae) woodlands and inferences regarding their fragmentation history. Bot J Linnean Soc. 2013;172(4):544–554.
  • Curtis H, Schnek A. Curtis Biología. 7th edition in Spanish. Buenos Aires-Argentina: Editorial Médica Panamericana; 2008. Chapter 18, Principios básicos de la genética de poblaciones; p. 355.
  • Ye Y, Li J, Wu C. Genetic diversity and population connectivity of the Asian green mussel Perna viridis in South China Sea, inferred from mitochondria DNA markers. Biochem Syst Ecol. 2015;61:470–476.
  • De Jong MA, Wahlberg N, Van Eijk M, et al. Mitochondrial DNA signature for range-wide populations of Bicyclus anynana suggests a rapid expansion from recent refugia. PLoS One. 2011;6(6):e21385.
  • Mao M, Liu HL. Genetic diversity of Trichomonas vaginalis clinical isolates from Henan province in central China. Pathog Glob Health. 2015;109(5):242–246.
  • Bonan GB, Doney SC. Climate, ecosystems, and planetary futures: the challenge to predict life in Earth system models. Science. 2018;359(6375):eaam8328.
  • Davis MB, Shaw RG. Range shifts and adaptive responses to quaternary climate change. Science. 2001;292(5517):673–679.
  • McLaughlin JF, Hellmann JJ, Boggs CL, et al. Climate change hastens population extinctions. Proceedings of the National Academy of Sciences. 2002;99( 9):6070–6607