Miguel Navascués

The complex history of the olive tree: from Late Quaternary diversification of Mediterranean lineages to primary domestication in the northern Levant

The location and timing of domestication of the olive tree, a key crop in Early Mediterranean societies, remain hotly debated. Here, we unravel the history of wild olives (oleasters), and then infer the primary origins of the domesticated olive. Phylogeography and Bayesian molecular dating analyses based on plastid genome profiling of 1263 oleasters and 534 cultivated genotypes reveal three main lineages of pre-Quaternary origin. Regional hotspots of plastid diversity, species distribution modelling and macrofossils support the existence of three long-term refugia; namely the Near East (including Cyprus), the Aegean area and the Strait of Gibraltar. These ancestral wild gene pools have provided the essential foundations for cultivated olive breeding. Comparison of the geographical pattern of plastid diversity between wild and cultivated olives indicates the cradle of first domestication in the northern Levant followed by dispersals across the Mediterranean basin in parallel with the expansion of civilizations and human exchanges in this part of the world.

Chloroplast microsatellites: measures of genetic diversity and the effect of homoplasy

Chloroplast microsatellites have been widely used in population genetic studies of conifers in recent years. However, their haplotype configurations suggest that they could have high levels of homoplasy, thus limiting the power of these molecular markers. A coalescent‐based computer simulation was used to explore the influence of homoplasy on measures of genetic diversity based on chloroplast microsatellites. The conditions of the simulation were defined to fit isolated populations originating from the colonization of one single haplotype into an area left available after a glacial retreat. Simulated data were compared with empirical data available from the literature for a species of Pinus that has expanded north after the Last Glacial Maximum. In the evaluation of genetic diversity, homoplasy was found to have little influence on Neis unbiased haplotype diversity (HE) while Goldsteins genetic distance estimates () were much more affected. The effect of the number of chloroplast microsatellite loci for evaluation of genetic diversity is also discussed.

Elevated substitution rate estimates from ancient DNA: model violation and bias of Bayesian methods

The increasing ability to extract and sequence DNA from noncontemporaneous tissue offers biologists the opportunity to analyse ancient DNA (aDNA) together with modern DNA (mDNA) to address the taxonomy of extinct species, evolutionary origins, historical phylogeography and biogeography. Perhaps more exciting are recent developments in coalescence‐based Bayesian inference that offer the potential to use temporal information from aDNA and mDNA for the estimation of substitution rates and divergence dates as an alternative to fossil and geological calibration. This comes at a time of growing interest in the possibility of time dependency for molecular rate estimates. In this study, we provide a critical assessment of Bayesian Markov chain Monte Carlo (MCMC) analysis for the estimation of substitution rate using simulated samples of aDNA and mDNA. We conclude that the current models and priors employed in Bayesian MCMC analysis of heterochronous mtDNA are susceptible to an upward bias in the estimation of substitution rates because of model misspecification when the data come from populations with less than simple demographic histories, including sudden short‐lived population bottlenecks or pronounced population structure. However, when model misspecification is only mild, then the 95% highest posterior density intervals provide adequate frequentist coverage of the true rates.

Chloroplast microsatellites reveal colonization and metapopulation dynamics in the Canary Island pine

Chloroplast microsatellites are becoming increasingly popular markers for population genetic studies in plants, but there has been little focus on their potential for demographic inference. In this work the utility of chloroplast microsatellites for the study of population expansions was explored. First, we investigated the power of mismatch distribution analysis and the FS test with coalescent simulations of different demographic scenarios. We then applied these methods to empirical data obtained for the Canary Island pine (Pinus canariensis). The results of the simulations showed that chloroplast microsatellites are sensitive to sudden population growth. The power of the FS test and accuracy of demographic parameter estimates, such as the time of expansion, were reduced proportionally to the level of homoplasy within the data. The analysis of Canary Island pine chloroplast microsatellite data indicated population expansions for almost all sample localities. Demographic expansions at the island level can be explained by the colonization of the archipelago by the pine, while population expansions of different ages in different localities within an island could be the result of local extinctions and recolonization dynamics. Comparable mitochondrial DNA sequence data from a parasite of P. canariensis, the weevil Brachyderes rugatus, supports this scenario, suggesting a key role for volcanism in the evolution of pine forest communities in the Canary Islands.

Mutation rate estimates for 110 Y-chromosome STRs combining population and father–son pair data

Y-chromosome microsatellites (Y-STRs) are typically used for kinship analysis and forensic identification, as well as for inferences on population history and evolution. All applications would greatly benefit from reliable locus-specific mutation rates, to improve forensic probability calculations and interpretations of diversity data. However, estimates of mutation rate from father–son transmissions are available for few loci and have large confidence intervals, because of the small number of meiosis usually observed. By contrast, population data exist for many more Y-STRs, holding unused information about their mutation rates. To incorporate single locus diversity information into Y-STR mutation rate estimation, we performed a meta-analysis using pedigree data for 80 loci and individual haplotypes for 110 loci, from 29 and 93 published studies, respectively. By means of logistic regression we found that relative genetic diversity, motif size and repeat structure explain the variance of observed rates of mutations from meiosis. This model allowed us to predict locus-specific mutation rates (mean predicted mutation rate 2.12 × 10−3, SD=1.58 × 10−3), including estimates for 30 loci lacking meiosis observations and 41 with a previous estimate of zero. These estimates are more accurate than meiosis-based estimates when a small number of meiosis is available. We argue that our methodological approach, by taking into account locus diversity, could be also adapted to estimate population or lineage-specific mutation rates. Such adjusted estimates would represent valuable information for selecting the most reliable markers for a wide range of applications.

Combining contemporary and ancient DNA in population genetic and phylogeographical studies

The analysis of ancient DNA in a population genetic or phylogeographical framework is an emerging field, as traditional analytical tools were largely developed for the purpose of analysing data sampled from a single time point. Markov chain Monte Carlo approaches have been successfully developed for the analysis of heterochronous sequence data from closed panmictic populations. However, attributing genetic differences between temporal samples to mutational events between time points requires the consideration of other factors that may also result in genetic differentiation. Geographical effects are an obvious factor for species exhibiting geographical structuring of genetic variation. The departure from a closed panmictic model require researchers to either exploit software developed for the analysis of isochronous data, take advantage of simulation approaches using algorithms developed for heterochronous data, or explore approximate Bayesian computation. Here, we review statistical approaches employed and available software for the joint analysis of ancient and modern DNA, and where appropriate we suggest how these may be further developed.

CpDNA‐based species identification and phylogeography: application to African tropical tree species

Despite the importance of the African tropical rainforests as a hotspot of biodiversity, their history and the processes that have structured their biodiversity are understood poorly. With respect to past demographic processes, new insights can be gained through characterizing the distribution of genetic diversity. However, few studies of this type have been conducted in Central Africa, where the identification of species in the field can be difficult. We examine here the distribution of chloroplast DNA (cpDNA) diversity in Lower Guinea in two tree species that are difficult to distinguish, Erythrophleum ivorense and Erythrophleum suaveolens (Fabaceae). By using a blind‐sampling approach and comparing molecular and morphological markers, we first identified retrospectively all sampled individuals and determined the limits of the distribution of each species. We then performed a phylogeographic study using the same genetic data set. The two species displayed essentially parapatric distributions that were correlated well with the rainfall gradient, which indicated different ecological requirements. In addition, a phylogeographic structure was found for E. suaveolens and, for both species, substantially higher levels of diversity and allelic endemism were observed in the south (Gabon) than in the north (Cameroon) of the Lower Guinea region. This finding indicated different histories of population demographics for the two species, which might reflect different responses to Quaternary climate changes. We suggest that a recent period of forest perturbation, which might have been caused by humans, favoured the spread of these two species and that their poor recruitment at present results from natural succession in their forest formations.

Estimating gametic introgression rates in a risk assessment context: a case study with Scots pine relicts.

The estimation of recent gene immigration is fundamental to a wide range of evolutionary and conservation studies. In a risk assessment context, gene flow estimation procedures are needed that are both accurate and readily amenable to formal evaluation of statistical uncertainty. However, genetic methods for estimating recent migration rates that are specific and have been thoroughly evaluated are scarce. Here we use an original and straightforward maximum-likelihood method to estimate recent uniparental gametic immigration from non-local plantations into an endangered population of the Iberian relict pine variety Pinus sylvestris var. nevadensis D. H. Christ. Our approach is not intended to ascertain population membership of individuals, but rather to obtain accurate immigration rate estimates with reliable confidence limits. We found very high (∼40%) pollen introgression at the seed-crop level into the Scots pine relict, and substantial (10–15%) male gametic introgression among naturally regenerated recruits. Using numerical simulation, we show that our method yields uniparental gametic immigration estimates that are expected to be virtually unbiased and usually accurate under our sampling conditions. Among four tested methods to estimate the confidence intervals for immigration estimates, the profile-likelihood method was the best, as it outperformed bootstrapping procedures and yielded coverage close to nominal limits under different sample sizes and migration rates. This study presents a method by which researchers can facilitate decision making within a gene flow risk assessment context.

Natural recovery of genetic diversity by gene flow in reforested areas of the endemic Canary Island pine, Pinus canariensis

The endemic pine, Pinus canariensis, forms one of the main forest ecosystems in the Canary Islands. In this archipelago, pine forest is a mosaic of natural stands (remnants of past forest overexploitation) and artificial stands planted from the 1940s. The genetic makeup of the artificially regenerated forest is of some concern. The use of reproductive material with uncontrolled origin or from a reduced number of parental trees may produce stands ill adapted to local conditions or unable to adapt in response to environmental change. The genetic diversity within a transect of reforested stands connecting two natural forest fragments has been studied with nuclear and chloroplast microsatellites. Little genetic differentiation and similar levels of genetic diversity to the surrounding natural stands were found for nuclear markers. However, chloroplast microsatellites presented lower haplotype diversity in reforested stands, and this may be a consequence of the lower effective population size of the chloroplast genome, meaning chloroplast markers have a higher sensitivity to bottlenecks. Understory natural regeneration within the reforestation was also analysed to study gene flow from natural forest into artificial stands. Estimates of immigration rate into artificially regenerated forest were high (0.68-0.75), producing a significant increase of genetic diversity (both in chloroplast and nuclear microsatellites), which indicates the capacity for genetic recovery for P. canariensis reforestations surrounded by larger natural stands.

Recent population decline and selection shape diversity of taxol-related genes.

Taxanes are defensive metabolites produced by Taxus species (yews) and used in anticancer therapies. Despite their medical interest, patterns of natural diversity in taxane‐related genes are unknown. We examined variation at five main genes of Taxus baccata in the Iberian Peninsula, a region where unique yew genetic resources are endangered. We looked at several gene features and applied complementary neutrality tests, including diversity/divergence tests, tests solely based on site frequency spectrum (SFS) and Zeng’s compound tests. To account for specific demography, microsatellite data were used to infer historical changes in population size based on an Approximate Bayesian Computation (ABC) approach. Polymorphism‐divergence tests pointed to positive selection for genes TBT and TAT and balancing selection for DBAT. In addition, neutrality tests based on SFS found that while a recent reduction in population size may explain most statistics’ values, selection may still be in action in genes TBT and DBAT, at least in some populations. Molecular signatures on taxol genes suggest the action of frequent selective waves with different direction or intensity, possibly related to varying adaptive pressures produced by the host–enemy co‐evolution on defence‐related genes. Such natural selection processes may have produced taxane variants still undiscovered.