The ancient tropical rainforest tree Symphonia globulifera L. f. (Clusiaceae) was not restricted to postulated Pleistocene refugia in Atlantic Equatorial Africa

Topography drives microgeographic adaptations of closely related species in two tropical tree species complexes

Closely related tree species that grow in sympatry are abundant in rainforests. However, little is known of the eco-evolutionary processes that govern their niches and local coexistence. We assessed genetic species delimitation in closely related sympatric species belonging to two Neotropical tree species complexes and investigated their genomic adaptation to a fine-scale topographic gradient with associated edaphic and hydrologic features. Combining LiDAR-derived topography, tree inventories, and single nucleotide polymorphisms (SNPs) from gene capture experiments, we explored genome-wide population genetic structure, covariation of environmental variables, and genotype-environment association to assess microgeographic adaptations to topography within the species complexes Symphonia (Clusiaceae), and Eschweilera (Lecythidaceae) with three species per complex and 385 and 257 individuals genotyped, respectively. Within species complexes, closely related tree species had different realized optima for topographic niches defined through the topographic wetness index or the relative elevation, and species displayed genetic signatures of adaptations to these niches. Symphonia species were genetically differentiated along water and nutrient distribution particularly in genes responding to water deprivation, whereas Eschweilera species were genetically differentiated according to soil chemistry. Our results suggest that varied topography represents a powerful driver of processes modulating tropical forest biodiversity with differential adaptations that stabilize local coexistence of closely related tree species.

WP, Sterkhova VV, Kusamba C, Rödel M, Penner J, Peterson AT, Brown RM. Rivers, not refugia, drove diversification in arboreal, sub-Saharan African snakes

The relative roles of rivers versus refugia in shaping the high levels of species diversity in tropical rainforests have been widely debated for decades. Only recently has it become possible to take an integrative approach to test predictions derived from these hypotheses using genomic sequencing and paleo-species distribution modeling. Herein, we tested the predictions of the classic river, refuge, and river-refuge hypotheses on diversification in the arboreal sub-Saharan African snake genus Toxicodryas. We used dated phylogeographic inferences, population clustering analyses, demographic model selection, and paleo-distribution modeling to conduct a phylogenomic and historical demographic analysis of this genus. Our results revealed significant population genetic structure within both Toxicodryas species, corresponding geographically to river barriers and divergence times from the mid-Miocene to Pliocene. Our demographic analyses supported the interpretation that rivers are indications of strong barriers to gene flow among populations since their divergence. Additionally, we found no support for a major contraction of suitable habitat during the last glacial maximum, allowing us to reject both the refuge and river-refuge hypotheses in favor of the river-barrier hypothesis. Based on conservative interpretations of our species delimitation analyses with the Sanger and ddRAD data sets, two new cryptic species are identified from east-central Africa. This study highlights the complexity of diversification dynamics in the African tropics and the advantages of integrative approaches to studying speciation in tropical regions.

Plastome phylogeography in two African rain forest legume trees reveals that Dahomey Gap populations originate from the Cameroon volcanic line

Paleo-environmental data show that the distribution of African rain forests was affected by Quaternary climate changes. In particular, the Dahomey Gap (DG) - a 200 km wide savanna corridor currently separating the West African and Central African rain forest blocks and containing relict rain forest fragments - was forested during the mid-Holocene and possibly during previous interglacial periods, whereas it was dominated by open vegetation (savanna) during glacial periods. Genetic signatures of past population fragmentation and demographic changes have been found in some African forest plant species using nuclear markers, but such events appear not to have been synchronous or shared across species. To better understand the colonization history of the DG by rain forest trees through seed dispersal, the plastid genomes of two widespread African forest legume trees, Anthonotha macrophylla and Distemonanthus benthamianus, were sequenced in 47 individuals for each species, providing unprecedented phylogenetic resolution of their maternal lineages (857 and 115 SNPs, respectively). Both species exhibit distinct lineages separating three regions: 1. Upper Guinea (UG, i.e. the West African forest block), 2. the area ranging from the DG to the Cameroon volcanic line (CVL), and 3. Lower Guinea (LG, the western part of the Central African forest block) where three lineages co-occur. In both species, the DG populations (including southern Nigeria west of Cross River) exhibit much lower genetic diversity than UG and LG populations, and their plastid lineages originate from the CVL, confirming the role of the CVL as an ancient forest refuge. Despite the similar phylogeographic structures displayed by A. macrophylla and D. benthamianus, molecular dating indicates very contrasting ages of lineage divergence (UG diverged from LG since c. 7 Ma and 0.7 Ma, respectively) and DG colonization (probably following the Mid Pleistocene Transition and the Last Glacial Maximum, respectively). The stability of forest refuge areas and repeated similar forest shrinking/expanding events during successive glacial periods might explain why similar phylogeographic patterns can be generated over contrasting timescales.

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.

A geographical traceability system for Merbau (Intsia palembanica Miq.), an important timber species from peninsular Malaysia

To inform product users about the origin of timber, the implementation of a traceability system is necessary for the forestry industry. In this study, we developed a comprehensive genetic database for the important tropical timber species Merbau, Intsia palembanica, to trace its geographic origin within peninsular Malaysia. A total of 1373 individual trees representing 39 geographically distinct populations of I. palembanica were sampled throughout peninsular Malaysia. We analyzed the samples using a combination of four chloroplast DNA (cpDNA) markers and 14 short tandem repeat (STR) markers to establish both cpDNA haplotype and STR allele frequency databases. A haplotype map was generated through cpDNA sequencing for population identification, resulting in six unique haplotypes based on 10 informative intraspecifically variable sites. Subsequently, an STR allele frequency database was developed from 14 STRs allowing individual identification. Bayesian cluster analysis divided the individuals into two genetic clusters corresponding to the northern and southern regions of peninsular Malaysia. Tests of conservativeness showed that the databases were conservative after the adjustment of the θ values to 0.2000 and 0.2900 for the northern (f = 0.0163) and southern (f = 0.0285) regions, respectively. Using self-assignment tests, we observed that individuals were correctly assigned to populations at rates of 40.54-94.12% and to the identified regions at rates of 79.80-80.62%. Both the cpDNA and STR markers appear to be useful for tracking Merbau timber originating from peninsular Malaysia. The use of these forensic tools in addition to the existing paper-based timber tracking system will help to verify the legality of the origin of I. palembanica and to combat illegal logging issues associated with the species.

A tale of two forests: ongoing aridification drives population decline and genetic diversity loss at continental scale in Afro-Macaronesian evergreen-forest archipelago endemics

BACKGROUND AND AIMS

Various studies and conservationist reports have warned about the contraction of the last subtropical Afro-Macaronesian forests. These relict vegetation zones have been restricted to a few oceanic and continental islands around the edges of Africa, due to aridification. Previous studies on relict species have generally focused on glacial effects on narrow endemics; however, little is known about the effects of aridification on the fates of previously widespread subtropical lineages.

METHODS

Nuclear microsatellites and ecological niche modelling were used to understand observed patterns of genetic diversity in two emblematic species, widely distributed in these ecosystems: Canarina eminii (a palaeoendemic of the eastern Afromontane forests) and Canarina canariensis (a palaeoendemic of the Canarian laurel forests). The software DIYABC was used to test alternative demographic scenarios and an ensemble method was employed to model potential distributions of the selected plants from the end of the deglaciation to the present.

KEY RESULTS

All the populations assessed experienced a strong and recent population decline, revealing that locally widespread endemisms may also be alarmingly threatened.

CONCLUSIONS

The detected extinction debt, as well as the extinction spiral to which these populations are subjected, demands urgent conservation measures for the unique, biodiversity-rich ecosystems that they inhabit.

Refining Species Traits in a Dynamic Vegetation Model to Project the Impacts of Climate Change on Tropical Trees in Central Africa

African tropical ecosystems and the services they provide to human society suffer from an increasing combined pressure of land use and climate change. How individual tropical tree species respond to climate change remains relatively unknown. In this study, we refined the species characterization in the CARAIB (CARbon Assimilation In the Biosphere) dynamic vegetation model by replacing plant functional type morpho-physiological traits by species-specific traits. We focus on 12 tropical tree species selected for their importance in both the plant community and human society. We used CARAIB to simulate the current species net primary productivity (NPP), biomass and potential distribution and their changes in the future. Our results indicate that the use of species-specific traits does not necessarily result in an increase of predicted current NPPs. The model projections for the end of the century highlight the large uncertainties in the future of African tropical species. Projected changes in species distribution vary greatly with the general circulation model (GCM) and, to a lesser extent, with the concentration pathway. The question about long-term plant response to increasing CO2 concentrations also leads to contrasting results. In absence of fertilization effect, species are exposed to climate change and might lose 25% of their current distribution under RCP8.5 (12.5% under RCP4.5), considering all the species and climatic scenarios. The vegetation model projects a mean biomass loss of −21.2% under RCP4.5 and −34.5% under RCP8.5. Potential range expansions, unpredictable due to migration limitations, are too limited for offsetting range contraction. By contrast, if the long-term species response to increasing [CO2] is positive, the range reduction is limited to 5%. However, despite a mean biomass increase of 12.2%, a positive CO2 feedback might not prevent tree dieback. Our analysis confirms that species will respond differently to new climatic and atmospheric conditions, which may induce new competition dynamics in the ecosystem and affect ecosystem services.

Approximate Bayesian computation analysis of EST-associated microsatellites indicates that the broadleaved evergreen tree Castanopsis sieboldii survived the Last Glacial Maximum in multiple refugia in Japan

Climatic changes have played major roles in plants' evolutionary history. Glacial oscillations have been particularly important, but some of their effects on plants' populations are poorly understood, including the numbers and locations of refugia in Asian warm temperate zones. In the present study, we investigated the demographic history of the broadleaved evergreen tree species Castanopsis sieboldii (Fagaceae) during the last glacial period in Japan. We used approximate Bayesian computation (ABC) for model comparison and parameter estimation for the demographic modeling using 27 EST-associated microsatellites. We also performed the species distribution modeling (SDM). The results strongly support a demographic scenario that the Ryukyu Islands and the western parts in the main islands (Kyushu and western Shikoku) were derived from separate refugia and the eastern parts in the main islands and the Japan Sea groups were diverged from the western parts prior to the coldest stage of the Last Glacial Maximum (LGM). Our data indicate that multiple refugia survived at least one in the Ryukyu Islands, and the other three regions of the western and eastern parts and around the Japan Sea of the main islands of Japan during the LGM. The SDM analysis also suggests the potential habitats under LGM climate conditions were mainly located along the Pacific Ocean side of the coastal region. Our ABC-based study helps efforts resolve the demographic history of a dominant species in warm temperate broadleaved forests during and after the last glacial period, which provides a basic model for future phylogeographical studies using this approach.

History of the fragmentation of the African rain forest in the Dahomey Gap: insight from the demographic history of Terminalia superba

Paleo-environmental reconstructions show that the distribution of tropical African rain forests was affected by Quaternary climate changes. They suggest that the Dahomey Gap (DG)-the savanna corridor that currently separates Upper Guinean (UG, West Africa) and Lower Guinean (LG, western Central Africa) rain forest blocks-was forested during the African Humid Holocene period (from at least 9 ka till 4.5 ka), and possibly during other interglacial periods, while an open vegetation developed in the DG under drier conditions, notably during glacial maxima. Nowadays, relics of semi-deciduous forests containing UG and LG forest species are still present within the DG. We used one of these species, the pioneer tree Terminalia superba (Combretaceae), to study past forest fragmentation in the DG and its impact on infraspecific biodiversity. A Bayesian clustering analysis of 299 individuals genotyped at 14 nuclear microsatellites revealed five parapatric genetic clusters (UG, DG, and three in LG) with low to moderate genetic differentiation (Fst from 0.02 to 0.24). Approximate Bayesian Computation analyses inferred a demographic bottleneck around the penultimate glacial period in all populations. They also supported an origin of the DG population by admixture of UG and LG populations around 54,000 (27,600-161,000) years BP, thus before the Last Glacial Maximum. These results contrast with those obtained on Distemonanthus benthamianus where the DG population seems to originate from the Humid Holocene period. We discuss these differences in light of the ecology of each species. Our results challenge the simplistic view linking population fragmentation/expansion with glacial/interglacial periods in African forest species.

Altitudinal gradients, biogeographic history and microhabitat adaptation affect fine-scale spatial genetic structure in African and Neotropical populations of an ancient tropical tree species

The analysis of fine-scale spatial genetic structure (FSGS) within populations can provide insights into eco-evolutionary processes. Restricted dispersal and locally occurring genetic drift are the primary causes for FSGS at equilibrium, as described in the isolation by distance (IBD) model. Beyond IBD expectations, spatial, environmental or historical factors can affect FSGS. We examined FSGS in seven African and Neotropical populations of the late-successional rain forest tree Symphonia globulifera L. f. (Clusiaceae) to discriminate the influence of drift-dispersal vs. landscape/ecological features and historical processes on FSGS. We used spatial principal component analysis and Bayesian clustering to assess spatial genetic heterogeneity at SSRs and examined its association with plastid DNA and habitat features. African populations (from Cameroon and São Tomé) displayed a stronger FSGS than Neotropical populations at both marker types (mean Sp = 0.025 vs. Sp = 0.008 at SSRs) and had a stronger spatial genetic heterogeneity. All three African populations occurred in pronounced altitudinal gradients, possibly restricting animal-mediated seed dispersal. Cyto-nuclear disequilibria in Cameroonian populations also suggested a legacy of biogeographic history to explain these genetic patterns. Conversely, Neotropical populations exhibited a weaker FSGS, which may reflect more efficient wide-ranging seed dispersal by Neotropical bats and other dispersers. The population from French Guiana displayed an association of plastid haplotypes with two morphotypes characterized by differential habitat preferences. Our results highlight the importance of the microenvironment for eco-evolutionary processes within persistent tropical tree populations.

Pleistocene refugia and genetic diversity patterns in West Africa: Insights from the liana Chasmanthera dependens (Menispermaceae)

Processes shaping the African Guineo-Congolian rain forest, especially in the West African part, are not well understood. Recent molecular studies, based mainly on forest tree species, confirmed the previously proposed division of the western African Guineo-Congolian rain forest into Upper Guinea (UG) and Lower Guinea (LG) separated by the Dahomey Gap (DG). Here we studied nine populations in the area of the DG and the borders of LG and UG of the widespread liana species, Chasmanthera dependens (Menispermaceae) by amplified fragment length polymorphism (AFLP), a chloroplast DNA sequence marker, and modelled the distribution based on current as well as paleoclimatic data (Holocene Climate Optimum, ca. 6 kyr BP and Last Glacial Maximum, ca. 22 kyr BP). Current population genetic structure and geographical pattern of cpDNA was related to present as well as historical modelled distributions. Results from this study show that past historical factors played an important role in shaping the distribution of C. dependens across West Africa. The Cameroon Volcanic Line seems to represent a barrier for gene flow in the present as well as in the past. Distribution modelling proposed refugia in the Dahomey Gap, supported also by higher genetic diversity. This is in contrast with the phylogeographic patterns observed in several rainforest tree species and could be explained by either diverging or more relaxed ecological requirements of this liana species.

Tropical rainforests that persisted: inferences from the Quaternary demographic history of eight tree species in the Guiana shield

How Quaternary climatic and geological disturbances influenced the composition of Neotropical forests is hotly debated. Rainfall and temperature changes during and/or immediately after the last glacial maximum (LGM) are thought to have strongly affected the geographical distribution and local abundance of tree species. The paucity of the fossil records in Neotropical forests prevents a direct reconstruction of such processes. To describe community-level historical trends in forest composition, we turned therefore to inferential methods based on the reconstruction of past demographic changes. In particular, we modelled the history of rainforests in the eastern Guiana Shield over a timescale of several thousand generations, through the application of approximate Bayesian computation and maximum-likelihood methods to diversity data at nuclear and chloroplast loci in eight species or subspecies of rainforest trees. Depending on the species and on the method applied, we detected population contraction, expansion or stability, with a general trend in favour of stability or expansion, with changes presumably having occurred during or after the LGM. These findings suggest that Guiana Shield rainforests have globally persisted, while expanding, through the Quaternary, but that different species have experienced different demographic events, with a trend towards the increase in frequency of light-demanding, disturbance-associated species.

Development of genomic tools in a widespread tropical tree, Symphonia globulifera L.f.: a new low-coverage draft genome, SNP and SSR markers

Population genetic studies in tropical plants are often challenging because of limited information on taxonomy, phylogenetic relationships and distribution ranges, scarce genomic information and logistic challenges in sampling. We describe a strategy to develop robust and widely applicable genetic markers based on a modest development of genomic resources in the ancient tropical tree species Symphonia globulifera L.f. (Clusiaceae), a keystone species in African and Neotropical rainforests. We provide the first low-coverage (11X) fragmented draft genome sequenced on an individual from Cameroon, covering 1.027 Gbp or 67.5% of the estimated genome size. Annotation of 565 scaffolds (7.57 Mbp) resulted in the prediction of 1046 putative genes (231 of them containing a complete open reading frame) and 1523 exact simple sequence repeats (SSRs, microsatellites). Aligning a published transcriptome of a French Guiana population against this draft genome produced 923 high-quality single nucleotide polymorphisms. We also preselected genic SSRs in silico that were conserved and polymorphic across a wide geographical range, thus reducing marker development tests on rare DNA samples. Of 23 SSRs tested, 19 amplified and 18 were successfully genotyped in four S. globulifera populations from South America (Brazil and French Guiana) and Africa (Cameroon and São Tomé island, F_ST  = 0.34). Most loci showed only population-specific deviations from Hardy-Weinberg proportions, pointing to local population effects (e.g. null alleles). The described genomic resources are valuable for evolutionary studies in Symphonia and for comparative studies in plants. The methods are especially interesting for widespread tropical or endangered taxa with limited DNA availability.

Variable gene dispersal conditions and spatial deforestation patterns can interact to affect tropical tree conservation outcomes

Tropical lowland rain forest (TLRF) biodiversity is under threat from anthropogenic factors including deforestation which creates forest fragments of different sizes that can further undergo various internal patterns of logging. Such interventions can modify previous equilibrium abundance and spatial distribution patterns of offspring recruitment and/or pollen dispersal. Little is known about how these aspects of deforestation and fragmentation might synergistically affect TLRF tree recovery demographics and population genetics in newly formed forest fragments. To investigate these TLRF anthropogenic disturbance processes we used the computer program NEWGARDEN (NG), which models spatially-explicit, individual-based plant populations, to simulate 10% deforestation in six different spatial logging patterns for the plant functional type of a long-lived TLRF canopy tree species. Further, each logging pattern was analyzed under nine varying patterns of offspring versus pollen dispersal distances that could have arisen post-fragmentation. Results indicated that gene dispersal condition (especially via offspring) had a greater effect on population growth and genetic diversity retention (explaining 98.5% and 88.8% of the variance respectively) than spatial logging pattern (0.2% and 4.7% respectively), with ‘Near’ distance dispersal maximizing population growth and genetic diversity relative to distant dispersal. Within logged regions of the fragment, deforestation patterns closer to fragment borders more often exhibited lower population recovery rates and founding genetic diversity retention relative to more centrally located logging. These results suggest newly isolated fragments have populations that are more sensitive to the way in which their offspring and pollen dispersers are affected than the spatial pattern in which subsequent logging occurs, and that large variation in the recovery rates of different TLRF tree species attributable to altered gene dispersal regimens will be a likely outcome of fragmentation. Conservation implications include possible manual interventions (manual manipulations of offspring dispersers and/or pollinators) in forest fragments to increase population recovery and genetic diversity retention.

Comparative LC-MS-based metabolite profiling of the ancient tropical rainforest tree Symphonia globulifera

In the last few years, several phytochemical studies have been undertaken on the tropical tree Symphonia globulifera leading to the isolation and characterisation of several compounds exhibiting antiparasitic activities against Plasmodium falciparum, Trypanosoma brucei and Leishmania donovani. The comparative LC-MS based metabolite profiling study conducted on the tree led to the identification of compounds originating from specific tissues. The results showed that renewable organs/tissues can be used as the starting material for the production of polycyclic poly-prenylated-acylphloroglucinols, therefore reducing impacts on biodiversity. This study also underlined the lack of knowledge on the secondary metabolites produced by S. globulifera since only a small number of the total detected features were putatively identified using the database of known compounds for the species.

Comparative phylogeography in rainforest trees from Lower Guinea, Africa

Comparative phylogeography is an effective approach to assess the evolutionary history of biological communities. We used comparative phylogeography in fourteen tree taxa from Lower Guinea (Atlantic Equatorial Africa) to test for congruence with two simple evolutionary scenarios based on physio-climatic features 1) the W-E environmental gradient and 2) the N-S seasonal inversion, which determine climatic and seasonality differences in the region. We sequenced the trnC-ycf6 plastid DNA region using a dual sampling strategy: fourteen taxa with small sample sizes (dataset 1, mean n = 16/taxon), to assess whether a strong general pattern of allele endemism and genetic differentiation emerged; and four taxonomically well-studied species with larger sample sizes (dataset 2, mean n = 109/species) to detect the presence of particular shared phylogeographic patterns. When grouping the samples into two alternative sets of two populations, W and E, vs. N and S, neither dataset exhibited a strong pattern of allelic endemism, suggesting that none of the considered regions consistently harboured older populations. Differentiation in dataset 1 was similarly strong between W and E as between N and S, with 3-5 significant F ST tests out of 14 tests in each scenario. Coalescent simulations indicated that, given the power of the data, this result probably reflects idiosyncratic histories of the taxa, or a weak common differentiation pattern (possibly with population substructure) undetectable across taxa in dataset 1. Dataset 2 identified a common genetic break separating the northern and southern populations of Greenwayodendron suaveolens subsp. suaveolens var. suaveolens, Milicia excelsa, Symphonia globulifera and Trichoscypha acuminata in Lower Guinea, in agreement with differentiation across the N-S seasonal inversion. Our work suggests that currently recognized tree taxa or suspected species complexes can contain strongly differentiated genetic lineages, which could lead to misinterpretation of phylogeographic patterns. Therefore the evolutionary processes of such taxa require further study in African tropical rainforests.