New genetic markers for Sapotaceae phylogenomics: More than 600 nuclear genes applicable from family to population levels.

Some tropical plant families, such as the Sapotaceae, have a complex taxonomy, which can be resolved using Next Generation Sequencing (NGS). For most groups however, methodological protocols are still missing. Here we identified 531 monocopy genes and 227 Short Tandem Repeats (STR) markers and tested them on Sapotaceae using target capture and NGS. The probes were designed using two genome skimming samples from Capurodendron delphinense and Bemangidia lowryi, both from the Tseboneae tribe, as well as the published Manilkara zapota transcriptome from the Sapotoideae tribe. We combined our probes with 261 additional ones previously published and designed for the entire angiosperm group. On a total of 792 low-copy genes, 638 showed no signs of paralogy and were used to build a phylogeny of the family with 231 individuals from all main lineages. A highly supported topology was obtained at high taxonomic ranks but also at the species level. This phylogeny revealed the existence of more than 20 putative new species. Single nucleotide polymorphisms (SNPs) extracted from the 638 genes were able to distinguish lineages within a species complex and to highlight geographical structuration. STR were recovered efficiently for the species used as reference (C. delphinense) but the recovery rate decreased dramatically with the phylogenetic distance to the focal species. Altogether, the new loci will help reaching a sound taxonomic understanding of the family Sapotaceae for which many circumscriptions and relationships are still debated, at the species, genus and tribe levels.

Plant NHX Antiporters: From Function to Biotechnological Application, with Case Study.

Salt stress is one of the major abiotic stresses that negatively affect crops worldwide. Plants have evolved a series of mechanisms to cope with the limitations imposed by salinity. Molecular mechanisms, including the upregulation of cation transporters such as the Na+/H+ antiporters, are one of the processes adopted by plants to survive in saline environments. NHX antiporters are involved in salt tolerance, development, cell expansion, growth performance and disease resistance of plants. They are integral membrane proteins belonging to the widely distributed CPA1 sub-group of monovalent cation/H+ antiporters and provide an important strategy for ionic homeostasis in plants under saline conditions. These antiporters are known to regulate the exchange of sodium and hydrogen ions across the membrane and are ubiquitous to all eukaryotic organisms. With the genomic approach, previous studies reported that a large number of proteins encoding Na+/H+ antiporter genes have been identified in many plant species and successfully introduced into desired species to create transgenic crops with enhanced tolerance to multiple stresses. In this review, we focus on plant antiporters and all the aspects from their structure, classification, function to their in silico analysis. On the other hand, we performed a genome-wide search to identify the predicted NHX genes in Argania spinosa L. We highlighted for the first time the presence of four putative NHX (AsNHX1-4) from the Argan tree genome, whose phylogenetic analysis revealed their classification in one distinct vacuolar cluster. The essential information of the four putative NHXs, such as gene structure, subcellular localization and transmembrane domains was analyzed.

Metapopulation Vicariance, Age of Island Taxa and Dispersal: A Case Study Using the Pacific Plant Genus Planchonella (Sapotaceae).

Oceanic islands originate from volcanism or tectonic activity without connections to continental landmasses, are colonized by organisms, and eventually vanish due to erosion and subsidence. Colonization of oceanic islands occurs through long-distance dispersals (LDDs) or metapopulation vicariance, the latter resulting in lineages being older than the islands they inhabit. If metapopulation vicariance is valid, island ages cannot be reliably used to provide maximum age constraints for molecular dating. We explore the relationships between the ages of members of a widespread plant genus (Planchonella, Sapotaceae) and their host islands across the Pacific to test various assumptions of dispersal and metapopulation vicariance. We sampled three nuclear DNA markers from 156 accessions representing some 100 Sapotaceae taxa, and analyzed these in BEAST with a relaxed clock to estimate divergence times and with a phylogeographic diffusion model to estimate range expansions over time. The phylogeny was calibrated with a secondary point (the root) and fossils from New Zealand. The dated phylogeny reveals that the ages of Planchonella species are, in most cases, consistent with the ages of the islands they inhabit. Planchonella is inferred to have originated in the Sahul Shelf region, to which it back-dispersed multiple times. Fiji has been an important source for range expansion in the Pacific for the past 23 myr. Our analyses reject metapopulation vicariance in all cases tested, including between oceanic islands, evolution of an endemic Fiji-Vanuatu flora, and westward rollback vicariance between Vanuatu and the Loyalty Islands. Repeated dispersal is the only mechanism able to explain the empirical data. The longest (8900 km) identified dispersal is between Palau in the Pacific and the Seychelles in the Indian Ocean, estimated at 2.2 Ma (0.4-4.8 Ma). The first split in a Hawaiian lineage (P. sandwicensis) matches the age of Necker Island (11.0 Ma), when its ancestor diverged into two species that are distinguished by purple and yellow fruits. Subsequent establishment across the Hawaiian archipelago supports, in part, progression rule colonization. In summary, we found no explanatory power in metapopulation vicariance and conclude that Planchonella has expanded its range across the Pacific by LDD. We contend that this will be seen in many other groups when analyzed in detail.

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.

Two tropical conifers show strong growth and water-use efficiency responses to altered CO2 concentration.

BACKGROUND AND AIMS: Conifers dominated wet lowland tropical forests 100 million years ago (MYA). With a few exceptions in the Podocarpaceae and Araucariaceae, conifers are now absent from this biome. This shift to angiosperm dominance also coincided with a large decline in atmospheric CO2 concentration (ca). We compared growth and physiological performance of two lowland tropical angiosperms and conifers at ca levels representing pre-industrial (280 ppm), ambient (400 ppm) and Eocene (800 ppm) conditions to explore how differences in ca affect the growth and water-use efficiency (WUE) of seedlings from these groups. METHODS: Two conifers (Araucaria heterophylla and Podocarpus guatemalensis) and two angiosperm trees (Tabebuia rosea and Chrysophyllum cainito) were grown in climate-controlled glasshouses in Panama. Growth, photosynthetic rates, nutrient uptake, and nutrient use and water-use efficiencies were measured. KEY RESULTS: Podocarpus seedlings showed a stronger (66 %) increase in relative growth rate with increasing ca relative to Araucaria (19 %) and the angiosperms (no growth enhancement). The response of Podocarpus is consistent with expectations for species with conservative growth traits and low mesophyll diffusion conductance. While previous work has shown limited stomatal response of conifers to ca, we found that the two conifers had significantly greater increases in leaf and whole-plant WUE than the angiosperms, reflecting increased photosynthetic rate and reduced stomatal conductance. Foliar nitrogen isotope ratios (δ15N) and soil nitrate concentrations indicated a preference in Podocarpus for ammonium over nitrate, which may impact nitrogen uptake relative to nitrate assimilators under high ca SIGNIFICANCE: Podocarps colonized tropical forests after angiosperms achieved dominance and are now restricted to infertile soils. Although limited to a single species, our data suggest that higher ca may have been favourable for podocarp colonization of tropical South America 60 MYA, while plasticity in photosynthetic capacity and WUE may help account for their continued persistence under large changes in ca since the Eocene.

Isolation and characterization of 8 microsatellite loci from Chrysophyllum gonocarpum (Sapotaceae).

Chrysophyllum gonocarpum is a tropical tree species that is very important in the recovery of heterogeneous forests and of degraded areas of permanent preservation. We identified microsatellite loci for C. gonocarpum to assess the genetic variability and the patterns of the population structure of the species. We isolated 8 microsatellite primers by using CT- and GT-enriched genomic libraries. We detected 2-4 alleles with 2.9 alleles per locus on average, by polymerase chain reaction. Test for cross-amplification showed that some loci were successfully amplified in 2 other Chrysophyllum species. The microsatellites can be used to assess the genetic diversity and population structure of C. gonocarpum. Some primer pairs can be amplified in C. marginatum and C. splendens.

Origins and close relatives of a semi-domesticated neotropical fruit tree: Chrysophyllum cainito (Sapotaceae).

PREMISE OF THE STUDY: Understanding patterns and processes associated with domestication has implications for crop development and agricultural biodiversity conservation. Semi-domesticated crops provide excellent opportunities to examine the interplay of natural and anthropogenic influences on plant evolution. The domestication process has not been thoroughly examined in many tropical perennial crop species. Chrysophyllum cainito (Sapotaceae), the star apple or caimito, is a semi-domesticated species widely cultivated for its edible fruits. It is known to be native to the neotropics, but the precise geographic origins of wild and cultivated forms are unresolved. METHODS: We used nuclear ribosomal ITS sequences to infer phylogenetic relationships among C. cainito and close relatives (section Chrysophyllum). We employed phylogeographic approaches using ITS and plastid sequence data to determine geographic origins and center(s) of domestication of caimito. KEY RESULTS: ITS data suggest a close relationship between C. cainito and C. argenteum. Plastid haplotype networks reveal several haplotypes unique to individual taxa but fail to resolve distinct lineages for either C. cainito or C. argenteum. Caimito populations from northern Mesoamerica and the Antilles exhibit a subset of the genetic diversity found in southern Mesoamerica. In Panama, cultivated caimito retains high levels of the diversity seen in wild populations. CONCLUSIONS: Chrysophyllum cainito is most closely related to a clade containing Central and South American C. argenteum, including subsp. panamense. We hypothesize that caimito is native to southern Mesoamerica and was domesticated from multiple wild populations in Panama. Subsequent migration into northern Mesoamerica and the Antilles was mediated by human cultivation.

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.

Ubiquitin promoter-terminator cassette promotes genetically stable expression of the taste-modifying protein miraculin in transgenic lettuce.

Lettuce is a commercially important leafy vegetable that is cultivated worldwide, and it is also a target crop for plant factories. In this study, lettuce was selected as an alternative platform for recombinant miraculin production because of its fast growth, agronomic value, and wide availability. The taste-modifying protein miraculin is a glycoprotein extracted from the red berries of the West African native shrub Richadella dulcifica. Because of its limited natural availability, many attempts have been made to produce this protein in suitable alternative hosts. We produced transgenic lettuce with miraculin gene driven either by the ubiquitin promoter/terminator cassette from lettuce or a 35S promoter/nos terminator cassette. Miraculin gene expression and miraculin accumulation in both cassettes were compared by quantitative real-time PCR analysis, Western blotting, and enzyme-linked immunosorbent assay. The expression level of the miraculin gene and protein in transgenic lettuce was higher and more genetically stable in the ubiquitin promoter/terminator cassette than in the 35S promoter/nos terminator cassette. These results demonstrated that the ubiquitin promoter/terminator cassette is an efficient platform for the genetically stable expression of the miraculin protein in lettuce and hence this platform is of benefit for recombinant miraculin production on a commercial scale.

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.

Phylogeny and generic limits in the Niemeyera complex of New Caledonian Sapotaceae: evidence of multiple origins of the anisomerous flower.

Traditional generic limits within the "Niemeyera complex" (Sapotaceae: Chrysophylloideae) in Australia and New Caledonia do not correspond to natural groups. We analyzed nuclear (ETS, ITS) and chloroplast (trnH-psbA, trnS-G) sequence data, and 42 morphological characters, using a near-complete taxon sampling. The resulting phylogeny provides a new generic framework where Leptostylis and Sebertia are monophyletic, Niemeyera is recognized as a small genus confined to Australia, and the circumscriptions of Achradotypus and Pycnandra require significant modification. This framework allows about 20 recently discovered species to be described in appropriate genera and assessed for their conservation status. Evolutionary changes in two widely used characters, anisomerous flowers and the number of stamens inserted opposite each corolla lobe, have occurred multiple times. There is no evidence that anisomery originated through hybridization as suggested in other groups. Instead, the two characters are closely coupled and often mutually exclusive. The diagnostic value of morphological characters varies considerably; for example, the presence, absence, and type of malpighiaceous hairs convey more phylogenetic information than many traditionally used features. Criteria and options for a new generic classification are discussed, and a reconstruction of the hypothesized ancestor that likely colonized New Caledonia in the Oligocene is presented.

Accounting for variation of substitution rates through time in Bayesian phylogeny reconstruction of Sapotoideae (Sapotaceae).

We used Bayesian phylogenetic analysis of 5 kb of chloroplast DNA data from 68 Sapotaceae species to clarify phylogenetic relationships within Sapotoideae, one of the two major clades within Sapotaceae. Variation in substitution rates through time was shown to be a very important aspect of molecular evolution for this data set. Relative rates tests indicated that changes in overall rate have taken place in several lineages during the history of the group and Bayes factors strongly supported a covarion model, which allows the rate of a site to vary over time, over commonly used models that only allow rates to vary across sites. Rate variation over time was actually found to be a more important model component than rate variation across sites. The covarion model was originally developed for coding gene sequences and has so far only been tested for this type of data. The fact that it performed so well with the present data set, consisting mainly of data from noncoding spacer regions, suggests that it deserves a wider consideration in model based phylogenetic inference. Repeatability of phylogenetic results was very difficult to obtain with the more parameter rich models, and analyses with identical settings often supported different topologies. Overparameterization may be the reason why the MCMC did not sample from the posterior distribution in these cases. The problem could, however, be overcome by using less parameter rich evolutionary models, and adjusting the MCMC settings. The phylogenetic results showed that two taxa, previously thought to belong in Sapotoideae, are not part of this group. Eberhardtia aurata is the sister of the two major Sapotaceae clades, Chrysophylloideae and Sapotoideae, and Neohemsleya usambarensis belongs in Chrysophylloideae. Within Sapotoideae two clades, Sideroxyleae and Sapoteae, were strongly supported. Bayesian analysis of the character history of some floral morphological traits showed that the ancestral type of flower in Sapotoideae may have been characterized by floral parts (sepals, petals, stamens, and staminodes) in single whorls of five, entire corolla lobes, and seeds with an adaxial hilum.

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

In this study we investigated the within- and between-population genetic variation using microsatellite markers and quantitative traits of the shea tree, Vitellaria paradoxa, an important agroforestry tree species of the Sudano-Sahelian region in Africa. Eleven populations were sampled across Mali and in northern Côte d'Ivoire. Leaf size and form and growth traits were measured in a progeny test at the nursery stage. Eight microsatellites were used to assess neutral genetic variation. Low levels of heterozygosity were recorded (1.6-3.0 alleles/locus; H(E) = 0.25-0.42) and the fixation index (F(IS) = -0.227-0.186) was not significantly different from zero suggesting that Hardy-Weinberg equilibrium is encountered in all populations sampled. Quantitative traits exhibited a strong genetic variation between populations and between families within populations. The degree of population differentiation of the quantitative traits (Q(ST) = 0.055-0.283, Q(STmean) = 0.189) strongly exceeds that in eight microsatellite loci (F(ST) = -0.011-0.142, F(STmean) = 0.047). Global and pairwise F(ST) values were very low and not significantly different from zero suggesting agroforestry practices are amplifying gene flow (Nm = 5.07). The population means for quantitative traits and the rainfall variable were not correlated, showing variation was not linked with this climatic cline. It is suggested that this marked differentiation for quantitative traits, independent of environmental clines and despite a high gene flow, is a result of local adaptation and human selection of shea trees. This process has induced high linkage disequilibrium between underlying loci of polygenic characters.

Temporal and spatial genetic structure in Vitellaria paradoxa (shea tree) in an agroforestry system in southern Mali.

Ten microsatellite loci were used to investigate the impact of human activity on the spatial and temporal genetic structure of Vitellaria paradoxa (Sapotaceae), a parkland tree species in agroforestry systems in southern Mali. Two stands (forest and fallow) and three cohorts (adults, juveniles and natural regeneration) in each stand were studied to: (i) compare their levels of genetic diversity (gene diversity, HE; allelic richness, Rs; and inbreeding, FIS); (ii) assess their genetic differentiation (FST); and (iii) compare their levels of spatial genetic structuring. Gene diversity parameters did not vary substantially among stands or cohorts, and tests for bottleneck events were nonsignificant. The inbreeding coefficients were not significantly different from zero in most cases (FIS = -0.025 in forest and 0.045 in fallow), suggesting that the species is probably outbreeding. There was a weak decrease in F(IS) with age, suggesting inbreeding depression. Differentiation of stands within each cohort was weak (FST = 0.026, 0.0005, 0.010 for adults, juveniles and regeneration, respectively), suggesting extensive gene flow. Cohorts within each stand were little differentiated (FST = -0.001 and 0.001 in forest and fallow, respectively). The spatial genetic structure was more pronounced in fallow than in forest where adults showed no spatial structuring. In conclusion, despite the huge influence of human activity on the life cycle of Vitellaria paradoxa growing in parkland systems, the impact on the pattern of genetic variation at microsatellite loci appears rather limited, possibly due to the buffering effect of extensive gene flow between unmanaged and managed populations.