Morphological diversity and phylogeography of the Georgian durmast oak (<Emphasis Type="Italic">Q. petraea</Emphasis> subsp. <Emphasis Type="Italic">iberica</Emphasis>) and related Caucasian oak species in Georgia (South Caucasus)

The Caucasus region is one of the 25 global biodiversity hotspots and constitutes a shelter area for Neogene relict species as well as a center of ongoing radiation. In order to elucidate the taxonomic identity, divergence patterns, and evolutionary history of the largely widespread Georgian durmast oak (Quercus petraea subsp. iberica), we examined leaf morphology and chloroplast DNA (cpDNA) (trnH-psbA, trnK-matK) sequence variation across its South Caucasian range. Six other oak taxa distributed throughout Georgia were included in the dataset and used for comparison. Evidence for differentiation in both sets of traits was found. Populations represented by different taxa from ecologically equivalent areas showed common morphological features and genetic structures. Molecular analysis clearly indicated the presence of two major haplotype lineages (West Caucasian vs. East Caucasian zonation type) and suggested a maternal lineage diversification of Q. petraea subsp. iberica in the Late Miocene, as a likely result of complex patterns associated with major orogenic and climatic changes. The Quaternary glacial oscillations resulted in a number of less common, derived haplotypes. Based on mismatch distribution analysis and neutrality tests, we found no evidence of demographic expansion for the populations from the West and East Caucasian zonation types. The two Caucasian provinces therefore acted as important shelter/diversification areas and as a lineage crossroad for the Georgian oaks. Close intra- and interspecific cpDNA relationships shared with other oaks from bordering countries support the relevant role played by the Colchis region as a primary refugium for the European temperate forest species.     

Genetic structure of Schrenck newt <Emphasis Type="Italic">Salamandrella schrenckii</Emphasis> populations by mitochondrial cytochrome <Emphasis Type="Italic">b</Emphasis> variation

The nucleotide sequence variation of the mitochondrial cytochrome b gene was studied in Schrenck newt Salamandrella schrenckii (Strauch, 1870) from populations of Primorye and the Khabarovsk region. Phylogenetic analysis revealed two haplotype clusters, southern cluster 1 and northern cluster 2, with a divergence of 3%. Analysis of the mtDNA and cytochrome b amino acid sequence variations made it possible to assume that the modern range of Schrenck newt was colonized from south Primorye northwards. In contrast to the southern cluster, the northern one demonstrated all the signs of demographic expansion (a unimodal distribution of pairwise nucleotide differences, specific results of tests for selective neutrality of mtDNA variation, and a good correspondence of genetic parameters to those expected from demographic expansion models).     

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 of the Asian lesser white-toothed shrew, <Emphasis Type="Italic">Crocidura shantungensis</Emphasis>, in East Asia: role of the Korean Peninsula as refugium for small mammals

Many peninsulas in the temperate zone played an important role as refugia of various flora and fauna, and the southern Korean Peninsula also served as a refugium for many small mammals in East Asia during the Pleistocene. The Asian lesser white-toothed shrew, Crocidura shantungensis, is a widely distributed species in East Asia, and is an appropriate model organism for exploring the role of the Korean Peninsula as a refugium of small mammals. Here, we investigated phylogenetic relationships and genetic diversity based on the entire sequence of the mitochondrial cytochrome b gene (1140 bp). A Bayesian tree for 98 haplotypes detected in 228 C. shantungensis specimens from East Asia revealed the presence of three major groups with at least 5 subgroups. Most haplotypes were distributed according to their geographic proximity. Pairwise F_ST’s and analysis of molecular variance (AMOVA) revealed a high degree of genetic differentiation and variance among regions as well as among populations within region, implying little gene flow among local populations. Genetic evidence from South Korean islands, Jeju-do Island of South Korea, and Taiwan leads us to reject the hypothesis of recent population expansion. We observed unique island-type genetic characteristics consistent with geographic isolation and resultant genetic drift. Phylogeographic inference, together with estimates of genetic differentiation and diversity, suggest that the southern most part the Korean Peninsula, including offshore islands, played an important role as a refugium for C. shantungensis during the Pleistocene. However, the presence of several refugia on the mainland of northeast Asia is also proposed.     

Taxonomic implications of the morphological and genetic variation of cultivated and domesticated populations of the <Emphasis Type="Italic">Agave angustifolia</Emphasis> complex (Agavoideae,Asparagaceae) in Oaxaca,Mexico

The Agave angustifolia complex, distributed from Mexico to Costa Rica, comprises four species and five varieties, including three species used for mescal production. The complex is represented in the Mexican state of Oaxaca by two wild taxa, A. angustifolia var. angustifolia and A. angustifolia var. rubescens, the cultivated form A. angustifolia “Espadín” and the partially cultivated species A. rodacantha. The aims of this study were to investigate the morphological and genetic variation of the A. angustifolia complex in the state of Oaxaca and to identify traits useful for taxonomic delimitation. Four wild and three cultivated populations of A. angustifolia from Oaxaca, one population of A. tequilana from Guanajuato and one population of A. angustifolia from Sonora were sampled for morphological, genetic and cytometric analyses. We showed that cultivated populations of A. angustifolia “Espadin,” A. rhodacantha and A. tequilana could be clearly differentiated from wild populations. Furthermore, the domesticated populations of A. angustifolia, known locally as “Espadin,” had a higher ploidy level and lower genetic variation than their related wild populations. The population of A. angustifolia from Sonora could be recognized as a different entity. Populations of A. rhodacantha need to be studied throughout their entire distribution area to further evaluate their taxonomic delimitation.     

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

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

Differential Rates of Male Genital Evolution in Sibling Species of <Emphasis Type="Italic">Drosophila</Emphasis>

Genital morphology in animals with internal fertilization is considered to be among the fastest evolving traits. Sexual selection is often proposed as the main driver of genital diversification but the exact selection mechanisms involved are usually unclear. In addition, the mechanisms operating may differ even between pairs of sibling species. We investigated patterns of male genital variation within and between natural populations of the cactophilic fly Drosophila koepferae ranging its entire geographic distribution and compared them with those previously observed in its sibling species, D. buzzatii. Using both mtDNA and nDNA markers we found that genital shape variation in D. koepferae is more restricted than expected for neutral evolution, suggesting the predominance of stabilizing selection. We also detected dissimilar patterns of divergence between populations of D. koepferae that were allopatric and sympatric with D. buzzatii. The constrained evolution inferred for D. koepferae’s genitalia clearly contrasts with the rapid divergence and higher morphological disparity observed in the populations of D. buzzatii. Finally, different possible scenarios of male genital evolution in each species and within the radiation of D. buzzatii cluster are discussed.     

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

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

Low Genetic Diversity and Low Gene Flow Corresponded to a Weak Genetic Structure of Ruddy-Breasted Crake (<Emphasis Type="Italic">Porzana fusca</Emphasis>) in China

The Ruddy-breasted Crake (Porzana fusca) is an extremely poorly known species. Although it is not listed as globally endangered, in recent years, with the interference of climate change and human activities, its habitat is rapidly disappearing and its populations have been shrinking. There are two different life history traits for Ruddy-breasted Crake in China, i.e., non-migratory population in the south and migratory population in the north of China. In this study, mitochondrial control sequences and microsatellite datasets of 88 individuals sampled from 8 sites were applied to analyze their genetic diversity, genetic differentiation, and genetic structure. Our results indicated that low genetic diversity and genetic differentiation exit in most populations. The neutrality test suggested significantly negative Fu’s Fs value, which, in combination with detection of the mismatch distribution, indicated that population expansion occurred in the interglacier approximately 98,000 years ago, and the time of the most recent common ancestor (TMRCA) was estimated to about 202,705 years ago. Gene flow analysis implied that the gene flow was low, but gene exchange was frequent among adjacent populations. Both phylogenetic and STRUCTURE analyses implied weak genetic structure. In general, the genetic diversity, gene flow, and genetic structure of Ruddy-breasted Crake were low.     

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.     

Intraspecific Variation in Mitogenomes of Five <Emphasis Type="Italic">Crassostrea</Emphasis> Species Provides Insight into Oyster Diversification and Speciation

A large number of Crassostrea oysters are found in Asia-Pacific. While analyses of interspecific variation have helped to establish historical relationships among these species, studies on intraspecific variation are necessary to understand their recent evolutionary history and current forces driving population biology. We resequenced 18 and analyzed 31 mitogenomes of five Crassostrea species from China: Crassostrea gigas, Crassostrea angulata, Crassostrea sikamea, Crassostrea ariakensis, and Crassostrea hongkongensis. Our analysis finds abundant insertions, deletions, and single-nucleotide polymorphisms in all species. Intraspecific variation varies greatly among species with polymorphic sites ranging from 54 to 293 and nucleotide diversity ranging from 0.00106 to 0.00683. In all measurements, C. hongkongensis that has the narrowest geographic distribution exhibits the least sequence diversity; C. ariakensis that has the widest distribution shows the highest diversity, and species with intermediate distribution show intermediate levels of diversity. Low sequence diversity in C. hongkongensis may reflect recent bottlenecks that are probably exacerbated by human transplantation. High diversity in C. ariakensis is likely due to divergence of northern and southern China populations that have been separated without gene flow. The significant differences in mitogenome diversity suggest that the five sister species of Crassostrea have experienced different evolutionary forces since their divergence. The recent divergence of two C. ariakensis populations and the C. gigas/angulata species complex provides evidence for continued diversification and speciation of Crassostrea species along China’s coast, which are shaped by unknown mechanisms in a north–south divide.     

Endangered species in small habitat patches can possess high genetic diversity: the case of the Tana River red colobus and mangabey

We used mtDNA sequence data from the Tana River red colobus and mangabey to determine how their population genetic structure was influenced by dispersal and habitat fragmentation. The colobus and mangabey are critically endangered primates endemic to gallery forests in eastern Kenya. The forests are a Pliocene–Pleistocene refugium that has recently undergone significant habitat loss and fragmentation due to human activities. We expected both primates to exhibit low levels of genetic diversity due to elevated genetic drift in their small populations, and to show a strong correspondence between genetic and geographic distance due to disruption of gene flow between forests by habitat fragmentation. Additionally, because mangabey females are philopatric, we expected their mtDNA variation to be homogeneous within forest patches but to be heterogeneous between patches. In contrast, colobus have a female-biased dispersal and so we expected their mtDNA variation to be homogeneous within and between forest patches. We found high levels of haplotype and nucleotide diversity as well as high levels of sequence divergence between haplotype groups in both species. The red colobus had significantly higher genetic variation than the mangabey did. Most of the genetic variation in both primates was found within forest fragments. Although both species showed strong inter-forest patch genetic structure we found no correspondence between genetic and geographic distances for the two primates. We attributed the high genetic diversity to recent high effective population size, and high sequence divergence and strong genetic structures to long-term habitat changes in the landscape.     

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.     

The Role of Geographical and Ecological Factors on Population Divergence of the Neotropical otter <Emphasis Type="Italic">Lontra longicaudis</Emphasis> (Carnivora,Mustelidae)

The geographic distribution of the populations of a species are influenced by the spatial structure of the ecosystems, the environmental factors and the presence of geographic barriers. The Neotropical otter, Lontra longicaudis, is widely distributed throughout the Americas, where a wide range of environmental conditions and geographical features could promote genetic and morphological variation on the three currently recognized subspecies. In this study, we combined phylogeographic, morphometric and environmental niche modelling analyses to examine whether: (1) genetic variation is associated with the presence of barriers to gene flow and/or hydrography; (2) genetic and morphologic variation are associated with environmental variation; and (3) the observed variation in L. longicaudis populations corresponds to the previously defined subspecies. We found strong phylogeographic structure between the northern (L. l. annectens) and the two-southern subspecies (L. l. longicaudis and L. l. enudris), and although shallower, we also detected genetic differentiation between the two South American subspecies. Such genetic differentiation corresponds to the hydrography and to the geographical barriers characteristic of the distributional area of the species. We found a correlation between the shape of the skull and mandible with the environmental variation through the distribution of the species, and we rejected the hypothesis of niche equivalency and similarity between the three identified genetic lineages, suggesting adaptations to different environmental conditions. Our results support that the variation in environmental conditions, in concert with geographical barriers to gene flow and hydrography, have led to population divergence of L. longicaudis along the Neotropics. These results have important taxonomic implications for the species and its conservation.     

Genetic structure analysis of mantis shrimp <Emphasis Type="Italic">Oratosquilla oratoria</Emphasis> based on mitochondrial DNA control region sequence

The genetic relationships of mantis shrimp Oratosquilla oratoria between the coastal waters of China and Japan were not well studied. To reveal the genetic differentiation and genetic structure among populations, we collected populations of mantis shrimp O. oratoria from the coastal waters of China and Japan to analyze the mtDNA control region variation. A total of 309 individuals of O. oratoria were collected from 13 localities (11 from China and 2 from Japan) and a segment of mitochondrial DNA control region was sequenced. Three hundred nine haplotypes were defined, yielding a very high haplotype diversity and nucleotide diversity. Two lineages of O. oratoria were revealed and displayed strong differences in the geographical distribution. In the coastal waters of China, the geographic distribution of the two lineages was completely separated by the Yangtze River estuary; however, the lineages showed geographic sympatry in two populations from Japan. Based on the lineage distribution, three groups were defined. There was no significant genetic differentiation among the populations within the three groups, indicating high gene flow within each group. Significant and negative values for Tajima D and Fu’s Fs tests, and mismatch distributions for two lineages indicated population expansion. The present result confirmed that the freshwater outflow from the Yangtze River formed a physical barrier and affected gene exchange. The different distribution patterns of the two lineages in coastal waters of China and Japan indicated that the larvae of O. oratoria were transferred from China to the coastal waters of Japan with a one-way gene flow.     

Genetic diversity of endangered orchid <Emphasis Type="Italic">Phaius australis</Emphasis> across a fragmented Australian landscape

Historical events such as colonisation, spatial distribution across different habitats, and contemporary processes, such as human-mediated habitat fragmentation can leave lasting imprints on the population genetics of a species. Orchids currently comprise 17% of threatened flora species in Australia (Environment Protection and Biodiversity Conservation Act 1999) due to the combination of fragmentation and illegal harvesting (Benwell in Recovery plan, swamp orchids Phaius australis, Phaius tancarvilliae, NSW National Parks and Wildlife Service, Sydney, 1994; Jones in A complete guide to native orchids of Australia including the island territories, 2nd edn, Reed Natural History, Sydney, 2006; DE in Phaius australis in species profile and threats database, Department of the Environment. http://www.environment.gov.au/sprat, 2015). The federally endangered Swamp Orchid Phaius australis has a disjunct distribution across an almost 2000 km latitudinal range along Australia’s east coast but it was estimated that 95% of the populations have been lost since European settlement (Benwell 1994). Phaius australis is endangered due to illegal collection and habitat loss that has resulted in limited connectivity between populations, in ecosystems that are vulnerable to climate change. Thus the genetic impacts of its history combined with more recent fragmentation may have impacts on its future viability especially in light of changing environmental conditions. Thirty-four populations were sampled from tropical north Queensland to the southern edge of the subtropics in New South Wales. Population genetics analysis was conducted using 13 polymorphic microsatellite markers developed for the species using NextGen sequencing. Spatial genetic patterns indicate post-colonisation divergence from the tropics southwards to its current climate niche limits. Genetic diversity is low across all populations (A?=?1.5, H _e  = 0.171), and there is little evidence of genetic differentiation between regions. Consistent with population genetic theory, the historic loss of populations has resulted in significantly lower genetic diversity in small populations compared to large (P, A, He; p?<?0.05). The viability and persistence of P. australis populations now and in a changing climate are discussed in the context of conservation priorities.     

黄河上游极边扁咽齿鱼种群遗传结构

极边扁咽齿鱼Platypharodon extremus是黄河上游的特有鱼类,近年来由于过度捕捞、环境变迁等因素,其资源量剧减,种群处于濒危状态。研究中采集了分布于黄河上游的3个种群（n=107）,基于线粒体DNA控制区695bp序列,共检测到101个可变位点,占总分析位点的14.5%,其中58个为简约信息位点。3种群共界定了87个单倍型,种群平均单倍型多样性h=0.995,核苷酸多样性π=0.0129。结果表明,极边扁咽齿鱼的种群遗传多样性水平较高。分子方差分析（AMOVA）表明,3群体间总遗传分化系数Fst=0.0528,群体间尚未有显著的遗传分化。系统发生树显示3种群的单倍型混合分布,各进化枝之间分歧度很低,没有形成明显的单倍型组,也没有显示出单倍型与地理位置的对应关系。3个群体共享1个单倍型（Hap₃1）,推测它们来自于共同的祖先。歧点分布和Fu’Fs中性检测显示极边扁咽齿鱼并未经历种群扩张。     

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.     

Population structure and gene flow among wild populations of the black muntjac (Muntiacus crinifrons) based on mitochondrial DNA control region sequences

The black muntjac (Muntiacus crinifrons) is a rare species endemic to China and its current distribution is confined to partial mountain ranges in eastern China. To assess the population structure and gene flow among the extant populations, we sequenced a portion of the mitochondrial control region (424 bp) for 47 samples collected from its current three large populations (Huangshan, Tiammushan and Suichang). A total of 18 unique haplotypes were defined based on 22 polymorphic sites. Of these haplotypes, 15 ones were population-specific and only one haplotype was shared among the three populations. Significant genetic differentiation was detected between Suichang and Huangshan populations (?(ST)=0.1677, P<0.001) or between Suichang and Tianmushan populations (?(ST)=0.2002, P<0.001), indicating that the Suichang population may be spatially structured from other two populations along matriline. The Mantel test revealed that this significant differentiation was not driven by geographic distance (P=0.14), implying that genetic divergence of Suichang population might result from recent human disturbances. Phylogenetic analyses suggested the mitochondrial control region haplotypes were split into two well divergent clades (Clade I and Clade II). Interestingly, the two distinct haplotype clades were found to coexist in Suichang area. The nested clade analysis revealed a significant phylogeographic structure among the black muntjac populations (total cladogram: chi2=18.68; P<0.001), which was inferred to result from past fragmentation followed by range expansion. The population expansion was supported by the analysis of mismatch distribution and the tests of neutrality. Therefore, we suggest that the coexistence of distinct haplotypes in Suichang population was induced by historical population expansion after fragmentation and that the current genetic differentiation should be attributed to the reduction of female-mediated gene flow due to recent habitat fragmentation and subsequent loss.