Characterization of highly variable (GA/CT) n microsatellites in the bur oak, Quercus macrocarpa

Genetic Variability and Admixture Zones in the Italian Populations of Turkey Oak (Quercus cerris L.)

The Turkey oak (Quercus cerris L.) is widely distributed in Italy, where it is the ecologically dominant oak on sandy and acidic soil. In this work, we analysed 23 natural populations by means of eight SSR (microsatellite) markers, to obtain the first synthetic map of genetic variability for this species and to study its dispersion during the Holocene, due to the possibility that at least one refugium during the Last Glacial Maximum was in Italy. The analyses showed a good amount of genetic variability together with fair differentiation between populations, as indicated by F_ST = 0.059. A Bayesian analysis of the amount of admixture among populations revealed the presence of four putative gene pools of origin and a rough subdivision of the populations according to their geographic location, as confirmed by the spatial analysis. No evidence for the existence of putative refugial populations was found; however, this study paves the way for the planning of conservation strategies also with regard to the relationship between Turkey oak and other oak species in Italy.

Population structure in Quercus suber L. revealed by nuclear microsatellite markers

Quercus suber L. is a sclerophyllous tree species native to the western Mediterranean, a region that is considered highly vulnerable to increased temperatures and severe dry conditions due to environmental changes. Understanding the population structure and demographics of Q. suber is essential in order to anticipate whether populations at greater risk and the species as a whole have the genetic background and reproductive dynamics to enable rapid adaptation. The genetic diversity of Q. suber has been subject to different studies using both chloroplast and nuclear data, but population structure patterns remain unclear. Here, we perform genetic analyses on Q. suber using 13 nuclear microsatellite markers, and analysed 17 distinct locations across the entire range of the species. Structure analyses revealed that Q. suber may contain three major genetic clusters that likely result from isolation in refugia combined with posterior admixture and putative introgression from other Quercus species. Our results show a more complex structure scenario than previously inferred for Q. suber using nuclear markers and suggest that different southern populations contain high levels of genetic variation that may contribute to the resilience of Q. suber in a context of environmental change and adaptive pressure.

Population and Landscape Genetics Provide Insights Into Species Conservation of Two Evergreen Oaks in Qinghai-Tibet Plateau and Adjacent Regions

The combination of population and landscape genetics can facilitate the understanding of conservation strategy under the changing climate. Here, we focused on the two most diverse and ecologically important evergreen oaks: Quercus aquifolioides and Quercus spinosa in Qinghai-Tibetan Plateau (QTP), which is considered as world's biodiversity hotspot. We genotyped 1,657 individuals of 106 populations at 15 nuclear microsatellite loci throughout the species distribution range. Spatial patterns of genetic diversity were identified by mapping the allelic richness (AR) and locally common alleles (LCA) according to the circular neighborhood methodology. Migration routes from QTP were detected by historical gene flow estimation. The response pattern of genetic variation to environmental gradient was assessed by the genotype-environment association (GEA) analysis. The overall genetic structure showed a high level of intra-species genetic divergence of a strong west-east pattern. The West-to-East migration route indicated the complex demographic history of two oak species. We found evidence of isolation by the environment in Q. aqu-East and Q. spi-West lineage but not in Q. aqu-West and Q. spi-East lineage. Furthermore, priority for conservation should be given to populations that retain higher spatial genetic diversity or isolated at the edge of the distribution range. Our findings indicate that knowledge of spatial diversity and migration route can provide valuable information for the conservation of existing populations. This study provides an important guide for species conservation for two oak species by the integration of population and landscape genetic methods.

Genetic Diversity and Differentiation of Pedunculate Oak (Quercus robur L.) Populations at the Southern Margin of Its Distribution Range—Implications for Conservation

Understanding intraspecific genetic variation is one of the principal requirements for the evaluation of tree species capacity to cope with intensive climatic changes, as well as designing long-term conservation programs. Herein, we evaluated the genetic diversity and genetic structure of seven pedunculate oak (Quercus robur L.) populations, located at the southern margin of its distribution range on the Balkan Peninsula (Serbia). The objective of the study was to propose future in situ conservation measures aimed at protection of pedunculate oak adaptive and neutral genetic diversity at the species rear-edge. Genetic diversity and structure were estimated using twelve highly polymorphic simple sequence repeat (SSR) markers. The mean expected heterozygosity (He) was 0.769, allelic richness (AR) 9.63, and private allelic richness (pAR) 0.79, indicating high genetic diversity in the studied populations. Genetic differentiation among the populations was low (Fst = 0.032). Structure analysis, the unweighted pair group method with arithmetic mean (UPGMA) showed the existence of two gene pools unrelated to the populations’ area of occurrence. Taking into consideration the results of the current study and previous conservation activities on the pedunculate oak in Serbia, as well as the importance of rear-edge populations in the long-term conservation of the species genetic diversity, we suggested establishing three additional gene conservation units for securing long-term sustainability of the species.

Unraveling the genetic diversity and structure of Quercus liaotungensis population through analysis of microsatellite markers

Background

Quercus liaotungensis Koidz. is an ecologically and economically important tree species widely distributed in Northern China. However, the effective assessment, utilization, and protection of Q. liaotungensis resources remain unexplored.

Methods

In total, 120 samples obtained from 12 Q. liaotungensis populations of Northern China were investigated for genetic diversity and structure using 19 simple sequence repeat (SSR) primer pairs.

Results

The total number of alleles detected was 293, the average number of effective allele (Ne) was 6.084, the genetic differentiation coefficient (Fst) was 0.033, and the mean observed heterozygosity (Ho) and expected heterozygosity (He) were 0.690 and 0.801, respectively. Moreover, analysis of molecular variance (AMOVA) showed a 5.5% genetic variation among 12 Q. liaotungensis populations, indicating that a high level of genetic diversity and a low degree of genetic differentiation among Q. liaotungensis populations. STRUCTURE and cluster analysis divided the 12 Q. liaotungensis populations into the following three subpopulations: Bashang Plateau subpopulation (SH), Liaodong Peninsula subpopulation (NC), and Loess Plateau subpopulation (other 10 populations). The cluster analysis based on 19 climatic factors was consistent with the genetic structure. A positive correlation was found between genetic distance and geographical distance (r = 0.638, p = 0.028) by the Mantel test, and two boundaries were found among the 12 Q. liaotungensis populations by the Barrier analysis, indicating that Q. liaotungensis populations existed isolated by geographical distance and physical barrier.

Conclusion

This study suggests that geographical isolation, physical barrier, climatic types, and natural hybridization promote the formation of genetic structures, which can contribute to future protection and genetic improvement of Q. liaotungensis.

Genetic, Morphological, and Environmental Differentiation of an Arid-Adapted Oak with a Disjunct Distribution

The patterns of genetic and morphological diversity of a widespread species can be influenced by environmental heterogeneity and the degree of connectivity across its geographic distribution. Here, we studied Quercus havardii Rydb., a uniquely adapted desert oak endemic to the Southwest region of the United States, using genetic, morphometric, and environmental datasets over various geographic scales to quantify differentiation and understand forces influencing population divergence. First, we quantified variation by analyzing 10 eastern and 13 western populations from the disjunct distribution of Q. havardii using 11 microsatellite loci, 17 morphological variables, and 19 bioclimatic variables. We then used regressions to examine local and regional correlations of climate with genetic variation. We found strong genetic, morphological and environmental differences corresponding with the large-scale disjunction of populations. Additionally, western populations had higher genetic diversity and lower relatedness than eastern populations. Levels of genetic variation in the eastern populations were found to be primarily associated with precipitation seasonality, while levels of genetic variation in western populations were associated with lower daily temperature fluctuations and higher winter precipitation. Finally, we found little to no observed environmental niche overlap between regions. Our results suggest that eastern and western populations likely represent two distinct taxonomic entities, each associated with a unique set of climatic variables potentially influencing local patterns of diversity.

Quercus species divergence is driven by natural selection on evolutionarily less integrated traits

Functional traits are organismal attributes that can respond to environmental cues, thereby providing important ecological functions. In addition, an organism's potential for adaptation is defined by the patterns of covariation among groups of functionally related traits. Whether an organism is evolutionarily constrained or has the potential for adaptation is based on the phenotypic integration or modularity of these traits. Here, we revisited leaf morphology in two European sympatric white oaks (Quercus petraea (Matt.) Liebl. and Quercus robur L.), sampling 2098 individuals, across much of their geographical distribution ranges. At the phenotypic level, leaf morphology traditionally encompasses discriminant attributes among different oak species. Here, we estimated in situ heritability, genetic correlation, and integration across such attributes. Also, we performed Selection Response Decomposition to test these traits for potential differences in oak species' evolutionary responses. Based on the uncovered functional units of traits (modules) in our study, the morphological module "leaf size gradient" was highlighted among functionally integrated traits. Equally, this module was defined in both oaks as being under "global regulation" in vegetative bud establishment and development. Lamina basal shape and intercalary veins' number were not, or, less integrated within the initially defined leaf functional unit, suggesting more than one module within the leaf traits' ensemble. Since these traits generally show the greatest species discriminatory power, they potentially underwent effective differential response to selection among oaks. Indeed, the selection of these traits could have driven the ecological preferences between the two sympatric oaks growing under different microclimates.

Genetic diversity of avocado from the southern highlands of Tanzania as revealed by microsatellite markers

BACKGROUND

Avocado is an important cash crop in Tanzania, however its genetic diversity is not thoroughly investigated. This study was undertaken to explore the genetic diversity of avocado in the southern highlands using microsatellite markers. A total of 226 local avocado trees originating from seeds were sampled in eight districts of the Mbeya, Njombe and Songwe regions. Each district was considered as a population. The diversity at 10 microsatellite loci was investigated.

RESULTS

A total of 167 alleles were detected across the 10 loci with an average of 16.7 ± 1.3 alleles per locus. The average expected and observed heterozygosity were 0.84 ± 0.02 and 0.65 ± 0.04, respectively. All but two loci showed a significant deviation from the Hardy-Weinberg principle. Analysis of molecular variance showed that about 6% of the variation was partitioned among the eight geographic populations. Population F_ST pairwise comparisons revealed lack of genetic differentiation for the seven of 28 population pairs tested. The principal components analysis (PCA) and hierarchical cluster analysis showed a mixing of avocado trees from different districts. The model-based STRUCTURE subdivided the trees samples into four major genetic clusters.

CONCLUSION

High diversity detected in the analysed avocado germplasm implies that this germplasm is a potentially valuable source of variable alleles that might be harnessed for genetic improvement of this crop in Tanzania. The mixing of avocado trees from different districts observed in the PCA and dendrogram points to strong gene flow among the avocado populations, which led to population admixture revealed in the STRUCTURE analysis. However, there is still significant differentiation among the tree populations from different districts that can be utilized in the avocado breeding program.

The draft genome sequence of cork oak

Cork oak (Quercus suber) is native to southwest Europe and northwest Africa where it plays a crucial environmental and economical role. To tackle the cork oak production and industrial challenges, advanced research is imperative but dependent on the availability of a sequenced genome. To address this, we produced the first draft version of the cork oak genome. We followed a de novo assembly strategy based on high-throughput sequence data, which generated a draft genome comprising 23,347 scaffolds and 953.3 Mb in size. A total of 79,752 genes and 83,814 transcripts were predicted, including 33,658 high-confidence genes. An InterPro signature assignment was detected for 69,218 transcripts, which represented 82.6% of the total. Validation studies demonstrated the genome assembly and annotation completeness and highlighted the usefulness of the draft genome for read mapping of high-throughput sequence data generated using different protocols. All data generated is available through the public databases where it was deposited, being therefore ready to use by the academic and industry communities working on cork oak and/or related species.

Genetic variation and structure of Ubame oak, Quercus phillyraeoides, in Japan revealed by chloroplast DNA and nuclear microsatellite markers

Genetic variation and structure of Ubame oak (Quercus phillyraeoides A. Gray), a shrub tree on the Pacific coast of Japan, were examined to elucidate historical population dynamics using five chloroplast DNA (cpDNA) marker regions and 11 nuclear microsatellite loci. Three cpDNA haplotypes (A, B and C) were identified in a screen of 41 populations across the entire distribution range in Japan. Haplotype A was the most prevalent and was found in the entire range. Haplotype B was locally restricted to the Kii Peninsula, while haplotype C was restricted to the southwestern part of Japan. These haplotypes corresponded with the nuclear genetic constitution revealed by microsatellite markers. Two genetically differentiated major groups were identified by STRUCTURE analysis applied to 536 individuals from 28 populations, and they mostly corresponded with the two major cpDNA haplotypes, A and C. These populations were further divided into three geographically identified groups: group 1 in the area including Kanto and Tokai regions, the Kii Peninsula and the Muroto-misaki Cape on the Pacific coast; group 2 in the Bungo-suido Channel area; and group 3 in southern Kyushu and Okinawa. Populations in the Seto Inland Sea were divided into two groups: one was included in group 1 and appears to have originated from last glacial maximum (LGM) refugia located in the Kii Peninsula, while the other was included in group 2 and appears to have originated from LGM refugia located in southern Kyushu. These groups can be considered as conservation units for the preservation of unique seashore ecosystems, or as a seed source to foster coastal protection forests and next-generation production forests. Considerable care should be taken to protect isolated populations that may be specialized to unique local environments, such as those on the islands of Koshikijima and Izenajima.

Species Boundaries Between Three Sympatric Oak Species: Quercus aliena, Q. dentata, and Q. variabilis at the Northern Edge of Their Distribution in China

Oaks are important timber trees with wide distributions in China, but few genetic studies have been conducted on a fine scale. In this study, we seek to investigate the genetic diversity and differentiation of three sympatric oak species (Quercus aliena Blume, Quercus dentata Thunb. ex Murray, and Quercus variabilis Blume) in their northern distribution in China using 17 bi-parentally inherited nSSRs markers and five maternally inherited chloroplast DNA (cpDNA) fragments. Both the cpDNA and the nSSRs show a high level of genetic differentiation between different oak sections. The chloroplast haplotypes are clustered into two lineages. Clear species boundaries are detected between Q. variabilis and either Q. aliena or Q. dentata. The sharing of chloroplast haplotype H1 between Q. aliena and Q. dentata suggests very recent speciation and incomplete lineage sorting or introgression of H1 from one species to another. The nSSRs data indicate a complete fixation of variation within sites for all three oak species, and that extensive gene flow occurs within species whereas only limited gene flow is detected between Q. aliena and Q. dentata and nearly no gene flow can be detected between Q. aliena and Q. variabilis and between Q. dentata and Q. variabilis. Prezygotic isolation may have contributed to the species boundaries of these three sympatric oak species.

Assessing Genetic Structure in Common but Ecologically Distinct Carnivores: The Stone Marten and Red Fox

The identification of populations and spatial genetic patterns is important for ecological and conservation research, and spatially explicit individual-based methods have been recognised as powerful tools in this context. Mammalian carnivores are intrinsically vulnerable to habitat fragmentation but not much is known about the genetic consequences of fragmentation in common species. Stone martens (Martes foina) and red foxes (Vulpes vulpes) share a widespread Palearctic distribution and are considered habitat generalists, but in the Iberian Peninsula stone martens tend to occur in higher quality habitats. We compared their genetic structure in Portugal to see if they are consistent with their differences in ecological plasticity, and also to illustrate an approach to explicitly delineate the spatial boundaries of consistently identified genetic units. We analysed microsatellite data using spatial Bayesian clustering methods (implemented in the software BAPS, GENELAND and TESS), a progressive partitioning approach and a multivariate technique (Spatial Principal Components Analysis-sPCA). Three consensus Bayesian clusters were identified for the stone marten. No consensus was achieved for the red fox, but one cluster was the most probable clustering solution. Progressive partitioning and sPCA suggested additional clusters in the stone marten but they were not consistent among methods and were geographically incoherent. The contrasting results between the two species are consistent with the literature reporting stricter ecological requirements of the stone marten in the Iberian Peninsula. The observed genetic structure in the stone marten may have been influenced by landscape features, particularly rivers, and fragmentation. We suggest that an approach based on a consensus clustering solution of multiple different algorithms may provide an objective and effective means to delineate potential boundaries of inferred subpopulations. sPCA and progressive partitioning offer further verification of possible population structure and may be useful for revealing cryptic spatial genetic patterns worth further investigation.

Landscape genetics and population structure in Valley Oak (Quercus lobata Née)

PREMISE OF THE STUDY

Although long-distance pollen movement is common in wind-pollinated trees, barriers to gene flow may occur in species that have discontinuous ranges or are confined to certain habitat types. We investigated the genetic structure of Quercus lobata Née populations throughout much of their range in California. We assessed the connectivity of populations and determined if barriers to gene flow occurred, and if so, if they corresponded to landscape features.

METHODS

We collected leaf samples from 270 trees from 12 stands of Quercus lobata and genotyped these trees using eight polymorphic microsatellite loci. Genetic structure and clustering was evaluated using genetic distance methods, Bayesian clustering approaches, and network analysis of spatial genetic structure.

KEY RESULTS

The southernmost population of Quercus lobata sampled from the Santa Monica area comprised a separate genetic cluster from the rest of the species, suggesting that Transverse Ranges such as the San Gabriel Mountains limit gene flow. Population differentiation among the other sites was small but significant. Network analysis reflected higher connectivity among populations along the Central Coast range, with few connections spanning the dry, low Central Valley.

CONCLUSIONS

While long distance pollen movement has been shown to be common in oaks, on larger spatial scales, topographic features such as mountain ranges and the large, flat Central Valley of California limit gene flow. Such landscape features explain gene flow patterns much better than geographic distance alone.

Cryopreservation of embryogenic tissues from mature holm oak trees

The development of a vitrification method for cryopreservation of embryogenic lines from mature holm oak (Quercus ilex L.) trees is reported. Globular embryogenic clusters of three embryogenic lines grown on gelled medium, and embryogenic clumps of one line collected from liquid cultures, were used as samples. The effect of both high-sucrose preculture and dehydration by incubation in the PVS2 solution for 30-90 min, on both survival and maintenance of the differentiation ability was evaluated in somatic embryo explants with and without immersion into liquid nitrogen. Growth recovery of the treated samples and ability to differentiate cotyledonary embryos largely depended on genotype. Overall, high growth recovery frequencies on gelled medium and increase of fresh weight in liquid medium were obtained in all the tested lines, also after freezing. However, the differentiation ability of the embryogenic lines was severely hampered following immersion into LN. Two of the embryogenic lines from gelled medium were able to recover the differentiation ability, one not. In the lines with reduced or no differentiation ability, variation in the microsatellite markers was observed when comparing samples taken prior to and after cryopreservation. The best results were achieved in the genotype Q8 in which 80% of explants grown on gelled medium differentiated into cotyledonary embryos following cryopreservation when they were precultured on medium with 0.3M sucrose and then incubated for 30 min in the PVS2 solution. Explants of the same genotype from liquid medium were unable to recover the differentiation ability. A 4-weeks storage period both in liquid nitrogen and in an ultra-low temperature freezer at -80°C was also evaluated with four embryogenic lines from gelled medium using the best vitrification treatment. Growth recovery frequencies of all lines from the two storage systems were very high, but their differentiation ability was completely lost.

The effect of drought stress on heterozygosity-fitness correlations in pedunculate oak (Quercus robur)

BACKGROUND AND AIMS

The interaction between forest fragmentation and predicted climate change may pose a serious threat to tree populations. In small and spatially isolated forest fragments, increased homozygosity may directly affect individual tree fitness through the expression of deleterious alleles. Climate change-induced drought stress may exacerbate these detrimental genetic consequences of forest fragmentation, as the fitness response to low levels of individual heterozygosity is generally thought to be stronger under environmental stress than under optimal conditions.

METHODS

To test this hypothesis, a greenhouse experiment was performed in which various transpiration and growth traits of 6-month-old seedlings of Quercus robur differing in multilocus heterozygosity (MLH) were recorded for 3 months under a well-watered and a drought stress treatment. Heterozygosity-fitness correlations (HFC) were examined by correlating the recorded traits of individual seedlings to their MLH and by studying their response to drought stress.

KEY RESULTS

Weak, but significant, effects of MLH on several fitness traits were obtained, which were stronger for transpiration variables than for the recorded growth traits. High atmospheric stress (measured as vapour pressure deficit) influenced the strength of the HFCs of the transpiration variables, whereas only a limited effect of the irrigation treatment on the HFCs was observed.

CONCLUSIONS

Under ongoing climate change, increased atmospheric stress in the future may strengthen the negative fitness responses of trees to low MLH. This indicates the necessity to maximize individual multilocus heterozygosity in forest tree breeding programmes.

High rates of gene flow by pollen and seed in oak populations across Europe

Gene flow is a key factor in the evolution of species, influencing effective population size, hybridisation and local adaptation. We analysed local gene flow in eight stands of white oak (mostly Quercus petraea and Q. robur, but also Q. pubescens and Q. faginea) distributed across Europe.

Adult trees within a given area in each stand were exhaustively sampled (range [239, 754], mean 423), mapped, and acorns were collected ([17,147], 51) from several mother trees ([3], [47], 23). Seedlings ([65,387], 178) were harvested and geo-referenced in six of the eight stands. Genetic information was obtained from screening distinct molecular markers spread across the genome, genotyping each tree, acorn or seedling. All samples were thus genotyped at 5–8 nuclear microsatellite loci. Fathers/parents were assigned to acorns and seedlings using likelihood methods. Mating success of male and female parents, pollen and seed dispersal curves, and also hybridisation rates were estimated in each stand and compared on a continental scale.

On average, the percentage of the wind-borne pollen from outside the stand was 60%, with large variation among stands (21–88%). Mean seed immigration into the stand was 40%, a high value for oaks that are generally considered to have limited seed dispersal. However, this estimate varied greatly among stands (20–66%). Gene flow was mostly intraspecific, with large variation, as some trees and stands showed particularly high rates of hybridisation.

Our results show that mating success was unevenly distributed among trees. The high levels of gene flow suggest that geographically remote oak stands are unlikely to be genetically isolated, questioning the static definition of gene reserves and seed stands.

Efficient mitigation of founder effects during the establishment of a leading-edge oak population

Numerous plant species are shifting their range polewards in response to ongoing climate change. Range shifts typically involve the repeated establishment and growth of leading-edge populations well ahead of the main species range. How these populations recover from founder events and associated diversity loss remains poorly understood. To help fill this gap, we exhaustively investigated a newly established population of holm oak (Quercus ilex) growing more than 30 km ahead of the nearest larger stands. Pedigree reconstructions showed that plants belong to two non-overlapping generations and that the whole population originates from only two founder trees. The four first-generation trees that have reached maturity showed disparate mating patterns despite being full-sibs. Long-distance pollen immigration was notable despite the strong isolation of the stand: 6 per cent gene flow events in acorns collected on the trees (n = 255), and as much as 27 per cent among their established offspring (n = 33). Our results show that isolated leading-edge populations of wind-pollinated forest trees can rapidly restore their genetic diversity through the interacting effects of efficient long-distance pollen flow and purging of inbred individuals during recruitment. They imply that range expansions of these species are primarily constrained by initial propagule arrival rather than by subsequent gene flow.

Molecular characterization and genetic structure of Quercus acutissima germplasm in China using microsatellites

Quercus acutissima is native to eastern Asia. It has a wide distribution in China and China is an important component in understanding the ecology and genetic structure of this species. Q. acutissima attained high economic value for hardwood product and can be managed as an energy tree species. To investigate the genetic variation of Q. acutissima provenances, 12 microsatellite primer pairs were used to analyze 672 trees sampled from 28 provenances of Q. acutissima in China. All of the tested microsatellite loci proved to be effective for the studied Q. acutissima provenances. The results revealed that allele numbers varied from 5 to 13 per locus, with an average of 8 alleles per locus. The mean observed heterozygosity and expected heterozygosity were 0.4927 and 0.7023, respectively. The relatedness of the provenances was studied using the arithmetic mean algorithm based on Nei's genetic distance and principal coordinates analysis. Interestingly, both approaches revealed two main groups: one consisted of the eastern Chinese provenances, and the other comprised of the western Chinese provenances. An analysis of molecular variance indicated that most genetic variation was contained within populations (84 %). The two microsatellite markers developed in this study may be employed for genetic characterization of other oak species. Considering the management or breeding programs of Q. acutissima provenances in China, we should treat each main group as a single gene resource.

Inter- and intra-specific variability in isoprene production and photosynthesis of Central European oak species

European deciduous oaks are closely related and are known for their strong emission of volatile isoprenoids. They are chemo-taxonomically diverse, but hybridise frequently. Four-year-old oak seedlings growing together in a model ecosystem facility under near-natural conditions were studied. The leaves were morphologically classified in the three oak species Quercus robur, Q. pubescens and Q. petraea (with four provenances each) and further investigated by a molecular-genetic approach. Q. robur was morphologically and genetically clearly different from Q. pubescens and Q. petraea, whereas Q. pubescens and Q. petraea individuals used in this study were morphologically and genetically more similar. There was a minor impact of among and within species variability on isoprene synthesis, isoprene emission and photosynthesis. Isoprene emission rates normalised to 25 °C leaf temperature ranged from 5.78 to 10.66 nmol m(-2)  s(-1) , whereas photosynthesis ranged from 12.8 to 17.6 μmol m(-2)  s(-1) . On cloudy days, among the provenances of each species, only net photosynthesis of the Q. robur provenance Hünenberg was reduced and isoprene synthase activity of the Q. pubescens provenance Promotogno increased. On sunny days, photosynthesis did not differ among the provenances. Over all provenances, gas exchange on cloudy days did not differ significantly from sunny days. In the combined data of cloudy and sunny days, no differences between the studied provenances and oak species were detected in isoprene emission and photosynthesis. Thus, isoprene emission and photosynthesis rates were remarkably stable among oak species and provenances. The results indicate that taxonomic differences in the studied oak species are not reflected in isoprene emission and photosynthesis, probably because of the high plasticity of gene expression resulting in high phenotypic flexibility.

Detecting the footprints of divergent selection in oaks with linked markers

Genome scans are increasingly used to study ecological speciation, providing a useful genome-wide perspective on divergent selection in the presence of gene flow. Here, we compare current approaches to detect footprints of divergent selection in closely related species. We analyzed 192 individuals from two interfertile European temperate oak species using 30 nuclear microsatellites from eight linkage groups. These markers present little intraspecific differentiation and can be used in combination to assign individual genotypes to species. We first show that different outlier detection tests give somewhat different results, possibly due to model constraints. Second, using linkage information for these markers, we further characterize the signature of divergent selection in the presence of gene flow. In particular, we show that recombination estimates for regions with outlier markers are lower than those for a control region, in line with a prediction from ecological speciation theory. Most importantly, we show that analyses at the haplotype level can distinguish between truly divergent (bi-directional) selection and positive selection in one of the two species, offering a new and improved method for characterizing the speciation process.

Comparative mapping in the Fagaceae and beyond with EST-SSRs

BACKGROUND

Genetic markers and linkage mapping are basic prerequisites for comparative genetic analyses, QTL detection and map-based cloning. A large number of mapping populations have been developed for oak, but few gene-based markers are available for constructing integrated genetic linkage maps and comparing gene order and QTL location across related species.

RESULTS

We developed a set of 573 expressed sequence tag-derived simple sequence repeats (EST-SSRs) and located 397 markers (EST-SSRs and genomic SSRs) on the 12 oak chromosomes (2n = 2x = 24) on the basis of Mendelian segregation patterns in 5 full-sib mapping pedigrees of two species: Quercus robur (pedunculate oak) and Quercus petraea (sessile oak). Consensus maps for the two species were constructed and aligned. They showed a high degree of macrosynteny between these two sympatric European oaks. We assessed the transferability of EST-SSRs to other Fagaceae genera and a subset of these markers was mapped in Castanea sativa, the European chestnut. Reasonably high levels of macrosynteny were observed between oak and chestnut. We also obtained diversity statistics for a subset of EST-SSRs, to support further population genetic analyses with gene-based markers. Finally, based on the orthologous relationships between the oak, Arabidopsis, grape, poplar, Medicago, and soybean genomes and the paralogous relationships between the 12 oak chromosomes, we propose an evolutionary scenario of the 12 oak chromosomes from the eudicot ancestral karyotype.

CONCLUSIONS

This study provides map locations for a large set of EST-SSRs in two oak species of recognized biological importance in natural ecosystems. This first step toward the construction of a gene-based linkage map will facilitate the assignment of future genome scaffolds to pseudo-chromosomes. This study also provides an indication of the potential utility of new gene-based markers for population genetics and comparative mapping within and beyond the Fagaceae.

Influence of environmental heterogeneity on genetic diversity and structure in an endemic southern Californian oak

Understanding how specific environmental factors shape gene flow while disentangling their importance relative to the effects of geographical isolation is a major question in evolutionary biology and a specific goal of landscape genetics. Here, we combine information from nuclear microsatellite markers and ecological niche modelling to study the association between climate and spatial genetic structure and variability in Engelmann oak (Quercus engelmannii), a wind-pollinated species with high potential for gene flow. We first test whether genetic diversity is associated with climatic niche suitability and stability since the Last Glacial Maximum (LGM). Second, we use causal modelling to analyse the potential influence of climatic factors (current and LGM niche suitability) and altitude in the observed patterns of genetic structure. We found that genetic diversity is negatively associated with local climatic stability since the LGM, which may be due to higher immigration rates in unstable patches during favourable climatic periods and/or temporally varying selection. Analyses of spatial genetic structure revealed the presence of three main genetic clusters, a pattern that is mainly driven by two highly differentiated populations located in the northern edge of the species distribution range. After controlling for geographic distance, causal modelling analyses showed that genetic relatedness decreases with the environmental divergence among sampling sites estimated as altitude and current and LGM niche suitability. Natural selection against nonlocal genotypes and/or asynchrony in reproductive phenology may explain this pattern. Overall, this study suggests that local environmental conditions can shape patterns of genetic structure and variability even in species with high potential for gene flow and relatively small distribution ranges.

The genetic structure of Quercus crispula in northeastern Japan as revealed by nuclear simple sequence repeat loci

Previous studies have reached different discussions about the genetic variation and genetic structure of Quercus crispula populations in northeastern Japan. This is a common oak species in Eastern Asia. Some studies have suggested that the populations in northeastern Japan were derived from those remaining in the southwest after the last glacial maximum (LGM), whilst other studies have found evidence that populations persisted in northeastern Japan during the LGM. Using seven highly polymorphic nuclear simple sequence repeat loci, we investigated the genetic structure of 16 Q. crispula populations along a latitudinal gradient in northeastern Japan (northern Honshu and Hokkaido), spanning about half of the species' biogeographic range in the country. Although the level of population differentiation was low (F (ST) = 0.021; [Formula: see text] = 0.090), two geographically differentiated clusters were detected by STRUCTURE analysis. The first cluster included most of the populations in Hokkaido, and may indicate continued survival throughout past glacial periods. We found a significant decrease in allelic richness with latitude, so the second cluster may represent an expansion of the lineage from Honshu during the post-glacial period. These results should enhance our understanding of historical north-south migrations of this species in northeastern Japan.

Genetic differences within natural and planted stands of Quercus petraea

Five sessile oak [Quercus petraea (Matt.) Liebl.] stands from the Czech Republic were studied to learn about the impact of different types of forest management regimes on the genetic differences among tree populations and on population structures. One population had not been markedly affected by human activity, two populations represented unplanted stands that were extensively managed for a long period of time using the coppice system, and two populations were planted stands. Approximately 100 trees from each stand were mapped and subsequently genotyped using 10 nuclear microsatellite loci. We determined the spatial genetic structure of each population and the genetic differentiation among the populations. We found that: (i) the populations were genetically differentiated, but the differences between the unplanted and planted stands were not markedly significant; (ii) the genetic differentiation among the populations depended on the geographical distribution of the populations; (iii) within unplanted stands, a strong spatial genetic structure was seen; and (iv) within planted stands, no spatial genetic structure was observed. Our findings implies that the analysis of spatial genetic structure of the sessile oak forest stand can help reveal and determine its origin.

Two highly validated multiplexes (12-plex and 8-plex) for species delimitation and parentage analysis in oaks (Quercus spp.)

Multiplex PCR is a fast and cost-effective technique allowing increased genotyping throughput of microsatellites. We developed two multiplexes for Quercus petraea and Q. robur, a 12-plex of EST-SSRs (eSSRs) and an 8-plex of genomic SSRs (gSSRs). We studied the origin of allele calling errors at the human reader and software levels. We showed that the robustness of allele identification can be improved by binning on raw peak sizes prior to genetic data analysis. We checked through simulation the power of these markers for species delimitation and hybrid detection. The resolution achieved with all 20 markers was greatly improved compared to that of previous studies based on a subset of the markers. Preliminary PCR tests suggest that these multiplexes might be useful to study other oak species as well. The strategy used for multiplex microsatellite development (from PCR conditions to the definition of allele calling rules) should be broadly applicable.

Pollen-mediated gene flow in isolated and continuous stands of bur oak, Quercus macrocarpa (Fagaceae)

PREMISE OF THE STUDY

Pollination patterns determine the reproductive neighborhood size of plants, the connectivity of populations, and the impacts of habitat fragmentation. We characterized pollination in three populations of Quercus macrocarpa occurring in a highly altered landscape in northeastern Illinois to determine whether isolated remnant stands were reproductively isolated. •

METHODS

We used microsatellites to genotype all adults and 787 acorns from two isolated savanna remnants and a stand in an old-growth forest. One isolated remnant occurred in a highly urbanized/industrialized landscape, and one occurred in an agricultural landscape. Parentage assignment was used to assess pollen-mediated gene flow. •

KEY RESULTS

Pollen donors from outside the study sites accounted for between 46% and 53% of paternities and did not differ significantly among sites, indicating that similar high levels of gene flow occurred at all three sites. Within stands, the mean pollination distance ranged from 42 to 70 meters, and when accounting for outside pollinations, mean pollination distances were well over 100 meters. Genetic diversity of incoming pollen was extremely high in all three stands. The number of effective pollen donors, N(ep), calculated from paternity assignment was higher than that estimated by an indirect correlated paternity approach. •

CONCLUSIONS

Our findings indicate that extremely isolated stands of oaks are unlikely to be genetically and reproductively isolated, and remnant stands may contribute to maintaining genetic connectivity in highly modified landscapes.

Exploring species limits in two closely related Chinese oaks

BACKGROUND

The species status of two closely related Chinese oaks, Quercus liaotungensis and Q. mongolica, has been called into question. The objective of this study was to investigate the species status and to estimate the degree of introgression between the two taxa using different approaches.

METHODOLOGY/PRINCIPAL FINDINGS

Using SSR (simple sequence repeat) and AFLP (amplified fragment length polymorphism) markers, we found that interspecific genetic differentiation is significant and higher than the differentiation among populations within taxa. Bayesian clusters, principal coordinate analysis and population genetic distance trees all classified the oaks into two main groups consistent with the morphological differentiation of the two taxa rather than with geographic locations using both types of markers. Nevertheless, a few individuals in Northeast China and many individuals in North China have hybrid ancestry according to Bayesian assignment. One SSR locus and five AFLPs are significant outliers against neutral expectations in the interspecific FST simulation analysis, suggesting a role for divergent selection in differentiating species.

MAIN CONCLUSIONS/SIGNIFICANCE

All results based on SSRS and AFLPS reached the same conclusion: Q. liaotungensis and Q. mongolica maintain distinct gene pools in most areas of sympatry. They should therefore be considered as discrete taxonomic units. Yet, the degree of introgression varies between the two species in different contact zones, which might be caused by different population history or by local environmental factors.

Natural hybridisation between kermes (Quercus coccifera L.) and holm oaks (Q. ilex L.) revealed by microsatellite markers

Hybridisation between species of the genus Quercus is a common phenomenon as a result of weak reproductive isolation mechanisms between phylogenetically close species that frequently co-occur in mixed stands. In this study, we use microsatellite markers to analyse introgression between kermes (Quercus coccifera L.) and holm (Q. ilex L.) oak, two closely related taxa that frequently dominate the landscape in extensive areas in the Mediterranean region. All tested microsatellites amplified and were polymorphic in both kermes and holm oaks. Bayesian admixture analyses showed a good correspondence between each species and one of the two inferred genetic clusters. Five sampled individuals were a priori tentatively identified as hybrids on the basis of intermediate morphological characteristics, and it was confirmed that they also presented mixed genotypes. However, we also detected different levels of genetic introgression among morphologically pure individuals, suggesting that successful backcrossing and/or reduced phenotypic expression of genetic variance in certain individuals may have resulted in strong convergence towards a single species phenotype.

Contrasting relations between diversity of candidate genes and variation of bud burst in natural and segregating populations of European oaks

Nucleotide diversity was assessed within nine candidate genes (in total 4.6 kb) for the time of bud burst in nine sessile oak (Quercus petraea) populations distributed in central and northern Europe. The sampled populations were selected on the basis of their contrasting time of bud burst observed in common garden experiments (provenance tests). The candidate genes were selected according to their expression profiles during the transition from quiescent to developing buds and/or their functional role in model plants. The overall nucleotide diversity was large (π(tot)=6.15 × 10(-3); π(silent)=11.2 × 10(-3)), but population differentiation was not larger than for microsatellites. No outlier single-nucleotide polymorphism (SNP), departing from neutral expectation, was found among the total of 125 SNPs. These results contrasted markedly with the significant associations that were observed between the candidate genes and bud burst in segregating populations. Quantitative trait loci (QTLs) for bud burst were identified for 13 year*site seasonal observations in a cloned mapping pedigree. Nineteen QTLs were detected, and QTLs located on linkage groups 2, 5 and 9 contributed repeatedly to more than 12% of the phenotypic variation of the trait. Eight genes were polymorphic in the two parents of the pedigree and could be mapped on the existing genetic map. Five of them located within the confidence intervals of QTLs for bud burst. Interestingly, four of them located within the three QTLs exhibiting the largest contributions to bud burst.

Patterns of contemporary hybridization inferred from paternity analysis in a four-oak-species forest

Background

Few studies address the issue of hybridization in a more than two-species context. The species-rich Quercus complex is one of the systems which can offer such an opportunity. To investigate the contemporary pattern of hybridization we sampled and genotyped 320 offspring from a natural mixed forest comprising four species of the European white oak complex: Quercus robur, Q. petraea, Q. pubescens, and Q. frainetto.

Results

A total of 165 offspring were assigned unambiguously to one of the pollen donors within the study plot. The minimum amount of effective pollen originating from outside the plot varied markedly among the seed parents, ranging from 0.18 to 0.87. The majority of the successful matings (64.1%) occurred between conspecific individuals indicating the existence of reproductive barriers between oak species. However, the isolation was not complete since we found strong evidence for both first-generation (8.4%) and later-generation hybrids (27.5%). Only two out of eight seed parents, belonging to Q. petraea and Q. robur, showed a high propensity to hybridize with Q. pubescens and Q. petraea, respectively. Significant structure of the effective pollen pools (Φ_pt = 0.069, P = 0.01) was detected in our sample. However, no support was found for the isolation by distance hypothesis. The proportion of hybrids was much higher (79%) in the seed generation when compared to the adult tree generation.

Conclusion

First-generation hybrids were observed only between three out of six possible species combinations. Hybrids between one pair of species preferred to mate with one of their parental species. The observation of first and later-generation hybrids in higher frequency in acorns than in adults might be explained by selection against hybrid genotypes, the history of this uneven-aged forest or past introgression between species.

Relevance of genetics for conservation policies: the case of Minorcan cork oaks

BACKGROUND AND AIMS

Marginal populations of widely distributed species can be of high conservation interest when they hold a significant or unique portion of the genetic diversity of the species. However, such genetic information is frequently lacking. Here the relevance of genetic surveys to develop efficient conservation strategies for such populations is illustrated using cork oak (Quercus suber) from Minorca (Balearic Islands, Spain) as a case study. Cork oak is highly endangered on the island, where no more than 67 individuals live in small, isolated stands in siliceous sites. As a consequence, it was recently granted protected status.

METHODS

Two Bayesian clustering approaches were used to analyse the genetic structure of the Minorcan population, on the basis of nuclear microsatellite data. The different groups within the island were also compared with additional island and continental populations surrounding Minorca.

KEY RESULTS

Very high genetic diversity was found, with values comparable with those observed in continental parts of the species' range. Furthermore, the Minorcan oak stands were highly differentiated from one another and were genetically related to different continental populations of France and Spain.

CONCLUSIONS

The high levels of genetic diversity and inter-stands differentiation make Minorcan cork oak eligible for specific conservation efforts. The relationship of Minorcan stands to different continental populations in France and Spain probably reflects multiple colonization events. However, discrepancy between chloroplast DNA- and nuclear DNA-based groups does not support a simple scenario of recent introduction. Gene exchanges between neighbouring cork oak stands and with holm oak have created specific and exceptional genetic combinations. They also constitute a wide range of potential genetic resources for research on adaptation to new environmental conditions. Conservation guidelines that take into account these findings are provided.

Natural hybridisation between Quercus petraea (Matt.) Liebl. and Quercus pubescens Willd. within an Italian stand as revealed by microsatellite fingerprinting

Interspecific gene flow is frequently reported in the genus Quercus. However, interfertile oak species often seem to remain distinct, even within areas of sympatry. This study employed molecular markers to verify, at a fine scale, the presence of interspecific gene flow in a natural population of Quercus petraea and Quercus pubescens. Within a delimited area of 6 ha, all adult trees belonging to the studied oak complex and seeds from a subsample of such trees were collected and analysed using molecular microsatellite markers. A low interspecific genetic differentiation and a high level of interspecific genetic admixture suggested past hybridisation. Paternity inference of seeds allowed the estimation of pollination frequencies from the three groups of pollen donors (Q. petraea, Q. pubescens, intermediate). We also assayed pollen viability and germinability of each species group. We observed natural hybridisation between Q. petraea and Q. pubescens, with a predominant component in the direction Q. petraea --> Q. pubescens: Q. pubescens displayed a higher level of heterospecific pollination by Q. petraea (25.8%) and intermediate morphotypes (14.7%), compared to Q. petraea acting as pollen receptor (with less than 5% heterospecific pollinations). Intermediate 'mother trees' were pollinated in similar proportions by Q. petraea (23.1%), Q. pubescens (37.8%) and intermediate morphotypes (39.1%). The asymmetrical introgression observed for the studied generation may be caused, among other factors, by the relative abundance of trees from each species group in the studied area.

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.

Detection of hybrids in nature: application to oaks (Quercus suber and Q. ilex)

Powerful and accurate detection of first-generation (F1) hybrids and backcrosses in nature is needed to achieve a better understanding of the function and dynamics of introgression. To document the frequency of ongoing interspecific gene exchange between two Mediterranean evergreen oaks, the cork oak (Quercus suber) and the holm oak (Q. ilex), we analyzed 1487 individuals originating from across the range of the two species using eight microsatellite loci and two Bayesian clustering approaches (implemented in the programs STRUCTURE and NEWHYBRIDS). Simulated data were used to assess the differences between the two clustering methods and to back up the choice of the threshold value for the posterior probability to discriminate admixed from pure individuals. We found that the use of STRUCTURE resulted in the highest power to detect hybrids, whereas NEWHYBRIDS provided the highest accuracy. Irrespective of the approach, the two species were clearly distinguished as independent genetic entities without any prior information. In contrast with previous reports, we found no evidence for unidirectional introgression. The overall hybridization rate was very low (<2% of introgressed individuals). Only two individuals were identified as F1 hybrids and five as early backcrosses. This work shows that the combined application of the two complementary Bayesian approaches and their systematic validation with simulations, fit for the case at hand, helps gain resolution in the identification of admixed individuals.

Utility of multilocus genotypes for taxon assignment in stands of closely related European white oaks from Switzerland

BACKGROUND AND AIMS

European white oaks (Quercus petraea, Q. pubescens, Q. robur) have long puzzled plant biologists owing to disputed species differentiation. Extensive hybridization or shared ancestry have been proposed as alternative hypotheses to explain why genetic differentiation between these oak species is low. Species delimitation is usually weak and often shows gradual transitions in leaf morphology. Hence, individual identification may be difficult, but remains a critical step for both scientific work and practical management.

METHODS

Multilocus genotype data (five nuclear microsatellites) were used from ten Swiss oak stands for taxon identification without a priori grouping of individuals or populations, using model-based Bayesian assignment tests.

KEY RESULTS

Three groups best structured the data, indicating that the taxonomical signal was stronger than the spatial signal. Most individuals showed high posterior probabilities for either of three genetic groups that were best circumscribed as taxonomical units. The assignment of a subset of trees, whose taxonomic status had been previously characterized in detail, supported this classification scheme.

CONCLUSIONS

Molecular-genetic assignment tests are useful in the identification of species status in critical taxon complexes such as the European white oaks. Such an approach is of practical importance for forest management, e.g. for stand certification or in seed trade to trace the origin of forest products.

Conserving the evolutionary potential of California valley oak (Quercus lobataNée): a multivariate genetic approach to conservation planning

California valley oak (Quercus lobata Née) is a seriously threatened endemic oak species in California and a keystone species for foothill oak ecosystems. Urban and agricultural development affects a significant fraction of the species' range and predicted climate change is likely to dislocate many current populations. Here, we explore spatial patterns of multivariate genotypes and genetic diversity throughout the range of valley oak to determine whether ongoing and future patterns of habitat loss could threaten the evolutionary potential of the species by eradicating populations of distinctive genetic composition. This manuscript will address three specific questions: (i) What is the spatial genetic structure of the chloroplast and nuclear genetic markers? (ii) What are the geographical trends in the distribution of chloroplast and nuclear genotypes? (iii) Is there any part of the species' range where allelic diversity in either the chloroplast or nuclear genomes is particularly high? We analysed six chloroplast and seven nuclear microsatellite genetic markers of individuals widespread across the valley oak range. We then used a multivariate approach correlating genetic markers and geographical variables through a canonical trend surface analysis, followed by GIS mapping of the significant axes. We visualized population allelic richness spatially with GIS tools to identify regions of high diversity. Our findings, based on the distribution of multivariate genotypes and allelic richness, identify areas with distinctive histories and genetic composition that should be given priority in reserve network design, especially because these areas also overlap with landscape change and little degree of protection. Thus, without a careful preservation plan, valuable evolutionary information will be lost for valley oak.

Evidence for hybridization and introgression within a species-rich oak (Quercus spp.) community

Background

Analysis of interspecific gene flow is crucial for the understanding of speciation processes and maintenance of species integrity. Oaks (genus Quercus, Fagaceae) are among the model species for the study of hybridization. Natural co-occurrence of four closely related oak species is a very rare case in the temperate forests of Europe. We used both morphological characters and genetic markers to characterize hybridization in a natural community situated in west-central Romania and which consists of Quercus robur, Q. petraea, Q. pubescens, and Q. frainetto, respectively.

Results

On the basis of pubescence and leaf morphological characters ~94% of the sampled individuals were assigned to pure species. Only 16 (~6%) individual trees exhibited intermediate morphologies or a combination of characters of different species. Four chloroplast DNA haplotypes were identified in the study area. The distribution of haplotypes within the white oak complex showed substantial differences among species. However, the most common haplotypes were present in all four species. Furthermore, based on a set of 7 isozyme and 6 microsatellite markers and using a Bayesian admixture analysis without any a priori information on morphology we found that four genetic clusters best fit the data. There was a very good correspondence of each species with one of the inferred genetic clusters. The estimated introgression level varied markedly between pairs of species ranging from 1.7% between Q. robur and Q. frainetto to 16.2% between Q. pubescens and Q. frainetto. Only nine individuals (3.4%) appeared to be first-generation hybrids.

Conclusion

Our data indicate that natural hybridization has occurred at relatively low rates. The different levels of gene flow among species might be explained by differences in flowering time and spatial position within the stand. In addition, a partial congruence between phenotypically and genetically intermediate individuals was found, suggesting that intermediate appearance does not necessarily mean hybridization. However, it appears that natural hybridization did not seriously affect the species identity in this area of sympatry.