Comparison of quantitative and molecular variation in agroforestry populations of the shea tree (Vitellaria paradoxa C.F. Gaertn) in Mali

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.

Directional and stabilizing selection shaped morphological, reproductive, and physiological traits of the invader Solidago canadensis

Trait evolution in invasive plant species is important because it can impact demographic parameters key to invasion success. Invasive plant species often show phenotypic clines along geographic and climatic gradients. However, the relative contributions of natural selection and neutral evolutionary processes to phenotypic trait variation among populations of invasive plants remain unclear. A common method to assess whether a trait has been shaped by natural selection or neutral evolutionary processes is to compare the geographical pattern for the trait of interest to the divergence in neutral genetic loci (i.e., Q ST -F ST comparisons). Subsequently, a redundancy analysis (RDA) can facilitate identification of putative agents of natural selection on the trait. Here, we employed both a Q ST -F ST comparisons approach and RDA to infer whether natural selection shaped traits of invasive populations of Solidago canadensis in China and identify the potential environmental drivers of natural selection. We addressed two questions: (1) Did natural selection drive phenotypic trait variation among S. canadensis populations? (2) Did climatic, latitudinal, longitudinal, and altitudinal gradients drive patterns of genetic variation among S. canadensis populations? We found significant directional selection for several morphological and reproductive traits (i.e., Q ST  > F ST) and stabilizing selection for physiological traits (i.e., Q ST  < F ST). The RDA showed that stem biomass of S. canadensis was strongly positively correlated with longitude, while leaf width ratio and specific leaf area were significantly positively correlated with the mean diurnal range. Stem biomass had a strong negative correlation with annual precipitation. Moreover, height of S. canadensis individuals was strongly positively correlated with altitude and precipitation of the wettest quarter. A longitudinal shift in precipitation seasonality likely selected for larger stem biomass in S. canadensis. Overall, these results suggest that longitudinal and altitudinal clines in climate exerted strong selection pressures that shaped the phenotypic traits of S. canadensis.

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.

Identification of genes involved in shea butter biosynthesis from Vitellaria paradoxa fruits through transcriptomics and functional heterologous expression

Shea tree (Vitellaria paradoxa) is one economically important plant species that mainly distributes in West Africa. Shea butter extracted from shea fruit kernels can be used as valuable products in the food and cosmetic industries. The most valuable composition in shea butter was one kind of triacylglycerol (TAG), 1,3-distearoyl-2-oleoyl-glycerol (SOS, C18:0–C18:1–C18:0). However, shea butter production is limited and little is known about the genetic information of shea tree. In this study, we tried to reveal genetic information of shea tree and identified shea TAG biosynthetic genes for future shea butter production in yeast cell factories. First, we measured lipid content, lipid composition, and TAG composition of seven shea fruits at different ripe stages. Then, we performed transcriptome analysis on two shea fruits containing obviously different levels of SOS and revealed a list of TAG biosynthetic genes potentially involved in TAG biosynthesis. In total, 4 glycerol-3-phosphate acyltransferase (GPAT) genes, 8 lysophospholipid acyltransferase (LPAT) genes, and 11 diacylglycerol acyltransferase (DGAT) genes in TAG biosynthetic pathway were predicted from the assembled transcriptome and 14 of them were cloned from shea fruit cDNA. Furthermore, the heterologous expression of these 14 potential GPAT, LPAT, and DGAT genes in Saccharomyces cerevisiae changed yeast fatty acid and lipid profiles, suggesting that they functioned in S. cerevisiae. Moreover, two shea DGAT genes, VpDGAT1 and VpDGAT7, were identified as functional DGATs in shea tree, showing they might be useful for shea butter (SOS) production in yeast cell factories.

Comparisons between Q(ST) and F(ST) --how wrong have we been?

The comparison between quantitative genetic divergence (Q(ST) ) and neutral genetic divergence (F(ST) ) among populations has become the standard test for historical signatures of selection on quantitative traits. However, when the mutation rate of neutral markers is relatively high in comparison with gene flow, estimates of F(ST) will decrease, resulting in upwardly biased comparisons of Q(ST) vs. F(ST) . Reviewing empirical studies, the difference between Q(ST) and F(ST) is positively related to marker heterozygosity. After refuting alternative explanations for this pattern, we conclude that marker mutation rate indeed has had a biasing effect on published Q(ST) -F(ST) comparisons. Hence, it is no longer clear that populations have commonly diverged in response to divergent selection. We present and discuss potential solutions to this bias. Comparing Q(ST) with recent indices of neutral divergence that statistically correct for marker heterozygosity (Hedrick's G'st and Jost's D) is not advised, because these indices are not theoretically equivalent to Q(ST) . One valid solution is to estimate F(ST) from neutral markers with mutation rates comparable to those of the loci underlying quantitative traits (e.g. SNPs). Q(ST) can also be compared to Φ(ST) (Phi(ST) ) of amova, as long as the genetic distance among allelic variants used to estimate Φ(ST) reflects evolutionary history: in that case, neutral divergence is independent of mutation rate. In contrast to their common usage in comparisons of Q(ST) and F(ST) , microsatellites typically have high mutation rates and do not evolve according to a simple evolutionary model, so are best avoided in Q(ST) -F(ST) comparisons.

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

While the genetic structure of many tree species in temperate, American and Asian regions is largely explained by climatic oscillations and subsequent habitat contractions and expansions, little is known about Africa. We investigated the genetic diversity and structure of shea tree (Vitellaria paradoxa,) in Western Africa, an economically important tree species in the Sudano-Sahelian zone. Eleven nuclear microsatellites (nuc) were used to genotype 673 trees selected in 38 populations. They revealed moderate to high within-population diversity: allelic richness ranged from R(nuc) = 3.99 to 5.63. This diversity was evenly distributed across West Africa. Populations were weakly differentiated (F(STnuc) = 0.085; P < 0.0001) and a pattern of isolation by distance was noted. No phylogeographic signal could be detected across the studied sample. Additionally, two chloroplast microsatellite loci, leading to 11 chlorotypes, were used to analyse a sub-set of 370 individuals. Some variation in chloroplast allelic richness among populations could be detected (R(cp) = 0.00 to 4.36), but these differences were not significant. No trend with latitude and longitude were observed. Differentiation was marked (G(STcp) = 0.553; P < 0.0001), but without a significant phylogeographical signal. Population expansion was detected considering the total population using approximate Bayesian computation (nuclear microsatellites) and mismatch distribution (chloroplast microsatellites) methods. This expansion signal and the isolation by distance pattern could be linked to the past climatic conditions in West Africa during the Pleistocene and Holocene which should have been favourable to shea tree development. In addition, human activities through agroforestry and domestication (started 10,000 bp) have probably enhanced gene flow and population expansion.

Past climate changes explain the phylogeography of Vitellaria paradoxa over Africa

The evolution of the savanna biome has been deeply marked by repeated contraction/expansion phases due to climate perturbations during the Quaternary period. In this study, we investigated the impact of the last glacial maximum (LGM) on the present genetic pattern of Vitellaria paradoxa (shea tree), a major African savanna tree. A range-wide sampling of the species enabled us to sample 374 individuals from 71 populations distributed throughout sub-Sahelian Africa. Trees were genotyped using 3 chloroplasts and 12 nuclear microsatellites, and were sequenced for 2 polymorphic chloroplast intergenic spacers. Analyses of genetic diversity and structure were based on frequency-based and Bayesian methods. Potential distributions of V. paradoxa at present, during the LGM and the last interglacial period, were examined using DIVA-GIS ecological niche modelling (ENM). Haplotypic and allelic richness varied significantly across the range according to chloroplast and nuclear microsatellites, which pointed to higher diversity in West Africa. A high but contrasted level of differentiation was revealed among populations with a clear phylogeographic signal, with both nuclear (F(ST) = 0.21; R(ST) = 0.28; R(ST) > R(ST) (permuted)) and chloroplast simple sequence repeats (SSRs) (G(ST) = 0.81; N(ST) = 0.90; N(ST) > N(ST) (permuted)). We identified a strong geographically related structure separating western and eastern populations, and a substructure in the eastern part of the area consistent with subspecies distinction. Using ENM, we deduced that perturbations during the LGM fragmented the potential eastern distribution of shea tree, but not its distribution in West Africa. Our main results suggest that climate variations are the major factor explaining the genetic pattern of V. paradoxa.

Adaptive divergence for a fitness-related trait among invasive Ambrosia artemisiifolia populations in France

The impact of natural selection on the adaptive divergence of invasive populations can be assessed by testing the null hypothesis that the extent of quantitative genetic differentiation (Q(ST) ) would be similar to that of neutral molecular differentiation (F(ST) ). Using eight microsatellite loci and a common garden approach, we compared Q(ST) and F(ST) among ten populations of an invasive species Ambrosia artemisiifolia (common ragweed) in France. In a common garden study with varying water and nutrient levels, we measured Q(ST) for five traits (height, total biomass, reproductive allocation, above- to belowground biomass ratio, and days to flowering). Although low F(ST) indicated weak genetic structure and strong gene flow among populations, we found significant diversifying selection (Q(ST) > F(ST) ) for reproductive allocation that may be closely related to fitness. It suggests that abiotic conditions may have exerted selection pressure on A. artemisiifolia populations to differentiate adaptively, such that populations at higher altitude or latitude evolved greater reproductive allocation. As previous studies indicate multiple introductions from various source populations of A. artemisiifolia in North America, our results suggest that the admixture of introduced populations may have increased genetic diversity and additive genetic variance, and in turn, promoted the rapid evolution and adaptation of this invasive species.

Spatial genetic structuring of baobab (Adansonia digitata, Malvaceae) in the traditional agroforestry systems of West Africa

This study evaluates the spatial genetic structure of baobab (Adansonia digitata) populations from West African agroforestry systems at different geographical scales using AFLP fingerprints. Eleven populations from four countries (Benin, Ghana, Burkina Faso, and Senegal) had comparable levels of genetic diversity, although the two populations in the extreme west (Senegal) had less diversity. Pairwise F(ST) ranged from 0.02 to 0.28 and increased with geographic distance, even at a regional scale. Gene pools detected by Bayesian clustering seem to be a byproduct of the isolation-by-distance pattern rather than representing actual discrete entities. The organization of genetic diversity appears to result essentially from spatially restricted gene flow, with some influences of human seed exchange. Despite the potential for relatively long-distance pollen and seed dispersal by bats within populations, statistically significant spatial genetic structuring within populations (SGS) was detected and gave a mean indirect estimate of neighborhood size of ca. 45. This study demonstrated that relatively high levels of genetic structuring are present in baobab at both large and within-population level, which was unexpected in regard to its dispersal by bats and the influence of human exchange of seeds. Implications of these results for the conservation of baobab populations are discussed.

Speciation with gene flow in the large white-headed gulls: does selection counterbalance introgression?

We investigated the role of selection in generating and maintaining species distinctness in spite of ongoing gene flow, using two zones of secondary contact between large gull species in Europe (Larus argentatus and Larus cachinnans) and North America (Larus glaucescens and Larus occidentalis). We used the pattern of neutral genetic differentiation at nine microsatellite loci (F(ST)) as an indicator of expected changes under neutral processes and compared it with phenotypic differentiation (P(ST)) for a large number of traits (size, plumage melanism and coloration of bare parts). Even assuming very low heritability, interspecific divergence between L. glaucescens and L. occidentalis in plumage melanism and orbital ring colour clearly exceeded neutral differentiation. Similarly, melanism of the central primaries was highly divergent between L. argentatus and L. cachinnans. Such divergence is unlikely to have arisen randomly and is therefore attributed to spatially varying selection. Variation in plumage melanism in both transects agrees with Gloger's rule, which suggests that latitude (and associated sun and humidity gradients) could be the selective pressure shaping differentiation in plumage melanism. We suggest that strong species differentiation in orbital ring colour results from sexual selection. We conclude that these large gull species, along with other recently diverged species that hybridize after coming into secondary contact, may differ only in restricted regions of the genome that are undergoing strong disruptive selection because of their phenotypic effects.

Comparative Molecular and Phytochemical Study of the Tree SpeciesSantalum austrocaledonicum (Santalaceae) Distributed in the New-Caledonian Archipelago

We have tried to elucidate the origin of phytochemical variation in trees by studying concomitantly the chemical and microsatellite variations in Santalum austrocaledonicum. Eight natural populations were sampled in the New-Caledonian archipelago, a total of 157 individuals being analyzed. The main components, as revealed by gas chromatography (GC), were alpha- and beta-santalol (as in other sandalwood species), although the level of (Z)-lanceol was particularly high. Most of the chemical variation was observed within populations (83.7%). With microsatellites, the variation between populations was more pronounced (32% of the total variation). Although the chemical variation between populations was small, we investigated the effects of genetic drift and migration by comparing the chemical- and molecular-differentiation patterns. The poor congruence between neighbor-joining trees, confirmed by the non-significant Mantel test between the molecular and chemical distance matrices (R=0.26, P=0.12), showed that genetic drift and migration are not the main evolutionary forces acting on chemical differentiation between populations. We could not find any effect of soil and rainfall conditions neither. Although the impact of drift and migration cannot be discounted in rationalizing between-population differentiation, the low variation among populations could result from a stabilizing selection caused by the same phytopathogen charge across the natural range.