Geological and climatic changes in quaternary shaped the evolutionary history of Calibrachoa heterophylla, an endemic South-Atlantic species of petunia
© Mäder et al.; licensee BioMed Central Ltd. 2013
Received: 27 February 2013
Accepted: 22 August 2013
Published: 29 August 2013
The glacial and interglacial cycles that characterized the Quaternary greatly affected the distribution and genetic diversity of plants. In the Neotropics, few phylogeographic studies have focused on coastal species outside of the Atlantic Rainforest. Climatic and sea level changes during the Quaternary played an important role in the evolutionary history of many organisms found in coastal regions. To contribute to a better understanding of plant evolution in this environment in Southern South America, we focused on Calibrachoa heterophylla (Solanaceae), an endemic and vulnerable wild petunia species from the South Atlantic Coastal Plain (SACP).
We assessed DNA sequences from two cpDNA intergenic spacers and analyzed them using a phylogeographic approach. The present phylogeographic study reveals the influence of complex geologic and climatic events on patterns of genetic diversification. The results indicate that C. heterophylla originated inland and subsequently colonized the SACP; the data show that the inland haplogroup is more ancient than the coastal one and that the inland was not affected by sea level changes in the Quaternary. The major diversification of C. heterophylla that occurred after 0.4 Myr was linked to sea level oscillations in the Quaternary, and any diversification that occurred before this time was obscured by marine transgressions that occurred before the coastal sand barrier’s formation. Results of the Bayesian skyline plot showed a recent population expansion detected in C. heterophylla seems to be related to an increase in temperature and humidity that occurred at the beginning of the Holocene.
The geographic clades have been formed when the coastal plain was deeply dissected by paleochannels and these correlate very well with the distributional limits of the clades. The four major sea transgressions formed a series of four sand barriers parallel to the coast that progressively increased the availability of coastal areas after the regressions and that may have promoted the geographic structuring of genetic diversity observed today. The recent population expansion for the entire species may be linked with the event of marine regression after the most recent sea transgression at ~5 kya.
The glacial and interglacial cycles that characterized the Pleistocene had major effects on the distribution and genetic diversity of plant species in the Northern Hemisphere caused by the advance and retreat of large ice sheets at higher latitudes [1, 2]. In South America, fully formed glaciers occurred only in the Andes, and the main climatic effects on most of the region were that closed vegetation types (forests) alternated with open formations (grasslands) between the glacial and interglacial periods . This might have resulted in more complex histories regarding the establishment and composition of current vegetation than simple extinctions of local biota and the formation of typical glacial refugia. In the Neotropics, most of the phylogeographic studies have focused on eastern South America [4, 5] and have investigated continental species from the Atlantic Rainforest e.g., , although a few have analyzed species from open areas such the grasslands located in the highlands and in the Pampas region of southern South America [7, 8].
After the Last Glacial Maximum (LGM) and at the beginning of the Holocene, approximately 11,700 years ago , the climate began to change. Evidence indicates a strong environmental change associated with the expansion of vegetation in southern South America due to a climate markedly warmer and moister than in the Pleistocene [10–12]. The accentuated moisture in this phase was also reported from the South Atlantic Coastal Plain (SACP) and adjacent area [13, 14]. According to these references, the pollen concentration of individual taxa increased dramatically at the beginning of the Holocene, providing evidence for rapid climatic change.
Several studies [5–8, 15, 16] have shown that demographic changes in South American plants are related to Quaternary climatic oscillations. Moreover, changes in sea level due to melting glaciers in the interglacial periods played an important role in the evolutionary history of many coastal region inhabitants . Studies of South American coastal regions are rare despite the potentially interesting evolutionary history of species in these areas, which have been strongly influenced by changes in sea level during glacial and interglacial periods . However, in the Northern Hemisphere, studies of coastal plants have been conducted more often e.g., [18, 19] and have found that changes in habitat availability and Quaternary sea level changes were responsible for the genetic structure of Asteraceae species in the Mediterranean Basin. Profound effects of Quaternary climatic cycling, tectonic shifts and changes in river dynamics, as well as sea level fluctuations, were highlighted as playing a role in shaping the genetic structure and phylogeographic patterns of plants .
After Quaternary sea level changes and the establishment of the current SACP, vegetation from the west could expand into this new environment and adapt to different conditions of salinity, soil composition, climate and hydric regime. Many plant species in the SACP originated in Andean–Pampean areas and invaded the subtropical region during the Pleistocene, especially from colder and higher elevated regions [23–25]. The central lowlands (Sedimentary depression; Figure 1) are considered an important gateway for the colonization of the SACP by Andean taxa. Small patches or long strips of gallery forests occur in the central lowlands along rivers between extensive areas of grassy vegetation . This mosaic of plant species was also influenced by the complex history and geography of the SACP region, which contains dunes, lakes, lagoons, and channels that could represent barriers or corridors to gene flow for many species that currently inhabit the area.
In the present work, we chose Calibrachoa heterophylla (Sendtn.) Wijsman (Solanaceae) as a model to study phylogeographic patterns in the SACP and thus better understand plant evolution in coastal environments in South America. Calibrachoa heterophylla is an endemic and vulnerable wild petunia species with a distribution range predominantly in the SACP, where it inhabits dunes and sandy grasslands. In addition to being found in the coastal region, this species also inhabits inland regions of the same longitude, including sandy soils that are distant (more than 500 Km) from the South Atlantic Coast. The ancestors of wild petunias have an Andean-Pampean origin, with a recent colonization of the present geographical range coupled with a fast and recent morphological radiation  and it was demonstrated  that much of the diversification of wild petunias occurred during the Pleistocene. The disjunct distribution and the recent evolutionary history of C. heterophylla make it an excellent model for understanding the Quaternary evolutionary processes of plants in the South Atlantic Coast.
Genetic markers can shed light on phylogeographical history of organisms . Chloroplast DNA (cpDNA) sequences are commonly used assess geographic structure in plants [29, 30] and they are especially informative in relatives of petunia [7, 31–33]. Plastid markers have some merits over nuclear markers for phylogeographical studies. These properties include the absence of recombination, low effective population size and the conservative mutation rate, allowing for primers to be designed over a wide taxonomic range [34, 35]. These classical cpDNA markers can answer perfectly well phylogeographic questions minimizing the effect of recent gene flow, which could mask historic biogeographic events .
The main aim of this work was to investigate how past climatic and geological changes in the SAPC could have affected the evolutionary history of C. heterophylla. Specifically we asked the following questions: (i) Did C. heterophylla colonize the SACP from inland populations growing in ancient soils that were not affected by sea level changes, or did the species arise in the SACP and then were brought or migrated inland later? (ii) Did the environmental and sea level changes during the Quaternary result in allopatric processes that isolated and shaped the genetic structure of populations? (iii) Did Holocene climate changes result in changes of population size in C. heterophylla? To answer these questions, we assessed the DNA sequences of two cpDNA intergenic spacers (trnH-psbA and trnS-trnG) that have been successfully employed in population studies of wild petunias [7, 31]. To relate the genetic variability of cpDNA markers to the geological and climatic history of the SACP during the Quaternary, several phylogeographic and phylogenetic analyses, including the dating of the formation of major clades, were performed.
Summary statistics for the datasets used
Summary statistics for each population
H14(27); H15(4); H16(1); H17(1)
South Coast (green)
0.324 ± 0.098
0.028 ± 0.031
H5(14); H7(2); H8(10)
Central Coast (blue)
0.579 ± 0.057
0.052 ± 0.047
Central Coast (blue)
H5(27); H11(1); H12(4), H13(1)
Central Coast (blue)
0.324 ± 0.098
0.033 ± 0.034
Central Coast (blue)
0.154 ± 0.126
0.013 ± 0.021
Central (blue) and North Coast (red)
0.508 ± 0.040
0.450 ± 0.249
H19(26); H20(1); H21(3)
Central (blue) and North Coast (red)
0.246 ± 0.098
0.156 ± 0.101
H2(20); H9(3); H10(1)
North Coast (red)
0.301 ± 0.112
0.030 ± 0.033
North Coast (red)
0.087 ± 0.078
0.007 ± 0.015
North Coast (red)
H4(3); H22(6); H23(1); H24(1)
0.673 ± 0.123
0.070 ± 0.060
H4(7); H25(1); H26(1)
0.417 ± 0.191
0.067 ± 0.060
0.879 ± 0.010
0.409 ± 0.220
Phylogeography and population structure
The degree of genetic structure of the populations was estimated by AMOVA, which indicated that 85.1% of the variation was among populations and 14.9% was within populations. SAMOVA analysis was used to infer the best grouping of populations on the basis of molecular variation. The configuration of four groups was considered the best by presenting the highest value of ΦCT (0.691; p < 0.001). Four groups were inferred: a group in the extreme south of the SACP (G1), a group in the central SACP (G2), a group in north of the SACP (G3) and a final group comprising inland populations (G4). Populations in each group are listed in Table 2.
A Mantel test showed a significant and moderate correlation between genetic and geographical distances between the populations (correlation coefficient = 0.39; p < 0.001). Spatial autocorrelation analysis with Alleles in Space (AIS), for the whole species, agrees with the above result, but showed a direct correlation between genetic and geographic distances up to approximately 230 km, after that the genetic divergence is maximum and shows no further increase (data not shown). Since the maximum distances between individuals within the groups are around the above value (~300 km), this suggests that, as expected, this isolation by distance pattern occurs mainly within groups. This is supported by the similar results of the spatial autocorrelation analyses for each of the two clusters with large sample sizes (G2 and G3).
Phylogenetic relationships and divergence times
The Bayesian analysis (Figure 4) presented the same four major geographically structured clades that were identified in the haplotype network analysis. As expected, given the low sequence divergence in this species, support values for most nodes were moderate to low, especially within the main clades. Haplogroups South, Central and North Coast formed well-supported clades. The estimated divergence time between C. heterophylla cpDNA haplotypes and the outgroup was 1.13 (0.53-2.26) Myr ago. The most ancient divergence within C. heterophylla was between the South Coast haplogroup and the others at approximately 1 Myr ago. The Central Coast haplogroup separated from the North Coast and Mainland haplogroups approximately 0.89 (0.37-1.68) Myr ago, and the latter two diverged approximately 0.85 (0.34-1.59) Myr ago (Figure 4). The major diversification inside the haplogroups occurred within the last 0.4 Myr. The times to the most recent common ancestors (TMRCA) for the four haplogroups were 0.23 (0.08-0.49) (North Coast), 0.38 (0.16-0.97) (Central Coast), 0.40 (0.18-1.08) (South Coast) and 0.77 (0.28-1.45) (Mainland) (see Figure 4).
Summary statistics and population growth rate for the four SAMOVA groups
0.324 ± 0.098
0.027 ± 0.031
0.607 ± 0.051
0.157 ± 0.099
0.647 ± 0.028
0.081 ± 0.061
0.690 ± 0.079
0.098 ± 0.072
0.879 ± 0.010
0.41 ± 0.22
Genetic diversity and population structure
Haplotype diversity in C. heterophylla (h = 0.879) was similar to that reported for Petunia axillaris (h = 0.831), P. altiplana (h = 0.801), P. bonjardinensis (h = 0.822) and P. guarapuavensis (h = 0.747), the species that showed the highest haplotype diversity in Petunia for the same plastidial intergenic spacers. However, the nucleotide diversity for C. heterophylla was higher than that of Petunia spp. e.g., [7, 31]. No other species of Calibrachoa has been studied so far, and thus, no comparisons within the genus can be made. The intergenic spacer trnS-trnG has proved to be an effective molecular marker for population studies in Calibrachoa, showing almost equivalent variability to trnH-psbA, in contrast to what was reported for Solanum (Solanaceae) .
Comparison of the haplotype network, phylogeny and geographical distribution of the haplotypes showed that C. heterophylla populations were genetically structured. This was confirmed by the results of AMOVA, in which most of the observed genetic diversity was found among populations. Similar results were also observed in Petunia exserta and P. axillaris and in seven highland Petunia species  using the same cpDNA molecular markers. The presence of haplotypes from different haplogroups in P6 and P7 is possibly the result of very recent contact because these populations are geographically adjacent between two distinct geographic clusters (G2 and G3). This highly structured distribution is not surprising when the biology of the species and the molecular markers used are taken into consideration. Chloroplasts are inherited maternally in Petunia and Calibrachoa species (J.R. Stehmann, Universidade Federal de Minas Gerais, personal communication) and mechanisms of long distance seed dispersal are absent in petunia species , which prevents cpDNA gene flow among populations. Studies of plants from coastal environments around the World also showed patterns of geographical structure that may be related to environmental instability mainly due to Quaternary sea level changes: South America , Europe [18, 19, 38], Asia , and North America .
Genetics structure vsgeology and climate
The only sample sites not affected by the Great Marine Transgression that occurred around 0.4 Myr ago were P12-14 (Figure 6a). The inland region (G4) may have acted as a refuge in relation to marine transgressions because the current populations in this region, despite the reduced sampling, showed the highest value of haplotype diversity among SAMOVA groups (Table 3). Moreover, the Bayesian chronogram indicated the Mainland haplogroup was the oldest (0.77 Myr), nearly twice the age of the remaining haplogroups (Figure 4). These results suggest that inland populations may have given rise to the current haplogroups found in the SACP by a colonization route from west to east, which was proposed for flora migration from Andean region to Atlantic Coastal in South America . In some regions the disappearance of geographical barriers may have favored contact among populations that were originally from different lineages, as in the case of populations P6 and P7.
The period of greatest diversification of C. heterophylla seems to have occurred in recent times, less than 0.4 Myr ago (Figure 4). A major transgression event occurred around this time where the sea level rose >20 m and most of the lowlands on the coast were submerged . This may have eliminated much of the diversity already present (Figure 6a). However, this and the three subsequent major transgressions formed a series of four sand barriers parallel to the coast that progressively increased the availability of coastal areas after the regressions and may have allowed for the genetic diversity observed today. Several studies have associated genetic diversity and historical biogeography of species distributed in eastern South America with influences of Quaternary climatic fluctuations . The same molecular markers showed that most of the genetic diversity of highlands petunia arose in a similar period .
After transgression/regression events, newly formed coastal habitats with the absence of competition provide ecological opportunities for species. These events allow the occurrence of evolutionary processes similar to those occurring on islands and/or mountaintops : rapid diversification driven by ecological opportunities. Similar results demonstrated a rapid diversification of freshwater fishes in the SACP and correlated this with intense environmental changes that were a consequence of Quaternary sea level fluctuations . Additionally, the Quaternary marine transgressions/regressions in the SACP were a major factor in Deutamys kempi (a rodent) diversification .
A similar pattern of population expansion was observed in each of our different analyses (BSP, neutrality tests, network diagram and phylogenetic tree), suggesting a general demographic pattern for the whole species. However, although all neutrality tests were negative (an indication of population expansion), most were not significant. This absence of significance may be explained by being conservative and not very powerful to detect older population size changes  compared with BSP. The star-like patterns of clades in the haplotype network and the relatively extensive geographical distribution of the majority of the derived haplotypes throughout coastal areas suggest that rapid colonization occurred in some regions. The recent expansion detected by the BSP (Figure 5) may be a sign of recent expansion linked to the Pleistocene-Holocene transition approximately 12 Kyr ago . At this time, environmental change occurred in the SACP due to marine transgressions/regressions, such as the formation of the Pleistocene Barrier III and Patos and Mirim Lagoons. The sea level was 70 m below what it is now , providing a large coastal plain that could be colonized by vegetation (Figure 6b). Moreover, with the beginning of the Holocene, the increase in temperature and humidity may have induced vegetation expansion [11, 12]. These environmental changes at the beginning of the Holocene are present in the geologic record from the SACP and adjacent areas in southern Brazil . Quaternary climatic fluctuations were likely strong enough to influence population sizes and, consequently, affect genetic diversification at shallow time scales and community dynamics at smaller scales . Global studies of plant species have also demonstrated population expansions in the same period and have correlated these results with Holocene environmental changes [53, 54].
Main insights and implications for conservation
Despite historical herbarium records, C. heterophylla was absent in Uruguay and in some Brazilian locations. This may be because of local extinctions that could have occurred naturally, but more likely, it is because of an anthropic influence, as the entire coastal area has undergone severe human-induced changes in recent decades, such as progressive urbanization, human occupation, development of cattle ranching and agriculture . To preserve the genetic diversity of C. heterophylla, specific strategies must be implemented. We propose that populations with higher genetic diversity belonging to each of the four haplogroups must be conservation priorities. Special attention should be given to the South Coast haplogroup represented in only one population (P1) and to P6 and P7 populations, the only ones with haplotypes of distinct haplogroups. The SACP deserves special attention mainly because it has been shaped by the confluence of taxa originating in tropical and temperate regions. These transitional regions are extremely important for the processes of diversification and speciation. The preservation of species that occur in these areas can be used as a model for biological responses to environmental or climatic changes .
The present phylogeographic study of C. heterophylla reveals the influence of complex geological and climatic events on patterns of diversification and distribution in this wild petunia species. The fact that the inland populations are grouped in a haplogroup found to be older than those from the coast and that the inland was not noticeably affected by Quaternary sea level changes suggest that C. heterophylla originated inland and subsequently colonized the SACP. The major diversification of C. heterophylla that occurred less than 0.4 Myr ago is likely to be related to Quaternary sea level oscillations, and earlier variability may have been lost as a consequence of marine transgressions that occurred before the formation of barriers I-IV. Furthermore, paleochannels may have been important for structuring the three haplogroups found. A recent population expansion detected by the BSP seems to be related to an increase in temperature and humidity in the SACP that occurred in the beginning of the Holocene.
The results highlight the substantial contribution of studies of species growing outside of the forests to our knowledge of past vegetation and climate dynamics in the Neotropics. However, more phylogeographic studies of SACP plant species are needed to determine if the phylogeographic history of C. heterophylla is representative of multiple taxa or if it is simply an isolated case. If the same or similar patterns were found in other endemic plants, it would suggest that a common mechanism was responsible for driving diversification in this area.
Even single locus phylogenies may not describe exactly the species’ evolutionary history, the cpDNA intergenic spacers trnG-trnS and psbA-trnH combined are successfully informative to draw the phylogeographic pattern in Calibrachoa species [32, 33] and related genus Petunia.
Sample collection, DNA extraction, amplification, and sequencing
Sampling details of Calibrachoa heterophylla populations used in the study
Santa Vitória do Palmar/RS
32° 59' 15''S/52° 43' 56''W
32° 31' 26''S/52° 32' 48''W
Mar Grosso, São José do Norte/RS
32° 02' 47''S/52° 00' 36''W
Bojuru, São José do Norte/RS
31° 40' 02''S/51° 25' 30''W
31° 06' 33''S/50° 54' 04''W
Barba Negra Island, Barra do Ribeiro/RS
30° 28' 31''S/51° 08' 46''W
Barra do Ribeiro/RS
30° 25' 13''S/51° 13' 30''W
Santo Antônio da Patrulha/RS
29° 53' 34''S/50° 25' 46''W
29° 25' 56''S/49° 47' 53''W
28° 27' 36''S/48° 45' 54''W
30° 55' 09''S/51° 29' 46''W
Cacequi River, Cacequi/RS
29° 53' 41''S/54° 51' 13''W
Santa Maria River, Cacequi/RS
29° 51' 17''S/54° 54' 31''W
São Francisco de Assis/RS
29° 34' 59''S/55° 06' 03''W
The trnH-psbA and trnS-trnG cpDNA intergenic spacers were amplified and sequenced using primers previously described [58, 59], respectively. Other plastid (psbB-psbH, trnL-trnF and trnS-trnM intergenic spacers) and nuclear (ITS and Leafy intron 2) loci were tested in representative samples, but no significant variation was found. PCR amplification and sequencing were conducted as described in other petunia species . Sequences were deposited in GenBank (accession numbers: JQ072006-JQ072022, trnH-psbA, and JQ082455-JQ082470, trnS-trnG).
For all collection localities, no specific permits were required because the land is not privately owned or protected and the field studies did not involve endangered or protected species. This work was conducted under MP 2.186-16 of the Brazilian Federal Government.
Sequence variation and population genetic structure
The sequences were aligned manually using GeneDoc . Because poly-A/T regions and small inversions were highly variable and homoplastic in preliminary phylogenetic analyses, they were not considered for further analyses, as described previously e.g., [61, 62]. The two cpDNA spacers were concatenated in all analyses. Contiguous insertion/deletion events longer than one base pair (bp) were treated as single mutation events .
Standard diversity indices including haplotype diversity (h), nucleotide diversity (π) , and analysis of molecular variance (AMOVA ) were obtained with Arlequin 22.214.171.124 . The AMOVAs were conducted using 1,000 permutations among collection sites and ΦST (pairwise differences). Mantel test between genetic and geographical distances and Spatial autocorrelation analyses (1,000 permutations) were performed with the program Alleles in Space 1.0 (AIS ). Close sites with small sample sizes were merged.
Historical barriers to gene flow were identified using Spatial Analysis of Molecular Variation in SAMOVA 1.0  with 1,000 permutations. This method identifies geographically homogeneous groups of populations that are maximally differentiated from each other, and it attempts to maximize the proportion of total genetic variation among groups of populations (ΦCT).
Phylogenetic analyses and divergence times
The evolutionary relationships between haplotypes were estimated using the Median-Joining method (ϵ = 0 ) implemented in the Network 4.6 program (available at http://www.fluxus-engineering.com). Sequences of Calibrachoa paranensis (Dusén) Wijsman, C. serrulata (L. B. Sm. & Downs) Stehmann & Semir, and C. elegans (Miers) Stehmann & Semir that belong to the same subgenus of C. heterophylla were also included in the analysis (Additional file 2).
Phylogeny and divergence times among haplotypes were estimated using a Bayesian approach with BEAST 1.6.1 . Two runs of 108 chains were conducted, sampling every 1,000 generations. The settings used were the Yule tree prior, the HKY substitution model with four gamma categories and the strict molecular clock. The substitution rate used was 2.8 × 10-9 (± 5.4 × 10-11) substitutions per site per year, or 0.56% per million years. This rate was chosen because it was previously estimated for Petunia for the same loci studied here  and it is similar to those found in other plant groups (0.22–0.58% per Myr ; 0.26–0.92% per Myr ). The best-fit model of sequence evolution was identified in Modeltest 3.7  using the Akaike Information Criterion . Tracer 1.5 (available at http://beast.bio.ed.ac.uk/Tracer) was used to check for convergence of the Markov chains and adequate effective sample sizes (>200) after the first 2 × 107 chains were deleted as “burn-in”. A maximum clade credibility tree was obtained using TreeAnnotator, part of the BEAST software package. Statistical support was determined by assessing the Bayesian posterior probabilities.
The demographic history of Calibrachoa heterophylla was investigated using several methods, such as Tajima’s D and Fu’s Fs neutrality tests calculated with Arlequin for the entire dataset and for each haplogroup (found in the network and phylogeny, see below). In addition, the Bayesian Skyline Plot (BSP ; implemented in BEAST and Tracer) was used to estimate the dynamics of changes in population size over time for the C. heterophylla dataset. The settings for the substitution model and substitution rate were the same as those used for the phylogenetic analysis.
Availability of supporting data
Collectively the group is interested in investigating evolutionary processes, plant speciation, conservation genetics, phylogeographic patterns and molecular systematics.
We are grateful to Cladinara Roberts Sarturi for help with the sequencing and João Renato Stehmann for taxonomic support. We also thank two anonymous reviewers for comments and suggestions that improved this manuscript. This project was supported by the Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq), Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES), Fundação de Amparo à Pesquisa do Estado do Rio Grande do Sul (FAPERGS), and Programa de Pós-Graduação em Genética e Biologia Molecular da Universidade Federal do Rio Grande do Sul (PPGBM-UFRGS).
- Comes HP, Kadereit JW: The effect of quaternary climatic changes on plant distribution and evolution. Trends Plant Sci. 1998, 3: 432-438. 10.1016/S1360-1385(98)01327-2.View ArticleGoogle Scholar
- Hewitt G: The genetic legacy of the quaternary ice ages. Nature. 2000, 405: 907-913. 10.1038/35016000.PubMedView ArticleGoogle Scholar
- Pennington RT, Lavin M, Oliveira-Filho A: Woody plant diversity, evolution, and ecology in the tropics: perspectives from seasonally dry tropical forests. Annu Rev Ecol Evol Syst. 2009, 40: 437-57. 10.1146/annurev.ecolsys.110308.120327.View ArticleGoogle Scholar
- Beheregaray LB: Twenty years of phylogeography: the state of the field and the challenges for the Southern Hemisphere. Mol Ecol. 2008, 17: 3754-3774.PubMedGoogle Scholar
- Turchetto-Zolet AC, Pinheiro F, Salgueiro F, Palma-Silva C: Phylogeographical patterns shed light on evolutionary process in South America. Mol Ecol. 2013, 22: 1193-1213. 10.1111/mec.12164.PubMedView ArticleGoogle Scholar
- Ribeiro RA, Lemos-Filho JP, Ramos ACS, Lovato MB: Phylogeography of the endangered rosewood Dalbergia nigra (Fabaceae): insights into the evolutionary history and conservation of the Brazilian Atlantic Forest. Heredity. 2011, 106: 46-57. 10.1038/hdy.2010.64.PubMed CentralPubMedView ArticleGoogle Scholar
- Lorenz-Lemke AP, Togni PD, Mäder G, Kriedt RA, Stehmann JR, Salzano FM, Bonatto SL, Freitas LB: Plant species diversification in eastern South American highland subtropical region: a phylogeographic perspective with native Petunia (Solanaceae). Mol Ecol. 2010, 19: 5240-5251. 10.1111/j.1365-294X.2010.04871.x.PubMedView ArticleGoogle Scholar
- Pinheiro F, Barros F, Palma-Silva C, Fay MF, Lexer C, Cozzolino S: Phylogeography and genetic differentiation along the distributional range of the orchid epidendrum fulgens: a neotropical coastal species not restricted to glacial refugia. J Biogeogr. 2011, 38: 1923-1935. 10.1111/j.1365-2699.2011.02539.x.View ArticleGoogle Scholar
- Gradstein FM, Ogg JG, Schmitz MD, Ogg GM: GM: The Geologic Time Scale 2012. 2012, Oxford: ElsevierGoogle Scholar
- Roth L, Lorscheitter ML: Palynology of a bog in Parque Nacional de Aparados da Serra, East Plateau of Rio Grande do Sul, Brazil. Quat S Am A. 1993, 8: 39-69.Google Scholar
- Behling H: South and Southeast Brazilian grasslands during late quaternary times: a synthesis. Palaeogeogr Palaeoclimatol Palaeoecol. 2002, 177: 19-27. 10.1016/S0031-0182(01)00349-2.View ArticleGoogle Scholar
- Leonhardt A, Lorscheitter ML: The last 25000 years in the eastern Plateau of Southern Brazil according to Alpes de São Francisco record. J South Am Earth Sci. 2010, 29: 454-463. 10.1016/j.jsames.2009.09.003.View ArticleGoogle Scholar
- Lorscheitter ML: Contribution to the Holocene history of Atlantic rain forest in the Rio Grande do Sul state, southern Brazil. Rev Mus Argent Cienc Nat. 2003, 5: 261-271.Google Scholar
- Leal MG, Lorscheitter ML: Plant succession in a forest on the Lower Northeast Slope of Serra Geral, Rio Grande do Sul, and Holocene palaeoenvironments, Southern Brazil. Acta Bot Bras. 2007, 21: 1-10.Google Scholar
- Ledru MP, Salatino MLF, Ceccantini G, Salatino A, Pinheiro F, Pintaud J-C: Regional assessment of the impact of climatic change on the distribution of a tropical conifer in the lowlands of South America. Divers Distrib. 2007, 13: 761-771. 10.1111/j.1472-4642.2007.00389.x.View ArticleGoogle Scholar
- Pessenda LCR, Gouveia SEM, Ribeiro AS, Oliveira PE, Aravena R: Late Pleistocene and Holocene vegetation changes in northeastern Brazil determined from carbon isotopes and charcoal records in soils. Palaeogeogr Palaeoclimatol Palaeoecol. 2010, 297: 597-608. 10.1016/j.palaeo.2010.09.008.View ArticleGoogle Scholar
- Westberg E, Kadereit JW: The influence of sea currents, past disruption of gene flow and species biology on the phylogeographical structure of coastal flowering plants. J Biogeogr. 2010, 36: 1398-1410.View ArticleGoogle Scholar
- Molins A, Mayol M, Rosselló JA: Phylogeographical structure in the coastal species Senecio rodriguezii (Asteraceae), a narrowly distributed endemic Mediterranean plant. J Biogeogr. 2009, 36: 1372-1383. 10.1111/j.1365-2699.2009.02108.x.View ArticleGoogle Scholar
- Ortiz MÁ, Tremetsberger K, Talavera S, Stuessy T, García-Castaño L: Population structure of Hypochaeris salzmanniana DC. (Asteraceae), an endemic species to the Atlantic coast on both sides of the Strait of Gibraltar, in relation to quaternary sea level changes. Mol Ecol. 2007, 16: 541-552. 10.1111/j.1365-294X.2006.03157.x.PubMedView ArticleGoogle Scholar
- Qiu YX, Fu CX, Comes HP: Plant molecular phylogeography in China and adjacent regions: Tracing the genetic imprints of quaternary climate and environmental change in the world’s most diverse temperate flora. Mol Phylogenet Evol. 2011, 59: 225-244. 10.1016/j.ympev.2011.01.012.PubMedView ArticleGoogle Scholar
- Tomazelli LJ, Villwock JA: Mapeamento geológico de planícies costeiras: o exemplo da costa do Rio Grande do Sul. Gravel. 2005, 3: 109-115.Google Scholar
- Villwock JA, Tomazelli LJ, Loss EL, Dehnhardt EA, Horn NO, Bachi FA, Dehnhardt BA: Geology of the Rio Grande do Sul Coastal Province. Quat S Am A. 1986, 4: 79-97.Google Scholar
- Safford HD: Brazilian Páramos I. An introduction to the physical environment and vegetation of the campos de altitude. J Biogeogr. 1999, 26: 693-712. 10.1046/j.1365-2699.1999.00313.x.View ArticleGoogle Scholar
- Stehmann JR, Semir J: New species of Calibrachoa and Petunia (Solanaceae) from subtropical South America. Festschrift for William G. Darcy: The Legacy of a Taxonomist. Edited by: Keating RC, Hollowell VC, Croat TB. 2005, Monographs in Systematic Botany from the Missouri Botanical Garden: Missouri Botanical Garden PressGoogle Scholar
- Safford HD: Brazilian Páramos IV. phytogeography of the campos de altitude. J Biogeogr. 2007, 34: 1701-1722. 10.1111/j.1365-2699.2007.01732.x.View ArticleGoogle Scholar
- Behling H, Pillar VD, Bauermann SG: Late Quaternary grassland (Campos), gallery forest, fire and climate dynamics, studied by pollen, charcoal and multivariate analysis of the São Francisco de Assis core in western Rio Grande do Sul (Southern Brazil). Rev Palaeobot Palynol. 2005, 133: 235-248. 10.1016/j.revpalbo.2004.10.004.View ArticleGoogle Scholar
- Kulcheski FR, Muschner VC, Lorenz-Lemke AP, Stehmann JR, Bonatto SL, Salzano FM, Freitas LB: Molecular phylogenetic analysis of Petunia Juss. (Solanaceae). Genetica. 2006, 126: 1-12.View ArticleGoogle Scholar
- Avise JC, Arnold J, Ball RM, Bermingham E, Lamb T, Neigel JE, Reeb CA, Saunders NC: Intraspecific phylogeography: the mitochondrial DNA bridge between population genetics and systematics. Ann Rev Ecol Syst. 1987, 18: 489-522.View ArticleGoogle Scholar
- Shaw J, Lickey EB, Beck JT, Farmer SB, Liu W, Miller J, Siripun KC, Winder CT, Schilling EE, Small RL: The tortoise and the hare II: relative utility of 21 noncoding chloroplast DNA sequences for phylogenetic analysis. Am J Bot. 2005, 92: 142-166. 10.3732/ajb.92.1.142.PubMedView ArticleGoogle Scholar
- Zuber D, Widmer A: Phylogeography and host race differentiation in the European mistletoe (Viscum album L.). Mol Ecol. 2009, 18: 1946-1962. 10.1111/j.1365-294X.2009.04168.x.PubMedView ArticleGoogle Scholar
- Lorenz-Lemke AP, Mäder G, Muschner VC, Stehmann JR, Bonatto Salzano FM, Freitas LB: Diversity and natural hybridization in a highly endemic species of Petunia (Solanaceae): a molecular and ecological analysis. Mol Ecol. 2006, 15: 4487-4497. 10.1111/j.1365-294X.2006.03100.x.PubMedView ArticleGoogle Scholar
- Fregonezi JN, Freitas LB, Bonatto SL, Semir J, Stehmann JR: Infrageneric classification of Calibrachoa (Solanaceae) based on morphological and molecular evidence. Taxon. 2012, 61: 120-130.Google Scholar
- Fregonezi JN, Turchetto C, Bonatto SL, Freitas LB: Biogeographical history and diversification of Petunia and Calibrachoa (Solanaceae) in the Neotropical Pampas grassland. Bot J Linn Soc. 2013, 171: 140-153. 10.1111/j.1095-8339.2012.01292.x.View ArticleGoogle Scholar
- Provan J, Soranzo N, Wilson NJ, Goldstein DB, Powell W: A low mutation rate for chloroplast microsatellites. Genetics. 1999, 153: 943-947.PubMed CentralPubMedGoogle Scholar
- Twyford AD, Kidner CA, Harrison N, Ennos RA: Population history and seed dispersal in widespread Central American Begonia species (Begoniaceae) inferred from plastome-derived microsatellite markers. Bot J Linn Soc. 2013, 171: 260-276. 10.1111/j.1095-8339.2012.01265.x.View ArticleGoogle Scholar
- Avise JC: Phylogeography: retrospect and prospect. J Biogeogr. 2009, 36: 3-15. 10.1111/j.1365-2699.2008.02032.x.View ArticleGoogle Scholar
- Derepas A, Dulieu H: Inheritance of the capacity to transfer plastids by pollen parent in Petunia hybrida. Hort J Hered. 1992, 83: 6-10.Google Scholar
- Clausing G, Vickers K, Kadereit JW: Historical biogeography in a linear system: genetic variation of Sea Rocket (Cakile maritima) and Sea Holly (Eryngium maritimum) along European coasts. Mol Ecol. 2000, 9: 1823-1833. 10.1046/j.1365-294x.2000.01083.x.PubMedView ArticleGoogle Scholar
- Franks SJ, Richards CL, Gonzales E, Cousins JE, Hamrick JL: Multi-scale genetic analysis of Uniola paniculata (Poaceae): a coastal species with a linear, fragmented distribution. Am J Bot. 2004, 91: 1345-1351. 10.3732/ajb.91.9.1345.PubMedView ArticleGoogle Scholar
- Weschenfelder J, Correa ICS, Aliotta S, Baitelli R: Paleochannels related to late Quaternary sea-level changes in Southern Brazil. Braz J Oceanogr. 2010, 58: 35-44. 10.1590/S1679-87592010000600005.View ArticleGoogle Scholar
- Pellegrino KCM, Rodrigues MT, Waite AN, Morando M, Yassuda YY, Sites JW: Phylogeography and species limits in the Gymnodactylus darwinii complex (Gekkonidae, Squamata): genetic structure coincides with river systems in the Brazilian Atlantic Forest. Biol J Linn Soc. 2005, 85: 13-26. 10.1111/j.1095-8312.2005.00472.x.View ArticleGoogle Scholar
- Cabanne GS, Santos FR, Miyaki CY: Phylogeography of Xiphorhynchus fuscus (Passeriformes, Dendrocolaptidae): vicariance and recent demographic expansion in southern Atlantic forest. Biol J Linn Soc Lond. 2007, 91: 73-84. 10.1111/j.1095-8312.2007.00775.x.View ArticleGoogle Scholar
- Thomé MTC, Zamudio KR, Giovanelli JGR, Haddad CFB, Baldissera FA, Alexandrino J: Phylogeography of endemic toads and post-Pliocene persistence of the Brazilian Atlantic Forest. Mol Phylogenet Evol. 2010, 55: 1018-1031. 10.1016/j.ympev.2010.02.003.PubMedView ArticleGoogle Scholar
- Wallace AR: On the monkeys of the Amazon. Proc Zool Soc Lond. 1852, 20: 107-110.Google Scholar
- Olson SL, Hearty PJ: A sustained +21 m sea-level highstand during MIS 11 (400 Ka): direct fossil and sedimentary evidence from Bermuda. Quat Sci Rev. 2009, 28: 271-285. 10.1016/j.quascirev.2008.11.001.View ArticleGoogle Scholar
- Werneck FP: The diversification of eastern South American open vegetation biomes: historical biogeography and perspectives. Quat Sci Rev. 2011, 30: 1630-1648. 10.1016/j.quascirev.2011.03.009.View ArticleGoogle Scholar
- Särkinen T, Pennington RT, Lavin M, Simon MF, Hughes CE: Evolutionary islands in the Andes: persistence and isolation explain high endemism in Andean dry tropical forests. J Biogeogr. 2012, 39: 884-900. 10.1111/j.1365-2699.2011.02644.x.View ArticleGoogle Scholar
- Beheregaray LB, Sunnucks P, Briscoe DA: A rapid fish radiation associated with the last sea level changes in Southern Brazil: the silverside Odontesthes perugiae complex. Proc R Soc Lond B. 2002, 269: 65-73. 10.1098/rspb.2001.1838.View ArticleGoogle Scholar
- Montes MA, Oliveira LFB, Bonatto SL, Callegari-Jacgues SM, Mattevi MS: DNA sequence analysis and the phylogeographical history of the rodent Deltamys kempi (Sigmodontinae, Cricetidae) on the Atlantic Coastal Plain of south of Brazil. J Evol Biol. 2008, 21: 1823-1835. 10.1111/j.1420-9101.2008.01586.x.PubMedView ArticleGoogle Scholar
- Fay JC, Wu CI: A human population bottleneck can account for the discordance between patterns of mitochondrial versus nuclear DNA variation. Mol Biol Evol. 1999, 16: 1003-1005. 10.1093/oxfordjournals.molbev.a026175.PubMedView ArticleGoogle Scholar
- Correa ICS: Evidence of sea-level fluctuation on the Rio Grande do Sul continental shelf, Brazil. Quat S Am A. 1986, 40: 237-249.Google Scholar
- Jansson R: Global patterns in endemism explained by past climatic change. Proc R Soc Lond B. 2003, 270: 583-590. 10.1098/rspb.2002.2283.View ArticleGoogle Scholar
- Giesecke T, Bennett KD: The Holocene spread of Picea abies (L.) Karst. in Fennoscandia and adjacent areas. J Biogeogr. 2004, 31: 1523-1548. 10.1111/j.1365-2699.2004.01095.x.View ArticleGoogle Scholar
- Zhang YH, Volis S, Sun H: Chloroplast phylogeny and phylogeography of Stellera chamaejasme on the Qinghai-Tibet Plateau and in adjacent regions. Mol Phylogenet Evol. 2010, 57: 1162-1172. 10.1016/j.ympev.2010.08.033.PubMedView ArticleGoogle Scholar
- Lopes CM, Freitas TRO: Human impact in naturally patched small populations: genetic structure and conservation of the burrowing rodent, tuco-tuco (Ctenomys lami). J Heredity. 2012, 103: 672-681. 10.1093/jhered/ess027.View ArticleGoogle Scholar
- Smith TB, Kark S, Schneider CJ, Wayne RK: Biodiversity hotspots and beyond: the need for preserving environmental transitions. Trends Ecol Evol. 2001, 16: 431-View ArticleGoogle Scholar
- Roy A, Frascaria N, MacKay J, Bousquet J: Segregating random amplified polymorphic DNAs (RAPDs) in Betula alleghaniensis. Theor Appl Genet. 1992, 85: 173-180.PubMedGoogle Scholar
- Sang T, Crawford DJ, Stuessy TF: Chloroplast DNA phylogeny, reticulate evolution, and biogeography of Paeonia (Paeoniaceae). Am J Bot. 1997, 84: 1120-1136. 10.2307/2446155.PubMedView ArticleGoogle Scholar
- Hamilton MB: Four primers pairs for the amplification of chloroplast intergenic regions with intraspecific variation. Mol Ecol. 1999, 8: 513-525.View ArticleGoogle Scholar
- Nicholas KB, Nicholas HBJ: GeneDoc: a tool for editing and annotating multiple sequences alignments. http://www.psc.edu/biomed/genedoc,
- Kelchner SA: The evolution of non-coding chloroplast DNA and its application in plant systematics. Ann Mo Bot Gard. 2000, 87: 482-498. 10.2307/2666142.View ArticleGoogle Scholar
- Kim K-J, Lee H-L: Widespread occurrence of small inversions in the chloroplast genomes of land plants. Mol Cells. 2005, 19: 104-113.PubMedGoogle Scholar
- Simmons M, Ochoterena H: Gaps as characters in sequence-based phylogenetic analyses. Syst Biol. 2000, 49: 369-381.PubMedView ArticleGoogle Scholar
- Nei M: Molecular evolutionary genetics. 1987, New York: Columbia University PressGoogle Scholar
- Excoffier L, Smouse PE, Quattro M: Analysis of molecular variance inferred from metric distances among DNA haplotypes: application to human mitochondrial DNA restriction data. Genetics. 1992, 131: 479-491.PubMed CentralPubMedGoogle Scholar
- 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: 564-567. 10.1111/j.1755-0998.2010.02847.x.PubMedView ArticleGoogle Scholar
- Miller M: Alleles in space: computer software for the joint analysis of interindividual spatial and genetic information. J Hered. 2005, 96: 722-724. 10.1093/jhered/esi119.PubMedView ArticleGoogle Scholar
- Dupanloup I, Schneider S, Excoffier L: A simulated annealing approach to define the genetic structure of populations. Mol Ecol. 2002, 11: 2571-2581. 10.1046/j.1365-294X.2002.01650.x.PubMedView ArticleGoogle Scholar
- Bandelt HJ, Forster P, Röhl A: Median-joining networks for inferring intraspecific phylogenies. Mol Biol Evol. 2005, 16: 37-48.View ArticleGoogle Scholar
- Drummond AJ, Suchard MA: Bayesian random local clocks or one rate to rule them all. BMC Biol. 2010, 8: 114-10.1186/1741-7007-8-114.PubMed CentralPubMedView ArticleGoogle Scholar
- Wolfe KH, Li WH, Sharp PM: Rates of nucleotide substitution vary greatly among plant mitochondrial, chloroplast, and nuclear DNAs. Proc Natl Acad Sci USA. 1987, 84: 9054-9058. 10.1073/pnas.84.24.9054.PubMed CentralPubMedView ArticleGoogle Scholar
- Alsos IG, Engelskjon T, Gielly L, Taberlet P, Brochmann C: Impact of ice ages on circumpolar molecular diversity insights from an ecological key species. Mol Ecol. 2005, 14: 2739-2753. 10.1111/j.1365-294X.2005.02621.x.PubMedView ArticleGoogle Scholar
- Posada D, Crandall KA: Modeltest: testing the model of DNA substitution. Bioinformatics. 1998, 14: 817-818. 10.1093/bioinformatics/14.9.817.PubMedView ArticleGoogle Scholar
- Kelchner SA, Thomas MA: Model use in phylogenetics: nine key questions. Trends Ecol Evol. 2007, 22: 87-94. 10.1016/j.tree.2006.10.004.PubMedView ArticleGoogle Scholar
- Tajima F: Statistical method for testing the neutral mutation hypothesis by DNA polymorphism. Genetics. 1989, 123: 585-595.PubMed CentralPubMedGoogle Scholar
- Fu YX: Statistical tests of neutrality of mutations against population growth, hitchhiking and background selection. Genetics. 1997, 147: 915-925.PubMed CentralPubMedGoogle Scholar
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 (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.