Phylogeographical structure inferred from cpDNA sequence variation of <Emphasis Type="Italic">Zygophyllum xanthoxylon</Emphasis> across north-west China

Zygophyllum xanthoxylon, a desert species, displaying a broad east–west continuous distribution pattern in arid Northwestern China, can be considered as a model species to investigate the biogeographical history of this region. We sequenced two chloroplast DNA spacers (psbK-psbI and rpl32-trnL) in 226 individuals from 31 populations to explore the phylogeographical structure. Median-joining network was constructed and analysis of AMOVA, SMOVA, neutrality tests and distribution analysis were used to examine genetic structure and potential range expansion. Using species distribution modeling, the geographical distribution of Z. xanthoxylon was modeled during the present and at the Last Glacial Maximum (LGM). Among 26 haplotypes, one was widely distributed, but most was restricted to either the eastern or western region. The populations with the highest levels of haplotype diversity were found in the Tianshan Mountains and its surroundings in the west, and the Helan Mountains and Alxa Plateau in the east. AMOVA and SAMOVA showed that over all populations, the species lacks phylogeographical structure, which is speculated to be the result of its specific biology. Neutrality tests and mismatch distribution analysis support past range expansions of the species. Comparing the current distribution to those cold and dry conditions in LGM, Z. xanthoxylon had a shrunken and more fragmented range during LGM. Based on the evidences from phylogeographical patterns, distribution of genetic variability, and paleodistribution modeling, Z. xanthoxylon is speculated most likely to have originated from the east and migrated westward via the Hexi Corridor.     

Phylogeography and historical demography of the orchid bee <Emphasis Type="Italic">Euglossa iopoecila</Emphasis>: signs of vicariant events associated to Quaternary climatic changes

The aim of this study was to investigate whether Pleistocene climatic instability influenced the phylogeographic structure and historical demography of an endemic Atlantic Forest (AF) orchid bee, Euglossa iopoecila Dressler, which shows two main patterns of integument colors over of its geographical distribution. We based our analysis on the concatenated sequence of four mtDNA segments belonging to genes 16S (357 bp), Cytb (651 bp) and COI (1206 bp), totaling 2234 bp. Samples of E. iopoecila populations were collected in 14 AF remnants along its geographic distribution. Median-Joining haplotype networks, SAMOVA and BAPS results indicated three lineages (southern, central and northern clusters) for E. iopoecila, with two important phylogeographic ruptures. We found higher genetic diversity among samples collected in the central region of the AF, which coincides with predicted areas of climatic stability, according to recent AF stability–extinction model. The demographic analysis suggests that only the southern cluster had undergone recent population expansion, which probably started after the last glacial maximum (LGM). Our data suggest that the differentiation observed in the three mitochondrial lineages of E. iopoecila is the result of past disconnections and multiple extinction/recolonization events involving climate fluctuations. In terms of conservation, we would emphasize the importance of considering: (1) the region of the central clade as the location of the highest genetic diversity of mtDNA of E. iopoecila populations; (2) the philopatric behavior of females that tends to restrict mtDNA gene flow in particular, with direct implications for the conservation of the total genetic diversity in euglossine populations.     

Genetic structure of Sakhalin spruce (<Emphasis Type="Italic">Picea glehnii</Emphasis>) in northern Japan and adjacent regions revealed by nuclear microsatellites and mitochondrial gene sequences

The genetic structure of Sakhalin spruce (Picea glehnii) was studied across the natural range of the species, including two small isolated populations in south Sakhalin and Hayachine, by using six microsatellite loci and maternally inherited mitochondrial gene sequences. We also analyzed P. jezoensis, a sympatric spruce in the range. Genetic diversity of P. glehnii was higher in central Hokkaido and the lowest in the Hayachine. Bayesian clustering and principal coordinate analysis by using the microsatellites indicated that the Hayachine was clearly distinct from other populations, implying that it had undergone strong genetic drift since the last glacial period. P. glehnii harbored four mitochondrial haplotypes, two of which were shared with P. jezoensis. One of the two was observed without geographical concentration, suggesting its derivation from ancestral polymorphism. Another was observed in south Sakhalin and in P. jezoensis across Sakhalin. The Bayesian clustering—by using four microsatellite loci, including P. jezoensis populations—indicated unambiguous species delimitation, but with possible admixture of P. jezoensis genes into P. glehnii in south Sakhalin, where P. glehnii is abundantly overwhelmed by P. jezoensis; this might explain the occurrence of introgression of the haplotype of P. jezoensis into P. glehnii.     

Demographic history and population genetic structure of <Emphasis Type="Italic">Hagenia abyssinica</Emphasis> (Rosaceae), a tropical tree endemic to the Ethiopian highlands and eastern African mountains

The distribution of genetic diversity among natural populations is significantly shaped by geographical and environmental heterogeneity. The key objectives of this study were to outline the population genetic structure and to investigate the effects of historical and current factors in shaping the population structure of an endemic tropical tree, Hagenia abyssinica. We used 11 polymorphic microsatellites to estimate genetic variability and evaluate gene flow among natural populations of H. abyssinica. Further, we employed ecological niche modeling approaches, to analyze the demographic history and map potential distributions of H. abyssinica during the Last Glacial Maximum and the present. Significant levels of genetic diversity (H _O = 0.477, H _E = 0.439) were observed among the sampled locations. High coefficient of genetic differentiation (F _ST = 0.32) and considerable genetic variation within the sampled locations (68.01%) were detected. Our results indicated the existence of three genetic groups with limited gene exchange and revealed positive correlations (r = 0.425, P < 0.05) between genetic diversity and geo-graphic distance. The ecological niche modeling (ENM) results support the existence of three distribution zones during the Last Glacial Maximum (LGM), with high probability of occurrence (0.8–1.0), and indicated slight distribution disturbances during and after the LGM. The fundamental patterns of genetic diversity and population structuring of H. abyssinica result from a combination of both environmental and geographical factors, including long-term isolation by distance and characteristic life history of this species. Our ENM results identified three zones that could have served as glacial refugia for this species and lay a foundation for further studies, outlining demographic histories and population structures of Afromontane species.     

Conservation genetics and geographic patterns of genetic variation of endangered shrub <Emphasis Type="Italic">Ammopiptanthus</Emphasis> (Fabaceae) in northwestern China

Phylogeographic patterns of Ammopiptanthus in northwestern China were examined with internal transcribed spacer (ITS) and three chloroplast intergenic spacers (trnH–psbA, trnL–trnF, and trnS–trnG). Two ITS genotypes (a–b) and 8 chloroplast haplotypes (A–H) were detected. Both ITS genotypes and chloroplast lineages were split in two geographic regions: western Xinjiang and the Alxa Desert. This lineage split was also supported by AMOVA analysis and the Mantel test. AMOVA showed that 89.81 % of variance in Ammopiptanthus occurred between the two geographic regions, and correlation between genetic distances and geographical distances was significant (r = 0.757, p < 0.0001). All populations in western Xinjiang shared haplotype A with high frequency, and range expansion was strongly supported by negative Fu’s F_S value, and mismatch distribution analysis, whereas populations in the Alxa Desert had higher genetic diversity and structure. We speculate that the cold and dry climate during the early Quaternary fragmented habitats of the species, limiting gene flow between regions, and interglacial periods most likely led to the range expansion in western Xinjiang. The low genetic diversity of Ammopiptanthus indicate a significant extinction risk, and protective measures should be taken immediately.     

Species radiation in the Alps: multiple range shifts caused diversification in Ringlet butterflies in the European high mountains

The distributions of European high mountain species are often characterised by small and geographically isolated populations and, in many cases, have highly complex biogeographic histories. The butterfly genus Erebia represents one of the best examples for small-scale diversification in the European high mountain systems and therefore to understand speciation processes and associated range dynamics of high mountain species. In this study, we analysed 17 polymorphic allozyme loci of 1731 individuals from 49 populations representing four species, one of which has three subspecies: Erebia nivalis; Erebia tyndarus; Erebia ottomana; and Erebia cassioides cassioides, Erebia cassioides arvernensis, and Erebia cassioides neleus. Samples were collected in the high mountain systems of Europe (i.e. Pyrenees, Massif Central, Alps, Apennines, Carpathians, Balkan high mountains). Genetic analyses supported all previously accepted species. However, the genetic differentiation within E. cassioides sensu lato into three geographically delimited groups is justifying species rank: E. arvernensis distributed in the Pyrenees, Massif Central and western Alps; E. cassioides sensu stricto in the eastern Alps and Apennines; and E. neleus in the Balkan mountains and the south-western Carpathians. While the differentiation between western Alps and Massif Central as well as eastern Alps and Apennines was low, the Pyrenees as well as the south-western Carpathians were significantly differentiated from the other regions within the respective taxon. In general, the differentiation among the populations of E. neleus was stronger than between populations of the other taxa. Within E. cassioides, we found a west-east gradient of genetic similarity over the eastern Alps. Based on the obtained genetic structures, we are able to delineate glacial refugia and interglacial range modifications. Based on the genetic structures and genetic diversity patterns, we conclude that, triggered by the glacial-interglacial cycles, repeated range modifications have taken place with subsequent differentiation and speciation in the region of the Alps and Balkans. Colonisations to Pyrenees (E. arvernensis pseudomurina, E. arvernensis pseudocarmenta), Massif Central (E. ottomana tardenota, E. a. arvernensis) and Apennines (E. cassioides majellana) appear to be recent and most probably not older than the last interglacial period.     

Potential Distribution of Fossil Xenarthrans in South America during the Late Pleistocene: co-Occurrence and Provincialism

Species distribution models (SDMs) are helpful for understanding actual and potential biogeographical traits of organisms. These models have recently started to be applied in the study of fossil xenarthrans. SDMs were generated for 15 South American late Pleistocene xenarthrans: eight Cingulata (Glyptodon clavipes, Doedicurus clavicaudatus, Panochthus tuberculatus, Neosclerocalyptus paskoensis, Pampatherium typum, Pampatherium humboldtii, Holmesina paulacoutoi, and Holmesina occidentalis) and seven Folivora (Glossotherium robustum, Lestodon armatus, Mylodon darwinii, Catonyx cuvieri, Catonyx (=Scelidodon) chilensis, Megatherium americanum, and Eremotherium laurillardi). Models were evaluated for three periods: the last interglacial (LIG), the last glacial maximum (LGM), and the Holocene climatic optimum (HCO). Co-occurrence records were studied based on the overlap of the potential distributions and compared with the available biome reconstructions of South America during the LGM to analyze species distribution patterns, ecological requirements, and possible interactions. Our results suggest the existence of provincialization within xenarthran megamammals grouped in at least three bioregions. Northern and southwestern taxa overlap in the Río de la Plata region where also some endemic taxa are found. We observed overlapping potential distributions but separated and continuous realized distributions between closely related xenarthrans suggesting competitive exclusion. A generalized reduction in potential habitats at the end of the Pleistocene was not obvious as some taxa show stable potential areas during HCO when comparing with LGM. Nonetheless, fragmentation of the most suitable areas due to climate variation and the impact of reduction in available land due to sea level changes cannot be ruled out as involved in the extinction.     

Spatial patterns and intraspecific diversity of the glacial relict legume species <Emphasis Type="Italic">Vavilovia formosa</Emphasis> (Stev.) Fed. in Eurasia

Vavilovia formosa is one of five genera in tribe Fabeae, (Fabaceae, Leguminosae) with close phylogenetic relationships to Pisum. It grows in subalpine and alpine levels in Armenia, Azerbaijan, Georgia, Iran, Iraq, Lebanon, Russia and Turkey and is recognized as an endangered and protected plant. This study was conducted to reveal its intraspecific variability, as well as to predict the past, extant and future species distribution range. We analysed 51 accessions with common phylogenetic markers (trnF-trnL, trnS-trnG, matK, rbcL, psbA-trnH and ITS). These represent in total up to 2551 bp of chloroplast and 664 bp of nuclear sequences per sample. Two populations from Turkey and Armenia were analysed for genetic diversity by AFLP. Leaf morphometry was conducted on 1457 leaflets from 43 specimens. Extracted bioclimatic parameters were used for niche-modelling approach. Analysis of cpDNA revealed two haplotypes, 12 samples from Armenia, Daghestan, Nakhichevan and Iran belonged to H1 group, while 39 samples of all Turkish and part of Armenian were in H2 group. The mean intrapopulation diversity based on AFLP was low (H _E = 0.088) indicating limited outcrossing rate. A significantly positive correlation between geographical latitude and leaf area (\(\rho\) = 0.527, p < 0.05) was found. Niche modelling has shown temporal variation of predicted occurrence across the projected time periods. Vavilovia formosa has suffered a range reduction following climate warming after last glacial maximum, which classify this species as cold-adapted among the Fabeae species as well as a glacial relict.     

Evidence of multiple introductions and genetic admixture of the Asian brush-clawed shore crab <Emphasis Type="Italic">Hemigrapsus takanoi</Emphasis> (Decapoda: Brachyura: Varunidae) along the Northern European coast

The Asian brush-clawed shore crab Hemigrapsus takanoi is a non-indigenous species along the Northern European coast. Although the history of range expansion of European H. takanoi has been well-documented, little is known about the genetic compositions of either the introduced European populations or the native Asian ones. We therefore collected H. takanoi broadly from their native Asian sites and introduced European ranges, and genotyped them by sequencing the mitochondrial 16S RNA gene and by analyzing nuclear microsatellite loci. Our results revealed that the H. takanoi Bay of Seine (France) populations consisted of a genetic admixture between populations in Japan and those in the Yellow Sea region. These French populations should be carefully monitored in the future, since the genetic admixture of multiple source populations may accelerate range expansion in non-indigenous organisms. Our results also suggested that shipping lines from East Asia were more probable vectors than historical juvenile oyster transportations from Japan for the foundation of present European H. takanoi populations. Interestingly, gene flow between populations in Japan and those in the Yellow Sea region (i.e., domestic invasion) was not observed despite the higher potential for artificial translocations via shipping lines in the native Asian range compared with those from Asia to Europe. The lack of domestic invasions implied that intra-specific priority effects of the resident H. takanoi populations played an important role in preventing the successful colonization of artificially-transferred individuals.     

The role of ecological factors in determining phylogeographic and population genetic structure of two sympatric island skinks (<Emphasis Type="Italic">Plestiodon kishinouyei</Emphasis> and <Emphasis Type="Italic">P. stimpsonii</Emphasis>)

We conducted comparative phylogeographic and population genetic analyses of Plestiodon kishinouyei and P. stimpsonii, two sympatric skinks endemic to islands in the southern Ryukyus, to explore different factors that have influenced population structure. Previous phylogenetic studies using partial mitochondrial DNA indicate similar divergence times from their respective closest relatives, suggesting that differences in population structure are driven by intrinsic attributes of either species rather than the common set of extrinsic factors that both presumably have been exposed to throughout their history. In this study, analysis of mtDNA sequences and microsatellite polymorphism demonstrate contrasting patterns of phylogeography and population structure: P. kishinouyei exhibits a lower genetic variability and lower genetic differentiation among islands than P. stimpsonii, consistent with recent population expansion. However, historical demographic analyses indicate that the relatively high genetic uniformity in P. kishinouyei is not attributable to recent expansion. We detected significant isolation-by-distance patterns among P. kishinouyei populations on the land bridge islands, but not among P. stimpsonii populations occurring on those same islands. Our results suggest that P. kishinouyei populations have maintained gene flows across islands until recently, probably via ephemeral Quaternary land bridges. The lower genetic variability in P. kishinouyei may also indicate smaller effective population sizes on average than that of P. stimpsonii. We interpret these differences as a consequence of ecological divergence between the two species, primarily in trophic level and habitat preference.     

Integrating genetics and suitability modelling to bolster climate change adaptation planning in Patagonian <Emphasis Type="Italic">Nothofagus</Emphasis> forests

We investigated the impact of past changes in habitat suitability on the current patterns of genetic diversity of two southern beeches (Nothofagus nervosa and Nothofagus obliqua) in their eastern fragmented range in Patagonian Argentina, and model likely future threats to their population genetic structure. Our goal was to develop a spatially-explicit strategy for guiding conservation and management interventions in light of climate change. We combined suitability modelling under current, past (Last Glacial Maximum ~ 21,000 bp), and future (2050s) climatic conditions with genetic characterization data based on chloroplast DNA, isozymes, and microsatellites. We show the complementary usefulness of the distribution of chloroplast haplotypes and locally common allelic richness calculated from microsatellite data for identifying the locations of putative glacial refugia. Our findings suggest that contemporary hotspots of genetic diversity correspond to convergence zones of different expansion routes, most likely as a consequence of admixture processes. Future suitability predictions suggest that climate change might differentially affect both species. All genetically most diverse populations of N. nervosa and several of N. obliqua are located in areas that may be most severely impacted by climate change, calling for forward-looking conservation interventions. We propose a practical spatially- explicit strategy to target conservation interventions distinguishing priority populations for (1) in situ conservation (hotspots of genetic diversity likely to remain suitable under climate change), (2) ex situ conservation in areas where high genetic diversity overlaps with high likelihood of drastic climate change, (3) vulnerable populations (areas expected to be negatively affected by climate change), and (4) potential expansion areas under climate change.     

Conservation genetics of redside dace (<Emphasis Type="Italic">Clinostomus elongatus</Emphasis>): phylogeography and contemporary spatial structure

Redside dace Clinostomus elongatus (Teleostei: Cyprinidae) is a species of conservation concern that is declining throughout its range as a result of habitat fragmentation, degradation and loss. We characterized the genetic structure and diversity of redside dace populations across the species range using mitochondrial and microsatellite data to inform conservation efforts and assess how historical and recent events have shaped genetic structure and diversity within and among populations. Phylogeographic structure among 28 redside dace populations throughout southern Ontario (Canada) and the United States was assessed by sequence analysis of the mitochondrial cytochrome b and ATPase 6 and 8 genes. Populations were also genotyped using ten microsatellite loci to examine genetic diversity within and among populations as well as contemporary spatial structuring. Mitochondrial DNA sequence data revealed three geographically distinct lineages, which were highly concordant with groupings identified by microsatellite analysis. The combined genetic data refute published glacial refugia hypotheses of a single Mississippian refugium or of two lineages associated with Mississippian and Atlantic refugia. Secondary contact between the two eastern groups was documented in the Allegheny River drainage and tributaries to Lake Ontario. With the exception of several allopatric populations within the Allegheny watershed, high genetic structuring among populations suggests their isolation, indicating that recovery efforts should be population-based.     

Population genetic analysis of Euro-Arctic polar cod <Emphasis Type="Italic">Boreogadus saida</Emphasis> suggests fjord and oceanic structuring

Polar cod, Boreogadus saida, is a key species in Arctic marine ecosystems; however, its genetic population structure is largely undescribed. The population genetic structure of 472 B. saida specimens among nine locations in the north-east Atlantic was revealed using 12 microsatellite loci. Pairwise F _ST comparisons showed significant population differentiation between B. saida sampled inside fjords in Svalbard and north-east Greenland, as compared to B. saida from the shelf. The observed genetic variation was not a function of isolation by distance, and it is speculated that B. saida populations inhabiting fjords may have become reproductively isolated from shelf-dwelling B. saida during the last post-glacial recolonization.     

Nuclear microsatellite markers reveal the low genetic structure of <Emphasis Type="Italic">Pinus mugo</Emphasis> Turra (dwarf mountain pine) populations in Europe

Twenty-one populations (555 individuals) covering the entire native range of Pinus mugo Turra (dwarf mountain pine) were investigated for genetic variation scored at 13 nuclear microsatellite markers (nSSRs). The main objective of the present study was to determine the genetic structure across the present distribution of the species and locate populations of different genetic compositions. Most of the genetic variation was observed within the populations (95%). The assignment of populations based on Bayesian clustering methods revealed that the Sudeten populations of P. mugo form a separate genetic cluster. These stands have likely been established through the founder effects of Alpine migrants. The distribution and level of SSR polymorphisms, along with no evidence of isolation by distance or phylogeographic structure, indicate that the present populations of P. mugo have diverged relatively recently and originate from a larger glacial distribution of the species. One peripheral stand from Italy had the lowest values of most calculated genetic variation indices. This stand could therefore be more susceptible to genetic drift and a negative impact of predicted environmental changes. We discuss our findings with respect to previously published results on the genetic and morphological variation of P. mugo and with consideration for the conservation genetics of the species.     

Genetic structure of the populations of <Emphasis Type="Italic">Dactylorhiza ochroleuca</Emphasis> and <Emphasis Type="Italic">D. incarnata</Emphasis> (Orchidaceae) in the area of their joint growth in Russia and Belarus

We carried out an allozyme analysis to investigate polymorphism and genetic structure of the populations of D. incarnata and D. ochroleuca in regions of their joint growth in Russia and Belarus. We found that D. ochroleuca individuals in the populations of the Urals and Siberia, which are distant fragments from the main range of the species, do not differ significantly from individuals within the main part of the area (Belarus) on the basis of the allelic composition of eight gene loci. We revealed that D. ochroleuca and D. incarnata are differentiated by different alleles of the GDH locus. Thus, we established a genetic marker suitable to distinguish these closely related taxa. In addition to the GDH locus, D. ochroleuca and D. incarnata in the places of their joint growth, differ in the allelic structure of the PGI and NADHD loci. D. incarnata from the Urals and Siberia were polymorphic for both loci, and individuals from Belarus were polymorphic for one locus (PGI). In contrast, all D. ochroleuca individuals growing in sympatric populations with polymorphic D. incarnata were homozygous for the same alleles. Thus, comparison of the genetic structure of D. ochroleuca and D. incarnata points to the existence of a genetic isolation and a functioning isolation mechanism even under conditions of their joint growth. We found that the GDH locus in D. incarnata is polymorphic only in populations which grow together with D. ochroleuca, with exception a few examples. Thus, we conclude that variability of the GDH locus in D. incarnata is associated with hybridization with D. ochroleuca.     

Molecular phylogeography and paleodistribution modeling of the boreal tree species <Emphasis Type="Italic">Ulmus lamellosa</Emphasis> (T.Wang et S. L. Chang) (Ulmaceae) in China

The uplift of mountains and climatic oscillations are important for understanding of the demographic history and genetic structure of species. We investigated the biogeographic history of the boreal tree species Ulmus lamellosa (Ulmaceae) in China, by using a combined phylogeographic and paleodistribution modeling approach. In this study, 14 populations of endangered U. lamellosa were analyzed by using chloroplast DNA (cpDNA) sequences. A high level of genetic differentiation (Φ _ST = 86.22%) among populations with a significant phylogeographic pattern (N _ST > G _ST, P < 0.05) was found in U. lamellosa. Ten haplotypes were detected by combining chloroplast DNA data, and haplotype 3 (H3) was found to be common and widespread. The intraspecific divergence of all U. lamellosa cpDNA haplotypes (9.27 Ma; 95% HPD 5.17–13.33 Ma) most probably began in the late Miocene. The pairwise difference among haplotypes and neutrality tests (Tajima’s D and Fu’s Fs statistic) indicated that populations of U. lamellosa, except group I, have not experienced recent sudden expansions. Multiple refuge areas were identified across the entire distribution ranges of U. lamellosa. The low level of gene flow (Nm = 0.14) among populations may have resulted from isolation resulting from distance and complex topography during climatic oscillations; this isolation was probably the major process that shaped the present distribution of haplotypes. These results support the hypothesis that U. lamellosa persisted in situ during glaciations and occupied multiple localized glacial refugia, contrary to the hypotheses of large-scale range contraction and long-distance southward migration.     

Phylogeography of endangered <Emphasis Type="Italic">Dendrobium moniliforme</Emphasis> in East Asia based on mitochondrial DNA sequence variations

Belonging to the genus Dendrobium of Orchidaceae, Dendrobium moniliforme is an endangered species with disjunct distribution in East Asia, possessing significant medicinal value. To investigate its phylogeography, this study compared sequences of two mitochondrial DNA (mtDNA) fragments (nad 1/b-c and nad 7/2-3) from 143 samples of 18 natural populations of D. moniliforme almost covering the entire range of the Sino-Japanese Floristic Region (SJFR) of East Asia. As a result, a total of 30 mtDNA haplotypes were identified in these populations which revealed high levels of haplotype diversity (H _d  = 0.8733) and total genetic diversity (H _T  = 0.8886). Additionally, G _ST value being significantly lower (0.451) than N _ST value (0.722) (P < 0.05) indicated the presence of strong phylogeographic structure in these populations, and the network of mtDNA haplotypes showed that all haplotypes were divided into two clades (A and B). Haplotype overlap was observed among several D. moniliforme population groups in mainland China, suggesting the occurrence of ongoing and/or historical gene flow among them. No common haplotypes were shared by D. moniliforme populations from mainland China and the CKJ Islands (representing Taiwan [China], South Korea, and Japan collectively), pointing to their allopatric evolution in the two regions. Moreover, mismatch distribution analysis and neutral test of mtDNA data rejected the population expansion model. According to the mtDNA-based results, we infer that multiple refugia for D. moniliforme existed in the SJFR of East Asia during the Quaternary glacial period.     

Population genetic diversity and geographical differentiation of MHC class II DAB genes in the vulnerable Chinese egret (<Emphasis Type="Italic">Egretta eulophotes</Emphasis>)

Major histocompatibility complex (MHC) genes are excellent markers for the study of adaptive genetic variation occurring over different geographical scales. The Chinese egret (Egretta eulophotes) is a vulnerable ardeid species with an estimated global population of 2600–3400 individuals. In this study, we sampled 172 individuals of this egret (approximately 6 % of the global population) from five natural populations that span the entire distribution range of this species in China. We examined their population genetic diversity and geographical differentiation at three MHC class II DAB genes by identifying eight exon 2 alleles at Egeu-DAB1, eight at Egeu-DAB2 and four at Egeu-DAB3. Allelic distributions at each of these three Egeu-DAB loci varied substantially within the five populations, while levels of genetic diversity varied slightly among the populations. Analysis of molecular variance showed low but significant genetic differentiation among five populations at all three Egeu-DAB loci (haplotype-based ?_ST: 0.029, 0.020 and 0.042; and distance-based ?_ST: 0.036, 0.027 and 0.043, respectively; all P < 0.01). The Mantel test suggested that this significant population genetic differentiation was likely due to an isolation-by-distance pattern of MHC evolution. However, the phylogenetic analyses and the Bayesian clustering analysis based on the three Egeu-DAB loci indicated that there was little geographical structuring of the genetic differentiation among five populations. These results provide fundamental population information for the conservation genetics of the vulnerable Chinese egret.     

Genetic differentiation among natural populations of the lizard complex <Emphasis Type="Italic">Darevskia raddei</Emphasis> as inferred from genome microsatellite marking

The article presents the genetic parameters of the populations of lizards of the Darevskia raddei complex (D. raddei nairensis and D. raddei raddei) and the populations of D. valentini calculated on the basis of the analysis of variability of 50 allelic variants of the three nuclear genome microsatellite-containing loci of 83 individuals. It was demonstrated that the F_st genetic distances between the populations of D. raddei nairensis and D. raddei raddei were not statistically significantly different from the F_st genetic distances between the populations of different species, D. raddei and D. valentini. At the same time, these distances were statistically significantly higher than the F_st distances between the populations belonging to one species within the genus Darevskia. These data suggest deep divergence between the populations of D. raddei raddei and D. raddei nairensis of the D. raddei complex and there arises the question on considering them as separate species.     

SSR-based identification of genetic groups within European populations of <Emphasis Type="Italic">Tuber aestivum</Emphasis> Vittad

Tuber species are ectomycorrhizal ascomycetes establishing relationships with different host trees and forming hypogeous fruiting bodies known as truffles. Among Tuber species, Tuber aestivum Vittad. has a wide distributional range being found naturally all over Europe. Here, we performed large-scale population genetic analyses in T. aestivum to (i) investigate its genetic diversity at the European scale, (ii) characterize its genetic structure and test for the presence of ecotypes and (iii) shed light into its demographic history. To reach these goals, 230 ascocarps from different populations were genotyped using 15 polymorphic simple sequence repeat markers. We identified 181 multilocus genotypes and four genetic groups which did not show a clear geographical separation; although, one of them was present exclusively in Southeast France, Italy and Spain. Fixation index values between pairs of genetic groups were generally high and ranged from 0.29 to 0.45. A significant deficit of heterozygosity indicated a population expansion instead of a recent population bottleneck, suggesting that T. aestivum is not endangered in Europe, not even in Mediterranean regions. Our study based on a large-scale population genetic analysis suggests that genetically distinct populations and likely ecotypes within T. aestivum are present. In turn, this study paves the way to future investigations aimed at addressing the biological and/or ecological factors that have concurred in shaping the population genetic structure of this species. Present results should also have implications for the truffle market since defining genetic markers are now possible at least for some specific T. aestivum genetic groups.