SNP Markers as a Successful Molecular Tool for Assessing Species Identity and Geographic Origin of Trees in the Economically Important South American Legume Genus Dipteryx.

Dipteryx timber has been heavily exploited in South America since 2000s due to the increasing international demand for hardwood. Developing tools for the genetic identification of Dipteryx species and their geographical origin can help to promote legal trading of timber. A collection of 800 individual trees, belonging to 6 different Dipteryx species, was genotyped based on 171 molecular markers. After the exclusion of markers out of Hardy-Weinberg equilibrium or with no polymorphism or low amplification, 83 nuclear, 29 chloroplast, 13 mitochondrial single nucleotide polymorphisms (SNPs), and 2 chloroplast and 5 mitochondrial INDELS remained. Six genetic groups were identified using Bayesian Structure analyses of the nuclear SNPs, which corresponded to the different Dipteryx species collected in the field. Seventeen highly informative markers were identified as suitable for species identification and obtained self-assignment success rates to species level of 78-96%. An additional set of 15 molecular markers was selected to determine the different genetic clusters found in Dipteryx odorata and Dipteryx ferrea, obtaining self-assignment success rates of 91-100%. The success to assign samples to the correct country of origin using all or only the informative markers improved when using the nearest neighbor approach (69-92%) compared to the Bayesian approach (33-80%). While nuclear and chloroplast SNPs were more suitable for differentiating the different Dipteryx species, mitochondrial SNPs were ideal for determining the genetic clusters of D. odorata and D. ferrea. These 32 selected SNPs will be invaluable genetic tools for the accurate identification of species and country of origin of Dipteryx timber.

Molecular evidence for three genetic species of Dipteryx in the Peruvian Amazon.

There is a high international demand for timber from the genus Dipteryx, or "shihuahuaco" as it is known in Peru. Developing tools that allow the identification and discrimination of Dipteryx species is therefore important for supporting management of natural populations and to underpin legal trade of its timber. The objective of this study was the molecular characterization of Dipteryx species in the Peruvian Amazonia. Two plastid regions (cpDNA: trnH-psbA and matK) were sequenced and 11 microsatellite markers (nDNA) were genotyped for 32 individuals identified as Dipteryx charapilla, D. micrantha morphotype 1 and D. micrantha morphotype 2. Using the concatenated sequences of the plastid genes, we identified ten haplotypes that were not shared between the species or between the D. micrantha morphotypes. Haplotypic diversity was greater in D. micrantha morphotype 2 and D. charapilla than in D. micrantha morphotype 1, which presented only one haplotype with a wide distribution in Peru. The microsatellites allowed the discrimination of the same three clades and identified diagnostic alleles for each clade. These results allowed us to demonstrate that the two morphotypes of D. micrantha are different at both the plastid and nuclear markers, which supports the existence of three genetically distinct species in Peru. This study provides information for the genetic discrimination of Dipteryx species and emphasises the importance of conserving the genetic variability of this genus in the Peruvian Amazonia.

Discovery and characterization of new microsatellite loci in Dipteryx alata Vogel (Fabaceae) using next-generation sequencing data.

The use of next-generation sequencing (NGS) technologies provides a great volume of genome sequence data even for non-model species. The development of microsatellite markers using these data is a relatively quick and easy process. Dipteryx alata Vogel (Fabaceae) is an arboreal species from the Cerrado biome and is considered an important plant genetic resource. Here, we report the development of microsatellite markers for D. alata using NGS data. DNA samples from four individuals were sequenced using the Illumina MiSeq platform and high-quality reads were assembled into contigs of the D. alata genome sequence. Microsatellite regions were identified using the IMEX webserver and primer pairs were designed using the Primer3 software. The amplification settings for each locus were optimized. Fluorescent-labeled primers were developed and used to genotype individuals derived from three natural populations of D. alata. Fifty-four microsatellite regions were identified, from which 27 were elected to primer design. Among the amplified loci, 11 were polymorphic, with the number of alleles ranging from 2 to 10. The expected heterozygosity under Hardy-Weinberg Equilibrium (HWE) per locus varied from 0.191 to 0.807. Genotype and allele frequencies for all loci agreed with those expected under HWE and linkage disequilibrium was not significant for all pairs of loci. The probabilities of exclusion of paternity and of combined identity were equal to 0.993 and 5.65 x 10-8, respectively. The markers developed in this study are useful to several types of population genetic studies with D. alata and, eventually, for closely related species.

Multi-objective optimization in systematic conservation planning and the representation of genetic variability among populations.

Biodiversity crises have led scientists to develop strategies for achieving conservation goals. The underlying principle of these strategies lies in systematic conservation planning (SCP), in which there are at least 2 conflicting objectives, making it a good candidate for multi-objective optimization. Although SCP is typically applied at the species level (or hierarchically higher), it can be used at lower hierarchical levels, such as using alleles as basic units for analysis, for conservation genetics. Here, we propose a method of SCP using a multi-objective approach. We used non-dominated sorting genetic algorithm II in order to identify the smallest set of local populations of Dipteryx alata (baru) (a Brazilian Cerrado species) for conservation, representing the known genetic diversity and using allele frequency information associated with heterozygosity and Hardy-Weinberg equilibrium. We worked in 3 variations for the problem. First, we reproduced a previous experiment, but using a multi-objective approach. We found that the smallest set of populations needed to represent all alleles under study was 7, corroborating the results of the previous study, but with more distinct solutions. In the 2nd and 3rd variations, we performed simultaneous optimization of 4 and 5 objectives, respectively. We found similar but refined results for 7 populations, and a larger portfolio considering intra-specific diversity and persistence with populations ranging from 8-22. This is the first study to apply multi-objective algorithms to an SCP problem using alleles at the population level as basic units for analysis.

orchard: Paternity program for autotetraploid species.

Advances in molecular marker technology have provided new opportunities to study the population genetics of polyploid taxa. Paternity analysis using microsatellite markers can be used in detection of gene flow between individuals and populations, in mating system analysis, to identify factors that influence fecundity and fertility, to identify behaviour of parent-offspring relationships and in the analysis of the reproductive success of different ecological groups. As there is no specific program for carrying out paternity analysis in tetraploid species, specialized software was designed for the assignment of paternity for autotetraploid species. orchard is a novel implementation of exclusion and likelihood statistics for carrying out paternity analysis of autotetraploids. First, the program performs an exclusion method, and then, a likelihood statistic is used with nonexcluded candidate fathers. Optional features include estimation of allele dosage of known mother trees and the estimation of pollen flow distances. orchard was tested using a data set of microsatellite data of Dipteryx odorata, a tetraploid Amazonian tree species.

Correlation between genetic diversity and environmental suitability: taking uncertainty from ecological niche models into account.

The hindcast of shifts in the geographical ranges of species as estimated by ecological niche modelling (ENM) has been coupled with phylogeographical patterns, allowing the inference of past processes that drove population differentiation and genetic variability. However, more recently, some studies have suggested that maps of environmental suitability estimated by ENM may be correlated to species' abundance, raising the possibility of using environmental suitability to infer processes related to population demographic dynamics and genetic variability. In both cases, one of the main problems is that there is a wide variation in ENM development methods and climatic models. In this study, we analyse the relationship between heterozygosity (He) and environmental suitability from multiple ENMs for 25 population estimates for Dipteryx alata, a widely distributed, endemic tree species of the Cerrado region of central Brazil. We propose a new approach for generating a statistical distribution of correlations under randomly generated ENM. The confidence intervals from these distributions indicate how model selection with different properties affects the ability to detect a correlation of interest (e.g. the correlation between He and suitability). Additionally, our approach allows us to explore which particular ensemble of ENMs produces the better result for finding an association between environmental suitability and He. Caution is necessary when choosing a method or a climatic data set for modelling geographical distributions, but the new approach proposed here provides a conservative way to evaluate the ability of ensembles to detect patterns of interest.

Impacts of selective logging on inbreeding and gene flow in two Amazonian timber species with contrasting ecological and reproductive characteristics.

Selective logging in Brazil allows for the removal of up to 90% of trees above 50 cm diameter of a given timber species, independent of a species' life history characteristics or how quickly it will recover. The genetic and demographic effects of selective logging on two Amazonian timber species (Dipteryx odorata Leguminosae, Jacaranda copaia Bignoniaceae) with contrasting ecological and reproductive characteristics were assessed in the same forest. Genetic diversity and gene flow were characterized by genotyping adults and seed sampled before and after logging, using hypervariable microsatellite markers. Overall, there were no short-term genetic impacts on the J. copaia population, with commercial application of current Brazilian forest management regulations. In contrast, for D. Odorata, selective logging showed a range of genetic impacts, with a 10% loss of alleles, and reductions in siring by pollen from trees within the 546-ha study area (23-11%) and in the number of pollen donors per progeny array (2.8-1.6), illustrating the importance of the surrounding landscape. Asynchrony in flowering between D. odorata trees led to trees with no breeding partners, which could limit the species reproduction and regeneration under current regulations. The results are summarized with other published studies from the same site and the implications for forest management discussed. The different types and levels of impacts associated with each species support the idea that ecological and genetic information by species, ecological guild or reproductive group is essential in helping to derive sustainable logging guidelines for tropical forests.

Long-term impacts of selective logging on two Amazonian tree species with contrasting ecological and reproductive characteristics: inferences from Eco-gene model simulations.

The impact of logging and subsequent recovery after logging is predicted to vary depending on specific life history traits of the logged species. The Eco-gene simulation model was used to evaluate the long-term impacts of selective logging over 300 years on two contrasting Brazilian Amazon tree species, Dipteryx odorata and Jacaranda copaia. D. odorata (Leguminosae), a slow growing climax tree, occurs at very low densities, whereas J. copaia (Bignoniaceae) is a fast growing pioneer tree that occurs at high densities. Microsatellite multilocus genotypes of the pre-logging populations were used as data inputs for the Eco-gene model and post-logging genetic data was used to verify the output from the simulations. Overall, under current Brazilian forest management regulations, there were neither short nor long-term impacts on J. copaia. By contrast, D. odorata cannot be sustainably logged under current regulations, a sustainable scenario was achieved by increasing the minimum cutting diameter at breast height from 50 to 100 cm over 30-year logging cycles. Genetic parameters were only slightly affected by selective logging, with reductions in the numbers of alleles and single genotypes. In the short term, the loss of alleles seen in J. copaia simulations was the same as in real data, whereas fewer alleles were lost in D. odorata simulations than in the field. The different impacts and periods of recovery for each species support the idea that ecological and genetic information are essential at species, ecological guild or reproductive group levels to help derive sustainable management scenarios for tropical forests.

A new eigenfunction spatial analysis describing population genetic structure.

Several methods of spatial analyses have been proposed to infer the relative importance of evolutionary processes on genetic population structure. Here we show how a new eigenfunction spatial analysis can be used to model spatial patterns in genetic data. Considering a sample of n local populations, the method starts by modeling the response variable (allele frequencies or phenotypic variation) against the eigenvectors sequentially extracted from a geographic distance matrix (n × n). The relationship between the coefficient of determination (R(2)) of the models and the cumulative eigenvalues, which we named the spatial signal-representation (SSR) curve, can be more efficient than Moran's I correlograms in describing different patterns. The SSR curve was also applied to simulated data (under distinct scenarios of population differentiation) and to analyze spatial patterns in alleles from microsatellite data for 25 local populations of Dipteryx alata, a tree species endemic to the Brazilian Cerrado. The SSR curves are consistent with previous phylogeographical patterns of the species, revealing combined effects of isolation-by-distance and range expansion. Our analyses demonstrate that the SSR curve is a useful exploratory tool for describing spatial patterns of genetic variability and for selecting spatial eigenvectors for models aiming to explain spatial responses to environmental variables and landscape features.

Demographic history and the low genetic diversity in Dipteryx alata (Fabaceae) from Brazilian Neotropical savannas.

Genetic effects of habitat fragmentation may be undetectable because they are generally a recent event in evolutionary time or because of confounding effects such as historical bottlenecks and historical changes in species' distribution. To assess the effects of demographic history on the genetic diversity and population structure in the Neotropical tree Dipteryx alata (Fabaceae), we used coalescence analyses coupled with ecological niche modeling to hindcast its distribution over the last 21 000 years. Twenty-five populations (644 individuals) were sampled and all individuals were genotyped using eight microsatellite loci. All populations presented low allelic richness and genetic diversity. The estimated effective population size was small in all populations and gene flow was negligible among most. We also found a significant signal of demographic reduction in most cases. Genetic differentiation among populations was significantly correlated with geographical distance. Allelic richness showed a spatial cline pattern in relation to the species' paleodistribution 21 kyr BP (thousand years before present), as expected under a range expansion model. Our results show strong evidences that genetic diversity in D. alata is the outcome of the historical changes in species distribution during the late Pleistocene. Because of this historically low effective population size and the low genetic diversity, recent fragmentation of the Cerrado biome may increase population differentiation, causing population decline and compromising long-term persistence.

Geographical patterns of turnover and nestedness-resultant components of allelic diversity among populations.

The analysis of geographical patterns in population divergence has always been a powerful way to infer microevolutionary processes involved in population differentiation, and several approaches have been used to investigate such patterns. Most frequently, multivariate spatial patterns of population differentiation are analyzed by computing pairwise genetic distances or F(ST) (or related statistics, such as ϕ(ST) from AMOVA), which are then correlated with geographical distances or landscape features. However, when calculating distances, especially based on presence-absence of alleles in local populations, there would be a confounding effect of allelic richness differences in the population differentiation. Moreover, the relative magnitude of these components and their spatial patterns can help identifying microevolutionary processes driving population differentiation. Here we show how recent methodological advances in ecological community analyses that allows partitioning dissimilarity into turnover (turnover) and richness differences, or nestedness-resultant dissimilarity, can be applied to allelic variation data, using an endemic Cerrado tree (Dipteryx alata) as a case study. Individuals from 15 local populations were genotyped for eight microsatellite loci, and pairwise dissimilarities were computed based on presence-absence of alleles. The turnover of alleles among populations represented 69 % of variation in dissimilarity, but only the richness difference component shows a clear spatial structure, appearing as a westward decrease of allelic richness. We show that decoupling richness difference and turnover components of allelic variation reveals more clearly how similarity among populations reflects geographical patterns in allelic diversity that can be interpreted in respect to historical range expansion in the species.

Geographic shifts in climatically suitable areas and loss of genetic variability in Dipteryx alata ("Baru" Tree; Fabaceae).

Many species are expected to suffer strong shifts in their geographic ranges due to climate changes in the next 50 years, with severe consequences for biodiversity patterns and population structure. We used here an ensemble forecast approach for obtaining species' range in which multiple species distribution models and climatic models were combined to model loss of genetic variability in Baru, Dipteryx alata (Fabaceae), an economically important Neotropical tree native to the Cerrado of Brazil. We estimated a series of genetic parameters (number of alleles per locus, expected heterozygosity under Hardy-Weinberg equilibrium and mutation-drift equilibrium) for this species based on eight microsatellite loci. We then recalculated these parameters assuming that local populations in areas of low future environmental suitability will go extinct. All genetic parameters remained approximately constant up to a 50% threshold of climatic suitability in the future; after this critical threshold there is an abrupt reduction in all parameters, although the magnitude of shift is only about 10% of current values, on average. Thus, despite the shifts in geographic range and climatically suitable areas towards southeastern Brazil, our analyses do not predict a strong loss of genetic diversity in D. alata because of the broad tolerance of this species, which ensures large future ranges, contrasting with other Cerrado species that have been analyzed in a similar manner.

Development of microsatellite markers for the neotropical tree species Dipteryx alata (Fabaceae).

PREMISE OF THE STUDY: Microsatellite markers were developed for the population genetic analyses of the neotropical tree Dipteryx alata (Fabaceae). METHODS AND RESULTS: Microsatellites were developed from a genomic shotgun library. Polymorphism at each microsatellite loci was analyzed based on 94 individuals from three populations. Eight loci amplified successfully and presented one to 10 alleles, and expected heterozygosities ranged from 0.097 to 0.862. Four loci also amplified in Pterodon emarginatus and presented similar polymorphism. CONCLUSION: The eight microsatellite primer pairs are potentially suitable for population genetic studies and successfully amplified in another Fabaceae species.

Pollen dispersal and genetic structure of the tropical tree Dipteryx panamensis in a fragmented Costa Rican landscape.

In the face of widespread deforestation, the conservation of rainforest trees relies increasingly on their ability to maintain reproductive processes in fragmented landscapes. Here, we analysed nine microsatellite loci for 218 adults and 325 progeny of the tree Dipteryx panamensis in Costa Rica. Pollen dispersal distances, genetic diversity, genetic structure and spatial autocorrelation were determined for populations in four habitats: continuous forest, forest fragments, pastures adjacent to fragments and isolated pastures. We predicted longer but less frequent pollen movements among increasingly isolated trees. This pattern would lead to lower outcrossing rates for pasture trees, as well as lower genetic diversity and increased structure and spatial autocorrelation among their progeny. Results generally followed these expectations, with the shortest pollen dispersal among continuous forest trees (240 m), moderate distances for fragment (343 m) and adjacent pasture (317 m) populations, and distances of up to 2.3 km in isolated pastures (mean: 557 m). Variance around pollen dispersal estimates also increased with fragmentation, suggesting altered pollination conditions. Outcrossing rates were lower for pasture trees and we found greater spatial autocorrelation and genetic structure among their progeny, as well as a trend towards lower heterozygosity. Paternal reproductive dominance, the pollen contributions from individual fathers, did not vary among habitats, but we did document asymmetric pollen flow between pasture and adjacent fragment populations. We conclude that long-distance pollen dispersal helps maintain gene flow for D. panamensis in this fragmented landscape, but pasture and isolated pasture populations are still at risk of long-term genetic erosion.

Landscape conservation genetics of Dipteryx alata ("baru" tree: Fabaceae) from Cerrado region of central Brazil.

In this paper random amplified polymorphic DNA (RAPD) was used to evaluate the degree of among-population differentiation and associated spatial patterns of genetic divergence for Dipteryx alata Vogel populations from Cerrado region of central Brazil, furnishing support for future programs of conservation of this species. We analyzed patterns of genetic and spatial population structure using 45 RAPD loci scored for 309 trees, sampled from five different regions with two populations each. Genetic structure analysis suggested that panmixia null hypothesis can be rejected, with significant among-population components of 15%. Hierarchical partition by Analysis of Molecular Variance (AMOVA) shows that 5% of genetic variation is within regions, whereas 10% of variation is among regions, and these results were confirmed by a Bayesian analyses on HICKORY. The Mantel correlogram revealed that this divergence is spatially structured, so that local populations situated at short geographic distances could not be considered independent units for conservation and management. However, genetic discontinuities among populations were found in the northwest and southeast parts of the study area, corresponding to regions of recent socio-economic expansion and high population density, respectively. Taking both geographic distances and genetic discontinuities into account it is possible to establish a group of population to be conserved, covering most of D. alata geographic distribution and congruent with previously established priority areas for conservation in the Cerrado region.