Potential refugia in Taiwan revealed by the phylogeographical study of Castanopsis carlesii Hayata (Fagaceae)

In this study, we examined spatial patterns of chloroplast DNA (cpDNA) variation in a total of 30 populations of Castanopsis carlesii Hayata (Fagaceae), a subtropical and temperate tree species, including 201 individuals sampled throughout Taiwan. By sequencing two cpDNA fragments using universal primers (the trnL intron and the trnV-trnM intergenic spacer), we found a total of 1663 bp and 21 polymorphic sites. These gave rise to a total of 28 cpDNA haplotypes. The level of differentiation among the populations studied was relatively high (GST = 0.723). Two ancestral haplotypes are widely distributed. The Central Mountain Ridge (CMR) of Taiwan represents an insurmountable barrier to the east-west gene flow of C. carlesii. Among the populations studied, three separated populations, at Lienhuachih, Fushan and Lichia, have high nucleotide diversity. Estimates of NST-GST for populations on both sides of the CMR indicate that no phylogeographical structure exists. According to the genealogical tree, number of rare haplotype and population genetic divergence, this study suggests that two potential refugia existed during the last glaciation: the first refugium was located in a region to the north of Hsuehshan Range (HR) and west of the CMR; the second refugium was located in south, especially southeastern Taiwan. In fact, the second refugium happens to be the same as that reported for Quercus glauca. A 'star-like' genealogy is characteristic when all haplotypes rapidly coalesce and is a general outcome of population expansion. The neutrality test and mismatch distribution also suggest demographic expansion recovering from a bottleneck.     

Spatial pattern of chloroplast DNA variation of Cyclobalanopsis glauca in Taiwan and East Asia

This study examined the spatial pattern of chloroplast DNA (cpDNA) variation in Cyclobalanopsis glauca (Thunb. ex Murray) Oerst. (Fagaceae) in 140 trees from Taiwan (25 populations), Japan (three), Ryukyus (two), Hong Kong (one) and Mainland China (one). By sequencing three cpDNA intergenic spacer fragments using universal primers (trnT-trnL, trnV-trnM, including the trnV intron, and petG-trnP), we found a total of 1,980 bp and 15 polymorphic sites. Among them, 12 sites were caused by point mutation, and three resulted from insertion. This gives rise to a total of 13 cpDNA haplotypes. The level of differentiation among the populations studied is relatively high (GST = 0.612). Two ancestral haplotypes (A and B) are distributed widely in East Asia. Interestingly, all the derived cpDNA variations are found only in Taiwan but not in other areas. The Central Mountain Ridge (CMR) of Taiwan creates an unsurpassed barrier to the east-west gene flow of C. glauca. Among the populations on the west of CMR, only three separated populations, Yangmingshan, Wushe and Chinshuiying, have high haplotype diversity, each consisting of sister haplotypes all mutated from the same ancestral haplotype. Thus, they have probably originated from de novo mutation after the last glaciation. This inference agrees with the observation that no spatial autocorrelation existed on the west side. Two unrelated dominant lineages on the east of the CMR (haplotypes D and F) showed significant spatial genetic structure. Estimate of NST - GST was -0.090 and differed significantly from zero. Thus at the local scale, the phylogeographical component of the genetic structure is significant on the east of the CMR. Accompanied by published palynological records of the last glaciation, this study suggests the possibility that these two types were colonized northward from the southeastern part of Taiwan. 'Star-like' genealogy is characterized, with all the haplotypes coalescing rapidly and as a general outcome of population expansion (Page & Holmes 1998). A neutrality test also suggested a demographic expansion recovered from a bottleneck. We therefore inferred that the southeastern part of Taiwan might be a potential refugium for C. glauca.     

Phylogeographic history of white spruce during the last glacial maximum: uncovering cryptic refugia

Although a recent study of white spruce using chloroplast DNA uncovered the presence of a glacial refuge in Alaska, chloroplast failed to provide information on the number or specific localities of refugia. Recent studies have demonstrated the utility of nuclear microsatellites to refine insights into postglacial histories. The greater relative rate of mutation may allow finer scale resolution of historic dynamics, including the number, location, and sizes of refugia. Genetic data were acquired from screening 6 microsatellite loci on approximately 14 trees from each of 22 populations located across the central and western boreal forests of Canada and Alaska. Our studies combining microsatellites with Bayesian analyses of population structure in white spruce support the phylogeographic patterns uncovered using chloroplast, separating Alaskan from non-Alaskan regions. Results also support the idea that north-central Alaska served as a glacial refugium during the last glacial maximum. Additionally, the relationship between the degree of genetic differentiation and geographic distance indicated that gene flow played a more important role in structuring non-Alaskan populations, whereas drift played a more important role in structuring Alaskan populations (R(ST)'s for non-Alaskan populations 0.029 ± 0.007 and 0.083 ± 0.012 for Alaskan populations). Microsatellites also substantiate the bidirectional patterns of gene flow previously uncovered using chloroplast DNA but indicate much greater movement and mixing. Results from our Bayesian analyses also suggest the existence of additional cryptic refugia. However, the locations have been obscured by high gene flow (R(ST) averaging 0.057 ± 0.004).     

Genetic evidence for glacial refugia of the temperate tree Eucryphia cordifolia (Cunoniaceae) in southern South America

? Premise of the study: The temperate forests of southern South America were greatly affected by glaciations. Previous studies have indicated that some cold-tolerant tree species were able to survive glacial periods in small, ice-free patches within glaciated areas in the Andes and in southern Patagonia. Here we asked whether populations of the mesothermic species Eucryphia cordifolia also were able to survive glaciations in these areas or only in unglaciated coastal areas. ? Methods: The chloroplast intergenic spacer trnV-ndhC was sequenced for 150 individuals from 22 locations. Genetic data were analyzed (standard indexes of genetic diversity, a haplotype network, and genetic differentiation) in a geographical context. ? Key results: Two of the nine haplotypes detected were widespread in high frequency across the entire range of the species. The highest levels of genetic diversity were found around 40°S, decreasing sharply northward and more moderately southward. No differences in genetic diversity were found between Andean and coastal populations. Notably, seven haplotypes were found in a small area of the Coast Range known as the Cordillera Pelada (40°S). The differentiation coefficients G(ST) and N(ST) revealed that most of the genetic variation detected was due to variation within populations. ? Conclusions: The low levels of population differentiation and the high genetic diversity found in the Cordillera Pelada suggest that this area was the main refugium for E. cordifolia during glaciations. Nevertheless, given the high levels of genetic diversity found in some Andean populations, we cannot discount that some local populations also survived the glaciation in the Andes.     

Phylogeography of the endangered Eryngium alpinum L. (Apiaceae) in the European Alps

We studied the phylogeography of Eryngium alpinum by sequencing two intergenic chloroplast spacers, trnH-psbA and trnS-trnG (1322 bp). The sampling design included 36 populations and 397 individuals spanning the entire distribution range of the species, from France to Bosnia. Twenty-one haplotypes were characterized and polymorphism was observed both within and among populations. Population differentiation was strong (F(ST) = 0.92) and largely explained by the distinction of five geographic regions: Southwestern, Western, Middle, Eastern Alps and Balkans (F(CT) = 0.62). Moreover, N(ST) was significantly higher than G(ST) (P < 0.05), showing the existence of a phylogeographic pattern. Six major lineages were recognized using samova and median-joining networks. One lineage, highly divergent from the other ones, was only found in the Balkans and probably persisted in situ during last glaciations. All other lineages might have survived in a Southwestern refugium (Mercantour) and colonized the entire Alpine arc (Southwestern, Western, Middle and Eastern Alps) through repeated colonization events at different time periods. This is the first empirical study suggesting Southern refugia for calcareous Alpine plants, although the existence of a secondary refugium in northern Italy/Austria is also suspected. We also observed recent haplotype diversification, especially in the Southwestern Alps.     

Late glacial history of the cold-adapted freshwater fish Cottus gobio,revealed by microsatellites

The distribution of genetic diversity at 10 highly polymorphic microsatellite loci within the European freshwater fish, Cottus gobio, L. was examined. The sampling range comprised a large geographical scale including lineages known to be highly divergent at both mitochondrial DNA (mtDNA) and allozymes. An analysis of genetic variability within populations showed that expected heterozygosity and allelic richness could be explained largely by current effective population sizes. Evidence was found, however, that historical processes predating the last major glaciation affected allelic richness. In addition to confirming the large-scale patterns from earlier studies, the microsatellite data revealed new insights into recent processes by analysing genetic structure within ancient lineages defined by mtDNA data. Stepwise mutation model (SMM) and nonSMM-based methods demonstrated a clear genetic structuring within the Northwestern European lineage comprising populations from Britain and Belgium, and within the Central European lineage populations from the rivers Danube, Elbe and Main. Supported by an analysis of genetic variability within populations these results showed that the bullhead populations most probably persisted throughout the last major glaciation within the British Isles and within the drainages of the rivers Elbe and Main. Such observations provide the first genetic evidence for a glacial refugium in such close proximity to the European glacial margins.     

Nuclear and chloroplast microsatellites reveal extreme population differentiation and limited gene flow in the Aegean endemic Brassica cretica (Brassicaceae)

Nuclear and chloroplast microsatellite markers were used to study population structure and gene flow among seven Cretan populations of the Aegean endemic plant species Brassica cretica (Brassicaceae). Both nuclear and chloroplast markers revealed exceptionally high levels of population differentiation (overall F(ST)=0.628 and 1.000, respectively) and relatively little within-population diversity (overall H(S)=0.211 and 0.000, respectively). Maximum-likelihood estimates of directional migration rates were low among all pairs of populations (average Nm=0.286). There was no evidence that differences in flower colour between populations had any influence on historical levels of gene flow. In addition, a haplotype network showed that all five chloroplast haplotypes found in the sample were closely related. Together, these results suggest that current patterns of diversification in B. cretica are mainly a result of genetic drift during the last half million years. The main conclusions from the present study are consistent with the prevailing hypothesis that plant diversification in the Aegean region is driven by random rather than adaptive differentiation among isolated populations.     

Contrasting evolutionary forces driving population structure at expressed sequence tag polymorphisms, allozymes and quantitative traits in white spruce

Patterns of variation in quantitative characters and genetic markers were compared among six regional populations of white spruce [Picea glauca (Moench) Voss]. Although some phenotypic characters were correlated with latitude (r = 0.791), longitude (r = -0.796) and precipitation during the growing season (r = 0.789), variability at genetic markers was not correlated with geographical or bioclimatic variables, and followed neutral expectations. Estimates of genetic diversity and population differentiation for 14 allozymes (translated regions of coding genes) were essentially indistinguishable from those observed for 11 expressed sequence tag polymorphisms (ESTPs) from untranslated regions of coding genes. Variation among populations for quantitative traits such as eighth year height (Q(ST) = 0.082), thirteenth year height (Q(ST) = 0.069), total wood density (Q(ST) = 0.102) and date of budset (Q(ST) = 0.246), was greater than for allozymes (G(ST) = 0.014) and ESTPs (G(ST) = 0.019). These trends suggest a strong adaptive response in quantitative traits, contrasting to allozymes and ESTPs where no selective response could be detected and where populations appeared to be essentially in a migration-drift equilibrium.     

Inferring the invasion history of coral berry Ardisia crenata from China to the USA using molecular markers

Genetic comparisons between native and invasive populations of a species can provide insights into its invasion history information, which is useful for guiding management and control strategies. The coral berry Ardisia crenata was introduced to Florida last century as a cultivated ornament plant, and has since spread widely throughout the southern regions of the USA. Previously, the genetic variation among 20 natural populations of A. crenata across its distribution center in southern China was quantified using seven microsatellite markers. Here we expand on that work by additionally sampling individuals from four other native populations in Taiwan and Japan, and from five invasive populations in the USA. We also examined the results from one chloroplast intergenic spacer region (trnF-trnL) in all 29 populations. Our aim is to identify the invasion source and subsequent history of the species?? spread throughout the southern USA. We observed lower genetic diversity in the invasive populations based on both microsatellite and chloroplast markers. Our data show that the invasive populations can be clustered with native populations in southeastern China, inferring this region as the geographic origin of A. crenata cultivars invading the USA. We further classified invasive individuals into invasive I and invasive II clusters. Nantou in Taiwan and Xihu in mainland China are the most closely related populations to those, which identify the former as potential sources for host-specific control agents. Our results, combined with the known introduction records, suggest that A. crenata was first multiply introduced into Florida and then secondarily colonized Louisiana and Texas from Florida.     

If FST does not measure neutral genetic differentiation,then comparing it with QST is misleading. Or is it?

The comparison between neutral genetic differentiation (F(ST) ) and quantitative genetic differentiation (Q(ST) ) is commonly used to test for signatures of selection in population divergence. However, there is an ongoing discussion about what F(ST) actually measures, even resulting in some alternative metrics to express neutral genetic differentiation. If there is a problem with F(ST) , this could have repercussions for its comparison with Q(ST) as well. We show that as the mutation rate of the neutral marker increases, F(ST) decreases: a higher within-population heterozygosity (He) yields a lower F(ST) value. However, the same is true for Q(ST) : a higher mutation rate for the underlying QTL also results in a lower Q(ST) estimate. The effect of mutation rate is equivalent in Q(ST) and F(ST) . Hence, the comparison between Q(ST) and F(ST) remains valid, if one uses neutral markers whose mutation rates are not too high compared to those of quantitative traits. Usage of highly variable neutral markers such as hypervariable microsatellites can lead to serious biases and the incorrect inference that divergent selection has acted on populations. Much of the discussion on F(ST) seems to stem from the misunderstanding that it measures the differentiation of populations, whereas it actually measures the fixation of alleles. In their capacity as measures of population differentiation, Hedrick's G'(ST) and Jost's D reach their maximum value of 1 when populations do not share alleles even when there remains variation within populations, which invalidates them for comparisons with Q(ST) .     

Phylogeographic structure of jack pine (Pinus banksiana; Pinaceae) supports the existence of a coastal glacial refugium in northeastern North America

? Premise of the study: The genetic structure of jack pine (Pinus banksiana Lamb.), a North American boreal conifer with large longitudinal distribution, was investigated to test for the possible existence of a genetically distinct lineage in the Maritimes region in northeastern North America, which could be indicative of a mid-latitude coastal refuge during the last glaciation. ? Methods: One maternally inherited mitochondrial DNA (mtDNA) minisatellite marker and four paternally inherited chloroplast DNA (cpDNA) microsatellite markers were used to assess the range-wide geographical structure of jack pine populations with particular focus on northeastern North America. ? Key results: The populations from the Maritimes region presented a unique mtDNA background characterized by very low diversity and the preponderance of a distinctive mitotype. The distribution of cpDNA diversity was not spatially structured, though three chlorotypes were restricted to the east. ? Conclusions: MtDNA data suggest that populations from the Maritimes region derive from a genetically depauperated north-coastal refugium. Contrastingly, the much higher geographical uniformity observed for cpDNA variation indicates that gene flow by pollen had been much more effective than seed gene flow at homogenizing population structure.     

Decoupled mitochondrial and chloroplast DNA population structure reveals Holocene collapse and population isolation in a threatened Mexican-endemic conifer

Chihuahua spruce (Picea chihuahuana Martínez) is a montane subtropical conifer endemic to the Sierra Madre Occidental in northwestern México. Range-wide variation was investigated using maternally inherited mitochondrial (mtDNA) and paternally inherited chloroplast (cpDNA) DNA markers. Among the 16 mtDNA regions analysed, only two mitotypes were detected, while the study of six cpDNA microsatellite markers revealed eight different chlorotypes. The average cpDNA diversity (H = 0.415) was low but much higher than that for mtDNA (H = 0). The distribution of mitotypes revealed two clear nonoverlapping areas (G(ST) = N(ST) = 1), one including northern populations and the second one including the southern and central stands, suggesting that these two regions may represent different ancestral populations. The cpDNA markers showed lower population differentiation (G(ST) = 0.362; R(ST) = 0.230), implying that the two ancestral populations continued to exchange pollen after their initial geographic separation. A lack of a phylogeographic structure was revealed by different spatial analyses of cpDNA (G(ST) > R(ST); and samova), and reduced cpDNA gene flow was noted among populations (Nm = 0.873). Some stands deviated significantly from the mutation-drift equilibrium, suggesting recent bottlenecks. Altogether, these various trends are consistent with the hypothesis of a population collapse during the Holocene warming and suggest that most of the modern P. chihuahuana populations are now effectively isolated with their genetic diversity essentially modelled by genetic drift. The conservation efforts should focus on most southern populations and on the northern and central stands exhibiting high levels of genetic diversity. Additional mtDNA sequence analysis confirmed that P. martinezii (Patterson) is not conspecific with P. chihuahuana, and thus deserves separate conservation efforts.     

Comparing Bayesian estimates of genetic differentiation of molecular markers and quantitative traits: an application to Pinus sylvestris

Comparison of the level of differentiation at neutral molecular markers (estimated as F(ST) or G(ST)) with the level of differentiation at quantitative traits (estimated as Q(ST)) has become a standard tool for inferring that there is differential selection between populations. We estimated Q(ST) of timing of bud set from a latitudinal cline of Pinus sylvestris with a Bayesian hierarchical variance component method utilizing the information on the pre-estimated population structure from neutral molecular markers. Unfortunately, the between-family variances differed substantially between populations that resulted in a bimodal posterior of Q(ST) that could not be compared in any sensible way with the unimodal posterior of the microsatellite F(ST). In order to avoid publishing studies with flawed Q(ST) estimates, we recommend that future studies should present heritability estimates for each trait and population. Moreover, to detect variance heterogeneity in frequentist methods (ANOVA and REML), it is of essential importance to check also that the residuals are normally distributed and do not follow any systematically deviating trends.     

Rapid divergence and postglacial colonization in western North American Steller's jays (Cyanocitta stelleri)

Post-Pleistocene avian colonization of deglaciated North America occurred from multiple refugia, including a coastal refugium in the northwest. The location of a Pacific Coastal refugium is controversial; however, multiple lines of evidence suggest that it was located near the Queen Charlotte Islands (also known as Haida Gwaii). The Queen Charlotte Islands contain a disproportionately large number of endemic plants and animals including the Steller's jay Cyanocitta stelleri carlottae. Using five highly variable microsatellite markers, we studied population structure among eight populations of Steller's jay (N = 150) from geographical areas representing three subspecies in western North America: C. s. carlottae, C. s. stelleri and C. s. annectens. Microsatellite analyses revealed genetic differentiation between each of the three subspecies, although more extensive sampling of additional C. s. annectens populations is needed to clarify the level of subspecies differentiation. High levels of population structure were found among C. s. stelleri populations with significant differences in all but two pairwise comparisons. A significant isolation by distance pattern was observed amongst populations in the Pacific Northwest and Alaska. In the C. s. carlottae population, there was evidence of reduced genetic variation, higher number of private alleles than northern C. s. stelleri populations and higher levels of divergence between Queen Charlotte Island and other populations. We were unable to reject the hypothesis that the Queen Charlotte Islands served as a refugium during the Pleistocene. Steller's jay may have colonized the Queen Charlotte Islands near the end of the last glaciation or persisted throughout the Pleistocene, and this subspecies may thus represent a glacial relic. The larger number of private alleles, despite reduced genetic variation, morphological distinctiveness and high divergence from other populations suggests that the Queen Charlotte Island colonization pre-dates that of the mainland. Furthermore, our results show rapid divergence in Steller's jay populations on the mainland following the retreat of the ice sheets.     

Influence of glaciation on divergence patterns of the endemic minnow, Zacco pachycephalus, in Taiwan

The genetic subdivision of the endemic minnow, Zacco pachycephalus, in 12 rivers of Taiwan was studied by allozyme electrophoresis. Among 26 loci surveyed, six were significantly differentiated among sites. Allozyme genotype frequencies within samples accord with Hardy-Weinberg expectations. A highly divergent genetic structure among sites (FST = 0.497) with low genetic variability within samples (H = 0.006-0.062, mean = 0.03) suggests that local populations originated from a small number of founders. UPGMA and Fitch trees derived from Nei's genetic distance and hierarchical analysis of FST values reveal that samples of Z. pachycephalus can be divided into northern, middle, and southern subgroups and that most genetic variation (87%) is distributed among rather than within groups. It is proposed that there were two isolated refugia during the most recent glacial period, and that populations of the southern group (Group I) originated from a southern refugium in Taiwan below the boundary of the Kaoping River, and other groups (Group II), including northern and middle subgroups, originated from a refugium located on the extinct land-bridge between Taiwan and mainland China. The postglacial division between northern and middle subgroups may be due to hydrological differences between these two regions.     

Adaptive divergence in moor frog (Rana arvalis) populations along an acidification gradient: inferences from Q(st) -F(st) correlations

Microevolutionary responses to spatial variation in the environment seem ubiquitous, but the relative role of selection and neutral processes in driving phenotypic diversification remain often unknown. The moor frog (Rana arvalis) shows strong phenotypic divergence along an acidification gradient in Sweden. We here used correlations among population pairwise estimates of quantitative trait (P(ST) or Q(ST) from common garden estimates of embryonic acid tolerance and larval life-history traits) and neutral genetic divergence (F(ST) from neutral microsatellite markers), as well as environmental differences (pond pH, predator density, and latitude), to test whether this phenotypic divergence is more likely due to divergent selection or neutral processes. We found that trait divergence was more strongly correlated with environmental differences than the neutral marker divergence, suggesting that divergent natural selection has driven phenotypic divergence along the acidification gradient. Moreover, pairwise P(ST) s of embryonic acid tolerance and Q(ST) s of metamorphic size were strongly correlated with breeding pond pH, whereas pairwise Q(ST) s of larval period and growth rate were more strongly correlated with geographic distance/latitude and predator density, respectively. We suggest that incorporating measurements of environmental variation into Q(ST) -F(ST) studies can improve our inferential power about the agents of natural selection in natural populations.     

Microsatellite variation and differentiation in African and non-African populations of Drosophila simulans

Drosophila simulans originated in sub-Saharan Africa or Madagascar and colonized the rest of the world after the last glaciation about 10 000 years ago. Consistent with this demographic history, sub-Saharan African populations have been shown to harbour higher levels of microsatellite and sequence variation than cosmopolitan populations. Nevertheless, only limited information is available on the population structure of D. simulans. Here, we analysed X-linked and autosomal microsatellite loci in four sub-Saharan African, one North African, one Israeli, and two European D. simulans populations. Bayesian clustering algorithms combined the North African, Israeli, and European populations into a single cosmopolitan group. The four sub-Saharan populations were split into two separate groups. Pairwise F(ST) analysis, however, indicated significant population differentiation between all eight populations surveyed. A significant signal for population reduction in cosmopolitan populations was found only for X-linked loci.     

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

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

Geographic patterns of genetic variation in nuclear and chloroplast genomes of two related oaks (<Emphasis Type="Italic">Quercus aliena</Emphasis> and <Emphasis Type="Italic">Q. serrata</Emphasis>) in Japan: implications for seed and seedling transfer

In this study, we assessed geographic patterns of genetic variations in nuclear and chloroplast genomes of two related native oaks in Japan, Quercus aliena and Q. serrata, in order to facilitate development of genetic guidelines for transfer of planting stocks for each species. A total of 12 populations of Q. aliena and 44 populations of Q. serrata were analyzed in this study. Genotyping of nuclear microsatellites in Q. aliena was done with only nine populations (n = 212) due to limited numbers of individuals in two populations, while all 12 populations (n = 89) were used in sequencing chloroplast DNA (cpDNA). In Q. serrata, 43 populations (n = 1032) were genotyped by nuclear microsatellite markers, while cpDNA of 44 populations (n = 350) was sequenced. As anticipated, geographic patterns detected in the variations of Q. aliena’s nuclear genome and its chloroplast haplotype distribution clearly distinguished northern and southern groups of populations. However, those of Q. serrata were inconsistent. The geographic distribution of its chloroplast haplotypes tends to show the predicted differentiation between northern and southern lineages, but geographic signals in the genetic structure of its nuclear microsatellites are weak. Therefore, treating northern and southern regions of Japan as genetically distinct transferrable zones for planting stocks is highly warranted for Q. aliena. For Q. serrata, the strong NE-SW geographic structure of cpDNA should be considered.     

Genetic isolation of fragmented populations is exacerbated by drift and selection

Reduced genetic variation at marker loci in small populations has been well documented, whereas the relationship between quantitative genetic variation and population size has attracted little empirical investigation. Here we demonstrate that both neutral and quantitative genetic variation are reduced in small populations of a fragmented plant metapopulation, and that both drift and selective change are enhanced in small populations. Measures of neutral genetic differentiation (F(ST)) and quantitative genetic differentiation (Q(ST)) in two traits were higher among small demes, and Q(ST) between small populations exceeded that expected from drift alone. This suggests that fragmented populations experience both enhanced genetic drift and divergent selection on phenotypic traits, and that drift affects variation in both neutral markers and quantitative traits. These results highlight the need to integrate natural selection into conservation genetic theory, and suggests that small populations may represent reservoirs of genetic variation adaptive within a wide range of environments.