Population genomics of the widespread African savannah trees Afzelia africana and Afzelia quanzensis reveals no significant past fragmentation of their distribution ranges.

PREMISE: Few studies have addressed the evolutionary history of tree species from African savannahs. Afzelia contains economically important timber species, including two species widely distributed in African savannahs: A. africana in the Sudanian region and A. quanzensis in the Zambezian region. We aimed to infer whether these species underwent range fragmentation and/or demographic changes, possibly reflecting how savannahs responded to Quaternary climate changes. METHODS: We characterized the genetic diversity and structure of these species across their distribution ranges using nuclear microsatellites (SSRs) and genotyping-by-sequencing (GBS) markers. Six SSR loci were genotyped in 241 A. africana and 113 A. quanzensis individuals, while 2800 high-quality single nucleotide polymorphisms (SNPs) were identified in 30 A. africana individuals. RESULTS: Both species appeared to be mainly outcrossing. The kinship between individuals decayed with the logarithm of the distance at similar rates across species and markers, leading to relatively small Sp statistics (0.0056 for SSR and 0.0054 for SNP in A. africana, 0.0075 for SSR in A. quanzensis). The patterns were consistent with isolation by distance expectations in the absence of large-scale geographic gradients. Bayesian clustering of SSR genotypes did not detect genetic clusters within species. In contrast, SNP data resolved intraspecific genetic clusters in A. africana, illustrating the higher resolving power of GBS. However, these clusters revealed low levels of differentiation and no clear geographical entities, so that they were interpreted as resulting from the isolation by distance pattern rather than from past population fragmentation. CONCLUSIONS: These results suggest that populations have remained connected throughout the large, continuous savannah landscapes. The absence of clear phylogeographic discontinuities, also found in a few other African savannah trees, indicates that their distribution ranges have not been significantly fragmented during the climatic oscillations of the Pleistocene, in contrast to patterns commonly found in African rainforest trees.

Highlighting convergent evolution in morphological traits in response to climatic gradient in African tropical tree species: The case of genus Guibourtia Benn.

Adaptive evolution is a major driver of organism diversification, but the links between phenotypic traits and environmental niche remain little documented in tropical trees. Moreover, trait-niche relationships are complex because a correlation between the traits and environmental niches displayed by a sample of species may result from (a) convergent evolution if different environmental conditions have selected different sets of traits, and/or (b) phylogenetic inertia if niche and morphological differences between species are simply function of their phylogenetic divergence, in which case the trait-niche correlation does not imply any direct causal link. Here, we aim to assess the respective roles of phylogenetic inertia and convergent evolution in shaping the differences of botanical traits and environmental niches among congeneric African tree species that evolved in different biomes.This issue was addressed with the tree genus Guibourtia Benn. (Leguminosae and Detarioideae), which contains 13 African species occupying various forest habitat types, from rain forest to dry woodlands, with different climate and soil conditions. To this end, we combined morphological data with ecological niche modelling and used a highly resolved plastid phylogeny of the 13 African Guibourtia species.First, we demonstrated phylogenetic signals in both morphological traits (Mantel test between phylogenetic and morphological distances between species: r = .24, p = .031) and environmental niches (Mantel test between phylogenetic and niche distances between species: r = .23, p = .025). Second, we found a significant correlation between morphology and niche, at least between some of their respective dimensions (Mantel's r = .32, p = .013), even after accounting for phylogenetic inertia (Phylogenetic Independent Contrast: r = .69, p = .018). This correlation occurred between some leaflet and flower traits and solar radiation, relative humidity, precipitations, and temperature range.Our results demonstrate the convergent evolution of some morphological traits in response to climatic factors in congeneric tree species and highlight the action of selective forces, along with neutral ones, in shaping the divergence between congeneric tropical plants.

A few north Appalachian populations are the source of European black locust.

The role of evolution in biological invasion studies is often overlooked. In order to evaluate the evolutionary mechanisms behind invasiveness, it is crucial to identify the source populations of the introduction. Studies in population genetics were carried out on Robinia pseudoacacia L., a North American tree which is now one of the worst invasive tree species in Europe. We realized large-scale sampling in both the invasive and native ranges: 63 populations were sampled and 818 individuals were genotyped using 113 SNPs. We identified clonal genotypes in each population and analyzed between and within range population structure, and then, we compared genetic diversity between ranges, enlarging the number of SNPs to mitigate the ascertainment bias. First, we demonstrated that European black locust was introduced from just a limited number of populations located in the Appalachian Mountains, which is in agreement with the historical documents briefly reviewed in this study. Within America, population structure reflected the effects of long-term processes, whereas in Europe it was largely impacted by human activities. Second, we showed that there is a genetic bottleneck between the ranges with a decrease in allelic richness and total number of alleles in Europe. Lastly, we found more clonality within European populations. Black locust became invasive in Europe despite being introduced from a reduced part of its native distribution. Our results suggest that human activity, such as breeding programs in Europe and the seed trade throughout the introduced range, had a major role in promoting invasion; therefore, the introduction of the missing American genetic cluster to Europe should be avoided.

Evolution in the Amphi-Atlantic tropical genus Guibourtia (Fabaceae, Detarioideae), combining NGS phylogeny and morphology.

Tropical rain forests support a remarkable diversity of tree species, questioning how and when this diversity arose. The genus Guibourtia (Fabaceae, Detarioideae), characterized by two South American and 13 African tree species growing in various tropical biomes, is an interesting model to address the role of biogeographic processes and adaptation to contrasted environments on species diversification. Combining whole plastid genome sequencing and morphological characters analysis, we studied the timing of speciation and diversification processes in Guibourtia through molecular dating and ancestral habitats reconstruction. All species except G. demeusei and G. copallifera appear monophyletic. Dispersal from Africa to America across the Atlantic Ocean is the most plausible hypothesis to explain the occurrence of Neotropical Guibourtia species, which diverged ca. 11.8 Ma from their closest African relatives. The diversification of the three main clades of African Guibourtia is concomitant to Miocene global climate changes, highlighting pre-Quaternary speciation events. These clades differ by their reproductive characters, which validates the three subgenera previously described: Pseudocopaiva, Guibourtia and Gorskia. Within most monophyletic species, plastid lineages start diverging from each other during the Pliocene or early Pleistocene, suggesting that these species already arose during this period. The multiple transitions between rain forests and dry forests/savannahs inferred here through the plastid phylogeny in each Guibourtia subgenus address thus new questions about the role of phylogenetic relationships in shaping ecological niche and morphological similarity among taxa.

Characterization of microsatellite markers in the African tropical tree species Guibourtia ehie (Fabaceae, Detarioideae).

PREMISE OF THE STUDY: Microsatellite primers (simple sequence repeats [SSRs]) were developed in Guibourtia ehie (Fabaceae, Detarioideae) to study population genetic structure and the history of African vegetation. METHODS AND RESULTS: We isolated 18 polymorphic SSRs from a nonenriched genomic library. This set of primer pairs was tested on four populations, and the results showed two to 16 alleles per locus with mean observed and expected heterozygosities of 0.27 ± 0.05 and 0.57 ± 0.05, respectively. Cross-amplification tests in 13 congeneric species were successful for the four taxa belonging to the subgenus Gorskia. CONCLUSIONS: This set of microsatellite markers will be useful to investigate the phylogeography and population genetics of G. ehie, a key representative of African semideciduous moist forests.

Evolution in African tropical trees displaying ploidy-habitat association: The genus Afzelia (Leguminosae).

Polyploidy has rarely been documented in rain forest trees but it has recently been found in African species of the genus Afzelia (Leguminosae), which is composed of four tetraploid rain forest species and two diploid dry forest species. The genus Afzelia thus provides an opportunity to examine how and when polyploidy and habitat shift occurred in Africa, and whether they are associated. In this study, we combined three plastid markers (psbA, trnL, ndhF), two nuclear markers (ribosomal ITS and the single-copy PEPC E7 gene), plastomes (obtained by High Throughput Sequencing) and morphological traits, with an extensive taxonomic and geographic sampling to explore the evolutionary history of Afzelia. Both nuclear DNA and morphological vegetative characters separated diploid from tetraploid lineages. Although the two African diploid species were well differentiated genetically and morphologically, the relationships among the tetraploid species were not resolved. In contrast to the nuclear markers, plastid markers revealed that one of the diploid species forms a well-supported clade with the tetraploids, suggesting historical hybridisation, possibly in relation with genome duplication (polyploidization) and habitat shift from dry to rain forests. Molecular dating based on fossil-anchored gene phylogenies indicates that extant Afzelia started diverging c. 14.5 or 20Ma while extant tetraploid species started diverging c. 7.0 or 9.4Ma according to plastid and nuclear DNA, respectively. Additional studies of tropical polyploid plants are needed to assess whether the ploidy-habitat association observed in African Afzelia would reflect a role of polyploidization in niche divergence in the tropics.

Revealing hidden species diversity in closely related species using nuclear SNPs, SSRs and DNA sequences - a case study in the tree genus Milicia.

BACKGROUND: Species delimitation in closely related plant taxa can be challenging because (i) reproductive barriers are not always congruent with morphological differentiation, (ii) use of plastid sequences might lead to misinterpretation, (iii) rare species might not be sampled. We revisited molecular-based species delimitation in the African genus Milicia, currently divided into M. regia (West Africa) and M. excelsa (from West to East Africa). We used 435 samples collected in West, Central and East Africa. We genotyped SNP and SSR loci to identify genetic clusters, and sequenced two plastid regions (psbA-trnH, trnC-ycf6) and a nuclear gene (At103) to confirm species' divergence and compare species delimitation methods. We also examined whether ecological niche differentiation was congruent with sampled genetic structure. RESULTS: West African M. regia, West African and East African M. excelsa samples constituted three well distinct genetic clusters according to SNPs and SSRs. In Central Africa, two genetic clusters were consistently inferred by both types of markers, while a few scattered samples, sympatric with the preceding clusters but exhibiting leaf traits of M. regia, were grouped with the West African M. regia cluster based on SNPs or formed a distinct cluster based on SSRs. SSR results were confirmed by sequence data from the nuclear region At103 which revealed three distinct 'Fields For Recombination' corresponding to (i) West African M. regia, (ii) Central African samples with leaf traits of M. regia, and (iii) all M. excelsa samples. None of the plastid sequences provide indication of distinct clades of the three species-like units. Niche modelling techniques yielded a significant correlation between niche overlap and genetic distance. CONCLUSIONS: Our genetic data suggest that three species of Milicia could be recognized. It is surprising that the occurrence of two species in Central Africa was not reported for this well-known timber tree. Globally, our work highlights the importance of collecting samples in a systematic way and the need for combining different nuclear markers when dealing with species complexes. Recognizing cryptic species is particularly crucial for economically exploited species because some hidden taxa might actually be endangered as they are merged with more abundant species.

Microsatellite development for the genus Guibourtia (Fabaceae, Caesalpinioideae) reveals diploid and polyploid species.

PREMISE OF THE STUDY: Nuclear microsatellites (nSSRs) were designed for Guibourtia tessmannii (Fabaceae, Caesalpinioideae), a highly exploited African timber tree, to study population genetic structure and gene flow. METHODS AND RESULTS: We developed 16 polymorphic nSSRs from a genomic library tested in three populations of G. tessmannii and two populations of G. coleosperma. These nSSRs display three to 14 alleles per locus (mean 8.94) in G. tessmannii. Cross-amplification tests in nine congeneric species demonstrated that the genus Guibourtia contains diploid and polyploid species. Flow cytometry results combined with nSSR profiles suggest that G. tessmannii is octoploid. CONCLUSIONS: nSSRs revealed that African Guibourtia species include both diploid and polyploid species. These markers will provide information on the mating system, patterns of gene flow, and genetic structure of African Guibourtia species.

Microsatellite development and flow cytometry in the African tree genus Afzelia (Fabaceae, Caesalpinioideae) reveal a polyploid complex.

PREMISE OF THE STUDY: Microsatellites were developed in the vulnerable African rainforest tree Afzelia bipindensis to investigate gene flow patterns. • METHODS AND RESULTS: Using 454 GS-FLX technique, 16 primer sets were identified and optimized, leading to 11 polymorphic and readable markers displaying each six to 25 alleles in a population. Up to four alleles per individual were found in each of the loci, without evidence of fixed heterozygosity, suggesting an autotetraploid genome. Cross-amplification succeeded for all loci in the African rainforest species A. pachyloba and A. bella, which appeared tetraploid, and for most loci in the African woodland species A. africana and A. quanzensis, which appeared diploid, but failed in the Asian species A. xylocarpa. Flow cytometry confirmed the suspected differences in ploidy. • CONCLUSIONS: African Afzelia species are diploid or tetraploid, a situation rarely documented in tropical trees. These newly developed microsatellites will help in the study of their mating system and gene flow patterns.

Soil charcoal to assess the impacts of past human disturbances on tropical forests.

The canopy of many central African forests is dominated by light-demanding tree species that do not regenerate well under themselves. The prevalence of these species might result from ancient slash-and-burn agricultural activities that created large openings, while a decline of these activities since the colonial period could explain their deficit of regeneration. To verify this hypothesis, we compared soil charcoal abundance, used as a proxy for past slash-and-burn agriculture, and tree species composition assessed on 208 rainforest 0.2 ha plots located in three areas from Southern Cameroon. Species were classified in regeneration guilds (pioneer, non-pioneer light-demanding, shade-bearer) and characterized by their wood-specific gravity, assumed to reflect light requirement. We tested the correlation between soil charcoal abundance and: (i) the relative abundance of each guild, (ii) each species and family abundance and (iii) mean wood-specific gravity. Charcoal was found in 83% of the plots, indicating frequent past forest fires. Radiocarbon dating revealed two periods of fires: "recent" charcoal were on average 300 years old (up to 860 BP, n = 16) and occurred in the uppermost 20 cm soil layer, while "ancient" charcoal were on average 1900 years old (range: 1500 to 2800 BP, n = 43, excluding one sample dated 9400 BP), and found in all soil layers. While we expected a positive correlation between the relative abundance of light-demanding species and charcoal abundance in the upper soil layer, overall there was no evidence that the current heterogeneity in tree species composition can be explained by charcoal abundance in any soil layer. The absence of signal supporting our hypothesis might result from (i) a relatively uniform impact of past slash-and-burn activities, (ii) pedoturbation processes bringing ancient charcoal to the upper soil layer, blurring the signal of centuries-old Human disturbances, or (iii) the prevalence of other environmental factors on species composition.

Phenological patterns in a natural population of a tropical timber tree species, Milicia excelsa (Moraceae): Evidence of isolation by time and its interaction with feeding strategies of dispersers.

PREMISE OF THE STUDY: Population genetic structuring over limited timescales is commonly viewed as a consequence of spatial constraints. Indirect approaches have recently revealed reproductive isolation resulting from flowering time (so-called isolation by time, IBT). Since phenological processes can be subject to selection, the persistence of flowering asynchrony may be due to opposing selective pressures during mating, dispersal, and regeneration phases. Our study aimed to investigate phenology, fruit handling by animals, and their interaction in a timber tree species, Milicia excelsa. METHODS: We analyzed phenological data collected over 6 years on 69 genotyped trees in a Cameroonian natural rainforest complemented by data from germination trials and field observations of dispersers. KEY RESULTS: Initiation of flowering was correlated with variation in temperature and relative humidity, but was also affected by genetic factors: pairwise differences in flowering time between nearby individuals correlated with kinship coefficient, and earliness of flowering remained stable over time. A decrease in mean seed production per fruit with increasing flowering time suggests selection against late bloomers. However, germination rate was not affected by seed collection date, and the main seed disperser, the bat Eidolon helvum, seemed to increase in abundance at the end of the reproductive season and preferred trees in open habitats where early and late bloomers are expected. CONCLUSIONS: The pairwise approach performs well in detecting IBT. The persistence of different mating pools in such a case may result from a trade off between selective forces during the mating and seed dispersal processes.

Forest refugia revisited: nSSRs and cpDNA sequences support historical isolation in a wide-spread African tree with high colonization capacity, Milicia excelsa (Moraceae).

The impact of the Pleistocene climate oscillations on the structure of biodiversity in tropical regions remains poorly understood. In this study, the forest refuge theory is examined at the molecular level in Milicia excelsa, a dioecious tree with a continuous range throughout tropical Africa. Eight nuclear microsatellites (nSSRs) and two sequences and one microsatellite from chloroplast DNA (cpDNA) showed a deep divide between samples from Benin and those from Lower Guinea. This suggests that these populations were isolated in separate geographical regions, probably for several glacial cycles of the Pleistocene, and that the nuclear gene pools were not homogenized despite M. excelsa's wind-pollination syndrome. The divide could also be related to seed dispersal patterns, which should be largely determined by the migration behaviour of M. excelsa's main seed disperser, the frugivorous bat Eidolon helvum. Within Lower Guinea, a north-south divide, observed with both marker types despite weak genetic structure (nSSRs: F(ST) = 0.035, cpDNA: G(ST) = 0.506), suggested the existence of separate Pleistocene refugia in Cameroon and the Gabon/Congo region. We inferred a pollen-to-seed dispersal distance ratio of c.1.8, consistent with wide-ranging gene dispersal by both wind and bats. Simulations in an Approximate Bayesian Computation framework suggested low nSSR and cpDNA mutation rates, but imprecise estimates of other demographic parameters, probably due to a substantial gene flow between the Lower Guinean gene pools. The decline of genetic diversity detected in some Gabonese populations could be a consequence of the relatively recent establishment of a closed canopy forest, which could negatively affect M. excelsa's reproductive system.