Temporal variation in pollen dispersal and breeding structure in a bee-pollinated Neotropical tree

Comparative population genomics in Tabebuia alliance shows evidence of adaptation in Neotropical tree species

The role of natural selection in shaping spatial patterns of genetic diversity in the Neotropics is still poorly understood. Here, we perform a genome scan with 24,751 probes targeting 11,026 loci in two Neotropical Bignoniaceae tree species: Handroanthus serratifolius from the seasonally dry tropical forest (SDTF) and Tabebuia aurea from savannas, and compared with the population genomics of H. impetiginosus from SDTF. OutFLANK detected 29 loci in 20 genes with selection signal in H. serratifolius and no loci in T. aurea. Using BayPass, we found evidence of selection in 335 loci in 312 genes in H. serratifolius, 101 loci in 92 genes in T. aurea, and 448 loci in 416 genes in H. impetiginosus. All approaches evidenced several genes affecting plant response to environmental stress and primary metabolic processes. The three species shared no SNPs with selection signal, but we found SNPs affecting the same gene in pair of species. Handroanthus serratifolius showed differences in allele frequencies at SNPs with selection signal among ecosystems, mainly between Caatinga/Cerrado and Atlantic Forest, while H. impetiginosus had one allele fixed across all populations, and T. aurea had similar allele frequency distribution among ecosystems and polymorphism across populations. Taken together, our results indicate that natural selection related to environmental stress shaped the spatial pattern of genetic diversity in the three species. However, the three species have different geographical distribution and niches, which may affect tolerances and adaption, and natural selection may lead to different signatures due to the differences in adaptive landscapes in different niches.

Patterns of genetic diversity, spatial genetic structure and gene flow in Campomanesia xanthocarpa: insights from SSR markers of different genomic origins

Campomanesia xanthocarpa (Mart.) O. Berg is a South American fruit tree species with important ecological and medicinal properties, which remnants are currently found mainly in isolated forest fragments. In this study, SSR markers from three different genomic origins (gene-linked, nuclear neutral, and organellar) were used to evaluate the patterns of genetic diversity, fine-scale spatial genetic structure and historical gene flow in fragmented forest formations of C. xanthocarpa from the Atlantic Forest in southern Brazil. Our results show that the forest fragments present moderate to high levels of genetic diversity in comparison to species presenting similar life traits, although a trend opposite to expected was observed concerning gene-linked and neutral SSR markers. The fine-scale spatial genetic structure revealed different patterns in short and large distance classes, with a distinct influence of gene-linked and neutral markers in driving the genetic structure in each distance class. The presence of an isolation-by-adaptation pattern implies the need for maintenance of the current remnants to assure the conservation of the private alleles. Finally, as the genetic diversity is found predominantly within forest fragments, programs of seed collection and/or genetic rescue should prioritize a larger number of individuals within each fragment, to increase the sampled diversity.

Landscape barriers to pollen and seed flow in the dioecious tropical tree Astronium fraxinifolium in Brazilian savannah

Gene flow studies provide information on gene exchange between populations, which is essential for developing genetic conservation strategies. Such analyses enable a better understanding of the life history and seed and pollen dispersal mechanisms of plant species. In this study, we investigate pollen and seed flow in a regenerant population of the pioneer species Astronium fraxinifolium in an area degraded during the construction of a hydroelectric dam. We mapped, sampled, sexed, and genotyped 386 individuals in the regenerant population (RP), as well as 128 adult trees located along two highways adjacent to the degraded area; one in Mato Grosso do Sul State (MS) and other in São Paulo State (SP). Parentage analyses was carried out for 370 individuals of the RP population, using as putative parents 348 individuals from RP and all 128 individuals sampled in MS and SP. Based on parentage analysis and eight microsatellite loci, our analyses revealed that for individuals of the RP with an identified father (pollen donor), 1.1% of the pollen was dispersed up to 532 m, while for those with an identified mother (seed donor), 0.5% of seeds were dispersed up to 4,782 m. However, a large proportion of pollen (76.5%) and seeds (57%) immigrated from trees outside the sampled populations. Pollen and seeds were dispersed through a pattern of isolation by distance. Genetic diversity was significantly similar between adults of both highway populations and individuals from RP, with significant levels of inbreeding detected only in RP. Our results demonstrate that the nearest trees contributed pollen and seeds for the recovery of the degraded area, indicating reproductive spatial isolation among the sampled populations due to the damming of the river. Such results help to understand the process of regeneration for A. fraxinifolium in regenerant populations to inform strategies for conservation and environmental recovery with this species.

Multi-Scale Landscape Influences on Genetic Diversity and Adaptive Traits in a Neotropical Savanna Tree

Changes in landscape structure can affect essential population ecological features, such as dispersal and recruitment, and thus genetic processes. Here, we analyze the effects of landscape metrics on adaptive quantitative traits variation, evolutionary potential, and on neutral genetic diversity in populations of the Neotropical savanna tree Tabebuia aurea. Using a multi-scale approach, we sampled five landscapes with two sites of savanna in each. To obtain neutral genetic variation, we genotyped 60 adult individuals from each site using 10 microsatellite loci. We measured seed size and mass. Seeds were grown in nursery in completely randomized experimental design and 17 traits were measured in seedlings to obtain the average, additive genetic variance (V a ) and coefficient of variation (CV a %), which measures evolvability, for each trait. We found that habitat loss increased genetic diversity (He) and allelic richness (AR), and decreased genetic differentiation among populations (F ST ), most likely due to longer dispersal distance of pollen in landscapes with lower density of flowering individuals. Habitat amount positively influenced seed size. Seeds of T. aurea are wind-dispersed and larger seeds may be dispersed to short distance, increasing genetic differentiation and decreasing genetic diversity and allelic richness. Evolvability (CV a %) in root length decreased with habitat amount. Savanna trees have higher root than shoot growth rate in the initial stages, allowing seedlings to obtain water from water tables. Landscapes with lower habitat amount may be more stressful for plant species, due to the lower plant density, edge effects and the negative impacts of agroecosystems. In these landscapes, larger roots may provide higher ability to obtain water, increasing survival and avoiding dying back because of fire. Despite the very recent agriculture expansion in Central Brazil, landscape changes are affecting neutral and adaptive variation in T. aurea. Several populations have low additive genetic variation for some traits and thus, may have limited evolvability, which may jeopardize species long-term persistence. The effect of habitat loss on highly variable neutral loci may only be detected after a certain threshold of population size is attained, that could become dangerously small masking important losses of heterozygosity endangering species conservation.

Endangered but genetically stable-Erythrophleum fordii within Feng Shui woodlands in suburbanized villages

Feng Shui woodlands are naturally or artificially formed green areas in southern China. They are precious for maintaining ecosystem balance in modern semiurban environments. However, they are generally small and geographically isolated from each other, and the status of genetic diversity of the plant species within them has been almost neglected. Therefore, we studied the genetic diversity of the endangered Erythrophleum fordii in eight Feng Shui woodlands (a total of 1,061 individuals) in Guangzhou, a large city in southern China, using microsatellites. For comparison, one population with 33 individuals sampled in a nature reserve was also studied. Although our results indicate that significant demographic declines occurred historically in E. fordii, such declines have not resulted in consistent reductions in genetic variation over generations in Feng Shui populations in the recent past, and the levels of genetic variation in these populations were higher than or comparable to the genetic variation of the population in the nature reserve. In addition, our parentage and paternity analyses indicated widespread and potential long‐distance pollen flow within one Feng Shui woodland, indicating the presence of an unbroken pollination network, which would at least partially alleviate the genetic erosion due to habitat fragmentation and the unequal gene contributions of E. fordii parents to their progenies when favorable recruitment habitats are absent under most of the parent trees. Overall, our results suggest that E. fordii in Feng Shui woodlands may not be driven to extinction in the near future. Nevertheless, uncontrolled fast urban development with a lack of awareness of Feng Shui woodlands will cause the local extinction of E. fordii, which has already happened in some Feng Shui woodlands.

Genome assembly of the Pink Ipê (Handroanthus impetiginosus, Bignoniaceae), a highly valued, ecologically keystone Neotropical timber forest tree

Background

Handroanthus impetiginosus (Mart. ex DC.) Mattos is a keystone Neotropical hardwood tree widely distributed in seasonally dry tropical forests of South and Mesoamerica. Regarded as the “new mahogany,” it is the second most expensive timber, the most logged species in Brazil, and currently under significant illegal trading pressure. The plant produces large amounts of quinoids, specialized metabolites with documented antitumorous and antibiotic effects. The development of genomic resources is needed to better understand and conserve the diversity of the species, to empower forensic identification of the origin of timber, and to identify genes for important metabolic compounds.

Findings

The genome assembly covers 503.7 Mb (N50 = 81 316 bp), 90.4% of the 557-Mbp genome, with 13 206 scaffolds. A repeat database with 1508 sequences was developed, allowing masking of ∼31% of the assembly. Depth of coverage indicated that consensus determination adequately removed haplotypes assembled separately due to the extensive heterozygosity of the species. Automatic gene prediction provided 31 688 structures and 35 479 messenger RNA transcripts, while external evidence supported a well-curated set of 28 603 high-confidence models (90% of total). Finally, we used the genomic sequence and the comprehensive gene content annotation to identify genes related to the production of specialized metabolites.

Conclusions

This genome assembly is the first well-curated resource for a Neotropical forest tree and the first one for a member of the Bignoniaceae family, opening exceptional opportunities to empower molecular, phytochemical, and breeding studies. This work should inspire the development of similar genomic resources for the largely neglected forest trees of the mega-diverse tropical biomes.

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.

Climatic changes can drive the loss of genetic diversity in a Neotropical savanna tree species

The high rates of future climatic changes, compared with the rates reported for past changes, may hamper species adaptation to new climates or the tracking of suitable conditions, resulting in significant loss of genetic diversity. Trees are dominant species in many biomes and because they are long-lived, they may not be able to cope with ongoing climatic changes. Here, we coupled ecological niche modelling (ENM) and genetic simulations to forecast the effects of climatic changes on the genetic diversity and the structure of genetic clusters. Genetic simulations were conditioned to climatic variables and restricted to plant dispersal and establishment. We used a Neotropical savanna tree as species model that shows a preference for hot and drier climates, but with low temperature seasonality. The ENM predicts a decreasing range size along the more severe future climatic scenario. Additionally, genetic diversity and allelic richness also decrease with range retraction and climatic genetic clusters are lost for both future scenarios, which will lead genetic variability to homogenize throughout the landscape. Besides, climatic genetic clusters will spatially reconfigure on the landscape following displacements of climatic conditions. Our findings indicate that climate change effects will challenge population adaptation to new environmental conditions because of the displacement of genetic ancestry clusters from their optimal conditions.

Mixed pollen load and late-acting self-incompatibility flexibility in Adenocalymma peregrinum (Miers) L.G. Lohmann (Bignonieae: Bignoniaceae)

Mixed cross and self-pollen load on the stigma (mixed pollination) of species with late-acting self-incompatibility system (LSI) can lead to self-fertilized seed production. This "cryptic self-fertility" may allow selfed seedling development in species otherwise largely self-sterile. Our aims were to check if mixed pollinations would lead to fruit set in LSI Adenocalymma peregrinum, and test for evidence of early-acting inbreeding depression in putative selfed seeds from mixed pollinations. Experimental pollinations were carried out in a natural population. Fruit and seed set from self-, cross and mixed pollinations were analysed. Further germination tests were carried out for the seeds obtained from treatments. Our results confirm self-incompatibility, and fruit set from cross-pollinations was three-fold that from mixed pollinations. This low fruit set in mixed pollinations is most likely due to a greater number of self- than cross-fertilized ovules, which promotes LSI action and pistil abortion. Likewise, higher percentage of empty seeds in surviving fruits from mixed pollinations compared with cross-pollinations is probably due to ovule discounting caused by self-fertilization. Moreover, germinability of seeds with developed embryos was lower in fruits from mixed than from cross-pollinations, and the non-viable seeds from mixed pollinations showed one-third of the mass of those from cross-pollinations. The great number of empty seeds, lower germinability, lower mass of non-viable seeds, and higher variation in seed mass distribution in mixed pollinations, strongly suggests early-acing inbreeding depression in putative selfed seeds. In this sense, LSI and inbreeding depression acting together probably constrain self-fertilized seedling establishment in A. peregrinum.

Relaxed random walk model coupled with ecological niche modeling unravel the dispersal dynamics of a Neotropical savanna tree species in the deeper Quaternary

Understanding the dispersal routes of Neotropical savanna tree species is an essential step to unravel the effects of past climate change on genetic patterns, species distribution and population demography. Here we reconstruct the demographic history and dispersal dynamics of the Neotropical savanna tree species Tabebuia aurea to understand the effects of Quaternary climate change on its current spatial patterns of genetic diversity. We sampled 285 individuals from 21 populations throughout Brazilian savannas and sequenced all individuals for three chloroplast intergenic spacers and ITS nrDNA. We analyzed data using a multi-model inference framework by coupling the relaxed random walk model (RRW), ecological niche modeling (ENM) and statistical phylogeography. The most recent common ancestor of T. aurea lineages dated from ~4.0 ± 2.5 Ma. T. aurea lineages cyclically dispersed from the West toward the Central-West Brazil, and from the Southeast toward the East and Northeast Brazil, following the paleodistribution dynamics shown by the ENMs through the last glacial cycle. A historical refugium through time may have allowed dispersal of lineages among populations of Central Brazil, overlapping with population expansion during interglacial periods and the diversification of new lineages. Range and population expansion through the Quaternary were, respectively, the most frequent prediction from ENMs and the most likely demographic scenario from coalescent simulations. Consistent phylogeographic patterns among multiple modeling inferences indicate a promising approach, allowing us to understand how cyclical climate changes through the Quaternary drove complex population dynamics and the current patterns of species distribution and genetic diversity.

Recovering species demographic history from multi-model inference: the case of a Neotropical savanna tree species

Background

Glaciations were recurrent throughout the Quaternary and potentially shaped species genetic structure worldwide by affecting population dynamics. Here, we implemented a multi-model inference approach to recover the distribution dynamics and demographic history of a Neotropical savanna tree, Tabebuia aurea (Bignoniaceae). Exploring different algorithms and paleoclimatic simulations, we used ecological niche modelling to generate alternative hypotheses of potential demographic changes through the last glacial cycle and estimated genetic parameters using coalescent modelling.

Results

Comparing predictions from demographic hypotheses with genetic parameters of modern populations, our findings revealed a likely scenario of population decline, with spatial displacement towards Northeast Brazil from the last glacial maximum to the mid-Holocene. Subsequently, populations expanded in response to the return of the climatically suitable conditions in Central-West Brazil. Nevertheless, a wide historical refugium across Central Brazil likely maintained large populations connected throughout time. The expected genetic signatures from such predicted distribution dynamics are also corroborated by spatial genetic structure observed in modern populations.

Conclusion

By exploring uncertainties inherent in multiple working hypotheses, we have shown that multi-model inference is a fruitful and efficient approach to recover the nature, timing and geographical context of the Tabebuia aurea population dynamic in response to the Quaternary climate changes.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-014-0213-0) contains supplementary material, which is available to authorized users.