Ancestry and divergence of subtropical montane forest isolates: molecular biogeography of the genus Abies (Pinaceae) in southern México and Guatemala

The genus Abies has a complex history in southern México and Guatemala. In this region, four closely related species, Abies flinckii , A. guatemalensis , A. hickelii , and A. religiosa , are distributed in fragmented and isolated montane populations. Range-wide genetic variation was investigated across species using cytoplasmic DNA markers with contrasted inheritance. Variation at two maternally inherited mitochondrial DNA markers was low. All species shared two of the nine mitotypes detected, while the remaining seven mitochondrial DNA types were restricted to a few isolated stands. Mitochondrial genetic differentiation across taxa was high ( G _ST = 0.933), it was not related to the taxonomic identity ( amova ; P  > 0.05) of the populations, and it was not phylogeographically structured ( G _ST ≈  N _ST). In contrast, variation at three paternally inherited chloroplast DNA microsatellites was high. Chloroplast genetic differentiation was lower ( G _ST = 0.402; R _ST = 0.547) than for mitochondrial DNA, but it was significantly related to taxonomy ( amova ; P  < 0.001), and exhibited a significant phylogeographical structure ( G _ST <  R _ST). Different analyses of population structure indicated that A. flinckii was the most divergent taxon, while the remaining three species formed a relatively homogeneous group. However, a small number of the populations of these three taxa, all located at the limits of their respective ranges or in the Transverse Volcanic Belt, diverged from this main cluster. These trends suggest that the Mesoamerican Abies share a recent common ancestor and that their divergence and speciation is mainly driven by genetic drift and isolation during the warm interglacial periods.     

Subtle population structuring within a highly vagile marine invertebrate, the veined squid Loligo forbesi, demonstrated with microsatellite DNA markers

Microsatellite DNA markers were applied for the first time in a population genetic study of a cephalopod and compared with previous estimates of genetic differentiation obtained using allozyme and mitochondrial DNA (mtDNA) markers. Levels of genetic variation detected with microsatellites were much higher than found with previous markers (mean number of alleles per locus=10.6, mean expected heterozygosity ( H _E)=0.79; allozyme H _E=0.08; mtDNA restriction fragment length polymorphism (RFLP) H _E=0.16). In agreement with previous studies, microsatellites demonstrated genetic uniformity across the population occupying the European shelf seas of the North East Atlantic, and extreme genetic differentiation of the Azores population ( R _ST/ F _ST=0.252/0.245; allozyme F _ST=0.536; mtDNA F _ST=0.789). In contrast to other markers, microsatellites detected more subtle, and significant, levels of differentiation between the populations of the North East Atlantic offshore banks (Rockall and Faroes) and the shelf population ( R _ST=0.048 and 0.057). Breakdown of extensive gene flow among these populations is indicated, with hydrographic (water depth) and hydrodynamic (isolating current regimes) factors suggested as possible barriers to migration. The demonstration of genetic subdivision in an abundant, highly mobile marine invertebrate has implications for the interpretation of dispersal and population dynamics, and consequent management, of such a commercially exploited species. Relative levels of differentiation indicated by the three different marker systems, and the use of measures of differentiation (assuming different mutation models), are discussed.     

Genetic variation at allozyme and RAPD loci in sessile oak Quercus petraea (Matt.) Liebl.: the role of history and geography

The nuclear genetic variation within and among 21 populations of sessile oak was estimated at 31 RAPD loci in conjunction with previous estimates of variation at eight allozyme loci. The aim of the study was to assess the relative role of isolation-by-distance and postglacial history on patterns of nuclear variation. Because of its small effective population size and maternal transmission, the chloroplast genome is a good marker of population history. Both kinds of nuclear variation (RAPD and allozyme) were therefore compared, first, to the geographical distances among populations and, secondly, to chloroplast DNA restriction polymorphism in the same populations. Multiple Mantel tests were used for this purpose. Although RAPDs revealed less genetic diversity than allozymes, levels of genetic differentiation ( G _ST) were identical. The standard genetic distance calculated at all RAPD loci was correlated with geographical distances but not with the genetic distance calculated from chloroplast DNA data. Conversely, allozyme variation was correlated with chloroplast DNA variation, but not with geography. Possibly, divergent selection at two allozyme loci during the glacial period could explain this pattern. Because of its greater number of loci assayed, RAPDs probably provided a less biased picture of the relative role of geography and history.     

Genetic Variation and Clonal Diversity of the Two Divergent Types of Clonal Populations of Leymus chinensis Tzvel on the Song Liao Steppe in the West of Northeastern China

The genetic variation and clonal diversity of two divergent types (grey-green and yellow-green) of clonal populations of Leymus chinensis Tzvel at 14 loci were compared. Total gene diversity (HT) and the coefficient of genetic differentiation (GST) were all higher for the yellow-green type (HT = 0.270; GST =0.186) than for the grey-green type (HT = 0.250; GST = 0.157) of L. chinensis. Rare alleles usually occurred as heterozygotes rather than homozygotes and significant deviations from Hardy-Weinberg equilibrium were found only at a few loci. This indicated that these two types of populations were mainly out-crossing. Clonal diversity, evenness of clones, and mean clone size were not significantly different between the two types. We found that differences between the clone size and genetic variation of the yellow-green type of populations occurred with different climate and habitat population groups. However, for the grey-green type of populations, these genetic variations decreased under conditions of different climate and habitat population groups.     

Allozyme Diversity in Populations of Cymbidium goeringii (Orchidaceae)

Abstract: Using 14 allozyme loci, we investigated levels of genetic diversity within populations, and degree of genetic divergence among 24 populations of Cymbidium goeringii (Orchidaceae) in Korea and Japan. Cymbidium goeringii maintains high levels of genetic diversity both at population (mean expected heterozygosity, H _e = 0.238) and species levels (0.260). Means of H _e found in 24 populations were not significantly different from each other. About 90 % of the total variation in the species is common to all populations (mean G _ST = 0.108). No unique allele was found in any population. The indirect estimate of gene flow based on the mean G _ST was high ( Nm = 2.06). Nei's genetic identities for pairs of populations had high values (mean = 0.974 [SD = 0.013]). The Mantel-Z test showed a significant correlation between genetic distance and geographic distance. However, the mean G _ST value between 17 populations in Korea and seven Japanese populations was relatively low (0.029), even though the land connection between the southern Korean peninsula and southern Japanese archipelagos has not existed since the middle Pleistocene. Large numbers of small seeds of C. goeringii might travel long distances by wind from populations to populations both in Korea and Japan, increasing genetic diversity within populations and maintaining low genetic differentiation among populations.     

Chloroplast DNA phylogeography of the argan tree of Morocco

Polymorphisms in the chloroplast genome of the argan tree (Sapotaceae), an endemic species of south-western Morocco, have been detected by restriction site studies of PCR-amplified fragments. A total of 12 chloroplast DNA (cpDNA) and two mitochondrial DNA (mtDNA) fragments were amplified and digested with a single restriction enzyme ( Hin fI). Polymorphisms were identified in six of the cpDNA fragments, whereas no mtDNA polymorphisms were detected in a survey of 95 individuals from 19 populations encompassing most of the natural range of the species. The cpDNA polymorphisms allowed the identification of 11 haplotypes. Two lineages, one in the south-east and the other in the north-west, divide the range of the argan tree into two distinct areas. The level of genetic differentiation measured at the haplotype level ( G _STc= 0.60) (i.e. with unordered haplotypes) was smaller than when phylogenetic relationships were taken into account ( N _STc= 0.71–0.74) (ordered haplotypes), indicating that population history must be considered in the study of the geographical distribution of cpDNA lineages in this species. If contrasted with the level of nuclear genetic differentiation measured in a previous study with isozymes ( G _STn= 0.25), the results indicate a relatively high level of gene flow by seeds, or conversely a relatively low level of gene flow by pollen, as compared with other tree species. Goats and camels could have played an important role in disseminating the fruits of this tree.     

Population divergence in the amphicarpic species Amphicarpaea edgeworthii Benth. (Fabaceae): microsatellite markers and leaf morphology

Comparative analyses of the genetic differentiation in microsatellite markers ( F _ST) and leaf morphology characters ( Q _ST) of Amphicarpaea edgeworthii Benth. were conducted to gain insight into the roles of random processes and natural selection in the population divergence. Simple sequence repeat analyses on 498 individuals of 19 natural populations demonstrate that a significant genetic differentiation occurs among populations (mean F _ST = 0.578), and A. edgeworthii is a highly self-fertilized species (mean selfing rate s  = 0.989). The distribution pattern of genetic diversity in this species shows that central populations possess high genetic diversity (e.g. population WL with H _E = 0.673 and population JG with H _E = 0.663), whereas peripheral ones have a low H _E as in population JD (0.011). The morphological divergence of leaf shape was estimated by the elliptical Fourier analysis on the data from 11 natural and four common garden populations. Leaf morphology analyses indicate the morphological divergence does not show strong correlation with the genetic differentiation ( R  = 0.260, P  = 0.069). By comparing the 95% confidence interval of Q _ST with that of F _ST, Q _ST values for five out of 12 quantitative traits are significantly higher than the average F _ST value over eight microsatellite loci. The comparison of F _ST and Q _ST suggests that two kinds of traits can be driven by different evolutionary forces, and the population divergence in leaf morphology is shaped by local selections.  © 2009 The Linnean Society of London, Biological Journal of the Linnean Society , 2009, 96 , 505–516.     

PATTERNS OF GENETIC ARCHITECTURE FOR LIFE-HISTORY TRAITS AND MOLECULAR MARKERS IN A SUBDIVIDED SPECIES

Abstract Understanding the utility and limitations of molecular markers for predicting the evolutionary potential of natural populations is important for both evolutionary and conservation genetics. To address this issue, the distribution of genetic variation for quantitative traits and molecular markers is estimated within and among 14 permanent lake populations of Daphnia pulicaria representing two regional groups from Oregon. Estimates of population subdivision for molecular and quantitative traits are concordant, with Q _ST generally exceeding G _ST. There is no evidence that microsatellites loci are less informative about subdivision for quantitative traits than are allozyme loci. Character-specific comparison of Q _ST and G _ST support divergent selection pressures among populations for the majority of life-history traits in both coast and mountain regions. The level of within-population variation for molecular markers is uninformative as to the genetic variation maintained for quantitative traits. In D. pulicaria , regional differences in the frequency of sex may contribute to variation in the maintenance of expressed within-population quantitative-genetic variation without substantially impacting diversity at the genic level. These data are compared to an identical dataset for 17 populations of the temporary-pond species, D. pulex .     

Inheritance of Anthracnose Resistance in the Common Bean Differential Cultivar G 2333 and Identification of a New Molecular Marker Linked to the Co-42 gene

The inheritance of anthracnose resistance of the common bean ( Phaseolus vulgaris L.) differential cultivar G 2333 to Colletotrichum lindemuthianum races 73 and 89 was studied in crosses with the susceptible cultivar Rudá. The segregation ratios of 15 : 1 in the F₂ and 3 : 1 in the backcrosses to Rudá indicate that for each of the races tested there are two independent resistance loci in G 2333. A random amplified polymorphic DNA (RAPD) molecular marker (OPH18_1200C) linked in resistance to race 73 was identified in a BC₃F_2:3 population derived from crosses between Rudá and G 2333. A RAPD molecular marker OPAS13_950C, previously identified as linked to gene Co-4² , was also amplified in this population. Co-segregation analyses showed that these two markers are located at 5.6 (OPH18_1200C) and 11.2 (OPAS13_950C) cM of the Co-4² gene. These markers were not present in BC₁F_2:3 plants resistant to race 89 indicating that this population carries a different resistance gene. DNA amplification of BC₁F_2:3 plants with RAPD molecular marker OPAB_450C, previously identified as linked to gene Co-5 , indicated that this gene is present in this population.     

The distribution of random amplified polymorphic DNA (RAPD) diversity amongst populations of Isotoma petraea (Lobeliaceae)

RAPDs were generated from plants of six populations of Isotoma petraea F. Muell. The species occurs on rock outcrops in southern and western Australia, with populations exhibiting different breeding systems, including complete autogamy, varying levels of outbreeding and complex hybridity. Non-metric multidimensional scaling (nMDS) analysis of the random amplified polymorphic DNA (RAPD) data set clearly resolved all populations. The Pigeon Rock population, which is home to both complex hybrid and structural homozygote plants, was divided into those two groups by the nMDS analysis. There was little diversity in highly autogamous populations, but levels were higher in the outbred Yackeyackine population. All complex hybrid populations and plants possessed numerous genetic system-specific RAPDs, some of which were shown to be held in fixed heterozygosity. Estimating G _ST using RAPDs has been problematical due to their dominance, and analytical methods usually rely on knowledge of the selfing rate or assume Hardy–Weinberg equilibrium. This assumption does not hold when populations exhibit fixed heterozygosity, and an alternative method, Shannon's Index, was used to partition genetic diversity. The distribution of genetic diversity fit expectations for an inbreeding species, with most of the variation (87.5%) occurring between populations. This compares to an average RAPD-based G _ST of 59.6% for inbreeding species generally and 15.5% for outbreeding species.     

POPULATION STRUCTURE IN AN INSHORE CETACEAN REVEALED BY MICROSATELLITE AND mtDNA ANALYSIS: BOTTLENOSE DOLPHINS (TURSIOPS SP.) IN SHARK BAY, WESTERN AUSTRALIA

We examined population substructure of bottlenose dolphins ( Tursiops sp). in Shark Bay, Western Australia, using 10 highly polymorphic microsatellite loci, and mitochondrial DNA (mtDNA). For microsatellite analysis, 302 different animals were sampled from seven localities throughout the bay. Analysis of genetic differentiation between sampling localities showed a significant correlation between the number of migrants ( Nm ) calculated from F _ST, R _ST and private alleles, and distance between localities–a pattern of isolation-by-distance. For mtDNA, 220 individuals from all seven localities were sequenced for a 351 base pair fragment of the control region, resulting in eight haplotypes, with two distinct clusters of haplotypes. Values of F _ST and (φ)_ST for mtDNA yielded statistically significant differences, mostly between localities that were not adjacent to each other, suggesting female gene flow over a scale larger than the sampled localities. We also observed a significant correlation between the number of female migrants calculated from F _ST and φ_ST and the distance of sampling localities. Our results indicate that dispersal in female dolphins in Shark Bay is more restricted than that of males.     

Conservation Genetics of an Endangered Herb, Hanabusaya asiatica (Campanulaceae)

Abstract: Hanabusaya asiatica (Nakai) Nakai (Campanulaceae), a bee- pollinated, perennial herb, is restricted to the mountainous regions of the eastern-central Korean peninsula. Allozyme analyses for 348 individuals assessed the levels of genetic diversity for five populations. Spatial autocorrelation statistics were also used to examine the spatial distribution of allozyme polymorphisms. The species maintains high levels of allozyme diversity ( H _eS = 0.217) and it exhibits low allozyme differentiation among populations ( G _ST = 0.132) compared with other endemics (mean H _e = 0.096, G _ST = 0.248). There is an apparent pattern of isolation by distance among populations. These results suggest that H. asiatica is at a genetic equilibrium. A considerable deficit in numbers of heterozygotes suggests mating among relatives in populations. At least three populations of H. asiatica should be sampled or conserved to capture or maintain > 99 % of the genetic diversity in the species as a whole. Within local populations, individuals are distributed in a structured, isolation by distance, manner. Approximate genetic patch width in the populations of H. asiatica examined is 5 - 8 m. For conservation purposes, it is suggested that, in general, the sampling of H. asiatica should be conducted at intervals in order to efficiently sample the genetic diversity across an entire population.     

Population genetics of a cyclostome species pair, river lamprey (Lampetra fluviatilis L.) and brook lamprey (Lampetra planeri Bloch)

Horizontal agarose gel electrophoresis of 24 allozyme loci in four species of Central European lampreys (321 Lampetra planeri , 83 L. fluviatilis , 11 Eudontomyzon mariae and nine Petromyzon marinus ) was used to study the 'paired species' L. fluviatilis and L. planeri . The genetic differentiation of the anadromous river lamprey ( L. fluviatilis ) from the stationary brook lamprey ( L. planeri ) was within the range of ingroup differentiation of the latter, but L. fluviatilis exhibited much greater population cohesion over a more extended geographic range: G _ST = 0.0537 versus G _ST = 0.3398, N _em = 4.402 versus N _em = 0.4856, mean genetic among-stock distances D = 0.0047 versus D = 0.0257. L. planeri populations coexisting geographically with L. fluviatilis in the Rhine and Elbe river systems were genetically more cohesive than L. planeri stocks from the Danubian basin where L. fluviatilis is absent. Danubian L. planeri populations exhibit a lower degree of heterozygosity than brook lampreys from the Rhine river system, but comprise deeper genetic lineages ( G _ST = 0.4629 versus G _ST = 0.2434), despite being sampled from a much more restricted area. Isolation-by-distance is observed for L. planeri from the Danubian but not from the Atlantic drainage basins. Transspecific gene flow between L. planeri from Atlantic drainage basins and the long-distant migrating L. fluviatilis is inferred, raising doubt on the validity of two separate biospecies. E. mariae and P. marinus are clearly differentiated from Lampetra spp. at several allozyme loci.     

Identification, cloning, and characterization of microsatellite DNA in Euglena gracilis

Microsatellite DNA has been developed into one of the most popular genetic markers. We have identified and cloned microsatellite loci in the genome of a free-living protozoan Euglena gracilis FACHB-848, using the random amplified microsatellites method (RAMS). The digoxigenin-labelled oligonucleotides (CT)10 and (GT)10 served as probes to detect complementary sequences in the randomly amplified polymorphic DNA (RAPD) fingerprints produced by means of Southern blotting. Subsequently, positive RAPD fragments were cloned. From a total of 31 RAPD primer profiles, eight microsatellite loci of E. gracilis were detected and characterized. Further, six sites (i.e. EGMS1, EGMS3, EGMS4, EGMS5, EGMS6, and EGMS7) showed polymorphisms. We found a GT or CT microsatellite every 10.5 kb in the genome of E. gracilis, and similar to animal genomes, the (GT)(n) motif was much more abundant than the (CT)(n) motif. These polymorphic microsatellite DNA will serve as advantageous molecular markers for studying the genetic diversity and molecular ecology of Euglena.     

Mitochondrial DNA analysis of the geographic structure of Indian scad mackerel in the Indo-Malay archipelago

Sequence analysis of a fragment of the cytochrome b gene in Decapterus russelli sampled from the Indo-Malay archipelago revealed two distinct mitochondrial clades (mean nucleotide divergence=2·2%) whose geographic distribution was heterogeneous (Nei's G _ST=0·416). This pointed to a complex pattern of genetic differentiation and demonstrated limited genetic exchange between populations in this highly mobile species.     

Evidence for new nuclear and mitochondrial genome organizations among high-frequency somatic embryogenesis-derived plants of allotetraploid Coffea arabica L. (Rubiaceae)

 The most important commercial species of coffee, Coffea arabica, which produces 73% of the world's coffee crop and almost all of the coffee in Latin America, is the only tetraploid (allotetraploid, 2n=4x=44) species known in the genus. High-frequency somatic embryogenesis, plant regeneration and plant recovery were achieved from leaf explants of a mature, elite plant of C. arabica cv. Cauvery (S-4347) using a two-step culture method. To assess the genetic integrity of the nuclear, mitochondrial and chloroplast genomes among the hardened regenerants, we employed multiple DNA markers (RFLP, RAPD, ISSR) for sampling various regions of the genome. Although the nuclear and mitochondrial genomes of the mother plant and five ramets derived from the mother ortet were similar in organization, this was not so in the somatic embryo-derived plants where both nuclear and mitochondrial genomes changed in different, characteristic ways and produced novel genome organizations. A total of 480 genetic loci, based on the data obtained from a total of 16 nuclear, mitochondrial and chloroplast gene probes, in combination with nine restriction enzyme digests, 38 RAPD and 17 SSR primers, were scored in 27 somatic embryo-derived plants and the single control. Among these, 44 loci were observed to be polymorphic. A relatively low level of polymorphism (4.36%) was found in the nuclear genome, while polymorphism in the mitochondrial genome (41%) was much higher. No polymorphism was detected in the chloroplast genome. The polymorphism in the mitochondrial genome was found in only 4 plants. Such selective polymorphism was not true for the nuclear genome. Thus, this in-depth and comprehensive study demonstrates, for the first time, the presence of subtle genetic variability and novel genome organizations in the commercially well-established somatic embryogenesis-derived plants of this important coffee species. Received: 2 July 1999 / Revision received: 1 February 2000 / Accepted: 17 February 2000     

DNA fragments of organellar origin in random amplified polymorphic DNA (RAPD) patterns of sugar beet (Beta vulgaris L.)

The technique of random amplified polymorphic DNA (RAPD) offers a broad range of applications in the investigation of plant genomes. A promising prospect is the use of RAPD products as genetic markers. We have investigated a possible organellar source of fragments in RAPD patterns of total DNA. Two nearly-isogenic lines of cytoplasmic male-sterile and male-fertile sugar beet (Beta vulgaris L.) were subjected to RAPD analysis with six different primers. Total, nuclear, mitochondrial (mt), and chloroplast (cp), DNA from each line were investigated. Reproducible DNA fingerprints could be obtained from both organellar DNAs. Differences in band patterns of mtDNA between cytoplasmic male-sterile and -fertile lines were observed with five out of six primers, whereas different cpDNA patterns were generated by one of the primers. Consequently, the RAPD technique can be used to discriminate between different cytoplasms. Clear evidence is provided for the organellar origin of fragments in genomic (total DNA) RAPD patterns. The consequences of these results for the interpretation of RAPD analyses are discussed.     

Scanning the genome for gene single nucleotide polymorphisms involved in adaptive population differentiation in white spruce

Conifers are characterized by a large genome size and a rapid decay of linkage disequilibrium, most often within gene limits. Genome scans based on noncoding markers are less likely to detect molecular adaptation linked to genes in these species. In this study, we assessed the effectiveness of a genome-wide single nucleotide polymorphism (SNP) scan focused on expressed genes in detecting local adaptation in a conifer species. Samples were collected from six natural populations of white spruce ( Picea glauca ) moderately differentiated for several quantitative characters. A total of 534 SNPs representing 345 expressed genes were analysed. Genes potentially under natural selection were identified by estimating the differentiation in SNP frequencies among populations ( F _ST) and identifying outliers, and by estimating local differentiation using a Bayesian approach. Both average expected heterozygosity and population differentiation estimates ( H _E = 0.270 and F _ST = 0.006) were comparable to those obtained with other genetic markers. Of all genes, 5.5% were identified as outliers with F _ST at the 95% confidence level, while 14% were identified as candidates for local adaptation with the Bayesian method. There was some overlap between the two gene sets. More than half of the candidate genes for local adaptation were specific to the warmest population, about 20% to the most arid population, and 15% to the coldest and most humid higher altitude population. These adaptive trends were consistent with the genes' putative functions and the divergence in quantitative traits noted among the populations. The results suggest that an approach separating the locus and population effects is useful to identify genes potentially under selection. These candidates are worth exploring in more details at the physiological and ecological levels.     

Genome- and species-specific markers and genome relationships of diploid perennial species in Triticeae based on RAPD analyses.

Eight different genomes (E, H, I, P, R, St, W, and Ns) represented by 22 diploid species of the tribe Triticeae were analyzed using the random amplified polymorphic DNA (RAPD) technique. The genome relationships were obtained based on 371 RAPD fragments produced with 30 primers. The four species of the genus Psathyrostachys (having various Ns genomes) were closely related. The genomes Ee and Eb had a similarly close relationship and were distinct from all other genomes analyzed. Genomes P, R, and St were grouped in one cluster and genomes H and I in another. Genome W had a distant relationship with all other genomes. These results agree with the conclusions from studies of chromosome pairing and isozyme and DNA sequence analyses. Twenty-nine and 11 RAPD fragments are considered to be genome- and species-specific markers, respectively. One to six genome-specific markers were identified for each genome. These RAPD markers are useful in studies of genome evolution, analysis of genome composition, and genome identification.     

Distinct geographical structure across species units evidenced by chloroplast DNA haplotypes and nuclear ribosomal ITS genotypes of Corylopsis (Hamamelidaceae) in the Japanese islands

The geographical distribution of chloroplast DNA (cpDNA) haplotypes and nuclear ribosomal internal transcribed spacer (nrITS) genotypes of Japanese Corylopsis (Hamamelidaceae), which consists of four species, was investigated. Two hundred and five individuals belonging to four species from 30 populations, covering the entire geographical range, were studied. Based on approximately 1108 bp of the three non-coding regions of cpDNA, nine haplotypes were detected, and each was distinguished from adjacent haplotypes by one substitution. Based on approximately 507-bp nrITS sequences, 47 genotypes were detected, for which three clades were identified in the phylogenetic analysis. There was inconsistency between the cpDNA haplotypes, nrITS genotypes, and classification of Corylopsis taxa, possibly because of incomplete lineage sorting or introgressive hybridization. The distribution of the haplotypes was highly structured geographically, and N _ST (0.893) was significantly greater than G _ST (0.819), implying that the current distribution of Corylopsis species was structured phylogeographically during Quaternary climatic oscillations. The haplotype composition and results of analysis of molecular variance showed that the populations in Hokuriku were highly divergent, suggesting that they are long-term persistent populations arising from refugia during the Quaternary climatic oscillations. Refugial populations in Chugoku and Shikoku may have lost genetic diversity because of a bottleneck resulting from a small population size, followed by post-glacial range expansion. Pre-existing refugia may have been so small that the subsequent range expansion replaced the pre-existing genetic structure.  © 2008 The Linnean Society of London, Botanical Journal of the Linnean Society , 2008, 157 , 501–518.