Demographic factors shaped diversity in the two gene pools of wild common bean Phaseolus vulgaris L.

Wild common bean (Phaseolus vulgaris L.) is distributed throughout the Americas from Mexico to northern Argentina. Within this range, the species is divided into two gene pools (Andean and Middle American) along a latitudinal gradient. The diversity of 24 wild common bean genotypes from throughout the geographic range of the species was described by using sequence data from 13 loci. An isolation–migration model was evaluated using a coalescent analysis to estimate multiple demographic parameters. Using a Bayesian approach, Andean and Middle American subpopulations with high percentage of parentages were observed. Over all loci, the Middle American gene pool was more diverse than the Andean gene pool (π_sil=0.0089 vs 0.0068). The two subpopulations were strongly genetically differentiated over all loci (F_st=0.29). It is estimated that the two current wild gene pools diverged from a common ancestor ∼111 000 years ago. Subsequently, each gene pool underwent a bottleneck immediately after divergence and lasted ∼40 000 years. The Middle American bottleneck population size was ∼46% of the ancestral population size, whereas the Andean was 26%. Continuous asymmetric gene flow was detected between the two gene pools with a larger number of migrants entering Middle American gene pool from the Andean gene pool. These results suggest that because of the complex population structure associated with the ancestral divergence, subsequent bottlenecks in each gene pool, gene pool-specific domestication and intense selection within each gene pool by breeders; association mapping would best be practised within each common bean gene pool.     

Phylogeography and genetic diversity of the widespread katydid Ducetia japonica (Thunberg, 1815) across China

Habitat fragmentation can lower migration rates and genetic connectivity among remaining populations of native species. Ducetia japonica is one of the most widespread katydids in China, but little is known about its genetic structure and phylogeographic distribution. We combined the five‐prime region of cytochrome c oxidase subunit I (COI‐5P), 11 newly developed microsatellite loci coupled with an ecological niche model (ENM) to examine the genetic diversity and population structure of D. japonica in China and beyond to Laos and Singapore. Both Bayesian inference (BI) and haplotype network methods revealed six mitochondrial COI‐5P lineages. The distribution of COI‐5P haplotypes may not demonstrate significant phylogeographic structure (N _ST > G _ST, p > .05). The STRUCTURE analysis based on microsatellite data also revealed six genetic clusters, but discordant with those obtained from COI‐5P haplotypes. For both COI‐5P and microsatellite data, Mantel tests revealed a significant positive correlation between geographic and genetic distances in mainland China. Bayesian skyline plot (BSP) analyses indicated that the population size of D. japonica''s three major mitochondrial COI‐5P lineages were seemingly not affected by last glacial maximum (LGM, 0.015–0.025 Mya). The ecological niche models showed that the current distribution of D. japonica was similar to the species’ distribution during the LGM period and only slightly extended in northern China. Further phylogeographic studies based on more extensive sampling are needed to identify specific locations of glacial refugia in northern China.     

AFLP markers suggest low population genetic differentiation of the black rockfish Sebastes schlegelii

AFLP markers were used to investigate the population genetic differentiation among populations of Sebastes schlegelii from China and Japan. A total of 405 bands were amplified for 180 individuals from 9 populations by 3 pairs of selective primers. S. schlegelii exhibited high Nei's gene diversity with an average value of 0.29 for all populations. No significant genealogical branches or clusters corresponding to sampling localities were detected by UPGMA tree. The results of AMOVA analysis and pairwise F_ST values showed no genetic divergence among different geographic populations. The result of STRUCTURE demonstrated all populations of S. schlegelii examined in the study shared one single gene pool.     

Genetic variation in cultivated Rheum tanguticum populations

To examine whether cultivation reduced genetic variation in the important Chinese medicinal plant Rheum tanguticum, the levels and distribution of genetic variation were investigated using ISSR markers. Fifty-eight R. tanguticum individuals from five cultivated populations were studied. Thirteen primers were used and a total of 320 DNA bands were scored. High levels of genetic diversity were detected in cultivated R. tanguticum (PPB = 82.19, H = 0.2498, H_B = 0.3231, I = 0.3812) and could be explained by the outcrossing system, as well as long-lived and human-mediated seed exchanges. Analysis of molecular variance (AMOVA) showed that more genetic variation was found within populations (76.1%) than among them (23.9%). This was supported by the coefficient of gene differentiation (G_st = 0.2742) and Bayesian analysis (θ_B = 0.1963). The Mantel test revealed no significant correlation between genetic and geographic distances among populations (r = 0.1176, p = 0.3686). UPGMA showed that the five cultivated populations were separated into three clusters, which was in good accordance with the results provided by the Bayesian software STRUCTURE (K = 3). A short domestication history and no artificial selection may be an effective way of maintaining and conserving the gene pools of wild R. tanguticum.     

Distribution and Differentiation of Wild,Feral, and Cultivated Populations of Perennial Upland Cotton (Gossypium hirsutum L.) in Mesoamerica and the Caribbean

Perennial forms of Gossypium hirsutum are classified under seven races. Five Mesoamerican races would have been derived from the wild race ‘yucatanense’ from northern Yucatán. ‘Marie-Galante’, the main race in the Caribbean, would have developed from introgression with G. barbadense. The racial status of coastal populations from the Caribbean has not been clearly defined. We combined Ecological Niche Modeling with an analysis of SSR marker diversity, to elucidate the relationships among cultivated, feral and wild populations of perennial cottons. Out of 954 records of occurrence in Mesoamerica and the Caribbean, 630 were classified into four categories cultivated, feral (disturbed and secondary habitats), wild/feral (protected habitats), and truly wild cotton (TWC) populations. The widely distributed three first categories cannot be differentiated on ecological grounds, indicating they mostly belong to the domesticated pool. In contrast, TWC are restricted to the driest and hottest littoral habitats, in northern Yucatán and in the Caribbean (from Venezuela to Florida), as confirmed by their climatic envelope in the factorial analysis. Extrapolating this TWC climatic model to South America and the Pacific Ocean points towards places where other wild representatives of tetraploid Gossypium species have been encountered. The genetic analysis sample comprised 42 TWC accessions from 12 sites and 68 feral accessions from 18 sites; at nine sites, wild and feral accessions were collected in close vicinity. Principal coordinate analysis, neighbor joining, and STRUCTURE consistently showed a primary divergence between TWC and feral cottons, and a secondary divergence separating ‘Marie-Galante’ from all other feral accessions. This strong genetic structure contrasts strikingly with the absence of geographic differentiation. Our results show that TWC populations of Mesoamerica and the Caribbean constitute a homogenous gene pool. Furthermore, the relatively low genetic divergence between the Mesoamerican and Caribbean domesticated pools supports the hypothesis of domestication of G. hirsutum in northern Yucatán.     

Phylogeographic and population genetic analyses reveal Pleistocene isolation followed by high gene flow in a wide ranging,but endangered,freshwater mussel

Freshwater organisms of North America have had their contemporary genetic structure shaped by vicariant events, especially Pleistocene glaciations. Life history traits promoting dispersal and gene flow continue to shape population genetic structure. Cumberlandia monodonta, a widespread but imperiled (IUCN listed as endangered) freshwater mussel, was examined to determine genetic diversity and population genetic structure throughout its range. Mitochondrial DNA sequences and microsatellite loci were used to measure genetic diversity and simulate demographic events during the Pleistocene using approximate Bayesian computation (ABC) to test explicit hypotheses explaining the evolutionary history of current populations. A phylogeny and molecular clock suggested past isolation created two mtDNA lineages during the Pleistocene that are now widespread. Two distinct groups were also detected with microsatellites. ABC simulations indicated the presence of two glacial refugia and post-glacial admixture of them followed by simultaneous dispersal throughout the current range of the species. The Ouachita population is distinct from others and has the lowest genetic diversity, indicating that this is a peripheral population of the species. Gene flow within this species has maintained high levels of genetic diversity in most populations; however, all populations have experienced fragmentation. Extirpation from the center of its range likely has isolated remaining populations due to the geographic distances among them.     

Spatial genetic structure of Aquilegia taxa endemic to the island of Sardinia

Background and Aims

The Mediterranean Basin is one of the most important regions for the Earth''s plant biodiversity; however, the scarcity of studies on fine scale patterns of genetic variation in this region is striking. Here, an assessment is made of the spatial genetic structure of all known locations of the three Sardinian endemic species of Aquilegia in order to determine the relative roles of gene flow and genetic drift as underlying evolutionary forces canalizing the divergence of Sardinian Aquilegia taxa, and to see if the spatial genetic structure found fits the current taxonomic differentiation of these taxa.

Methods

DNA from 89 individuals from all known locations of Aquilegia across Sardinia was analysed by means of amplified fragment length polymorphism (AFLP) markers. Both principal co-ordinates analysis (PCoA) and Bayesian clustering analyses were used to determine the spatial genetic structure irrespective of any taxonomic affiliation. Historical effects of gene flow and genetic drift were assessed by checking for the existence of isolation-by-distance patterns.

Key Results

STRUCTURE and PCoA analyses revealed a pattern of genetic variation geographically structured into four spatial genetic groups. No migration–drift equilibrium was detected for Aquilegia in Sardinia, when analysed either as a whole or in individual groups. The scenario approached a Case III pattern sensu Hutchinson and Templeton, which is associated with extreme isolation conditions where genetic drift has historically played a dominant role over gene flow.

Conclusions

The pattern of genetic variation of Sardinian taxa of Aquilegia indicates that genetic drift has been historically more influential than gene flow on population structure of Sardinian species of Aquilegia. Limited seed dispersal and divergent selection imposed by habitat conditions have been probably the main causes reinforcing post-Pleistocene geographical isolation of Aquilegia populations. The spatial genetic structure found here is not fully compatible with current taxonomic affiliations of Sardinian Aquilegia taxa. This is probably a consequence of the uncoupling between morphological and genetic patterns of differentiation frequently found in recently radiated taxa.     

The ancient tropical rainforest tree Symphonia globulifera L. f. (Clusiaceae) was not restricted to postulated Pleistocene refugia in Atlantic Equatorial Africa

Understanding the history of forests and their species'' demographic responses to past disturbances is important for predicting impacts of future environmental changes. Tropical rainforests of the Guineo-Congolian region in Central Africa are believed to have survived the Pleistocene glacial periods in a few major refugia, essentially centred on mountainous regions close to the Atlantic Ocean. We tested this hypothesis by investigating the phylogeographic structure of a widespread, ancient rainforest tree species, Symphonia globulifera L. f. (Clusiaceae), using plastid DNA sequences (chloroplast DNA [cpDNA], psbA-trnH intergenic spacer) and nuclear microsatellites (simple sequence repeats, SSRs). SSRs identified four gene pools located in Benin, West Cameroon, South Cameroon and Gabon, and São Tomé. This structure was also apparent at cpDNA. Approximate Bayesian Computation detected recent bottlenecks approximately dated to the last glacial maximum in Benin, West Cameroon and São Tomé, and an older bottleneck in South Cameroon and Gabon, suggesting a genetic effect of Pleistocene cycles of forest contraction. CpDNA haplotype distribution indicated wide-ranging long-term persistence of S. globulifera both inside and outside of postulated forest refugia. Pollen flow was four times greater than that of seed in South Cameroon and Gabon, which probably enabled rapid population recovery after bottlenecks. Furthermore, our study suggested ecotypic differentiation—coastal or swamp vs terra firme—in S. globulifera. Comparison with other tree phylogeographic studies in Central Africa highlighted the relevance of species-specific responses to environmental change in forest trees.     

Multiple refugia and barriers explain the phylogeography of the Valais shrew,Sorex antinorii (Mammalia: Soricomorpha)

The aim of the present study was to investigate the genetic structure of the Valais shrew (Sorex antinorii) by a combined phylogeographical and landscape genetic approach, and thereby to infer the locations of glacial refugia and establish the influence of geographical barriers. We sequenced part of the mitochondrial cytochrome b (cyt b) gene of 179 individuals of S. antinorii sampled across the entire species' range. Six specimens attributed to S. arunchi were included in the analysis. The phylogeographical pattern was assessed by Bayesian molecular phylogenetic reconstruction, population genetic analyses, and a species distribution modelling (SDM)‐based hindcasting approach. We also used landscape genetics (including isolation‐by‐resistance) to infer the determinants of current intra‐specific genetic structure. The phylogeographical analysis revealed shallow divergence among haplotypes and no clear substructure within S. antinorii. The starlike structure of the median‐joining network is consistent with population expansion from a single refugium, probably located in the Apennines. Long branches observed on the same network also suggest that another refugium may have existed in the north‐eastern part of Italy. This result is consistent with SDM, which also suggests several habitable areas for S. antinorii in the Italian peninsula during the LGM. Therefore S. antinorii appears to have occupied disconnected glacial refugia in the Italian peninsula, supporting previous data for other species showing multiple refugia within southern refugial areas. By coupling genetic analyses and SDM, we were able to infer how past climatic suitability contributed to genetic divergence of populations. The genetic differentiation shown in the present study does not support the specific status of S. arunchi. © 2012 The Linnean Society of London, Biological Journal of the Linnean Society, 2012, 105 , 864–880.     

Northern glacial refugia for the pygmy shrew Sorex minutus in Europe revealed by phylogeographic analyses and species distribution modelling

The southern European peninsulas (Iberian, Italian and Balkan) are traditionally recognized as glacial refugia from where many species colonized central and northern Europe after the Last Glacial Maximum (LGM). However, evidence that some species had more northerly refugia is accumulating from phylogeographic, palaeontological and palynological studies, and more recently from species distribution modelling (SDM), but further studies are needed to test the idea of northern refugia in Europe. Here, we take a rarely implemented multidisciplinary approach to assess if the pygmy shrew Sorex minutus, a widespread Eurasian mammal species, had northern refugia during the LGM, and if these influenced its postglacial geographic distribution. First, we evaluated the phylogeographic and population expansion patterns using mtDNA sequence data from 123 pygmy shrews. Then, we used SDM to predict present and past (LGM) potential distributions using two different training data sets, two different algorithms (Maxent and GARP) and climate reconstructions for the LGM with two different general circulation models. An LGM distribution in the southern peninsulas was predicted by the SDM approaches, in line with the occurrence of lineages of S. minutus in these areas. The phylogeographic analyses also indicated a widespread and strictly northern‐central European lineage, not derived from southern peninsulas, and with a postglacial population expansion signature. This was consistent with the SDM predictions of suitable LGM conditions for S. minutus occurring across central and eastern Europe, from unglaciated parts of the British Isles to much of the eastern European Plain. Hence, S. minutus likely persisted in parts of central and eastern Europe during the LGM, from where it colonized other northern areas during the late‐glacial and postglacial periods. Our results provide new insights into the glacial and postglacial colonization history of the European mammal fauna, notably supporting glacial refugia further north than traditionally recognized.     

Microclimatic differentiation of gene pools in the Lobaria pulmonaria symbiosis in a primeval forest landscape

Population genetics of the tree‐colonizing lichen Lobaria pulmonaria were studied in the largest primeval beech forest of Europe, covering 10 000 ha. During an intensive survey of the area, we collected 1522 thallus fragments originating from 483 trees, which were genotyped with eight mycobiont‐ and 14 photobiont‐specific microsatellite markers. The mycobiont and photobiont of L. pulmonaria were found to consist of two distinct gene pools, which are co‐existing within small areas of 3–180 ha in a homogeneous beech forest. The small‐scale distribution pattern of the symbiotic gene pools show habitat partitioning of lineages associated with either floodplains or mountain forests. Using approximate Bayesian computation (ABC), we dated the divergence of the two fungal gene pools of L. pulmonaria as the Early Pleistocene. Both fungal gene pools survived the Pleistocene glacial cycles in the Carpathians, although possibly in climatically different refugia. Fungal diversification prior to these cycles and the selection of photobionts with different altitudinal distributions explain the current sympatric, but ecologically differentiated habitat partitioning of L. pulmonaria. In addition, the habitat preferences of the mycobiont are determined by other factors and are rather independent of those of the photobiont at the landscape level. The distinct gene pools should be considered evolutionarily significant units and deserve specific conservation priorities in the future, for example gene pool A, which is a Pliocene relict.     

Use of multiple markers demonstrates a cryptic western refugium and postglacial colonisation routes of Atlantic salmon (Salmo salar L.) in northwest Europe

Glacial and postglacial processes are known to be important determinants of contemporary population structuring for many species. In Europe, refugia in the Italian, Balkan and Iberian peninsulas are believed to be the main sources of species colonising northern Europe after the glacial retreat; however, there is increasing evidence of small, cryptic refugia existing north of these for many cold-tolerant species. This study examined the glacial history of Atlantic salmon in western Europe using two independent classes of molecular markers, microsatellites (nuclear) and mitochondrial DNA variation. Alongside the well-documented refuge in the Iberian Peninsula, evidence for a cryptic refuge in northwest France is also presented. Critically, methods utilised to estimate divergence times between the refugia indicated that salmon in these two regions had diverged a long time before the last glacial maximum; coalescence analysis (as implemented in the program IMa2) estimated divergence times at around 60 000 years before present. Through the examination of haplotype frequencies, previously glaciated areas of northwest Europe, that is, Britain and Ireland, appear to have been colonised from salmon expanding out of both refugia, with the southwest of England being the primary contact zone and exhibiting the highest genetic diversity.     

Genetic diversity and structure of remnant <Emphasis Type="Italic">Magnolia stellata</Emphasis> populations affected by anthropogenic pressures and a conservation strategy for maintaining their current genetic diversity

Understanding the factors shaping rare species’ current genetic diversity and structure, particularly the impact of recent anthropogenic pressures, is important in order to develop appropriate conservation strategies based on robust predictions. Thus, we have genotyped all 585 surviving individuals of Magnolia stellata from six remnant populations and seven isolated tree sites in northern Mie Prefecture, Japan, using nuclear and chloroplast microsatellites. Three genetic clusters were detected by STRUCTURE analysis, with an oldest divergence time between pairs within 25 generations according to coalescent analysis. We attribute this recent divergence to recent anthropogenic environmental changes. Evidence of only one significant recent migration event between pairs of the six populations was detected, indicating that most of the remnant populations are isolated now. Their future genetic status was predicted using Monte Carlo simulations, under four scenarios. It declined more than twice as rapidly in a scenario assuming variations in fecundity among both female and male parents than in a scenario assuming no fecundity variations, but strongly improved in a scenario including promotion of migrations between adjacent pairs of populations. These predictions indicate that sexual reproduction of the species should be promoted by providing more suitable habitats and migration between populations should be enhanced by restoring remnant isolated tree sites and extinct populations. In addition, all the remnant populations should be conserved because they host current genetic variation that may be important for coping with future climate change, and they could provide important stepping-stones for gene flow.     

Exploring Germplasm Diversity to Understand the Domestication Process in Cicer spp. Using SNP and DArT Markers

To estimate genetic diversity within and between 10 interfertile Cicer species (94 genotypes) from the primary, secondary and tertiary gene pool, we analysed 5,257 DArT markers and 651 KASPar SNP markers. Based on successful allele calling in the tertiary gene pool, 2,763 DArT and 624 SNP markers that are polymorphic between genotypes from the gene pools were analyzed further. STRUCTURE analyses were consistent with 3 cultivated populations, representing kabuli, desi and pea-shaped seed types, with substantial admixture among these groups, while two wild populations were observed using DArT markers. AMOVA was used to partition variance among hierarchical sets of landraces and wild species at both the geographical and species level, with 61% of the variation found between species, and 39% within species. Molecular variance among the wild species was high (39%) compared to the variation present in cultivated material (10%). Observed heterozygosity was higher in wild species than the cultivated species for each linkage group. Our results support the Fertile Crescent both as the center of domestication and diversification of chickpea. The collection used in the present study covers all the three regions of historical chickpea cultivation, with the highest diversity in the Fertile Crescent region. Shared alleles between different gene pools suggest the possibility of gene flow among these species or incomplete lineage sorting and could indicate complicated patterns of divergence and fusion of wild chickpea taxa in the past.     

The importance of Anatolian mountains as the cradle of global diversity in Arabis alpina, a key arctic-alpine species

Background and Aims

Anatolia is a biologically diverse, but phylogeographically under-explored region. It is described as either a centre of origin and long-term Pleistocene refugium, or as a centre for genetic amalgamation, fed from distinct neighbouring refugia. These contrasting hypotheses are tested through a global phylogeographic analysis of the arctic–alpine herb, Arabis alpina.

Methods

Herbarium and field collections were used to sample comprehensively the entire global range, with special focus on Anatolia and Levant. Sequence variation in the chloroplast DNA trnL-trnF region was examined in 483 accessions. A haplotype genealogy was constructed and phylogeographic methods, demographic analysis and divergence time estimations were used to identify the centres of diversity and to infer colonization history.

Key Results

Fifty-seven haplotypes were recovered, belonging to three haplogroups with non-overlapping distributions in (1) North America/Europe/northern Africa, (2) the Caucuses/Iranian Plateau/Arabian Peninsula and (3) Ethiopia–eastern Africa. All haplogroups occur within Anatolia, and all intermediate haplotypes linking the three haplogroups are endemic to central Anatolia and Levant, where haplotypic and nucleotide diversities exceeded all other regions. The local pattern of haplotype distribution strongly resembles the global pattern, and the haplotypes began to diverge approx. 2·7 Mya, coinciding with the climate cooling of the early Middle Pleistocene.

Conclusions

The phylogeographic structure of Arabis alpina is consistent with Anatolia being the cradle of origin for global genetic diversification. The highly structured landscape in combination with the Pleistocene climate fluctuations has created a network of mountain refugia and the accumulation of spatially arranged genotypes. This local Pleistocene population history has subsequently left a genetic imprint at the global scale, through four range expansions from the Anatolian diversity centre into Europe, the Near East, Arabia and Africa. Hence this study also illustrates the importance of sampling and scaling effects when translating global from local diversity patterns during phylogeographic analyses.     

Subspecies divergence and pronounced phylogenetic incongruence in the East-Asia-endemic shrub Magnolia sieboldii

Background and AimsThe biogeographic patterns of the East-Asia-endemic shrub Magnolia sieboldii, in which the range of the subsp. sieboldii is interposed with the disjunct distribution of subsp. japonica, implies a complex evolutionary history, involving rapid speciation and hybridization. Here, we aim to reveal the evolutionary and phylogeographic histories of the species with a particular focus on the time of subspecies divergence, the hypothesis of secondary hybridization and the Pleistocene survival of each subspecies, using a combination of genetic analyses and ecological niche modelling.MethodsGenetic variation, genetic structures and phylogenetic relationships were elucidated based on nuclear low-copy genes, chloroplast DNA, and nuclear simple sequence repeats (SSRs). A scenario selection analysis and divergence time estimation were performed using coalescent simulation in DIYABC and *BEAST. Ecological niche modelling and a test of niche differentiation were performed using Maxent and ENMTools.Key ResultsAll marker types showed deep, but pronouncedly incongruent, west–east genetic divergences, with the subspecies being delineated only by the nuclear low-copy genes. Phylogenetic tree topologies suggested that ancient hybridization and introgression were likely to have occurred; however, this scenario did not receive significant support in the DIYABC analysis. The subspecies differentiated their niches, but both showed a dependence on high humidity and were predicted to have persisted during the last glacial cycle by maintaining a stable latitudinal distribution via migration to lower altitudes.ConclusionsWe found a deep genetic divergence and a pronounced phylogenetic incongruence among the two subspecies of M. sieboldii, which may have been driven by major paleogeographic and paleoclimatic events that have occurred since the Neogene in East Asia, including global cooling, climate oscillations and the formation of land bridges. Both subspecies were, however, considered to persist in situ in stable climatic conditions during the late Pleistocene.     

Genetic variation and geographical divergence in Elymus nutans Griseb. (Poaceae: Triticeae) from West China

Elymus nutans Griseb. as an important forage grass and gene pools for improving cereal crops, widely distributes in west China. However, little is known about its genetic and geographical patterns. Inter-simple sequence repeat (ISSR) markers were applied to assess genetic diversity and geographical divergence among 63 E. nutans accessions from west China. The cluster analysis separated the accessions into several groups with geographical origins. The molecular variance analysis (AMOVA) showed that the proportion of variance explained by within- and among-geographical group diversity was 43.2% and 56.8%, respectively. Based on pairwise genetic distances (Φ_ST) between geographical groups, the relationships were congruent with the cluster of accessions. The distinct geographical divergence of E. nutans was revealed between Qinghai-Tibet Plateau and Xinjiang. The ecogeographical conditions such as climate and mountain ranges and elevation behaved as the crucial factors for genetic divergence. Furthermore, the study also indicated that Qinghai-Tibet Plateau might be the diversity differentiation center of E. nutans. The result will facilitate the breeding program and germplasm collection and conservation.     

Interglacial refugia preserved high genetic diversity of the Chinese mole shrew in the mountains of southwest China

The mountains of southwest China (MSC) harbor extremely high species diversity; however, the mechanism behind this diversity is unknown. We investigated to what degree the topography and climate change shaped the genetic diversity and diversification in these mountains, and we also sought to identify the locations of microrefugia areas in these mountains. For these purposes, we sampled extensively to estimate the intraspecific phylogenetic pattern of the Chinese mole shrew (Anourosorex squamipes) in southwest China throughout its range of distribution. Two mitochondrial genes, namely, cytochrome b (CYT B) and NADH dehydrogenase subunit 2 (ND2), from 383 archived specimens from 43 localities were determined for phylogeographic and demographic analyses. We used the continuous-diffusion phylogeographic model, extensive Bayesian skyline plot species distribution modeling (SDM) and approximate Bayesian computation (ABC) to explore the changes in population size and distribution through time of the species. Two phylogenetic clades were identified, and significantly higher genetic diversity was preserved in the southern subregion of the mountains. The results of the SDM, continuous-diffusion phylogeographic model, extensive Bayesian skyline plot and ABC analyses were congruent and supported that the Last Interglacial Maximum (LIG) was an unfavorable period for the mole shrews because of a high degree of seasonality; A. squamipes survived in isolated interglacial refugia mainly located in the southern subregion during the LIG and rapidly expanded during the last glacial period. These results furnished the first evidence for major Pleistocene interglacial refugia and a latitudinal effect in southwest China, and the results shedding light on the higher level of species richness in the southern subregion.     

Evidence for Introduction Bottleneck and Extensive Inter-Gene Pool (Mesoamerica x Andes) Hybridization in the European Common Bean (Phaseolus vulgaris L.) Germplasm

Common bean diversity within and between Mesoamerican and Andean gene pools was compared in 89 landraces from America and 256 landraces from Europe, to elucidate the effects of bottleneck of introduction and selection for adaptation during the expansion of common bean (Phaseolus vulgaris L.) in Europe. Thirteen highly polymorphic nuclear microsatellite markers (nuSSRs) were used to complement chloroplast microsatellite (cpSSRs) and nuclear markers (phaseolin and Pv-shatterproof1) data from previous studies. To verify the extent of the introduction bottleneck, inter-gene pool hybrids were distinguished from “pure” accessions. Hybrids were identified on the basis of recombination of gene pool specific cpSSR, phaseolin and Pv-shatterproof1 markers with a Bayesian assignments based on nuSSRs, and with STRUCTURE admixture analysis. More hybrids were detected than previously, and their frequency was almost four times larger in Europe (40.2%) than in America (12.3%). The genetic bottleneck following the introduction into Europe was not evidenced in the analysis including all the accessions, but it was significant when estimated only with “pure” accessions, and five times larger for Mesoamerican than for Andean germplasm. The extensive inter-gene pool hybridization generated a large amount of genotypic diversity that mitigated the effects of the bottleneck that occurred when common bean was introduced in Europe. The implication for evolution and the advantages for common bean breeding are discussed.