Genetic structure of gall oak (<Emphasis Type="Italic">Quercus infectoria</Emphasis>) characterized by nuclear and chloroplast SSR markers

Quercus infectoria, commonly known as gall oak, is a small shrub found in Iran. Unfortunately, it is subjected to genetic erosion, and so, its conservation and evaluation are desirable. Thus, in the current research, 16 microsatellite primer pairs (seven nuclear simple sequence repeats (nSSRs) and nine chloroplast simple sequence repeats (cpSSRs)) were used in an attempt to assess the genetic diversity of 121 individuals of Q. infectoria belonging to 11 populations from three provinces in northern Zagros forests of Iran. In total, 69 alleles of nSSR and 18 alleles of cpSSR were detected among the individuals. The results of the overall analysis of molecular variance based on nSSRs indicated that 89.00% of the variation was due to differences within populations and 11.00% occurred among populations, while according to cpSSRs, 94.00% of the variation resided among populations, and only 6.00% could be attributed to variation within populations. A higher genetic differentiation of Q. infectoria populations was found according to cpSSR data in comparison to nSSR data. Cophenetic correlation coefficient values were statistically insignificant between nSSR and cpSSR data. The unweighted pair group method with arithmetic mean and Bayesian cluster analyses grouped the studied individuals into two main clusters based on both nSSR and cpSSR data. nSSR data could not completely clustered individuals next each other according to their geographical collection area. Information detailed by nSSR loci revealed that north-Zagros gall oak preserves average levels of genetic diversity at the species level, high level of within-population genetic diversity, and moderate level of genetic variation among populations. The present results provide valuable data for in situ or ex situ conservation and utilization of the studied germplasm.     

Development and Characterization of Genomic and Gene-Based Microsatellite Markers in North American Red Oak Species

Oaks (Quercus: Fagaceae) are ecological and economic keystones of many forested ecosystems but effective genetic management strategies are hindered by high levels of phenotypic plasticity within species and frequent hybridization among them. These same features, however, make oak communities suited for the study of speciation, hybridization, and genetic adaptation. Efforts to develop new and to adapt existing genomic resources to less-studied members of this genus should not only improve oak conservation and management but also aid the study of fundamental evolutionary processes. Here, we present a suite of 27 highly polymorphic simple sequence repeat (SSR) markers tested in four North American red oak (Quercus section Lobatae) species: Q. rubra, Q. ellipsoidalis, Q. coccinea, and Q. velutina. Five markers are genomic SSRs (gSSRs) — four novel and one previously transferred from Q. petraea — and 22 are gene-based SSRs derived from Q. robur and Q. petraea expressed sequence tags (EST-SSRs). Overall, levels of polymorphism detected with these primer pairs were high, with gene diversity (H _e) averaging 0.66 across all loci in natural populations. In addition, we show that EST-SSR markers may have the potential to detect divergent selection at stress-resistance candidate genes among closely related oak species.     

Genetic mapping of EST-derived simple sequence repeats (EST-SSRs) to identify QTL for leaf morphological characters in a Quercus robur full-sib family

The availability of genomic resources such as expressed sequence tag-derived simple sequence repeat (EST-SSR) markers in adaptive genes with high transferability across related species allows the construction of genetic maps and the comparison of genome structure and quantitative trait loci (QTL) positions. In the present study, genetic linkage maps were constructed for both parents of a Quercus robur × Q. robur ssp. slavonica full-sib pedigree. A total of 182 markers (61 AFLPs, 23 nuclear SSRs, 98 EST-SSRs) and 172 markers (49 AFLPs, 21 nSSRs, 101 EST-SSRs, 1 isozyme) were mapped on the female and male linkage maps, respectively. The total map length and average marker spacing were 1,038 and 5.7 cM for the female map and 998.5 and 5.8 cM for the male map. A total of 68 nuclear SSRs and EST-SSRs segregating in both parents allowed to define homologous linkage groups (LG) between both parental maps. QTL for leaf morphological traits were mapped on all 12 LG at a chromosome-wide level and on 6 LG at a genome-wide level. The phenotypic effects explained by each single QTL ranged from 4.0 % for leaf area to 15.8 % for the number of intercalary veins. QTL clusters for leaf characters that discriminate between Q. robur and Quercus petraea were mapped reproducibly on three LG, and some putative candidate genes among potentially many others were identified on LG3 and LG5. Genetic linkage maps based on EST-SSRs can be valuable tools for the identification of genes involved in adaptive trait variation and for comparative mapping.     

Primers for 52 polymorphic regions in the Quercus rubra chloroplast, 47 of which amplify across 11 tracheophyte clades

Postglacial migration studies in Quercus rubra L. (northern red oak) are hampered by low levels of population differentiation in the widely used universal chloroplast primers. We sequenced the large single copy (LSC) regions of the Q. rubra and Quercus ellipsoidalis chloroplasts to enable us to query additional regions for future studies on migration and speciation. Using 454 sequencing of long-range PCR amplicons and Sanger sequencing for gap closure, we report 65 coding sequences from Q. rubra and 59 from Q. ellipsoidalis. Comparison of our de novo assembly of the LSC region sequence for Q. rubra to Q. rubra chloroplast sequence (NCBI Reference Sequence: NC_020152.1) from a different tree revealed 106 polymorphisms, all within intergenic regions, that can serve as tools for postglacial migration studies and taxonomic studies within the Lobatae. Sequence alignment for the 59 complete coding regions in common for theQ. rubrachloroplast reference sequence, our Q. rubra sequence and our Q. ellipsoidalis sequence revealed no sequence polymorphisms and no indels. We also report the 52 primer pairs we used for gap closure, including 53 new primer pairs not previously reported. We tested these 52 primer pairs against 11 species representing the Tracheophyta and detected 47 that produced amplicons in all 11 species. The new universal primers we have identified provide additional tools for resolving the taxonomic relationships among the congeneric taxa of forest trees in the temperate and subtropical forests of the Northern Hemisphere.     

Effect of a red oak species gradient on genetic structure and diversity of Quercus castanea (Fagaceae) in Mexico

Incipient reproductive barriers are a common characteristic of oak species. Disruption of these barriers promotes changes in diversity and genetic structure of the species involved. Quercus castanea is a red oak with wide geographic distribution in Mexico, which presents atypically high morphological variability when it occurs in sympatry with other red oak species, suggesting that hybridization may explain the observed variation. We tested if the genetic structure and diversity levels of Q. castanea are related to the number of red oak species growing in sympatry. In total, 14 microsatellite (SSRs) primers (six nSSRs and eight cpSSRs) were used in 120 Q. castanea individuals (20/site) belonging to six populations, where the number of red oak species associated varied from zero to five. Results showed a positive and significant relationship between the genetic diversity of Q. castanea and the number of red oak species growing in sympatry, regardless of the marker type or the parameter of genetic diversity analyzed. Also, we found a higher genetic differentiation of Q. castanea populations using cpSSRs in comparison with nSSRs. Our results suggest that temperate forests with high red oaks species richness co-dominated by Q. castanea promote the increase in this species genetic diversity. From a conservation perspective, high genetic diversity levels of foundation species such as Q. castanea may have positive cascade effects extending to other species in the community.     

Genetic variation and differentiation in populations of Japanese emperor oak <Emphasis Type="Italic">Quercus dentata</Emphasis> Thunb. and Mongolian oak <Emphasis Type="Italic">Quercus mongolica</Emphasis> Fisch. ex Ledeb.) in the south of the Russian Far East

Allozyme variation of Japanese emperor oak Quercus dentata Thunb. and Mongolian oak Quercus mongolica Fisch. ex Ledeb.) was examined in 11 populations of these species from southern Primorye. Analysis of 18 loci showed that in these populations, 66.7% of genes of Mongolian oak are polymorphic, the number of alleles per locus being 2.28. In Japanese emperor oak, these parameters were respectively 59.8 and 2.36%. The observed heterozygosity in both species did not differ from the expected values, constituting 0.198 in Japanese emperor oak and 0.161 in Mongolian oak. The Nei’s distance between Quercus mongolica and Q. dentata was 0.053. The main contribution to the species differentiation is made by loci Lap-1, Idh-1, Pgm-1, and Ndh-1, although none of them is diagnostic. The Japanese emperor oak and Mongolian oak populations formed separate clusters in an UPGMA dendrogram. Taking into account morphological and ecological isolation of the species, we suggest that the differences in gene frequencies may reflect selection acting in favor of adaptation to different growth conditions and counteracting between-species hybridization.     

Genetic Differentiation of Pedunculate Oak Quercus robur L. in the European Part of Russia Based on RAPD Markers

Using randomly amplified polymorphic DNA markers (RAPD), genetic variation and differentiation in four populations of pedunculate oak Quercus robur L. were examined. The populations occupy a large part of the Quercus robur range in the European Russia (Voronezh and Novgorod oblasts; Republics of Mordovia and Bashkortostan). With each of six random primers (A02, A09, A17, B01, B08, B11), 96 DNA samples were analyzed by PCR. In all, 48 putative polymorphic RAPD loci were detected. We failed to reveal population-specific DNA fragments for any primer although the frequencies of 14 fragments were significantly different among populations. The oak populations studied exhibited high variability: 73–90% of genes were polymorphic and the effective allele number was about 1.4. The total genetic variation varied from 0.202 (Vor) to 0.245 (Nov), which corresponded to the estimates for populations of this species from Central and Western Europe. The populations examined showed low among-population differentiation (G _ST = 0.098); gene flow N _e m was 4.61. The proportion of among-population variation of the RAPD loci studied accounted for 7% of the total variability; more than 93% of the total variability was explained by individual and within-population variation.     

Hybridization and divergence in multi‐species oak (Quercus) communities

Oaks (Quercus: Fagaceae) commonly interbreed yet retain their morphological, genetic and ecological distinctiveness. Post‐zygotic isolation mechanisms, such as ecologically dependent selection on adaptive loci, may therefore limit introgression. To test this hypothesis, we quantified hybridization and genetic divergence across the contact zone of four red oaks (Quercus section Lobatae) in the Great Lakes region of North America using a suite of 259 amplified fragment length polymorphisms and 27 genic and genomic microsatellite markers. First, we identified hybrids using genetic structure analysis and confirmed the reliability of our assignments via simulations. Then, we identified candidate loci for species maintenance with three complementary tests for selection and obtained partial gene sequences linked to an outlier locus and three other loci. We detected evidence of recent hybridization among all species and considerable gene flow between Q. ellipsoidalis and Q. velutina. Overall, c. 20% of Q. velutina had recent ancestry from Q. ellipsoidalis, whereas nearly 30% of Q. ellipsoidalis had a Q. velutina ancestor. Most loci were negligibly to weakly differentiated among species, but two gene‐linked microsatellites deviated significantly from neutral expectations in multiple, complementary outlier tests. Both outlier loci were located in the same 15‐cM bin on an existing Q. robur linkage map, a region under divergent selection in other oak species. Adaptive loci in this highly differentiated genomic region may contribute to ecological divergence among species and limit introgression.     

Genetic differentiation of pedunculate oak Quercus robur L. in the European part of Russia based on RAPD markers

Using randomly amplified polymorphic DNA markers (RAPD), genetic variation and differentiation in four populations of pedunculate oak Quercus robur L. were examined. The populations occupy a large part of the Quercus robur range in the European Russia (Voronezh and Novgorod oblasts; Republics of Mordovia and Bashkortostan). With each of six random primers (A02, A09, A17, B01, B08, B11), 96 DNA samples were analyzed by PCR. In all, 48 putative polymorphic RAPD loci were detected. We failed to reveal population-specific DNA fragments for any primer although the frequencies of 14 fragments were significantly different among populations. The oak populations studied exhibited high variability: 73-90% of genes were polymorphic and the effective allele number was about 1.4. The total genetic variation varied from 0.202 (Vor) to 0.245 (Nov), which corresponded to the estimates for populations of this species from Central and Western Europe. The populations examined showed low among-population differentiation (GST = 0.098); gene flow Nem was 4.61. The proportion of among-population variation of the RAPD loci studied accounted for 7% of the total variability; more than 93% of the total variability was explained by individual and within--population variation.     

Pyrosequencing of the northern red oak (Quercus rubra L.) chloroplast genome reveals high quality polymorphisms for population management

Given the low intraspecific chloroplast diversity detected in northern red oak (Quercus rubra L.), more powerful genetic tools are necessary to accurately characterize Q. rubra chloroplast diversity and structure. We report the sequencing, assembly, and annotation of the chloroplast genome of northern red oak via pyrosequencing and a combination of de novo and reference-guided assembly (RGA). Chloroplast DNA from 16 individuals was separated into four MID-tagged pools for a Genome Sequencer 20 quarter-run (Roche Life Sciences, Indianapolis, IN, USA). A four-step assembly method was used to generate the Q. rubra chloroplast consensus sequence: (1) reads were assembled de novo into contigs, (2) de novo contigs were aligned to a reference genome and merged to produce a consensus sequence, (3) the consensus sequence was aligned to the reference sequence and gaps between contigs were filled with reference sequence to generate a "pseudoreference", and (4) reads were mapped to the pseudoreference using RGA to generate the draft chloroplast genome. One hundred percent of the pseudoreference sequence was covered with a minimum coverage of 2× and an average coverage of 43.75×. The 161,304-bp Q. rubra chloroplast genome draft sequence contained 137 genes and one rps19 pseudogene. The sequence was compared to that of Quercus robur and Q. nigra with 951 and 186 insertion/deletion or SNP polymorphisms detected, respectively. A total of 51 intraspecific polymorphisms were detected among four northern red oak individuals. The fully sequenced and annotated Q. rubra chloroplast genome containing locations of interspecific and intraspecific polymorphisms will be essential for studying population differentiation, phylogeography, and evolutionary history of this species as well as meeting management goals such as monitoring reintroduced populations, tracking wood products, and certifying seed lots and forests.     

A range wide geographic pattern of genetic diversity and population structure of Castanea mollissima populations inferred from nuclear and chloroplast microsatellites

Understanding geographical pattern of genetic diversity and population structure is of great importance for formulating conservation and utilization strategies. In this study, we investigated the genetic diversity and population structure of 28 natural populations of Castanea mollissima in China using eight nuclear and six chloroplast microsatellite makers (nSSRs and cpSSRs). Populations from central China harbored the highest genetic diversity at both nSSR and cpSSR markers (nSSR: H _E?=?0.705; cpSSR: H?=?0.461). The standardized measure of genetic differentiation estimated as G′ _ST was 0.447 for nSSR and 0.803 for cpSSR, respectively. The G′_ST-based pollen to seed flow ratio is 3.043, indicating that pollen flow is not extensive among C. mollissima populations. No obvious population genetic structure by geographical locations was found by STRUCTURE analysis based on nSSR data, and similarly, no signal of phylogeographic structure was detected for cpSSR analysis. Five boundaries defining zones of maximum genetic differences within the network of the C. mollissima populations were found, and the locations of those barriers were consistent with those of four mountains, i.e., Daloushan Mountain, Dabashan Mountain, Wushan Mountain, and Qingliangfeng Mountain, indicating that those mountains might act as genetic barriers obstructing the genetic exchange among natural C. mollissima populations. These results provide valuable baseline data for conservation and utilization of this species.     

Extent of nuclear genome sharing among white oak species (Quercus L. subgen. Lepidobalanus (Endl.) Oerst.) in Slovakia estimated by allozymes

European white oak species (Quercus subg. Lepidobalanus) are known to share a substantial part of the nuclear and cytoplasmic genome. In spite of this fact, they have preserved their integrity. We analyzed genetic variation at thirteen isozyme loci in nine Quercus robur, ten Q. petraea, two Q. pubescens and two mixed (two-species) populations in Slovakia to assess the extent of allele sharing. Moreover, a detailed morphometric study based on thirteen leaf traits was done in the mixed stands. A Bayesian analysis of population structure in the set of pure populations revealed the existence of three clusters, which coincided with the supposed taxonomical identity of populations. However, the distribution of posterior probabilities of data was bimodal with a second peak at six groups, what indicates a possible substructure in the gene pools of the investigated oak species. Moreover, gene admixture was close to 50% in all investigated populations. Morphometric analysis of the two mixed stands revealed a clear morphological separation of the studied oak species, sharply contrasting with a high degree of gene admixture. Differentiation is not uniformly distributed over the genome. High levels of differentiation (F _ST > 0.10) were observed at only three out of thirteen loci. The hypotheses explaining extensive allele sharing (shared ancestral polymorphisms vs. high levels of interspecific gene flow) are briefly discussed.     

High genetic variation among closely related red oak (<Emphasis Type="Italic">Quercus rubra</Emphasis>) populations in an ecosystem under metal stress: analysis of gene regulation

Metal toxicity is a major abiotic stressor of plants. It has been established that changes in genetic variation occur very rapidly in plants in response to environmental stressors such as increased levels of metals. Quercus rubra (red oak) is a pioneer species in mining regions contaminated with metals in Northern Ontario (Canada). The objectives of the study were to (1) determine the level of genetic variation in Q. rubra populations from mining damaged ecosystems using RAPD marker system and (2) assess the level of gene expression of candidate genes for nickel resistance. Total gene diversity (H_T) and the mean gene diversity among populations (H_S) were 0.22 and 0.19, respectively. The percent of polymorphic loci within populations was high ranging from 61 % (Capreol) to 72 % (Daisy Lake) despite a high level of gene flow (2.4). The population differentiation (G_ST) value was low (0.17). No significant difference was found among the contaminated and reference sites for all the genetic parameters estimated. Hence, all the Q. rubra populations from the metal-contaminated and damaged ecosystems are genetically sustainable. Moreover, this study reveals that all populations were genetically closely related with genetic distance values varying from 0.17 to 0.35. A zinc finger protein of Arabidopsis thaliana (ZAT11) gene involved in nickel resistance was differentially expressed in samples analyzed. There was a 120 times higher of ZAT11 expression in samples from metal contaminated areas of Wahnapitae Dam compared to other metal contaminated and uncontaminated sites, but no association between soil metal levels and expression of ZAT11 was established.     

Geographic structure of genetic and phenotypic variation in the hybrid zone between Quercus affinis and Q. laurina in Mexico

Analyzing the structure of hybrid zones is important for inferring their origin, dynamics and evolutionary significance. We examined the geographic structure of phenotypic and genetic variation in the contact zone between two Mexican red oaks, Quercus affinis and Q. laurina. A total of 105 individuals from seven populations were sampled along a 600‐km latitudinal gradient representing the distribution area of the two species and their contact zone. Individuals were genotyped for nine nuclear and four chloroplast DNA microsatellite loci (ncSSR and cpSSR, respectively), and characterized for several leaf and acorn traits. The cpSSR data revealed extensive haplotype sharing among populations of the two species, while a Bayesian assignment analysis based on ncSSRs identified two main genetic groups, each corresponding to one of the species, and two populations in the contact zone showing evidence of admixture. The proportion of genetic ancestry in the populations was strongly associated with latitude and showed a pattern of variation with the shape of a narrow sigmoidal cline. The variation in three of the seven phenotypic traits was partially congruent with molecular variation, while the other traits did not conform to a geographic cline but instead were correlated with environmental variables. In conclusion, the hybrid zone between the two oak species has some of the characteristics of a tension zone, but heterogeneous variation across traits suggests differential introgression and the action of extrinsic selection.     

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₁ 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.Key words: cpDNA PCR-RFLPs, nuclear microsatellite (nSSR) variation, hybridization, interspecies genetic introgression, Quercus suber, Quercus ilex     

A set of primers for the amplification of chloroplast microsatellites in Quercus

The increase in demand for the certification of oak seed lots, as well as control of the geographical origin of oak wood, has led us to develop powerful genetic markers permitting us to discriminate among provenance regions. With the aim of detecting new chloroplast variants, we have identified 17 potential cpSSRs motifs from available oak sequences and tested their variability among French oak populations. Six loci were polymorphic at the intraspecific level in Quercus petraea and Q. robur. Moreover, conservation of the primer pairs was checked on a set of 21 forest tree species and they were all shown to work well on several Quercus species, and even within Fagacaea.     

DNA methylation of Quercus robur L. plumules following cryo-pretreatment and cryopreservation

Pedunculate oak (Quercus robur) is an ecologically and economically important forest tree species which produces seeds that are classified as recalcitrant. Thus, cryopreservation of seed meristems is a method for long-term preservation of this germplasm in gene banks. During cryopreservation, many factors, such as desiccation, cryoprotection and cooling/rewarming, can induce stress in the frozen meristems. In this study, in vitro survival and the global DNA methylation level of plumules after cryoprotection, desiccation and cryostorage was evaluated. Results indicated that both desiccation and storage in liquid nitrogen have negligible influence on DNA methylation status of Q. robur plumules. These findings support the cryopreservation of plumules as an appropriate method for conservation of Q. robur germplasm.