Molecular data reveal isolation by distance and past population expansion for the shea tree (Vitellaria paradoxa C.F. Gaertn) in West Africa

An innovative optimized protocol for high-quality genomic DNA extraction from recalcitrant Shea tree (Vitellaria paradoxa, C.F. Gaertn) plant and its suitability for downstream applications

BACKGROUND

It is not always easy to find a universal protocol for the extraction of genomic DNA (gDNA) from plants. Extraction of gDNA from plants such as shea with a lot of polysaccharides in their leaves is done in two steps: a first step to remove the polysaccharides and a second step for the extraction of the gDNA. In this work, we designed a protocol for extracting high-quality gDNA from shea tree and demonstrate its suitability for downstream molecular applications.

METHODS

Fifty milligrams of leaf and root tissues were used to test the efficiency of our protocol. The quantity of gDNA was measured with the NanoDrop spectrometer and the quality was checked on agarose gel. Its suitability for use in downstream applications was tested with restriction enzymes, SSRs and RAPD polymerase chain reactions and Sanger sequencing.

RESULTS

The average yield of gDNA was 5.17; 3.96; 2.71 and 2.41 µg for dry leaves, dry roots, fresh leaves and fresh roots respectively per 100 mg of tissue. Variance analysis of the yield showed significant difference between all tissue types. Leaf gDNA quality was better compared to root gDNA at the absorbance ratio A260/280 and A260/230. The minimum amplifiable concentration of leaf gDNA was 1 pg/µl while root gDNA remained amplifiable at 10 pg/µl. Genomic DNA obtained was also suitable for sequencing.

CONCLUSION

This protocol provides an efficient, convenient and cost effective DNA extraction method suitable for use in various vitellaria paradoxa genomic studies.

Climate change will disproportionally affect the most genetically diverse lineages of a widespread African tree species

Global climate change is proceeding at an alarming rate with major ecological and genetic consequences for biodiversity, particularly in drylands. The response of species to climate change may differ between intraspecific genetic groups, with major implications for conservation. We used molecular data from 10 nuclear and two chloroplast genomes to identify phylogeographic groups within 746 individuals from 29 populations of Senegalia senegal, a savannah tree species in sub-Saharan Africa. Three phylogroups are identified corresponding to Sudano-Sahelian, Zambezian and Southern African biogeographic regions in West, East and Southern Africa. Genetic diversity was highest in Southern and Zambesian and lowest in the Sudano-Sahelian phylogroups. Using species distribution modeling, we infer highly divergent future distributions of the phylogroups under three climate change scenarios. Climate change will lead to severe reductions of distribution area of the genetically diverse Zambezian (- 41-- 54%) and Southern (- 63-- 82%) phylogroups, but to an increase for the genetically depauperate Sudano-Sahelian (+ 7- + 26%) phylogroups. This study improves our understanding of the impact of climate change on the future distribution of this species. This knowledge is particularly useful for biodiversity management as the conservation of genetic resources needs to be considered in complementary strategies of in-situ conservation and assisted migration.

Genomic Resources to Guide Improvement of the Shea Tree

A defining component of agroforestry parklands across Sahelo-Sudanian Africa (SSA), the shea tree (Vitellaria paradoxa) is central to sustaining local livelihoods and the farming environments of rural communities. Despite its economic and cultural value, however, not to mention the ecological roles it plays as a dominant parkland species, shea remains semi-domesticated with virtually no history of systematic genetic improvement. In truth, shea's extended juvenile period makes traditional breeding approaches untenable; but the opportunity for genome-assisted breeding is immense, provided the foundational resources are available. Here we report the development and public release of such resources. Using the FALCON-Phase workflow, 162.6 Gb of long-read PacBio sequence data were assembled into a 658.7 Mbp, chromosome-scale reference genome annotated with 38,505 coding genes. Whole genome duplication (WGD) analysis based on this gene space revealed clear signatures of two ancient WGD events in shea's evolutionary past, one prior to the Astrid-Rosid divergence (116-126 Mya) and the other at the root of the order Ericales (65-90 Mya). In a first genome-wide look at the suite of fatty acid (FA) biosynthesis genes that likely govern stearin content, the primary determinant of shea butter quality, relatively high copy numbers of six key enzymes were found (KASI, KASIII, FATB, FAD2, FAD3, and FAX2), some likely originating in shea's more recent WGD event. To help translate these findings into practical tools for characterization, selection, and genome-wide association studies (GWAS), resequencing data from a shea diversity panel was used to develop a database of more than 3.5 million functionally annotated, physically anchored SNPs. Two smaller, more curated sets of suggested SNPs, one for GWAS (104,211 SNPs) and the other targeting FA biosynthesis genes (90 SNPs), are also presented. With these resources, the hope is to support national programs across the shea belt in the strategic, genome-enabled conservation and long-term improvement of the shea tree for SSA.

Fine-scale spatial genetic structure, mating, and gene dispersal patterns in Parkia biglobosa populations with different levels of habitat fragmentation

PREMISE

A good understanding of genetic variation and gene dispersal in tree populations is crucial for their sustainable management, particularly in a context of rapid environmental changes. West African Sudanian savannahs are being fragmented and degraded, partly due to expansion of crop cultivation and monocultures that reduce tree density and may impact pollinators. The population dynamics of important indigenous trees could also be affected. We investigated the influence of habitat fragmentation on patterns of genetic diversity and gene dispersal of a key Sudanian agroforestry tree species, Parkia biglobosa.

METHODS

Using 10 highly polymorphic nuclear microsatellites, we genotyped 2475 samples from reproductive trees, seedlings, and embryos in four tree populations presenting different levels of habitat fragmentation.

RESULTS

Parkia biglobosa presented similar high genetic diversity across the four populations studied. Genetic diversity and inbreeding were similar between adults and embryo cohorts. In all four populations, the selfing rate was less than 1%. The effective number of pollen donors per tree was high (N_EP ~ 18-22), as was the pollen immigration rate (from 34 to 74%). Pollen dispersal was characterized by a fat-tailed distribution with mean estimates exceeding 200 m. In three populations, stem diameter had a pronounced effect on male reproductive success. Here, the highest male reproductive success was observed in trees with a diameter at breast height between 60 and 75 cm.

CONCLUSIONS

At the scale analyzed, fragmentation does not seem to pose limitations to gene flow in any of the sites investigated, regardless of the landscape configuration associated with the different tree stands. The study provides useful insights on the reproductive biology of an important tree species in the West African savannahs.

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.

Gene pool sharing and genetic bottleneck effects in subpopulations of Eschweilera ovata (Cambess.) Mart. ex Miers (Lecythidaceae) in the Atlantic Forest of southern Bahia, Brazil

Forest loss and fragmentation are the main threats to the maintenance of the Atlantic Forest, an important global biodiversity hotspot. Because of the current critical level of deforestation, ecological corridors are needed to facilitate species dispersion and gene flow among fragments. This study was conducted to investigate the genetic variability and gene pool sharing of Eschweilera ovata in five forest remnants in southern Bahia, Brazil using nuclear simple sequence repeat (nSSR) and plastid simple sequence repeat (cpSSR) microsatellite markers. cpSSR marker analysis revealed the domains of four haplotypes, showing that 80% of the individuals had only four maternal origins, reflecting a founder effect and/or genetic bottleneck. The results of cpSSR and nSSR analyses indicated moderate genetic diversity, particularly in conservation units with full protection, which showed the best parameters of all areas evaluated. Another indication of the susceptibility of these populations to forest loss and fragmentation was the strong genetic bottleneck observed. In contrast, genetic structure analyses (FST and discriminant analysis of principal components) revealed gene pool sharing between the subpopulations, which may reflect the historical gene flow that occurred before forest fragmentation.

Phylogeography of African Locust Bean (Parkia biglobosa) Reveals Genetic Divergence and Spatially Structured Populations in West and Central Africa

The evolutionary history of African savannah tree species is crucial for the management of their genetic resources. In this study, we investigated the phylogeography of Parkia biglobosa and its modeled distribution under past and present climate conditions. This tree species is very valued and widespread in West Africa, providing edible and medicinal products. A large sample of 1610 individuals from 84 populations, distributed across 12 countries in Western and Central Africa, were genotyped using 8 nuclear microsatellites. Individual-based assignments clearly distinguished 3 genetic clusters, extreme West Africa (EWA), center of West Africa (CWA), and Central Africa (CA). Overall, estimates of genetic diversity were moderate to high, with lower values for populations in EWA (allelic richness after rarefaction [A_R] = 6.4, expected heterozygosity [H_E] = 0.78, and observed heterozygosity [H_O] = 0.7) and CA (A_R = 5.9, H_E = 0.67, and H_O = 0.61) compared with populations in CWA (A_R = 7.3, H_E = 0.79, and H_O = 0.75). The overall population differentiation was found to be moderate (F_ST = 0.09). A highly significant isolation by distance pattern was detected, with a marked phylogeographic signature suggesting possible effects of past climate and geographic barriers to migration. Modeling the potential distribution of the species showed a contraction during the last glaciations followed by expansion events. The exploratory approximate Bayesian computation conducted suggests a best-supported scenario in which the cluster CWA traced back to the ancestral populations and a first split between EWA and CWA took place about 160000 years before present (BP), then a second split divided CA and CWA, about 100000 years BP. However, our genetic data do not enable us to conclusively distinguish among a few alternative possible scenarios.

Comparative phylogeography of Aedes mosquitoes and the role of past climatic change for evolution within Africa

The study of demographic processes involved in species diversification and evolution ultimately provides explanations for the complex distribution of biodiversity on earth, indicates regions important for the maintenance and generation of biodiversity, and identifies biological units important for conservation or medical consequence. African and forest biota have both received relatively little attention with regard to understanding their diversification, although one possible mechanism is that this has been driven by historical climate change. To investigate this, we implemented a standard population genetics approach along with Approximate Bayesian Computation, using sequence data from two exon-primed intron-crossing (EPIC) nuclear loci and mitochondrial cytochrome oxidase subunit I, to investigate the evolutionary history of five medically important and inherently forest dependent mosquito species of the genus Aedes. By testing different demographic hypotheses, we show that Aedes bromeliae and Aedes lilii fit the same model of lineage diversification, admixture, expansion, and recent population structure previously inferred for Aedes aegypti. In addition, analyses of population structure show that Aedes africanus has undergone lineage diversification and expansion while Aedes hansfordi has been impacted by population expansion within Uganda. This congruence in evolutionary history is likely to relate to historical climate-driven habitat change within Africa during the late Pleistocene and Holocene epoch. We find differences in the population structure of mosquitoes from Tanzania and Uganda compared to Benin and Uganda which could relate to differences in the historical connectivity of forests across the continent. Our findings emphasize the importance of recent climate change in the evolution of African forest biota.

Genetic structure and hierarchical population divergence history of Acer mono var. mono in South and Northeast China

Knowledge of the genetic structure and evolutionary history of tree species across their ranges is essential for the development of effective conservation and forest management strategies. Acer mono var. mono, an economically and ecologically important maple species, is extensively distributed in Northeast China (NE), whereas it has a scattered and patchy distribution in South China (SC). In this study, the genetic structure and demographic history of 56 natural populations of A. mono var. mono were evaluated using seven nuclear microsatellite markers. Neighbor-joining tree and STRUCTURE analysis clearly separated populations into NE and SC groups with two admixed-like populations. Allelic richness significantly decreased with increasing latitude within the NE group while both allelic richness and expected heterozygosity showed significant positive correlation with latitude within the SC group. Especially in the NE region, previous studies in Quercus mongolica and Fraxinus mandshurica have also detected reductions in genetic diversity with increases in latitude, suggesting this pattern may be common for tree species in this region, probably due to expansion from single refugium following the last glacial maximum (LGM). Approximate Bayesian Computation-based analysis revealed two major features of hierarchical population divergence in the species' evolutionary history. Recent divergence between the NE group and the admixed-like group corresponded to the LGM period and ancient divergence of SC groups took place during mid-late Pleistocene period. The level of genetic differentiation was moderate (FST  = 0.073; G'ST  = 0.278) among all populations, but significantly higher in the SC group than the NE group, mirroring the species' more scattered distribution in SC. Conservation measures for this species are proposed, taking into account the genetic structure and past demographic history identified in this study.

Divergent pattern of nuclear genetic diversity across the range of the Afromontane Prunus africana mirrors variable climate of African highlands

BACKGROUND AND AIMS

Afromontane forest ecosystems share a high similarity of plant and animal biodiversity, although they occur mainly on isolated mountain massifs throughout the continent. This resemblance has long provoked questions on former wider distribution of Afromontane forests. In this study Prunus africana (one of the character trees of Afromontane forests) is used as a model for understanding the biogeography of this vegetation zone.

METHODS

Thirty natural populations from nine African countries covering a large part of Afromontane regions were analysed using six nuclear microsatellites. Standard population genetic analysis as well as Bayesian and maximum likelihood models were used to infer genetic diversity, population differentiation, barriers to gene flow, and recent and all migration among populations.

KEY RESULTS

Prunus africana exhibits strong divergence among five main Afromontane regions: West Africa, East Africa west of the Eastern Rift Valley (ERV), East Africa east of the ERV, southern Africa and Madagascar. The strongest divergence was evident between Madagascar and continental Africa. Populations from West Africa showed high similarity with East African populations west of the ERV, whereas populations east of the ERV are closely related to populations of southern Africa, respectively.

CONCLUSIONS

The observed patterns indicate divergent population history across the continent most likely associated to Pleistocene changes in climatic conditions. The high genetic similarity between populations of West Africa with population of East Africa west of the ERV is in agreement with faunistic and floristic patterns and provides further evidence for a historical migration route. Contrasting estimates of recent and historical gene flow indicate a shift of the main barrier to gene flow from the Lake Victoria basin to the ERV, highlighting the dynamic environmental and evolutionary history of the region.