Effect of geographical discontinuity on genetic variation in Quercus ilex L. (holm oak). Evidence from enzyme polymorphism

Genome skimming and microsatellite analysis reveal contrasting patterns of genetic diversity in a rare sandhill endemic (Erysimum teretifolium, Brassicaceae)

Barriers between islands often inhibit gene flow creating patterns of isolation by distance. In island species, the majority of genetic diversity should be distributed among isolated populations. However, a self-incompatible mating system leads to higher genetic variation within populations and very little between-population subdivision. We examine these two contrasting predictions in Erysimum teretifolium, a rare self-incompatible plant endemic to island-like sandhill habitats in Santa Cruz County, California. We used genome skimming and nuclear microsatellites to assess the distribution of genetic diversity within and among eight of the 13 remaining populations. Phylogenetic analyses of the chloroplast genomes revealed a deep separation of three of the eight populations. The nuclear ribosomal DNA cistron showed no genetic subdivision. Nuclear microsatellites suggest 83% of genetic variation resides within populations. Despite this, 18 of 28 between-population comparisons exhibited significant population structure (mean F_ST = 0.153). No isolation by distance existed among all populations, however when one outlier population was removed from the analysis due to uncertain provenance, significant isolation by distance emerged (r² = 0.5611, p = 0.005). Population census size did not correlate with allelic richness as predicted on islands. Bayesian population assignment detected six genetic groupings with substantial admixture. Unique genetic clusters were concentrated at the periphery of the species’ range. Since the overall distribution of nuclear genetic diversity reflects E. tereifolium’s self-incompatible mating system, the vast majority of genetic variation could be sampled within any individual population. Yet, the chloroplast genome results suggest a deep split and some of the nuclear microsatellite analyses indicate some island-like patterns of genetic diversity. Restoration efforts intending to maximize genetic variation should include representatives from both lineages of the chloroplast genome and, for maximum nuclear genetic diversity, should include representatives of the smaller, peripheral populations.

Morphological and physiological divergences within Quercus ilex support the existence of different ecotypes depending on climatic dryness

BACKGROUND AND AIMS

Several studies show apparently contradictory findings about the functional convergence within the Mediterranean woody flora. In this context, this study evaluates the variability of functional traits within holm oak (Quercus ilex) to elucidate whether provenances corresponding to different morphotypes represent different ecotypes locally adapted to the prevaling stress levels.

METHODS

Several morphological and physiological traits were measured at leaf and shoot levels in 9-year-old seedlings of seven Q. ilex provenances including all recognized morphotypes. Plants were grown in a common garden for 9 years under the same environmental conditions to avoid possible biases due to site-specific characteristics.

KEY RESULTS

Leaf morphometry clearly separates holm oak provenances into 'ilex' (more elongated leaves with low vein density) and 'rotundifolia' (short and rounded leaves with high vein density) morphotypes. Moreover, these morphotypes represent two consistent and very contrasting functional types in response to dry climates, mainly in terms of leaf area, major vein density, leaf specific conductivity, resistance to drought-induced cavitation and turgor loss point.

CONCLUSIONS

The 'ilex' and 'rotundifolia' morphotypes correspond to different ecotypes as inferred from their contrasting functional traits. To the best of our knowledge, this is the first time that the combined use of morphological and physiological traits has provided support for the concept of these two holm oak morphotypes being regarded as two different species.

Field evidence of colonisation by Holm Oak, at the northern margin of its distribution range, during the Anthropocene period

A major unknown in the context of current climate change is the extent to which populations of slowly migrating species, such as trees, will track shifting climates. Niche modelling generally predicts substantial northward shifts of suitable habitats. There is therefore an urgent need for field-based forest observations to corroborate these extensive model simulations. We used forest inventory data providing presence/absence information from just over a century (1880-2010) for a Mediterranean species (Quercus ilex) in forests located at the northern edge of its distribution. The main goals of the study were (i) to investigate whether this species has actually spread into new areas during the Anthropocene period and (ii) to provide a direct estimation of tree migration rate. We show that Q. ilex has colonised substantial new areas over the last century. However, the maximum rate of colonisation by this species (22 to 57 m/year) was much slower than predicted by the models and necessary to follow changes in habitat suitability since 1880. Our results suggest that the rates of tree dispersion and establishment may also be too low to track shifts in bioclimatic envelopes in the future. The inclusion of contemporary, rather than historical, migration rates into models should improve our understanding of the response of species to climate change.

Ancient and current gene flow between two distantly related Mediterranean oak species, Quercus suber and Q. ilex

BACKGROUND AND AIMS

Quercus suber and Q. ilex are distantly related and their distributions partially overlap. They hybridize occasionally, but the complete replacement of Q. suber chloroplast DNA (cpDNA) by that of Q. ilex was identified in two specific geographical areas. The objective of this study was to determine whether the contrasting situation reflected current or recent geographical interspecies gene flow variation or was the result of ancient introgression.

METHODS

cpDNA PCR-RFLPs (restriction fragment length polymorphisms) and variation at ten nuclear microsatellite loci were analysed in populations of each species, in 16 morphologically intermediate individuals and the progeny of several of them. Interspecies nuclear introgression was based on individual admixture rates using a Bayesian approach with no a priori species assignment, and on a maximum-likelihood (ML) method, using allele frequencies in the allopatric populations of each species as controls. Gene flow was compared specifically between populations located within and outside the specific areas.

KEY RESULTS

High interspecies nuclear genetic differentiation was observed, with twice the number of alleles in Q. ilex than in Q. suber. According to Bayesian assignment, approx. 1 % of individuals had a high probability of being F(1) hybrids, and bidirectional nuclear introgression affected approx. 4 % of individuals in each species. Hybrid and introgressed individuals were identified predominantly in mixed stands and may have a recent origin. Higher proportions including allospecific genes recovered from past hybridization were obtained using the ML method. Similar rates of hybridization and of nuclear introgression, partially independent of cpDNA interspecies transfer suggestive of gene filtering, were obtained in the populations located within and outside the areas of complete cpDNA replacement.

CONCLUSIONS

The results did not provide evidence for geographical variation in interspecies gene flow. In contrast, historical introgression is supported by palynological records and constitutes the more reliable origin of cpDNA replacement in specific regions.

Contrasting nuclear and cytoplasmic exchanges between phylogenetically distant oak species (Quercus suber L. and Q. ilex L.) in Southern France: inferring crosses and dynamics

Gene flow is particularly frequent in the genus Quercus (oaks), especially between closely related species. We focus here on Quercus ilex and the cork-producing Quercus suber, which occasionally hybridize although they are phylogenetically markedly separated. Morphological observations were combined with both allozymic and chloroplastic diagnostic markers to characterize hybridization and introgression and to infer their dynamics in two French regions (French Catalonia and Provence), which are separated by several hundred kilometres. Some hybrids were found in both regions, indicating recent hybridization events. As expected from previous studies, most hybrids resulted from female symbol Q. ilex x male symbol Q. suber crosses, but our data showed that the reciprocal cross is also possible. Partial independence between nuclear and chloroplastic introgression was observed in the two species. Nuclear introgression was limited in both species and both regions, with no preferred direction. In Provence, chloroplastic introgression was very rare in both species. Conversely, all Q. suber individuals from French Catalonia were introgressed by Q. ilex chlorotypes. This might be explained by introgression in the Iberian Peninsula antedating the first occurrence of the two species in French Catalonia. We also observed a new chlorotype that was created locally, and was exchanged between the two species. However, the two species still remain genetically differentiated. The dynamics and complexity of exchanges and the factors determining them (including human management of Q. suber) are discussed.

Variation in the holm oak leaf proteome at different plant developmental stages, between provenances and in response to drought stress

Major proteins of the holm oak leaf proteome have been previously identified using a combination of 2-DE, MS analysis and BLAST similarity search (Jorge et al., Proteomics 2005, 5, 222-234). That study, conducted with field samples from mature trees, revealed the existence of a great variability in the 2-DE protein map, with qualitative as well as quantitative changes, both analytical and biological. A similar study has been carried out with 2-year-old seedlings to analyze and study: (i) changes in the 2-DE protein profile at different tree developmental stages; (ii) the 2-DE protein map variability between three different Spanish provenances; and (iii) variations in the 2-DE protein profile in response to drought stress. Although the protein profile of leaves from seedlings and mature trees was fairly similar, the biological variance found was lower in the former. In the present study, new proteins have been identified. At least four different protein spots differentiated Spanish provenances, two of them identified as an ATP synthase alpha chain, and a 2,3-bisphosphoglycerate-independent phosphoglycerate mutase. Fourteen different protein spots were qualitatively variable between well-watered and drought-stressed seedlings, with some of them corresponding to enzymes of carbohydrate and protein metabolism. Data presented indicated the mobilization of storage proteins and carbohydrates, as well as photosynthesis inhibition under drought conditions.

Endemic North African Quercus afares Pomel originates from hybridisation between two genetically very distant oak species (Q. suber L. and Q. canariensis Willd.): evidence from nuclear and cytoplasmic markers

Hybridisation is a potent force in plant evolution, although there are few reported examples of stabilised species that have been created through homoploid hybridisation. We focus here on Quercus afares, an endemic North African species that combines morphological, physiological and ecological traits of both Q. suber and Q. canariensis, two phylogenetically distant species. These two species are sympatric with Q. afares over most of its distribution. We studied two Q. afares populations (one from Algeria and one from Tunisia), as well as several populations of both Q. suber and Q. canariensis sampled both within and outside areas where these species overlap with Q. afares. A genetic analysis was conducted using both nuclear (allozymes) and chloroplastic markers, which shows that Q. afares originates from a Q. suber x Q. canariensis hybridisation. At most loci, Q. afares predominantly possesses alleles from Q. suber, suggesting that the initial cross between Q. suber and Q. canariensis was followed by backcrossing with Q. suber. Other hypotheses that can account for this result, including genetic drift, gene silencing, gene conversion and selection, are discussed. A single Q. suber chlorotype was detected, and all Q. afares individuals displayed this chlorotype, indicating that Q. suber was the maternal parent. Q. afares is genetically, morphologically and ecologically differentiated from its parental species, and can therefore be considered as a stabilised hybrid species.

Chloroplast DNA variation of Quercus rubra L. in North America and comparison with other Fagaceae

Quercus rubra is one of the most important timber and ornamental tree species from eastern North America. It is a widespread species growing under variable ecological conditions. Chloroplast DNA variation was studied by PCR-RFLP (polymerase chain reaction-restriction fragment length polymorphism) in 290 individuals from 66 populations sampled throughout the natural range. A total of 12 haplotypes were detected, with one found in 75% of the trees. Population differentiation is relatively low (G(ST) = 0.46), even when similarities between haplotypes are taken into account (N(ST) = 0.50), pointing to a weak phylogeographical structure. Furthermore, no spatial structure of genetic diversity could be detected. The genetic differentiation increased northwards, reflecting the postglacial history of Q. rubra. The unusual aspect of this study was the low level of chloroplast DNA genetic differentiation in Q. rubra compared to that typically observed in other oak species. Palynological evidence indicates that during the last glacial maximum, Q. rubra had one major distribution range with populations located relatively far to the north, resulting in only modest movement northwards when climate improved, whereas European white oaks were largely restricted to the southern European peninsulas and experienced extensive movements during the postglacial period. The contrasted geographical features and levels of tree species richness of both continents might further explain why congeneric species sharing similar life history traits have genetic structures that are so different.

High variability of chloroplast DNA in three Mediterranean evergreen oaks indicates complex evolutionary history

Chloroplast DNA variation was studied in three evergreen Quercus species (Q. suber L., Q. ilex L. and Q. coccifera L.) from the Western Mediterranean Basin using PCR-RFLP. We studied five primer pair/enzyme combinations, four of them previously used in other European Quercus, obtaining a large number of haplotypes (81) grouped in three main types (suber type, ilex-coccifera I type and ilex-coccifera II type). Such level of haplotype diversity is higher than previously reported for the genus. Remarkable differences in haplotype richness between species have been found. Q. ilex and Q. coccifera usually share the same haplotypes, while a number of Q. suber populations possesses variants of the ilex-coccifera I type. This fact is interpreted as a result of genetic introgression between Q. suber and Q. ilex. Reproductive factors that could determine this exchange are discussed, as well as the influence of different species histories on the present structure of evergreen Quercus in the Western Mediterranean Basin.

Genetic variation and differentiation of populations within the Quercus affinis – Quercus laurina (Fagaceae) complex analyzed with RAPD markers

The population genetics of two hybridizing Mexican red oaks, Quercus affinis Schweid. and Quercus laurina Humb. & Bonpl., was investigated with 54 randomly amplified polymorphic DNA (RAPD) markers scored in 415 individuals from 16 populations representing the distribution area of the two species and a probable secondary hybrid zone. Genetic relationships among populations, depicted in a unweighted pair group method with arithmetic averaging (UPGMA) dendrogram, were largely incongruent with the morphological classification of populations as Q. affinis-like or Q. laurina-like that was obtained in previous studies. In contrast, the two main population clusters in the UPGMA dendrogram corresponded to the location of populations in two distinct geographical areas: southwestern and northeastern. A Mantel test confirmed a significant association between geographic and genetic distances among populations. Analyses of molecular variance (AMOVA) indicated that most genetic variation is contained within populations (84%), while 10.5% (P < 0.0001) is among populations, and 5.1% (P = 0.007) is between the two morphological groups. Differentiation between the southwestern and northeastern geographical groups (as recognized by the UPGMA), was 7.8% (P < 0.0001). The incongruence between genetic and phenotypic patterns suggests that introgression of neutral markers has been considerable between the two species in the hybrid zone, while morphological differentiation has remained comparatively stable.Key words: hybridization, population genetics, Quercus, RAPD markers.

Phylogeographical variation of chloroplast DNA in holm oak (Quercus ilex L.)

Variation in the lengths of restriction fragments (RFLPs) of the whole chloroplast DNA molecule was studied in 174 populations of Quercus ilex L. sampled over the entire distribution of this evergreen and mainly Mediterranean oak species. By using five endonucleases, 323 distinct fragments were obtained. From the 29 and 17 cpDNA changes identified as site and length mutations, respectively, 25 distinct chlorotypes were distinguished, mapped and treated cladistically with a parsimony analysis, using as an outgroup Q. alnifolia Poech, a closely related evergreen oak species endemic to Cyprus where Q. ilex does not grow. The predominant role of Q. ilex as maternal parent in hybridization with other species was reflected by the occurrence of a single very specific lineage of related chlorotypes, the most ancestral and recent ones being located in the southeastern and in the northwestern parts of the species' geographical distribution, respectively. The lineage was constituted of two clusters of chlorotypes observed in the 'ilex' morphotyped populations of the Balkan and Italian Peninsulas (including the contiguous French Riviera), respectively. A third cluster was divided into two subclusters identified in the 'rotundifolia' morphotyped populations of North Africa, and of Iberia and the adjacent French regions, respectively. Postglacial colonization probably started from three distinct southerly refugia located in each of the three European peninsulas, and a contact area between the Italian and the Iberian migration routes was identified in the Rhône valley (France). Chlorotypes identical or related to those of the Iberian cluster were identified in the populations from Catalonia and the French Languedoc region, which showed intermediate morphotypes, and in the French Atlantic populations which possessed the 'ilex' morphotype, suggesting the occurrence of adaptive morphological changes in the northern part of the species' distribution.

Allozyme characterisation of four Mediterranean evergreen oak species

Polymorphism at 11 enzyme loci was used to compare the four evergreen oak species of the Mediterranean Basin and to establish their taxonomical relationships. As a comparison, two evergreen oak species from Tibet, located in the primary distribution area of evergreen oaks, were analysed at the same loci. Cladistic analyses led to the separation of a cluster of four species, namely the three Mediterranean Q. ilex, Q. coccifera and Q. alnifolia and the Tibetan Q. aquifolioides. The other Tibetan species, Q. semicarpifolia, and Q. suber, from the Western Mediterranean Basin, were very distinct genetically from the four other species. These results were not consistent with previous taxonomical treatment, based on morphology alone, which classified Q. ilex and Q. coccifera in "ilex group" (=subgenus Sclerophyllodrys O. Schwartz), and the other four species in "cerris group" (subgenus Cerris, according to Schwartz). Allozyme variation has thus provided useful information to clarify the taxonomy of evergreen oaks.

Random amplified polymorphic DNA diversity of marginal and central populations in Pinus contorta subsp. latifolia

Fifteen populations of lodgepole pine (Pinus contorta subsp. latifolia) were surveyed for diversity across 52 random amplified polymorphic DNAs (RAPDs). The objective was to compare single-locus and multilocus structures in four marginal, three intermediate, and eight central populations. Single-locus estimates indicated average observed and expected heterozygosity to be 0.19 and 0.17, respectively. When these estimates were split into population categories, a clear trend of increasing diversity was detected in the direction of marginal to central populations. F-statistics indicated an excess of heterozygotes, with FIS ranging from -0.08 for marginal populations to -0.15 for central populations and averaging -0.12 over 15 populations. The estimates of FST decreased towards the margins of the species range, indicating increased population differentiation. Forty-nine of 52 RAPDs tested neutral in the Ewens-Watterson analysis. Multilocus analysis showed significant two-locus and high-order gametic disequilibria in all 15 populations. The most prominent components of the two-locus analysis were the variance of disequilibrium (VD, 46.2%) and the multilocus Wahlund effect (31.9%). This high value for VD indicated that founder effects could explain much of the observed multilocus associations. When analyzed by population categories, the VD showed a decreasing trend indicating that variation due to founder effects was more prominent in marginal populations. The two-locus Wahlund effect (WC) that is characteristic of strong population subdivision was highest in the central populations. This indicated significant levels of gene flow between populations with different allelic combinations.Key words: multilocus genetic structure, central and marginal populations, RAPD, Pinus contorta subsp. latifolia.

Population differentiation in mediterranean plants: insights into colonization history and the evolution and conservation of endemic species

Colonization and isolation are critical events in the evolutionary dynamics of plant populations. In this paper I review how spatial population structure of genetic markers provides insights into the evolutionary significance of episodes of colonization and isolation in the Mediterranean flora. I use as themes to structure my review the following topics: spatial structure induced by historical associations among populations of widespread species; population differentiation in relation to the evolution of closely related species with disjunct distributions; the potential effect of founder events during colonization on character evolution; and the conservation implications of spatial population structure. My review illustrates that the Mediterranean flora is full of examples that provide key insights into such evolutionary and conservation issues.