Reproductive and genetic consequences of extreme isolation in Salix herbacea L. at the rear edge of its distribution.

BACKGROUND AND AIMS: At the rear edge of the distribution of species, extreme isolation and small population size influence the genetic diversity and differentiation of plant populations. This may be particularly true for Arctic-alpine species in mid-latitude mountains, but exactly how peripherality has shaped their genetic and reproductive characteristics is poorly investigated. The present study, focused on Salix herbacea, aims at providing new insights into the causes behind ongoing demographic dynamics and their consequences for peripheral populations of Arctic-alpine species. METHODS: We performed a whole-population, highly detailed sampling of the only two S. herbacea populations in the northern Apennines, comparing their clonal and genetic diversity, sex ratio and spatial genetic structure with a reference population from the Alps. After inspecting ~1800 grid intersections in the three populations, 563 ramets were genotyped at 11 nuclear microsatellite markers (nSSRs). Past demography and mating patterns of Apennine populations were investigated to elucidate the possible causes of altered reproductive dynamics. KEY RESULTS: Apennine populations, which experienced a Holocene bottleneck and are highly differentiated (FST = 0.15), had lower clonal and genetic diversity compared with the alpine population (RMLG = 1 and HE = 0.71), with the smaller population exhibiting the lowest diversity (RMLG = 0.03 and HE = 0.24). An unbalanced sex ratio was found in the larger (63 F:37 M) and the smaller (99 F:1 M) Apennine population. Both were characterized by the presence of extremely large clones (up to 2500 m2), which, however, did not play a dominant role in local reproductive dynamics. CONCLUSIONS: Under conditions of extreme isolation and progressive size reduction, S. herbacea has experienced an alteration of genetic characteristics produced by the prevalence of clonal growth over sexual reproduction. However, our results showed that the larger Apennine population has maintained levels of sexual reproduction enough to counteract a dramatic loss of genetic and clonal diversity.

Linking dendroecology and association genetics in natural populations: Stress responses archived in tree rings associate with SNP genotypes in silver fir (Abies alba Mill.).

Genetic association studies in forest trees would greatly benefit from information on the response of trees to environmental stressors over time, which can be provided by dendroecological analysis. Here, we jointly analysed dendroecological and genetic data of surviving silver fir trees to explore the genetic basis of their response to the iconic stress episode of the 1970s and 1980s that led to large-scale forest dieback in Central Europe and has been attributed to air pollution. Specifically, we derived dendrophenotypic measures from 190 trees in the Bavarian Forest that characterize the resistance, resilience and recovery during this growth depression, and in the drought year in 1976. By focusing on relative growth changes of trees and by standardizing the dendrophenotypes within stands, we accounted for variation introduced by micro- and macroscale environmental differences. We associated the dendrophenotypes with single nucleotide polymorphisms (SNPs) in candidate genes using general linear models (GLMs) and the machine learning algorithm random forest with subsequent feature selection. Most trees at our study sites experienced a severe growth decline from 1974 until the mid-1980s with minimum values during the drought year. Fifteen genes were associated with the dendrophenotypes, including genes linked to photosynthesis and drought stress. With our study, we show that dendrophenotypes can be a powerful resource for genetic association studies that permit to account for micro- and macroenvironmental variation when data are derived from natural populations. We call for a wider collaboration of dendroecologists and forest geneticists to integrate individual tree-level dendrophenotypes in genetic association studies.

High-density SNP assay development for genetic analysis in maritime pine (Pinus pinaster).

Maritime pine provides essential ecosystem services in the south-western Mediterranean basin, where it covers around 4 million ha. Its scattered distribution over a range of environmental conditions makes it an ideal forest tree species for studies of local adaptation and evolutionary responses to climatic change. Highly multiplexed single nucleotide polymorphism (SNP) genotyping arrays are increasingly used to study genetic variation in living organisms and for practical applications in plant and animal breeding and genetic resource conservation. We developed a 9k Illumina Infinium SNP array and genotyped maritime pine trees from (i) a three-generation inbred (F2) pedigree, (ii) the French breeding population and (iii) natural populations from Portugal and the French Atlantic coast. A large proportion of the exploitable SNPs (2052/8410, i.e. 24.4%) segregated in the mapping population and could be mapped, providing the densest ever gene-based linkage map for this species. Based on 5016 SNPs, natural and breeding populations from the French gene pool exhibited similar level of genetic diversity. Population genetics and structure analyses based on 3981 SNP markers common to the Portuguese and French gene pools revealed high levels of differentiation, leading to the identification of a set of highly differentiated SNPs that could be used for seed provenance certification. Finally, we discuss how the validated SNPs could facilitate the identification of ecologically and economically relevant genes in this species, improving our understanding of the demography and selective forces shaping its natural genetic diversity, and providing support for new breeding strategies.

Genomic exploration and molecular marker development in a large and complex conifer genome using RADseq and mRNAseq.

We combined restriction site associated DNA sequencing (RADseq) using a hypomethylation-sensitive enzyme and messenger RNA sequencing (mRNAseq) to develop molecular markers for the 16 gigabase genome of Cedrus atlantica, a conifer tree species. With each method, Illumina(®) reads from one individual were used to generate de novo assemblies. SNPs from the RADseq data set were detected in a panel of one single individual and three pools of three individuals each. We developed a flexible script to estimate the ascertainment bias in SNP detection considering the pooling and sampling effects on the probability of not detecting an existing polymorphism. Gene Ontology (GO) and transposable element (TE) search analyses were applied to both data sets. The RADseq and the mRNAseq assemblies represented 0.1% and 0.6% of the genome, respectively. Genome complexity reduction resulted in 17% of the RADseq contigs potentially coding for proteins. This rate was doubled in the mRNAseq data set, suggesting that RADseq also explores noncoding low-repeat regions. The two methods gave very similar GO-slim profiles. As expected, the two assemblies were poor in TE-like sequences (<4% of contigs length). We identified 17,348 single nucleotide polymorphisms (SNPs) in the RADseq data set and 5,714 simple sequence repeats (SSRs) in the transcriptome. A subset of 282 SNPs was validated using the Fluidigm genotyping technology, giving a conversion rate of 50.4%, falling within the expected range for conifers. Increasing sample size had the greatest effect for ascertainment bias reduction. These results validated the utility of the RADseq approach for highly complex genomes such as conifers.

Estimating exotic gene flow into native pine stands: zygotic vs. gametic components.

Monitoring contemporary gene flow from widespread exotic plantations is becoming an important problem in forest conservation genetics. In plants, where both seed and pollen disperse, three components of exotic gene flow with potentially unequal consequences should be, but have not been, explicitly distinguished: zygotic, male gametic and female gametic. Building on a previous model for estimating contemporary rates of zygotic and male gametic gene flow among plant populations, we present here an approach that additionally estimates the third (female gametic) gene flow component, based on a combination of uni- and biparentally inherited markers. Using this method and a combined set of chloroplast and nuclear microsatellites, we estimate gene flow rates from exotic plantations into two Iberian relict stands of maritime pine (Pinus pinaster) and Scots pine (Pinus sylvestris). Results show neither zygotic nor female gametic gene flow but moderate (6-8%) male gametic introgression for both species, implying significant dispersal of pollen, but not of seeds, from exotic plantations into native stands shortly after introduced trees reached reproductive maturity. Numerical simulation results suggest that the model yields reasonably accurate estimates for our empirical data sets, especially for larger samples. We discuss conservation management implications of observed levels of exposure to nonlocal genes and identify research needs to determine potentially associated hazards. Our approach should be useful for plant ecologists and ecosystem managers interested in the vectors of contemporary genetic connectivity among discrete plant populations.

First insights into the transcriptome and development of new genomic tools of a widespread circum-Mediterranean tree species, Pinus halepensis Mill.

Aleppo pine (Pinus halepensis Mill.) is a relevant conifer species for studying adaptive responses to drought and fire regimes in the Mediterranean region. In this study, we performed Illumina next-generation sequencing of two phenotypically divergent Aleppo pine accessions with the aims of (i) characterizing the transcriptome through Illumina RNA-Seq on trees phenotypically divergent for adaptive traits linked to fire adaptation and drought, (ii) performing a functional annotation of the assembled transcriptome, (iii) identifying genes with accelerated evolutionary rates, (iv) studying the expression levels of the annotated genes and (v) developing gene-based markers for population genomic and association genetic studies. The assembled transcriptome consisted of 48,629 contigs and covered about 54.6 Mbp. The comparison of Aleppo pine transcripts to Picea sitchensis protein-coding sequences resulted in the detection of 34,014 SNPs across species, with a Ka /Ks average value of 0.216, suggesting that the majority of the assembled genes are under negative selection. Several genes were differentially expressed across the two pine accessions with contrasted phenotypes, including a glutathione-s-transferase, a cellulose synthase and a cobra-like protein. A large number of new markers (3334 amplifiable SSRs and 28,236 SNPs) have been identified which should facilitate future population genomics and association genetics in this species. A 384-SNP Oligo Pool Assay for genotyping with the Illumina VeraCode technology has been designed which showed an high overall SNP conversion rate (76.6%). Our results showed that Illumina next-generation sequencing is a valuable technology to obtain an extensive overview on whole transcriptomes of nonmodel species with large genomes.

Nuclear microsatellites for Pinus pinea (Pinaceae), a genetically depauperate tree, and their transferability to P. halepensis.

PREMISE OF THE STUDY: Pinus pinea is one of the few widespread plant species that are also genetically depauperate. It is also an important commercial species with high market value seeds. A deeper knowledge of the existing population genetic variation was needed. METHODS AND RESULTS: Twelve nuclear microsatellites were isolated from genomic and cDNA sequences and screened for variability in 729 individuals from 33 natural populations. Low level of genetic variability was confirmed with average expected heterozygosity of 0.11. Hardy-Weinberg equilibrium expectations were not met in only ∼10% of the possible locus/population combinations. All loci were in linkage equilibrium, and the frequency of null alleles was very low (≤1% in 332 out of 396 locus/population combinations). Nine out of the 12 microsatellites were successfully transferred to P. halepensis. CONCLUSIONS: Despite low polymorphism, these new markers will be useful to resolve population structure and hold potential for seed origin identification and traceability.

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

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

Comparison of pollen gene flow among four European beech (Fagus sylvatica L.) populations characterized by different management regimes.

The study of the dispersal capability of a species can provide essential information for the management and conservation of its genetic variability. Comparison of gene flow rates among populations characterized by different management and evolutionary histories allows one to decipher the role of factors such as isolation and tree density on gene movements. We used two paternity analysis approaches and different strategies to handle the possible presence of genotyping errors to obtain robust estimates of pollen flow in four European beech (Fagus sylvatica L.) populations from Austria and France. In each country one of the two plots is located in an unmanaged forest; the other plots are managed with a shelterwood system and inside a colonization area (in Austria and France, respectively). The two paternity analysis approaches provided almost identical estimates of gene flow. In general, we found high pollen immigration (~75% of pollen from outside), with the exception of the plot from a highly isolated forest remnant (~50%). In the two unmanaged plots, the average within-population pollen dispersal distances (from 80 to 184 m) were higher than previously estimated for beech. From the comparison between the Austrian managed and unmanaged plots, that are only 500 m apart, we found no evidence that either gene flow or reproductive success distributions were significantly altered by forest management. The investigated phenotypic traits (crown area, height, diameter and flowering phenology) were not significantly related with male reproductive success. Shelterwood seems to have an effect on the distribution of within-population pollen dispersal distances. In the managed plot, pollen dispersal distances were shorter, possibly because adult tree density is three-fold (163 versus 57 trees per hectare) with respect to the unmanaged one.

Genetic evidence for a Janzen-Connell recruitment pattern in reproductive offspring of Pinus halepensis trees.

Effective seed dispersal, combining both dispersal and postdispersal (establishment) processes, determines population dynamics and colonization ability in plants. According to the Janzen-Connell (JC) model, high mortality near the mother plant shifts the offspring establishment distribution farther away from the mother plant relative to the seed dispersal distribution. Yet, extending this prediction to the distribution of mature (reproductive) offspring remains a challenge for long-living plants. To address this challenge, we selected an isolated natural Aleppo pine (Pinus halepensis) population in Mt. Pithulim (Israel), which expanded from five ancestor trees in the beginning of the 20th century into ∼2000 trees today. Using nine microsatellite markers, we assigned parents to trees established during the early stages of population expansion. To elucidate the effect of the distance from the mother plant on postdispersal survival, we compared the effective seed dispersal kernel, based on the distribution of mother-offspring distances, with the seed dispersal kernel, based on simulations of a mechanistic wind dispersal model. We found that the mode of the effective dispersal kernel is shifted farther away than the mode of the seed dispersal kernel, reflecting increased survival with increasing distance from the mother plant. The parentage analysis demonstrated a highly skewed reproductive success and a strong directionality in effective dispersal corresponding to the wind regime. We thus provide compelling evidence that JC effects act also on offspring that become reproductive and persist as adults for many decades, a key requirement in assessing the role of postdispersal processes in shaping population and community dynamics.

Chloroplast DNA phylogeography of the shrub Cistus ladanifer L. (Cistaceae) in the highly diverse Western Mediterranean region.

This study investigated the phylogeographic structure of Cistus ladanifer, in order to locate its Quaternary refugia, reconstruct its recolonisation patterns and assess the role of geographical features (mountain ranges, rivers and the Strait of Gibraltar) as barriers to its seed flow and expansion through the Western Mediterranean. Thirty-eight populations were screened for length variation of polymorphic chloroplast simple sequence repeats (cpSSRs). Statistical analyses included estimation of haplotypic diversity, hierarchical analysis of molecular variation (amova) and fixation indices. Mantel tests, SAMOVA and BARRIER analyses were applied to evaluate the geographical partitioning of genetic diversity across the entire species range. Pollen data from bibliography were used to complement molecular inferences. Chlorotype diversity within populations was similar throughout the natural range of C. ladanifer (mean haplotypic diversity=0.32). High differentiation among populations was estimated (G(ST)=0.60). Our data suggest that the barriers of the Strait of Gibraltar and the Betic ranges may have favoured the divergence during glacial periods of four different lineages of populations inferred with SAMOVA. The main northward colonisation of in the Iberian Peninsula occurred from refugia in southwest Iberia. This process may have been influenced by human activities (forest clearance, livestock grazing and even commerce) in the Iberian Peninsula. In contrast, populations in the Betic area have conserved a specific haplotype.

The Strait of Gibraltar as a major biogeographic barrier in Mediterranean conifers: a comparative phylogeographic survey.

The Strait of Gibraltar (SG) is reputed for being both a bridge and a geographic barrier to biological exchanges between Europe and Africa. Major genetic breaks associated with this strait have been identified in various taxa, but it is unknown whether these disjunctions have been produced simultaneously or by independent biogeographic processes. Here, the genetic structure of five conifers distributed on both sides of the SG was investigated using mitochondrial (nad1 b/c, nad5-1, nad5-4 and nad7-1) and chloroplast (Pt1254, Pt15169, Pt30204, Pt36480, Pt71936 and Pt87268) DNA markers. The distribution of genetic variation was partially congruent between types of markers within the same species. Across taxa, there was a significant overlapping between the SG and the genetic breaks detected, especially for the four Tertiary species surveyed (Abies pinsapo complex, Pinus nigra, Pinus pinaster and Taxus baccata). For most of these taxa, the divergence of populations across the SG could date back to long before the Pleistocene glaciations. However, their strongly different cpDNA G(ST) and R(ST) values point out that they have had dissimilar population histories, which might include contrasting amounts of pollen-driven gene flow since their initial establishment in the region. The fifth species, Pinus halepensis, was genetically depauperated and homogenous on both sides of the SG. A further analysis of nuclear DNA sequences with coalescent-based isolation with migration models suggests a Pleistocene divergence of P. halepensis populations across the SG, which is in sharp contrast with the pre-Pleistocene divergence dates obtained for P. pinaster. Altogether, these results indicate that the genetic breaks observed across this putative biogeographical barrier have been produced by independent evolutionary processes related to the biological history of each individual species instead of a common vicariant phenomenon.

Spatial genetic structure of Taxus baccata L. in the western Mediterranean Basin: past and present limits to gene movement over a broad geographic scale.

English yew (Taxus baccata L., Taxaceae), a Tertiary relict, provides a seminal example of a widespread albeit locally endangered (often close to extinction) tree species. In order to gain detailed insights into the evolutionary dynamics of the species on a broad geographical scale, over 1000 trees from 91 populations of English yew in the western Mediterranean were analyzed using seven nuclear microsatellite markers. Our results revealed contrasting patterns of genetic structure at different spatial scales: genetic variation was highly structured at the local scale, while only a low proportion of the observed variation was attributed to regional differences. We also found a geographic gradient of decreasing diversity and increasing population divergence from northwest (central Europe and northern Iberian Peninsula) to southeast (Mediterranean Iberia and North Africa). The patterns revealed in this study probably reflect the combined effects of Quaternary climatic changes and recent impact of human activities, and potentially also more ancient events dating back to the Tertiary. Both climatic and anthropogenic factors seem to have conducted to a long history of population isolation, which may have contributed significantly to enhance population divergence through restricted gene flow and genetic drift.

Current range characteristics of Swiss stone pine (Pinus cembra L.) along the Carpathians revealed by chloroplast SSR markers.

We investigated the diversity pattern of nine Swiss stone pine (Pinus cembra L.) populations along the Carpathian range including the High Tatras, by using six chloroplast DNA microsatellites (cpSSR). Our aim was to detect genetically distinct regions by clustering of populations, and to tackle possible historical colonization routes. Our analysis referred to an investigated geographical range with the two most distant populations situated at about 500 air km. We found that the most diverse populations are situated at the two edges of the investigated part, in the Retezat Mts. (South Carpathians) and the High Tatras, and diversity decreases towards the populations of the Eastern Carpathians. Hierarchical clustering and NMDS revealed that the populations of the South Carpathians with the Tatras form a distinct cluster, significantly separated from those of the Eastern Carpathians. Moreover, based on the most variable chloroplast microsatellites, the four populations of the two range edges are not significantly different. Our results, supported also by palynological and late glacial macrofossil evidences, indicate refugial territories within the Retezat Mts. that conserved rich haplotype composition. From this refugial territory Pinus cembra might have colonized the Eastern Carpathians, and this was accompanied by a gradual decrease in population diversity. Populations of the High Tatras might have had the same role in the colonizing events of the Carpathians, as positive correlation was detected among populations lying from each other at a distance of 280 km, the maximum distance between neighbouring populations.

Molecular evidence for the natural production of homozygous Cupressus sempervirens L. lines by Cupressus dupreziana seed trees.

Paternal apomixis was recently reported in the endangered Mediterranean cypress, Cupressus dupreziana. This species acts as a surrogate mother for the development of all-paternal embryos from pollen grains. C. dupreziana production of Cupressus sempervirens haploid or diploid seedlings from C. sempervirens pollen was also demonstrated. The haploid progeny was derived from the embryogenic development of haploid gametes, but the origin of the diploid progeny remained unknown. To determine the ontogenic origin of the diploid C. sempervirens progeny, we analyzed the heterozygozity of 63 diploid all-paternal C. sempervirens seedlings using highly variable co-dominant nuclear microsatellite markers. The bi-parental inheritance of the markers was checked in C. sempervirens controlled crosses. A high level of polymorphism was observed among the diploid all-paternal trees. All but three individuals exhibited single-band profiles as expected for homozygotes, which may originate from natural diploidization of a C. sempervirens haploid embryo or from the fusion of two male gametes produced by the same C. sempervirens microgametophyte. The three heterozygous seedlings must be derived from the fusion of male gametes produced by two different C. sempervirens microgametophytes. These findings offer a unique opportunity in conifers to produce homozygous lines, highly valuable for genetic analyses or breeding.

Spatial genetic structure in continuous and fragmented populations of Pinus pinaster Aiton.

Habitat fragmentation, i.e., the reduction of populations into small isolated remnants, is expected to increase spatial genetic structure (SGS) in plant populations through nonrandom mating, lower population densities and potential aggregation of reproductive individuals. We investigated the effects of population size reduction and genetic isolation on SGS in maritime pine (Pinus pinaster Aiton) using a combined experimental and simulation approach. Maritime pine is a wind-pollinated conifer which has a scattered distribution in the Iberian Peninsula as a result of forest fires and habitat fragmentation. Five highly polymorphic nuclear microsatellites were genotyped in a total of 394 individuals from two population pairs from the Iberian Peninsula, formed by one continuous and one fragmented population each. In agreement with predictions, SGS was significant and stronger in fragments (Sp = 0.020 and Sp = 0.026) than in continuous populations, where significant SGS was detected for one population only (Sp = 0.010). Simulations suggested that under fat-tailed dispersal, small population size is a stronger determinant of SGS than genetic isolation, while under normal dispersal, genetic isolation has a stronger effect. SGS was always stronger in real populations than in simulations, except if unrealistically narrow dispersal and/or high variance of reproductive success were modelled (even when accounting for potential overestimation of SGS in real populations as a result of short-distance sampling). This suggests that factors such as nonrandom mating or selection not considered in the simulations were additionally operating on SGS in Iberian maritime pine populations.

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.

Is Cupressus sempervirens native in Italy? An answer from genetic and palaeobotanical data.

This study represents the first large-scale analysis using nuclear molecular markers to assess genetic diversity and structure of Cupressus sempervirens L.. Genetic and fossil data were combined to infer the possible role of human activity and evolutionary history in shaping the diversity of cypress populations. We analysed 30 populations with six polymorphic nuclear microsatellite markers. Dramatic reductions in heterozygosity and allelic richness were observed from east to west across the species range. Structure analysis assigned individuals to two main groups separating central Mediterranean and eastern populations. The two main groups could be further divided into five subgroups which showed the following geographical distributions: Turkey with the Greek islands Rhodes and Samos, Greece (Crete), Southern Italy, Northern Italy, Tunisia with Central Italy. This pattern of genetic structure is also supported by SAMOVA and Barrier analyses. Palaeobotanical data indicated that Cupressus was present in Italy in the Pliocene, Pleistocene and Holocene. Furthermore, our molecular survey showed that Italian cypress populations experienced bottlenecks that resulted in reduced genetic diversity and allelic richness and greater genetic differentiation. Recent colonization or introduction may also have influenced levels of diversity detected in the Italian populations, as most individuals found in this range today have multilocus genotypes that are also present in the eastern range of the species. The data reveal a new interpretation of the history of cypress distribution characterized by ancient eastern populations (Turkey and Greek islands) and a mosaic of recently introduced trees and remnants of ancient, depauperate populations in the central Mediterranean range.

Combined analysis of nuclear and mitochondrial markers provide new insight into the genetic structure of North European Picea abies.

Norway spruce of northern Europe expanded at the end of the last glacial out of one refugium in Russia. To provide a detailed insight into how the genetic structure in the northern European lineage of this species has been shaped by postglacial migration, recurrent pollen flow and marginality, we here compare variation at seven highly variable nuclear microsatellite loci in 37 populations (1715 trees) with mitochondrial DNA variation. Microsatellite diversity was high (H(E)=0.640) and genetic differentiation was low (F(ST)=0.029). The microsatellite structure supported a mitochondrial DNA (mtDNA)-based hypothesis of two migration routes out of a single Russian refugium; one northwestern over Finland to northern Scandinavia, and one southwestern across the Baltic Sea into southern Scandinavia. Microsatellite diversity was maintained along the southwestern migration routes, whereas a significant decrease was observed towards the north. In contrast, the mtDNA diversity suggested higher amounts of historical gene flow towards the north than along the southwestern migration route. This suggests that the loss of nuclear diversity after postglacial colonization has been efficiently replenished by pollen-mediated gene flow in the south. Towards the north, smaller effective population size because of more limited seed and pollen production may have caused decreased nuclear diversity and increased inbreeding, reflecting the ecological marginality of the species in the north.

Introgression from modern hybrid varieties into landrace populations of maize (Zea mays ssp. mays L.) in central Italy.

Landraces are domesticated local plant varieties that did not experience a deliberate and intensive selection during a formal breeding programme. In Europe, maize landraces are still cultivated, particularly in marginal areas where traditional farming is often practiced. Here, we have studied the evolution of flint maize landraces from central Italy over 50 years of on-farm cultivation, when dent hybrid varieties were introduced and their use was widespread. We have compared an 'old' collection, obtained during the 1950s, before the introduction of hybrids, and a recent collection of maize landraces. For comparison, a sample of maize landraces from north Italy, and of improved germplasm, including hybrids and inbred lines were also used. A total of 296 genotypes were analysed using 21 microsatellites. Our results show that the maize landraces collected in the last 5-10 years have evolved directly from the flint landrace gene pool cultivated in central Italy before the introduction of modern hybrids. The population structure, diversity and linkage disequilibrium analyses indicate a significant amount of introgression from hybrid varieties into the recent landrace populations. No evidence of genetic erosion of the maize landraces was seen, suggesting that in situ conservation of landraces is an efficient strategy for preserving genetic diversity. Finally, the level of introgression detected was very variable among recent landraces, with most of them showing a low level of introgression; this suggests that coexistence between different types of agriculture is possible, with the adoption of correct practices that are aimed at avoiding introgression from undesired genetic sources.