Phylogeography of Eriotheca species complex: insights into the origin and range expansion of apomictic and polyploid trees in Neotropical Savannas

Polyploidy and whole genome duplication are major evolutionary drivers in plants. Climate variations during the Pleistocene have influenced distribution and range expansion worldwide. Similar trends have been reported for Cerrado plants, but no attempt has been made to link phylogeography with ploidy and breeding changes. Thus, we aimed to (i) assess ploidy and genome size of Eriotheca estevesiae Carv.-Sobr., and compare it with E. pubescens (Mart.) Schott & Endl. (Both included into the Eriotheca Stellate Trichome Species Complex - ESTSC). (ii) Subsequently, we investigated their phylogeography to see whether genetic structure and range expansion trends were similar to those previously described for the Cerrado biome. Finally (iii), we discuss whether ESTSC phylogeographic patterns could be associated with geographic parthenogenesis processes. Common cytogenetic techniques and flow cytometry were used to confirm chromosome number and genome size of E. estevesiae. We used three cpDNA regions to analyse 14 ESTSC Cerrado populations, for which we also obtained ploidy level and breeding information. We investigated haplotype diversity, population structure and tested neutrality, aiming to reconstruct phylogeographic scenarios. We found three ploidy levels and eight cpDNA haplotypes in ESTSC, one shared by most populations. Haplotype and ploidy distribution corroborated that E. pubescens, the widely distributed polyploid and apomictic species, may have originated from northern diploid and probably sexual E. estevesiae. Matrilinear cpDNA links support the idea that apomixis and polyploidy in ESTSC may have allowed range expansion during the Pleistocene, in a process analogous to the geographic parthenogenesis described elsewhere.

Stomatal size, ploidy and polyembryony in Eriotheca Stellate Trichome Species Complex (Bombacoideae - Malvaceae) in the Cerrados of Brazil

Geographic parthenogenesis, range expansion of apomictic plants after climate changes, has been described for Northern Hemisphere gametophytic apomicts. But similar trends have been observed for sporophytic apomicts of Cerrado, the savannas in Brazil. Eriotheca pubescens is a common Cerrado tree, an agamic complex of either hexaploid/polyembryonic apomicts or tetraploid/monoembryonic sexual individuals. Some populations have been described as a new species, Eriotheca estevesiae, all included in the Eriotheca Stellate Trichome Species Complex (ESTSC). Since breeding systems and ploidy are clearly associated with polyembryony and stomatal size, we used these ancillary features to map the reproductive and ploidy level traits of E. pubescens and E. estevesiae. Leaves and seeds were collected from individuals of 19 populations. Seeds were evaluated for the presence of polyembryony and leaves for stomatal measurements. Eight populations were monoembryonic while another eight were polyembryonic and for other three, the embryonic pattern was not readily verified. E. pubescens polyembryonic and hexaploid populations formed a homogeneous group, but monoembryonic plants were more variable. E. estevesiae populations were monoembryonic with smaller stomata. In contrast, some E. pubescens monoembryonic populations further south presented larger stomata. Despite these outliers, possibly mixed populations, stomatal size and embryonic pattern differed from northern to southern populations. Embryonic pattern and stomatal size indicated that northernmost populations of Eriotheca STSC (E. estevesiae) are diploid and sexual. Southernmost populations, mostly polyembryonic and with large stomata, are hexaploid and apomictic. This is in agreement with geographic parthenogenesis and range expansion of apomictic lineages to southern habitats available after the last glacial maximum.

Long-distance pollen and seed dispersal and inbreeding depression in Hymenaea stigonocarpa (Fabaceae: Caesalpinioideae) in the Brazilian savannah

Hymenaea stigonocarpa is a neotropical tree that is economically important due to its high-quality wood; however, because it has been exploited extensively, it is currently considered threatened. Microsatellite loci were used to investigate the pollen and seed dispersal, mating patterns, spatial genetic structure (SGS), genetic diversity, and inbreeding depression in H. stigonocarpa adults, juveniles, and open-pollinated seeds, which were sampled from isolated trees in a pasture and trees within a forest fragment in the Brazilian savannah. We found that the species presented a mixed mating system, with population and individual variations in the outcrossing rate (0.53-1.0). The studied populations were not genetically isolated due to pollen and seed flow between the studied populations and between the populations and individuals located outside of the study area. Pollen and seed dispersal occurred over long distances (>8 km); however, the dispersal patterns were isolated by distance, with a high frequency of mating occurring between near-neighbor trees and seeds dispersed near the parent trees. The correlated mating for individual seed trees was higher within than among fruits, indicating that fruits present a high proportion of full-sibs. Genetic diversity and SGS were similar among the populations, but offspring showed evidence of inbreeding, mainly originating from mating among related trees, which suggests inbreeding depression between the seed and adult stages. Selfing resulted in a higher inbreeding depression than mating among relatives, as assessed through survival and height. As the populations are not genetically isolated, both are important targets for in situ conservation to maintain their genetic diversity; for ex situ conservation, seeds can be collected from at least 78 trees in both populations separated by at least 250 m.

Genetic and Historical Colonization Analyses of an Endemic Savanna Tree, Qualea grandiflora, Reveal Ancient Connections Between Amazonian Savannas and Cerrado Core

The evolutionary processes underlying the high diversity and endemism in the Cerrado, the most extensive Neotropical savanna, remain unclear, including the factors promoting the presence and evolution of savanna enclaves in the Amazon forest. In this study, we investigated the effects of past climate changes on genetic diversity, dynamics of species range and the historical connections between the savanna enclaves and Cerrado core for Qualea grandiflora, a tree species widely distributed in the biome. Totally, 40 populations distributed in the Cerrado core and Amazon savannas were analyzed using chloroplast and nuclear DNA sequences. We used phylogeographic, coalescent and ecological niche modeling approaches. Genetic data revealed a phylogeographic structure shaped by Pleistocene climatic oscillations. An eastern-western split in the Cerrado core was observed. The central portion of the Cerrado core harbored most of the sampled diversity for cpDNA. Ecological niche models predicted the presence of a large historical refuge in this region and multiple small refuges in peripheral areas. Relaxed Random Walk (RRW) models indicated the ancestral population in the north-western border of the central portion of the Cerrado core and cyclical dynamics of colonization related to Pleistocene climatic oscillations. Central and western ancient connections between Cerrado core and Amazonian savannas were observed. No evidence of connections among the Amazonian savannas was detected. Our study highlights the importance of Pleistocene climatic oscillations for structuring the genetic diversity of Q. grandiflora and complex evolutionary history of ecotonal areas in the Cerrado. Our results do not support the recent replacement of a large area in the Amazon forest by savanna vegetation. The Amazonian savannas appear to be fragmented and isolated from each other, evolving independently a long ago.

Assessing the Ability of Chloroplast and Nuclear DNA Gene Markers to Verify the Geographic Origin of Jatoba (Hymenaea courbaril L.) Timber

Deforestation-reinforced by illegal logging-is a serious problem in many tropical regions and causes pervasive environmental and economic damage. Existing laws that intend to reduce illegal logging need efficient, fraud resistant control methods. We developed a genetic reference database for Jatoba (Hymenaea courbaril), an important, high value timber species from the Neotropics. The data set can be used for controls on declarations of wood origin. Samples from 308 Hymenaea trees from 12 locations in Brazil, Bolivia, Peru, and French Guiana have been collected and genotyped on 10 nuclear microsatellites (nSSRs), 13 chloroplast SNPs (cpSNP), and 1 chloroplast indel marker. The chloroplast gene markers have been developed using Illumina DNA sequencing. Bayesian cluster analysis divided the individuals based on the nSSRs into 8 genetic groups. Using self-assignment tests, the power of the genetic reference database to judge on declarations on the location has been tested for 3 different assignment methods. We observed a strong genetic differentiation among locations leading to high and reliable self-assignment rates for the locations between 50% to 100% (average of 88%). Although all 3 assignment methods came up with similar mean self-assignment rates, there were differences for some locations linked to the level of genetic diversity, differentiation, and heterozygosity. Our results show that the nuclear and chloroplast gene markers are effective to be used for a genetic certification system and can provide national and international authorities with a robust tool to confirm legality of timber.

Pervasive migration across rainforest and sandy coastal plain Aechmea nudicaulis (Bromeliaceae) populations despite contrasting environmental conditions

Understanding the colonization of extreme marginal habitats and the relative roles of space and environment in maintaining peripheral populations remains challenging. Here, we leverage a system of pairs of rainforest and sandy coastal plain communities that allow us to decouple spatial and environmental effects in the population structure and migration rates of the bromeliad Aechmea nudicaulis. Structure and gene flow between populations were estimated from Bayesian clustering and coalescent-based migration models applied to chloroplast sequence and nuclear microsatellite data. Contrary to our initial expectation, the sharp environmental gradient between rainforest and sandy plains does not seem to have affected the colonization and migration dynamics in A. nudicaulis. Our analyses uncover pervasive gene flow between neighbouring habitats in both chloroplast and nuclear data despite the striking differences in environmental conditions. This result is consistent with a scenario of repeated colonization of the sandy coastal plains from forest populations through seed dispersal, as well as the maintenance of gene flow between habitats through pollination. We also recovered a broad north/south population structure that has been found in other Atlantic rainforest groups and possibly reflects older phylogeographic dynamics.

Adulteration and Contamination of Commercial Sap of Hymenaea Species

The Hymenaea stigonocarpa and Hymenaea martiana species, commonly known as “jatobá,” produce a sap which is extracted by perforation of the trunk and is commonly used in folk medicine as a tonic. For this study, the authenticity of commercial samples of jatobá was verified by the identification of the main compounds and multivariate analysis and contamination by microbial presence analysis. The acute toxicity of the authentic jatobá sap was also evaluated. The metabolites composition and multivariate analysis revealed that none of the commercial samples were authentic. In the microbiological contamination analysis, five of the six commercial samples showed positive cultures within the range of 1,700–100,000 CFU/mL and the authentic sap produced no signs of toxicity, and from a histological point of view, there was the maintenance of tissue integrity. In brief, the commercial samples were deemed inappropriate for consumption and represent a danger to the population.

Climatic drivers of leaf traits and genetic divergence in the tree Annona crassiflora: a broad spatial survey in the Brazilian savannas

The Cerrado is the largest South American savanna and encompasses substantial species diversity and environmental variation. Nevertheless, little is known regarding the influence of the environment on population divergence of Cerrado species. Here, we searched for climatic drivers of genetic (nuclear microsatellites) and leaf trait divergence in Annona crassiflora, a widespread tree in the Cerrado. The sampling encompassed all phytogeographic provinces of the continuous area of the Cerrado and included 397 individuals belonging to 21 populations. Populations showed substantial genetic and leaf trait divergence across the species' range. Our data revealed three spatially defined genetic groups (eastern, western and southern) and two morphologically distinct groups (eastern and western only). The east-west split in both the morphological and genetic data closely mirrors previously described phylogeographic patterns of Cerrado species. Generalized linear mixed effects models and multiple regression analyses revealed several climatic factors associated with both genetic and leaf trait divergence among populations of A. crassiflora. Isolation by environment (IBE) was mainly due to temperature seasonality and precipitation of the warmest quarter. Populations that experienced lower precipitation summers and hotter winters had heavier leaves and lower specific leaf area. The southwestern area of the Cerrado had the highest genetic diversity of A. crassiflora, suggesting that this region may have been climatically stable. Overall, we demonstrate that a combination of current climate and past climatic changes have shaped the population divergence and spatial structure of A. crassiflora. However, the genetic structure of A. crassiflora reflects the biogeographic history of the species more strongly than leaf traits, which are more related to current climate.

Morphological analyses suggest a new taxonomic circumscription for Hymenaea courbaril L. (Leguminosae, Caesalpinioideae)

Hymenaea is a genus of the Resin-producing Clade of the tribe Detarieae (Leguminosae: Caesalpinioideae) with 14 species. Hymenaea courbaril is the most widespread species of the genus, ranging from southern Mexico to southeastern Brazil. As currently circumscribed, Hymenaea courbaril is a polytypic species with six varieties: var. altissima, var. courbaril, var. longifolia, var. stilbocarpa, var. subsessilis, and var. villosa. These varieties are distinguishable mostly by traits related to leaflet shape and indumentation, and calyx indumentation. We carried out morphometric analyses of 14 quantitative (continuous) leaf characters in order to assess the taxonomy of Hymenaea courbaril under the Unified Species Concept framework. Cluster analysis used the Unweighted Pair Group Method with Arithmetic Mean (UPGMA) based on Bray-Curtis dissimilarity matrices. Principal Component Analyses (PCA) were carried out based on the same morphometric matrix. Two sets of Analyses of Similarity and Non Parametric Multivariate Analysis of Variance were carried out to evaluate statistical support (1) for the major groups recovered using UPGMA and PCA, and (2) for the varieties. All analyses recovered three major groups coincident with (1) var. altissima, (2) var. longifolia, and (3) all other varieties. These results, together with geographical and habitat information, were taken as evidence of three separate metapopulation lineages recognized here as three distinct species. Nomenclatural adjustments, including reclassifying formerly misapplied types, are proposed.

Rock outcrop orchids reveal the genetic connectivity and diversity of inselbergs of northeastern Brazil

BACKGROUND

Because of their fragmented nature, inselberg species are interesting biological models for studying the genetic consequences of disjoint populations. Inselbergs are commonly compared with oceanic islands, as most of them display a marked ecological isolation from the surrounding area. The isolation of these rock outcrops is reflected in the high number of recorded endemic species and the strong floristic differences between individual inselbergs and adjacent habitats. We examined the genetic connectivity of orchids Epidendrum cinnabarinum and E. secundum adapted to Neotropical inselbergs of northeastern Brazil. Our goals were to identify major genetic divergences or disjunctions across the range of the species and to investigate potential demographic and evolutionary mechanisms leading to lineage divergence in Neotropical mountain ecosystems.

RESULTS

Based on plastid markers, high genetic differentiation was found for E. cinnabarinum (FST = 0.644) and E. secundum (FST = 0.636). Haplotypes were not geographically structured in either taxon, suggesting that restricted gene flow and genetic drift may be significant factors influencing the diversification of these inselberg populations. Moreover, strong differentiation was found between populations over short spatial scales, indicating substantial periods of isolation among populations. For E. secundum, nuclear markers indicated higher gene flow by pollen than by seeds.

CONCLUSIONS

The comparative approach adopted in this study contributed to the elucidation of patterns in both species. Our results confirm the ancient and highly isolated nature of inselberg populations. Both species showed similar patterns of genetic diversity and structure, highlighting the importance of seed-restricted gene flow and genetic drift as drivers of plant diversification in terrestrial islands such as inselbergs.

Geographical genetic structuring and phenotypic variation in the Vellozia hirsuta (Velloziaceae) ochlospecies complex

PREMISE OF THE STUDY

Vellozia hirsuta forms a complex presenting wide morphological and anatomical variation, resulting in five specific names and 14 morpho-anatomical patterns occurring in disjunct populations. We carried out a phylogeographical study to investigate the existence of correlation among the genetic and morphological patterns within this complex, and to determine whether it is composed of various species or should be treated as an ochlospecies, a species having widely polymorphic and weakly polytypic complex variation, with morphological characteristics varying independently.

METHODS

We carried out phylogeographical analyses using cpDNA rpl32F-trnL intergenic region.

KEY RESULTS

We found 20 haplotypes in 23 populations sampled. The populations are genetically structured (Φ(ST) = 0.818) into four phylogeographical groups demonstrating geographical structuring but with no correlation with morpho-anatomical patterns. Our analyses do not support recognizing any of the species now synonymized under Vellozia hirsuta. The northern populations were the most genetically differentiated and could be considered a distinct taxon, as they are also morphologically different.

CONCLUSIONS

It is recommended that Vellozia hirsuta be considered a single enormously variable species. The patterns of variation within V. hirsuta probably are related to climatic changes that occurred during the Pleistocene Epoch in tropical Brazil when reductions in forest cover favored the expansion of V. hirsuta populations into extensive lowland areas. The expansion of forest cover at the end of the glaciations would have again restricted the occurrence of campos rupestres vegetation to high elevations, which constitute the current centers of diversity of this species.

Reproductive biology and species geographical distribution in the Melastomataceae: a survey based on New World taxa

BACKGROUND AND AIMS

Apomictic plants are less dependent on pollinator services and able to occupy more diverse habitats than sexual species. However, such assumptions are based on temperate species, and comparable evaluation for species-rich Neotropical taxa is lacking. In this context, the Melastomataceae is a predominantly Neotropical angiosperm family with many apomictic species, which is common in the Campos Rupestres, endemism-rich vegetation on rocky outcrops in central Brazil. In this study, the breeding system of some Campo Rupestre Melastomataceae was evaluated, and breeding system studies for New World species were surveyed to test the hypothesis that apomixis is associated with wide distributions, whilst sexual species have more restricted areas.

METHODS

The breeding systems of 20 Campo Rupestre Melastomataceae were studied using hand pollinations and pollen-tube growth analysis. In addition, breeding system information was compiled for 124 New World species of Melastomataceae with either wide (>1000 km) or restricted distributions.

KEY RESULTS

Most (80 %) of the Campo Rupestre species studied were self-compatible. Self-incompatibility in Microlicia viminalis was associated with pollen-tube arrest in the style, as described for other Melastomataceae, but most self-incompatible species analysed showed pollen-tube growth to the ovary irrespective of pollination treatment. Apomictic species showed lower pollen viability and were less frequent among the Campo Rupestre plants. Among the New World species compiled, 43 were apomictic and 77 sexual (24 self-incompatible and 53 self-compatible). Most apomictic (86 %) and self-incompatible species (71 %) presented wide distributions, whilst restricted distributions predominate only among the self-compatible ones (53 %).

CONCLUSIONS

Self-compatibility and dependence on biotic pollination were characteristic of Campo Rupestre and narrowly distributed New World Melastomataceae species, whilst apomictics are widely distributed. This is, to a certain extent, similar to the geographical parthenogenesis pattern of temperate apomictics.

Large-scale phylogeography of the disjunct Neotropical tree species Schizolobium parahyba (Fabaceae-Caesalpinioideae)

Neotropical rainforests exhibit high levels of endemism and diversity. Although the evolutionary genetics of plant diversification has garnered increased interest, phylogeographic studies of widely distributed species remain scarce. Here we describe chloroplast and nuclear variation patterns in Schizolobium parahyba (Fabaceae), a widespread tree in Neotropical rainforests that harbor two varieties with a disjunct distribution. Chloroplast and nuclear sequence analyses yielded 21 and 4 haplotypes, respectively. Two genetic diversity centers that correlate with the two known varieties were identified: the Southeastern Atlantic forest and the Amazonian basin. In contrast, the populations from southern and northeastern Atlantic forests and Andean-Central American forests exhibited low levels of genetic diversity and divergent haplotypes, likely related to historical processes that impact the flora and fauna in these regions, such as a founder's effect after dispersion and demographic expansion. Phylogeographic and demographic patterns suggest that episodes of genetic isolation and dispersal events have shaped the evolutionary history for this species, and different patterns have guided the evolution of S. parahyba. Moreover, the results of this study suggest that the dry corridor formed by Cerrado and Caatinga ecoregions and the Andean uplift acted as barriers to this species' gene flow, a picture that may be generalized to most of the plant biodiversity tropical woodlands and forests. These results also reinforce the importance of evaluating multiple genetic markers for a more comprehensive understanding of population structure and history. Our results provide insight into the conservation efforts and ongoing work on the genetics of population divergence and speciation in these Neotropical rainforests.

Low but structured chloroplast diversity in Atherosperma moschatum (Atherospermataceae) suggests bottlenecks in response to the Pleistocene glacials

BACKGROUND AND AIMS

The cool temperate rainforests of Australia were much reduced in range during the cold and dry glacial periods, although genetic evidence indicates that two key rainforest species, Nothofagus cunninghamii and Tasmannia lanceolata, survived within multiple locations and underwent only local range expansions at the end of the Last Glacial. To better understand the glacial response of a co-occurring but wind-dispersed and less cold-tolerant rainforest tree species, Atherosperma moschatum, a chloroplast phylogeographic study was undertaken.

METHODS

A total of 3294 bp of chloroplast DNA sequence was obtained for 155 samples collected from across the species' range.

KEY RESULTS

The distribution of six haplotypes observed in A. moschatum was geographically structured with an inferred ancestral haplotype restricted to Tasmania, while three non-overlapping and endemic haplotypes were found on the mainland of south-eastern Australia. Last glacial refugia for A. moschatum are likely to have occurred in at least one location in western Tasmania and in Victoria and within at least two locations in the Great Dividing Range of New South Wales. Nucleotide diversity of A. moschatum was lower (π = 0·00021) than either N. cunninghamii (0·00101) or T. lanceolata (0·00073), and was amongst the lowest recorded for any tree species.

CONCLUSIONS

This study provides evidence for past bottlenecks having impacted the chloroplast diversity of A. moschatum as a result of the species narrower climatic niche during glacials. This hypothesis is supported by the star-like haplotype network and similar estimated rates of chloroplast DNA substitution for A. moschatum and the two more cold tolerant and co-occurring species that have higher chloroplast diversity, N. cunninghamii and T. lanceolata.

Past climate changes explain the phylogeography of Vitellaria paradoxa over Africa

The evolution of the savanna biome has been deeply marked by repeated contraction/expansion phases due to climate perturbations during the Quaternary period. In this study, we investigated the impact of the last glacial maximum (LGM) on the present genetic pattern of Vitellaria paradoxa (shea tree), a major African savanna tree. A range-wide sampling of the species enabled us to sample 374 individuals from 71 populations distributed throughout sub-Sahelian Africa. Trees were genotyped using 3 chloroplasts and 12 nuclear microsatellites, and were sequenced for 2 polymorphic chloroplast intergenic spacers. Analyses of genetic diversity and structure were based on frequency-based and Bayesian methods. Potential distributions of V. paradoxa at present, during the LGM and the last interglacial period, were examined using DIVA-GIS ecological niche modelling (ENM). Haplotypic and allelic richness varied significantly across the range according to chloroplast and nuclear microsatellites, which pointed to higher diversity in West Africa. A high but contrasted level of differentiation was revealed among populations with a clear phylogeographic signal, with both nuclear (F(ST) = 0.21; R(ST) = 0.28; R(ST) > R(ST) (permuted)) and chloroplast simple sequence repeats (SSRs) (G(ST) = 0.81; N(ST) = 0.90; N(ST) > N(ST) (permuted)). We identified a strong geographically related structure separating western and eastern populations, and a substructure in the eastern part of the area consistent with subspecies distinction. Using ENM, we deduced that perturbations during the LGM fragmented the potential eastern distribution of shea tree, but not its distribution in West Africa. Our main results suggest that climate variations are the major factor explaining the genetic pattern of V. paradoxa.

Phylogeography of the endangered rosewood Dalbergia nigra (Fabaceae): insights into the evolutionary history and conservation of the Brazilian Atlantic Forest

The Brazilian rosewood (Dalbergia nigra) is an endangered tree endemic to the central Brazilian Atlantic Forest, one of the world's most threatened biomes. The population diversity, phylogeographic structure and demographic history of this species were investigated using the variation in the chloroplast DNA (cpDNA) sequences of 185 individuals from 19 populations along the geographical range of the species. Fifteen haplotypes were detected in the analysis of 1297 bp from two non-coding sequences, trnV-trnM and trnL. We identified a strong genetic structure (F(ST)=0.62, P<0.0001), with a latitudinal separation into three phylogeographic groups. The two northernmost groups showed evidence of having maintained historically larger populations than the southernmost group. Estimates of divergence times between these groups pointed to vicariance events in the Middle Pleistocene (ca. 350,000-780,000 years ago). The recurrence of past climatic changes in the central part of the Atlantic forest, with cycles of forest expansion and contraction, may have led to repeated vicariance events, resulting in the genetic differentiation of these groups. Based on comparisons among the populations of large reserves and small, disturbed fragments of the same phylogeographic group, we also found evidence of recent anthropogenic effects on genetic diversity. The results were also analysed with the aim of contributing to the conservation of D. nigra. We suggest that the three phylogeographic groups could be considered as three distinct management units. Based on the genetic diversity and uniqueness of the populations, we also indicate priority areas for conservation.