BACKGROUND
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  <cit> .
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  alternated with open formations  between the glacial and interglacial periods  <cit> .
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  <cit>  and have investigated continental species from the atlantic rainforest e.g.,  <cit> , 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  <cit> .
after the last glacial maximum  and at the beginning of the holocene, approximately  <dig>  years ago  <cit> , 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  <cit> .
the accentuated moisture in this phase was also reported from the south atlantic coastal plain  and adjacent area  <cit> .
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  <cit>  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  <cit> .
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  <cit> .
however, in the northern hemisphere, studies of coastal plants have been conducted more often e.g.,  <cit>  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  <cit> .
the south atlantic coastal plain is the largest coastal plain in south america, covering approximately  <dig>  km <dig>  from 28° north to 34° south.
it is a flat area of lowlands occupied mostly by large systems of coastal lakes.
it extends ne-sw for approximately  <dig> km and has an average width of 60– <dig> km in the south central portion and 15– <dig> km in the northern portion .
the main factors that formed the sacp as it is today were oscillatory glaciation cycles during the pleistocene, which promoted sea level transgressions and regressions  <cit> .
during the late quaternary, four sand barriers parallel to the coast were formed and identified in the sacp; three of these barriers originated during the pleistocene, and one originated during the holocene.
they were associated with the limits of the marine transgressions of the four major interglacial periods in the last  <dig> myr  that occurred in the marine isotope stages   <dig>   <dig>   <dig>  and  <dig>  which are referred to as barrier-lagoon systems i to iv, respectively .
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  <cit> .
the central lowlands  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  <cit> .
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  wijsman  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  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  <cit>  and it was demonstrated  <cit>  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  <cit> .
chloroplast dna  sequences are commonly used assess geographic structure in plants  <cit>  and they are especially informative in relatives of petunia  <cit> .
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  <cit> .
these classical cpdna markers can answer perfectly well phylogeographic questions minimizing the effect of recent gene flow, which could mask historic biogeographic events  <cit> .
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:  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?
did the environmental and sea level changes during the quaternary result in allopatric processes that isolated and shaped the genetic structure of populations?
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  that have been successfully employed in population studies of wild petunias  <cit> .
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.
RESULTS
genetic diversity
the mean number of individuals genotyped per population was  <dig> .
the representative lengths of the alignment for trnh-psba and trns-trng were 456 bp and 784 bp, respectively.
overall,  <dig> sites were informative,  <dig> in each region, resulting in a value of  <dig>  ±  <dig>  for nucleotide diversity  and a value of  <dig>  ±  <dig>  for haplotype diversity .
moreover,  <dig> different haplotypes were identified from concatenated spacers .
the marker trnh-psba had diversity indices slightly higher than trns-trng, though the amount of variability is similar .
however, intra-population diversity varied widely throughout the populations sampled.
five populations were monomorphic and most others exhibited relatively low nucleotide diversity.
p <dig> and p <dig> populations had the highest values of nucleotide diversity.
the summary statistics obtained for the markers and populations are shown in tables  <dig> and  <dig>  respectively.
π, nucleotide diversity; h, haplotype diversity; sd, standard deviation;
* duplication of 5 bp; tv, transversion; ts, transition.
n, number of samplings; h, haplotype diversity; π, nucleotide diversity; sd, standard deviation.
phylogeography and population structure
the degree of genetic structure of the populations was estimated by amova, which indicated that  <dig> % of the variation was among populations and  <dig> % 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 .
four groups were inferred: a group in the extreme south of the sacp , a group in the central sacp , a group in north of the sacp  and a final group comprising inland populations .
populations in each group are listed in table  <dig>
a mantel test showed a significant and moderate correlation between genetic and geographical distances between the populations .
spatial autocorrelation analysis with alleles in space , 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 .
since the maximum distances between individuals within the groups are around the above value , 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 .
the haplotype network  resulted in four geographically structured haplogroups that were concordant with the samova results and are depicted in figure 3a.
the four haplotypes of the southernmost population  formed one haplogroup .
the haplotypes from the other coastal populations formed the haplogroups central coast  and north coast , respectively .
the remaining haplotypes belonging to the populations p <dig>  p <dig>  p <dig>  and p <dig> were grouped into a mainland haplogroup , considering that the last three were collected outside of the sacp.
the central haplotype h <dig>  whose grouping may be considered ambiguous in the network, was grouped with its respective clades according to the bayesian phylogeny .
the p <dig> and p <dig> populations were the only ones with haplotypes belonging to different haplogroups in the network.
the h <dig> haplotype that belongs to the north coast haplogroup  was also sampled in p <dig> and the h <dig> haplotype that belongs to the central coast haplogroup  was sampled in p <dig> .
a preeminent feature of the network was that all four clades have a star-like pattern, in which infrequent haplotypes are derived from a common central haplotype with only one or two substitutions.
it is also noteworthy that the four haplotypes shared among populations  were internal to the network haplotypes.
phylogenetic relationships and divergence times
the bayesian analysis  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  <dig>   myr ago.
the most ancient divergence within c. heterophylla was between the south coast haplogroup and the others at approximately  <dig> myr ago.
the central coast haplogroup separated from the north coast and mainland haplogroups approximately  <dig>   myr ago, and the latter two diverged approximately  <dig>   myr ago .
the major diversification inside the haplogroups occurred within the last  <dig>  myr.
the times to the most recent common ancestors  for the four haplogroups were  <dig>   ,  <dig>   ,  <dig>    and  <dig>    .
intraspecific demography
the neutrality tests for each group of populations and for the whole species were mostly non-significant .
however, tajima’s d and fu’s fs were negative for all results but significant only for samova groups g <dig> and g <dig> fs .
the bayesian skyline plot for the species  suggested that c. heterophylla underwent an increase in effective population size fairly recently, approximately  <dig> kyr ago.
the very small number of haplotypes, which was most likely due to the sample size, precludes the use of this methodology for each clade separately.
π, nucleotide diversity; h, haplotype diversity; sd, standard deviation; n, sampling size;.
§, significant values .
discussion
genetic diversity and population structure
haplotype diversity in c. heterophylla  was similar to that reported for petunia axillaris , p. altiplana , p. bonjardinensis  and p. guarapuavensis , 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.,  <cit> .
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   <cit> .
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 <cit>  and in seven highland petunia species  <cit>  using the same cpdna molecular markers.
the presence of haplotypes from different haplogroups in p <dig> and p <dig> is possibly the result of very recent contact because these populations are geographically adjacent between two distinct geographic clusters .
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 <cit>  and calibrachoa species  and mechanisms of long distance seed dispersal are absent in petunia species  <cit> , 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  <cit> , europe  <cit> , asia  <cit> , and north america  <cit> .
genetics structure vs geology and climate
the bayesian chronogram indicated that the diversification that gave rise to the four geographic haplogroups occurred between  <dig> and  <dig>  myr ago , before the beginning of the contemporary configuration of the sacp .
in this period, the sacp and the adjacent continental shelf that were exposed during regression periods  were deeply dissected by paleochannels that carried water from inland basins to the atlantic ocean, before the formation of the patos and mirim lagoons .
the paleochannels were interspersed by pre-cambrian geological formations located in the central-north region of the sacp  and in the southern portion .
in the area of the present day distribution of c. heterophylla, two main paleochannels  correlate very well with the distributional limits of the three haplogroups .
although there were several transgression/regression events between ~ <dig> and  <dig>  myr ago, the transgressions were less intense than the most recent ones .
these data suggest a scenario where the paleochannels may have been the main barriers responsible for the geographically restricted range of the three haplogroups.
the presence of riverine barriers has been suggested to play a significant role in delimiting the geographical distribution of species and generating phylogeographic breaks for species of fauna  <cit>  and flora  <cit>  on the brazilian atlantic coast.
because few phylogeographic studies have been conducted in the sacp, we report here for the first time the result of rivers acting as barriers to gene flow in a species from this region.
this corroborates the proposal  <cit>  of riverine barriers acting as a limiting factor for the geographical distribution of species and contributing to the diversification among populations.
the only sample sites not affected by the great marine transgression that occurred around  <dig>  myr ago were p12- <dig> .
the inland region  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 .
moreover, the bayesian chronogram indicated the mainland haplogroup was the oldest , nearly twice the age of the remaining haplogroups .
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  <cit> .
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 p <dig> and p <dig>
the period of greatest diversification of c. heterophylla seems to have occurred in recent times, less than  <dig>  myr ago .
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  <cit> .
this may have eliminated much of the diversity already present .
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  <cit> .
the same molecular markers showed that most of the genetic diversity of highlands petunia arose in a similar period  <cit> .
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  <cit> : 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  <cit> .
additionally, the quaternary marine transgressions/regressions in the sacp were a major factor in deutamys kempi  diversification  <cit> .
demography
a similar pattern of population expansion was observed in each of our different analyses , suggesting a general demographic pattern for the whole species.
however, although all neutrality tests were negative , most were not significant.
this absence of significance may be explained by being conservative and not very powerful to detect older population size changes  <cit>  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  may be a sign of recent expansion linked to the pleistocene-holocene transition approximately  <dig> kyr ago  <cit> .
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  <cit> , providing a large coastal plain that could be colonized by vegetation .
moreover, with the beginning of the holocene, the increase in temperature and humidity may have induced vegetation expansion  <cit> .
these environmental changes at the beginning of the holocene are present in the geologic record from the sacp and adjacent areas in southern brazil  <cit> .
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  <cit> .
global studies of plant species have also demonstrated population expansions in the same period and have correlated these results with holocene environmental changes  <cit> .
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  <cit> .
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  and to p <dig> and p <dig> 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  <cit> .
CONCLUSIONS
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  <dig>  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  <cit>  and related genus petunia <cit> .

 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. 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 alternated with open formations 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 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. after the last glacial maximum and at the beginning of the holocene, approximately  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. the accentuated moisture in this phase was also reported from the south atlantic coastal plain and adjacent area. 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 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., 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. the south atlantic coastal plain is the largest coastal plain in south america, covering approximately km from 28° north to 34° south. it is a flat area of lowlands occupied mostly by large systems of coastal lakes. it extends ne-sw for approximately km and has an average width of 60– km in the south central portion and 15– km in the northern portion. the main factors that formed the sacp as it is today were oscillatory glaciation cycles during the pleistocene, which promoted sea level transgressions and regressions. during the late quaternary, four sand barriers parallel to the coast were formed and identified in the sacp; three of these barriers originated during the pleistocene, and one originated during the holocene. they were associated with the limits of the marine transgressions of the four major interglacial periods in the last myr that occurred in the marine isotope stages and which are referred to as barrier-lagoon systems i to iv, respectively. 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. the central lowlands 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 wijsman 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 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 sequences are commonly used assess geographic structure in plants and they are especially informative in relatives of petunia. 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. 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: 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? did the environmental and sea level changes during the quaternary result in allopatric processes that isolated and shaped the genetic structure of populations? 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 that have been successfully employed in population studies of wild petunias. 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.  genetic diversity the mean number of individuals genotyped per population was. the representative lengths of the alignment for trnh-psba and trns-trng were 456 bp and 784 bp, respectively. overall, sites were informative, in each region, resulting in a value of  ±  for nucleotide diversity and a value of  ±  for haplotype diversity. moreover, different haplotypes were identified from concatenated spacers. the marker trnh-psba had diversity indices slightly higher than trns-trng, though the amount of variability is similar. however, intra-population diversity varied widely throughout the populations sampled. five populations were monomorphic and most others exhibited relatively low nucleotide diversity. p and p populations had the highest values of nucleotide diversity. the summary statistics obtained for the markers and populations are shown in tables  and respectively. π, nucleotide diversity; h, haplotype diversity; sd, standard deviation; * duplication of 5 bp; tv, transversion; ts, transition. n, number of samplings; h, haplotype diversity; π, nucleotide diversity; sd, standard deviation. phylogeography and population structure the degree of genetic structure of the populations was estimated by amova, which indicated that % of the variation was among populations and % 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. four groups were inferred: a group in the extreme south of the sacp, a group in the central sacp, a group in north of the sacp and a final group comprising inland populations. populations in each group are listed in table  a mantel test showed a significant and moderate correlation between genetic and geographical distances between the populations. spatial autocorrelation analysis with alleles in space, 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. since the maximum distances between individuals within the groups are around the above value, 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. the haplotype network resulted in four geographically structured haplogroups that were concordant with the samova results and are depicted in figure 3a. the four haplotypes of the southernmost population formed one haplogroup. the haplotypes from the other coastal populations formed the haplogroups central coast and north coast, respectively. the remaining haplotypes belonging to the populations p p p and p were grouped into a mainland haplogroup, considering that the last three were collected outside of the sacp. the central haplotype h whose grouping may be considered ambiguous in the network, was grouped with its respective clades according to the bayesian phylogeny. the p and p populations were the only ones with haplotypes belonging to different haplogroups in the network. the h haplotype that belongs to the north coast haplogroup was also sampled in p and the h haplotype that belongs to the central coast haplogroup was sampled in p. a preeminent feature of the network was that all four clades have a star-like pattern, in which infrequent haplotypes are derived from a common central haplotype with only one or two substitutions. it is also noteworthy that the four haplotypes shared among populations were internal to the network haplotypes. phylogenetic relationships and divergence times the bayesian analysis 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 myr ago. the most ancient divergence within c. heterophylla was between the south coast haplogroup and the others at approximately myr ago. the central coast haplogroup separated from the north coast and mainland haplogroups approximately myr ago, and the latter two diverged approximately myr ago. the major diversification inside the haplogroups occurred within the last myr. the times to the most recent common ancestors for the four haplogroups were,, and. intraspecific demography the neutrality tests for each group of populations and for the whole species were mostly non-significant. however, tajima’s d and fu’s fs were negative for all results but significant only for samova groups g and g fs. the bayesian skyline plot for the species suggested that c. heterophylla underwent an increase in effective population size fairly recently, approximately kyr ago. the very small number of haplotypes, which was most likely due to the sample size, precludes the use of this methodology for each clade separately. π, nucleotide diversity; h, haplotype diversity; sd, standard deviation; n, sampling size;. §, significant values. discussion genetic diversity and population structure haplotype diversity in c. heterophylla was similar to that reported for petunia axillaris, p. altiplana, p. bonjardinensis and p. guarapuavensis, 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.,. 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. 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 p and p is possibly the result of very recent contact because these populations are geographically adjacent between two distinct geographic clusters. 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 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, asia, and north america. genetics structure vs geology and climate the bayesian chronogram indicated that the diversification that gave rise to the four geographic haplogroups occurred between and myr ago, before the beginning of the contemporary configuration of the sacp. in this period, the sacp and the adjacent continental shelf that were exposed during regression periods were deeply dissected by paleochannels that carried water from inland basins to the atlantic ocean, before the formation of the patos and mirim lagoons. the paleochannels were interspersed by pre-cambrian geological formations located in the central-north region of the sacp and in the southern portion. in the area of the present day distribution of c. heterophylla, two main paleochannels correlate very well with the distributional limits of the three haplogroups. although there were several transgression/regression events between ~ and myr ago, the transgressions were less intense than the most recent ones. these data suggest a scenario where the paleochannels may have been the main barriers responsible for the geographically restricted range of the three haplogroups. the presence of riverine barriers has been suggested to play a significant role in delimiting the geographical distribution of species and generating phylogeographic breaks for species of fauna and flora on the brazilian atlantic coast. because few phylogeographic studies have been conducted in the sacp, we report here for the first time the result of rivers acting as barriers to gene flow in a species from this region. this corroborates the proposal of riverine barriers acting as a limiting factor for the geographical distribution of species and contributing to the diversification among populations. the only sample sites not affected by the great marine transgression that occurred around myr ago were p12-. the inland region 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. moreover, the bayesian chronogram indicated the mainland haplogroup was the oldest, nearly twice the age of the remaining haplogroups. 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 p and p the period of greatest diversification of c. heterophylla seems to have occurred in recent times, less than myr ago. 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. 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 diversification. demography a similar pattern of population expansion was observed in each of our different analyses, suggesting a general demographic pattern for the whole species. however, although all neutrality tests were negative, 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 may be a sign of recent expansion linked to the pleistocene-holocene transition approximately 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. moreover, with the beginning of the holocene, the increase in temperature and humidity may have induced vegetation expansion. 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. 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 and to p and p 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 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 and related genus petunia. 
climatic and sea level changes during the quaternary played an important role in the evolutionary history of many organisms found in coastal regions. we assessed dna sequences from two cpdna intergenic spacers and analyzed them using a phylogeographic approach. 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 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. to contribute to a better understanding of plant evolution in this environment in southern south america, we focused on calibrachoa heterophylla, an endemic and vulnerable wild petunia species from the south atlantic coastal plain. 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 glacial and interglacial cycles that characterized the quaternary greatly affected the distribution and genetic diversity of plants. the present phylogeographic study reveals the influence of complex geologic and climatic events on patterns of genetic diversification. in the neotropics, few phylogeographic studies have focused on coastal species outside of the atlantic rainforest. the recent population expansion for the entire species may be linked with the event of marine regression after the most recent sea transgression at ~ kya. 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 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, an endemic and vulnerable wild petunia species from the south atlantic coastal plain.  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 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 ~ kya. south-atlantic coastal plainphylogeographypleistocenequaternarygenetic diversityclimatic changespetunia 