Chloroplast phylogeography of a desert shrub,Calligonum calliphysa (Calligonum,Polygonaceae), in arid Northwest China

In order to investigate the influence of Pleistocene climate oscillations and paleogeographic events on the evolutionary history of xerophytic plants in arid Northwest China, the phylogeography of a well adapted desert shrub, Calligonum calliphysa (Calligonum, Polygonaceae), was studied. A total of 11 natural populations and 90 individuals were sampled. Two cpDNA intergenic spacer regions (rpl32-trnL and ycf6-psbM) were sequenced and 13 haplotypes were identified. Significant genetic differentiation was found among populations and groups, suggesting that short distance seed transfer and geographic isolation have restricted gene flow. Based on SAMOVA analysis, the 11 populations were divided into 3 clades with no haplotypes were shared between these. The time of divergence within C. calliphysa was estimated at between 2.36 and 0.18 Ma, during the early to middle Pleistocene, and for Clade 3, in the range of 0.88–0.29 Ma. The largest number of populations (7) and haplotypes (9) were found in Clade 3. We found fragmentation of genetic variation, with most unique haplotypes among populations located at the edges of the Gurbantunggut Desert, which may have occurred because vegetation there was sensitive to habitat variation caused by climate change, and we detected that a demographic expansion event followed the expansion of this desert.     

Genetic diversity and population differentiation of Capparis spinosa (Capparaceae) in Northwestern China

To explore the influence of vicariance on differentiation patterns of taxa in arid regions, we systematically assessed the genetic diversity and variability of Capparis spinosa, a typical xerophyte that is widely distributed in the Tianshan Mountains and adjacent areas. In total, 300 individuals from 25 populations were sampled, and 14 haplotypes were identified using two cpDNA sequences (rpS12-rpL20 and ndhF). A high level of total genetic diversity (H_T = 0.859) was detected, and this was attributed to the extensive distribution range, which included numerous large populations. The SAMOVA results suggested that the 25 populations were clustered into 4 major geographical groups; a similar divergence trend was found by constructing a BEAST phylogenetic tree and a network diagram. The AMOVA results revealed that significant genetic differentiation occurred among groups. Our results indicated a considerable correlation between genetic divergence and geographical distribution. Isolation due to complex mountain and desert geography might limit gene exchange between disjunct populations, resulting in high differentiation between geographical groups.     

Genetic structure and demographic history of Lymantria dispar (Linnaeus, 1758) (Lepidoptera: Erebidae) in its area of origin and adjacent areas

We analyzed the population genetic structure and demographic history of 20 Lymantria dispar populations from Far East Asia using microsatellite loci and mitochondrial genes. In the microsatellite analysis, the genetic distances based on pairwise F_ST values ranged from 0.0087 to 0.1171. A NeighborNet network based on pairwise F_ST genetic distances showed that the 20 regional populations were divided into five groups. Bayesian clustering analysis (K = 3) demonstrated the same groupings. The populations in the Korean Peninsula and adjacent regions, in particular, showed a mixed genetic pattern. In the mitochondrial genetic analysis based on 98 haplotypes, the median‐joining network exhibited a star shape that was focused on three high‐frequency haplotypes (Haplotype 1: central Korea and adjacent regions, Group 1; Haplotype 37: southern Korea, Group 2; and Haplotype 90: Hokkaido area, Group 3) connected by low‐frequency haplotypes. The mismatch distribution dividing the three groups was unimodal. In the neutral test, Tajima's D and Fu's FS tests were negative. We can thus infer that the Far East Asian populations of L. dispar underwent a sudden population expansion. Based on the age expansion parameter, the expansion time was inferred to be approximately 53,652 years before present (ybp) for Group 1, approximately 65,043 ybp for Group 2, and approximately 76,086 ybp for Group 3. We propose that the mixed genetic pattern of the inland populations of Far East Asia is due to these expansions and that the inland populations of the region should be treated as valid subspecies that are distinguishable from other subspecies by genetic traits.     

Genetic structure of Salsola junatovii (Chenopodiaceae) in the northern edge of the Taklimakan desert and conservational implications

Salsola junatovii Botsch. (Chenopodiaceae), a desert shrub endemic to the northern edge of the Taklimakan desert in China, is listed as an endangered species because of its narrow habitat requirements and a decline in numbers of individuals. In order to provide information on conservation strategies, we investigated the genetic patterns and evolutionary history of this species. Two cpDNA spacers, (psbA-trnH and psbK-psbI), and the nrITS sequence were employed in 116 individuals from 15 populations across the whole range of the species. A total of 7 haplotypes and 12 ribotypes were detected. Based on separate BEAST analyses, we suggest that intraspecific differentiation of S. junatovii occurred at 2.284–0.264 Ma for cpDNA haplotypes, and 1.304–0.1 Ma for nrITS ribotypes, both in the Pleistocene. Only mismatch analysis and neutrality tests for nrITS sequences over all populations were indicative of demographic expansion. Inconsistency between results from chloroplast and nuclear sequences was detected, apparently due to a stronger seed-mediated gene flow among populations than pollen-mediated flow under the severe desert conditions. Based on genetic structure results, one population, KMS, with the highest number of unique genotypes, should be a conservation priority; another nine populations containing unique genotypes were suggested for ex situ conservation.     

Remarkable genetic divergence of Gymnodiptychus dybowskii between south and north of Tianshan Mountain in northwest China

Gymnodiptychus dybowskii is endemic to Xinjiang, China and has been locally listed as protected animals. To investigate its genetic diversity and structure, specimens were collected from six localities in Yili River system and Kaidu River. Fragments of 1092bp Cyt b gene were sequenced for 116 individuals. A total of 21 haplotypes were found in all samples, and no haplotype was shared between Yili River system and Kaidu River population. Sequence comparisons revealed 123 variable sites, with eight singleton sites and 115 parsimony informative sites. For all the populations examined, the haplotype diversity (h) was 0.8298 ± 0.0226, nucleotide diversity (π) was 0.2521 ± 0.1202, and average number of pairwise nucleotide differences (k) was 275.3369 ± 118.5660. AMOVA analysis showed that the differences were significant for total populations except for Yili River system populations. The pairwise F_st values revealed same conclusion with AMOVA analysis: Kaidu River population was divergent from Yili River system populations. The genetic distance between two groups was 0.108 and the divergence time was estimated at 5.4–6.6 Ma, the uplift of Tianshan Mountain might have separated them and resulted in the genetic differentiation. The neutrality test and mismatch analysis indicated that both two groups of G. dybowskii had went through population expansion, the expansion time of Yili River system and Kaidu River population was estimated at 0.5859–0.7146 Ma and 0.5151–0.6282 Ma, respectively. The climate changes of Qinghai-Tibetan Plateau might have influenced the demographic history of G. dybowskii.     

Genetic diversity and population structure of Schizopygopsis younghusbandi Regan in the Yarlung Tsangpo River inferred from mitochondrial DNA sequence analysis

The genetic diversity and population structure of Schizopygopsis younghusbandi from six sites in the middle reach of the Yarlung Tsangpo River were examined based on mitochondrial DNA cytochrome b (Cyt b) gene and control region (D-loop) sequences. 1141 bp of complete Cyt b sequences and 737 bp of partial D-loop sequences for 153 individuals of S. younghusbandi were obtained by using PCR amplification and sequencing. The results showed that S. younghusbandi populations had high haplotype diversity and low nucleotide diversity, and the population genetic diversity of this species was at a low level. Analysis of molecular variance revealed that most genetic variation occurred within populations, and that genetic differentiation was at low or moderate levels. The network of haplotypes based on Cyt b gene showed that there were two groups within the examined populations. Neutrality tests and mismatch distribution analysis suggested that this species might have undergone a population expansion and the expansion time was estimated to be 0.25–0.46 Ma BP. All these results would be crucial to establish scientific strategies for effective conservation and sustainable exploitation of wild resources of S. younghusbandi.     

Genetic diversity and population structure analysis of Qinghai-Tibetan plateau schizothoracine fish (Gymnocypris dobula) based on mtDNA D-loop sequences

Gymnocypris dobula is a commercially important fishery species and mostly distributed in the freshwater of Tibetan plateau. In this study, genetic diversity and intraspecific population differentiation were examined by using mitochondrial DNA D-loop sequences in 97 individuals sampled from three localities (Pali, Lasa and Yanghu) in the Tibetan plateau. Two hundred and fourteen polymorphic sites and 50 haplotypes were defined among the three localities. Genetic diversity analysis showed that the highest genetic diversity level was found in Pali population. Phylogenetic relationships analysis results indicated that closer phylogenetic relationships were found between the Yanghu and Lasa populations. Genetic population differentiation analysis indicated that the majority of variation (84.91%) was attributed to variations among populations and the largest differentiation was found between Pali and Lasa localities (F_st = 0.874, Nm = 0.036). Moreover, the historical demographic events were assessed by implementing the mismatch distribution analysis, Fajima's D test and Fu's Fs test. The results indicated that the Pali population had undergone a demographic expansion, possibly within the last 0.163 MYA (Million Years Ago). Our study firstly identified the population genetic structure of the G. dobula, which could be helpful for artificial breeding, fishery stock identification and resource conservation for this species.     

Species‐level phylogeographical history of the endemic species Calligonum roborovskii and its close relatives in Calligonum section Medusa (Polygonaceae) in arid north‐western China

Quaternary climatic oscillations appear to have influenced the genetic diversity and evolutionary history of arid‐adapted plants. To understand the processes involved and reveal evolutionary relationships, haplotypes were examined from Calligonum roborovskii, an endemic species occurring in the arid zones across the desert regions of north‐western China, and seven other species also from Calligonum section Medusa, including C. gobicum, C. mongolicum and the narrow endemic species C. ebi‐nuricum, C. pumilum, C. taklimakanense, C. trifarium and C. yengisaricum. Forty‐three haplotypes were identified in 422 individuals from 51 natural populations, from variation of two plastid DNA intergenic spacers (rpl32‐trnL and ycf6‐psbM). A high level of total genetic diversity was found across species for which more than two populations were examined, including C. gobicum, C. mongolicum, C. pumilum and C. roborovskii. A distinct isolation‐by‐distance pattern in each of these species was suggested by the Mantel test, indicating that restricted gene flow caused high genetic differentiation among populations. Three haplotypes were shared by two or three species each, but the other 40 haplotypes were species‐specific. The 43 haplotypes split into three major clades, but not species‐specific lineages; most of the Calligonum species were not reciprocally monophyletic, probably due to incomplete lineage sorting or introgression. The identified haplotypes were dated to 1.97 Mya (95% highest posterior density: 2.95–0.99 Mya) and diverged until the late Pleistocene, possibly linked to aridification and enlargement of deserts caused by climate changes. Variation of desert habitats during the Pleistocene might play a key role in causing the divergence.     

Genetic population structure of miiuy croaker (Miichthys miiuy) in the Yellow and East China Seas base on mitochondrial COI sequences

The miiuy croaker (Miichthys miiuy) is an economically important fish species in China. However, information about population structure in this species is limited. In this study, the genetic diversity and population structure of four populations of Miichthys miiuy were investigated using a mitochondrial DNA (mtDNA) cytochrome c oxidase subunit I (COI) segment of 599 bp. A total of 35 polymorphic sites, and 31 haplotypes were defined in 109 individuals. High levels of haplotype diversity (0.731 ± 0.047) and low nucleotide diversity (0.00278 ± 0.0019) were observed in miiuy croaker populations. The results of AMOVA and pairwise fixation index (F_ST) analysis did not indicate any significant genetic differentiation among populations obtained from four sites (P > 0.05). In addition, no significant genealogical branches or clusters were recognized in the NJ tree. Analysis of neutral evolution and mismatch distribution suggested that the miiuy croaker might have experienced a recent population expansion.     

Genetic population structure of the swimming crab, Portunus trituberculatus in the East China Sea based on mtDNA 16S rRNA sequences

Portunus trituberculatus is a commercially important fishery species. In this study, intraspecific variation was examined by using mitochondrial DNA 16S rRNA gene in 213 individual crabs sampled from six localities along the coast of the East China Sea. Twenty-two polymorphic sites defined 25 distinct haplotypes, revealing a moderately high haplotype diversity and relatively low sequence divergence among the six localities. An excess of within population unique haplotypes at most sample locations were detected, which might influence genetic structure of the swimming crab populations. Neither neighbor-joining tree nor minimum spanning network (MSN) based on the haplotype data indicated distinct patterns of phylogeographic structure among the 25 haplotypes. Analyses of molecular variance (AMOVA) and F_ST statistics supported the hypothesis that population samples from the East China Sea were genetically nonhomogenous, indicating that gene flow might be restricted across those regions, despite the high potential of dispersal. In addition, tests of neutral evolution and analysis of mismatch distribution suggested that P. trituberculatus might have undergone a population expansion, possibly within the last 127,000 and 429,000 years. Our study unraveled the extant population genetic structure of the P. trituberculatus, and addressed the related fishery management issues including artificial breeding, fishery stock identification and conservation.     

A range wide geographic pattern of genetic diversity and population structure of Hexinia polydichotoma (Asteraceae) in Tarim Basin and adjacent areas

In order to investigate the genetic diversity and influence of climate oscillations on evolutionary processes of organisms in northwest China, we selected Hexinia polydichotoma, a species endemic to China, and examined the geographic pattern of genetic variation in its entire cover range, Tarim Basin and adjacent areas. In the study, 22 chloroplast (cp) haplotypes were identified based on two cpDNA sequences (trnH–psbA and ycf6–psbM), and five ITS sequence variations were found. Shown in the cp haplotype network, the two common cp haplotypes mainly distributed along the northern and southern rims of the basin respectively and intersected in the center of the basin, whereas in the ITS network, ITS genotype 1 was widespread across the whole distribution area, and rare genotypes were concentrated in the western rim of the basin. Genetic variation primarily occurred among populations and SAMOVA groups. Fragmented desert habitat may have caused gene flow barrier among populations or groups far from each other, leading to significant genetic differentiation at these levels. It appears that expansion and contraction cycles of river systems and oases during the middle Pleistocene was the source of the observed fragmentation. Geographic range expansion in H. polydichotoma was supported by the significant value for Tajima's D, Fu's F_S, and by a unimodal mismatch distribution. It was possible that the enlargement of the Taklimakan Desert during the middle Pleistocene may have provided appropriate conditions for the range dispersal. We identified western rim of the basin as the center of genetic diversity of H. polydichotoma based on the present dataset.     

Diversification and vicariance of desert plants: Evidence inferred from chloroplast DNA sequence variation of Lagochilus ilicifolius (Lamiaceae)

To document arid Northern China as a diversification center for desert plants, and to better understand the mechanisms of desert taxa diversification, we used five cpDNA spacers (trnL-trnF, rps16, psbA-trnH, psbK-psbI and trnS-trnG) to investigate Lagochilus ilicifolius in all distribution areas, including Northern China, Mongolia and Russia. Phylogenetic analyses showed that L. ilicifolius comprises two distinctive lineages, one distributed in China, and another in Mongolia–Russia. Our data confirmed that arid Northern China, is a distinctive area with many endemic genera. Biogeographic inferences, based on a Bayesian uncorrelated lognormal model together with molecular dating, suggested that the main diversification within the species occurred in the Pleistocene (ca. 1.38–0.3 Ma), resulting from the transition of the climate of Eurasia to a dry-cold pattern as well as the desertification caused the vicariance of desert plants.     

Phylogeography of Chrysanthemum indicum L. (Compositae) in China based on trnL-F sequences

Using chloroplast DNA (cpDNA) trnL-F sequences, we studied the phylogeographic pattern and demographic history of Chrysanthemum indicum with 81 individuals from 27 populations representing the partial geographic range in China. Fifteen haplotypes were identified based on nucleotide variation. Relatively high levels of haplotype diversity (Hd = 0.783) and low levels of nucleotide diversity (Pi = 0.00521) were detected in C. indicum. An analysis of molecular variance (AMOVA) indicated that most of the genetic variation (59.9%, P < 0.001) was partitioned within populations. Genetic differentiation among populations was high (F_ST = 0.35), suggesting that there is little gene flow among populations. Nested clade analysis showed that 15 haplotypes formed twelve 1-step clades, five 2-step clades and two 3-step clades, with four clades showing significant geographical association. We concluded that contiguous range expansion and restricted gene flow/dispersal but with some long distance dispersal seemed to have shaped the present-day distribution of the haplotypes of C. indicum in China.     

Multiple Pleistocene refugia and Holocene range expansion of an abundant southwestern American desert plant species (Melampodium leucanthum,Asteraceae)

Pleistocene climatic fluctuations had major impacts on desert biota in southwestern North America. During cooler and wetter periods, drought‐adapted species were isolated into refugia, in contrast to expansion of their ranges during the massive aridification in the Holocene. Here, we use Melampodium leucanthum (Asteraceae), a species of the North American desert and semi‐desert regions, to investigate the impact of major aridification in southwestern North America on phylogeography and evolution in a widespread and abundant drought‐adapted plant species. The evidence for three separate Pleistocene refugia at different time levels suggests that this species responded to the Quaternary climatic oscillations in a cyclic manner. In the Holocene, once differentiated lineages came into secondary contact and intermixed, but these range expansions did not follow the eastwardly progressing aridification, but instead occurred independently out of separate Pleistocene refugia. As found in other desert biota, the Continental Divide has acted as a major migration barrier for M. leucanthum since the Pleistocene. Despite being geographically restricted to the eastern part of the species’ distribution, autotetraploids in M. leucanthum originated multiple times and do not form a genetically cohesive group.     

Genetic structure of loggerhead turtle (Caretta caretta) populations in Turkey

This study investigates genetic structure of loggerhead turtle (Caretta caretta) for conservation purposes using both mitochondrial and nuclear DNA markers. To evaluate overall genetic structure of C. caretta at the Turkish nesting beaches and to compare it with previously published data from other nesting sites of the Mediterranean Sea, we have studied 256 hatchlings from 18 different nesting beaches. Seven distinct haplotypes were detected, of which three have previously been reported from the Mediterranean and one from the Atlantic. The remaining three haplotypes are described for the first time for C. caretta in the present study. Distribution of these haplotypes among the nesting sites showed a significant genetic structuring, indicating that females are philopatric and that gene flow among populations is restricted. Both mtDNA and microsatellite analyses determined genetic structuring (mtDNA: γ_st = 0.214, p < 0.01; nDNA F_st = 0.0004 p < 0.05) among nesting aggregates of C. caretta throughout the study area and enabled the detection of the different haplotypes to inform conservation strategies.     

Low mtDNA Cytb diversity and shallow population structure of Eleutheronema tetradactylum in the East China Sea and the South China Sea

Eleutheronema tetradactylum is an economically important fish species in China water. To investigate the genetic diversity and describe population structure of it, an 1151 base pair (bp) fragment of the mitochondrial DNA Cytb sequence was analyzed in 120 individuals from four populations in the East China Sea and the South China Sea. A total of 16 haplotypes were defined by 24 variable nucleotide sites. High level of haplotype diversity and low nucleotide diversity were observed in all populations. The results of AMOVA detected that 89.44% of the genetic variation occurred within populations. Significant genetic differentiations were detected among populations (0.05097, P < 0.05), but no large-scale regional differences were detected. Analysis of neutral evolution and mismatch distribution suggested no recent population expansion happened. The present results provided new information for genetic assessment, fishery management and conservation of this species.     

Population structure and genetic variation in the genus Dipteronia Oliv. (Aceraceae) endemic to China as revealed by cpSSR analysis

The genus Dipteronia Oliv. endemic to central and southern China consists of two species, D. sinensis Oliv. and D. dyerana Henry, both of them are rare and endangered. In the present study, genetic variation in 17 populations of D. sinensis and four populations of D. dyerana was estimated using ten polymorphic chloroplast microsatellite loci (cpSSR). Forty-nine chloroplast haplotypes were identified from 204 individuals analyzed. Thirty-nine haplotypes were found in D. sinensis, while ten in D. dyerana. No haplotype was shared between the species. AMOVA analysis revealed that the majority of the genetic variation was partitioned among populations within D. sinensis (F _ST = 0.7980) and D. dyerana (F _ST = 0.654). Cluster analysis grouped the 21 populations into two groups according to their species delimitation. The populations of D. sinensis were further divided into two subgroups corresponding to their geographical distributions. Correlation analysis revealed significant correlation between geographical distance and genetic distance of these populations (r = 0.326, p < 0.01 for D. sinensis vs. r = 0.777, p < 0.05 for D. dyerana). Genetic structure of these populations and a calculated pollen-to-seed flow ratio of (3.2:1 vs. 0.6:1) within the species suggested that little gene flow has occurred among the populations over an extended period. Thus, it implies that the genus Dipteronia might have experienced a genetic bottleneck and limited expansion during its evolutionary history.     

Genetic diversity and population structure of miiuy croaker populations in East China Sea revealed by the mitochondrial DNA control region sequence

The genetic diversity and population structure of the miiuy croaker (Miichthys miiuy) were investigated. Seventy-six individuals were sampled from six localities of the East China Sea. Genetic variation in DNA sequences was examined from the 5′-end of the mitochondrial DNA control region. High levels of haplotype diversity (h = 0.99333) in the control region were detected, indicating a high level of genetic diversity. A total of 81 polymorphic sites were found, and 65 haplotypes were defined. The results of AMOVA analysis indicated that 97.9% of the genetic variation contained within populations and 2.14% occurred among populations. No significant genealogical branches or clusters were recognized on the NJ tree. Tests of neutral evolution and analysis of mismatch distribution suggested that miiuy croaker might have undergone a population expansion. The knowledge on genetic diversity and genetic structure will be crucial to establish appropriate fishery management stocks for the species.     

Multiple introgression events and range shifts in Schizocodon (Diapensiaceae) during the Pleistocene

Range shifts during the Pleistocene shaped the unique phylogeographical structures of numerous species. Accompanying species migration, sister taxa may have experienced multiple introgression events. Here, we report the signature of introgression events in multiple areas in Schizocodon, herbs endemic to Japan, using amplified fragment length polymorphism (AFLP) fingerprinting and plastid DNA haplotyping in 48 populations. Although the present distributions of S. soldanelloides and S. ilicifolius are mainly allopatric, the species share plastid DNA haplotypes in each region (north‐eastern, north‐central, south‐central and south‐western Japan); in contrast, the specific groups were highly supported by AFLP analyses. These results support the occurrence of multiple introgression events in Schizocodon. Notably, the disjunct plastid haplotypes found only in S. ilicifolius var. intercedens suggest complete plastid DNA replacement at local areas from S. soldanelloides into S. ilicifolius var. ilicifolius. Furthermore, we found that S. soldanelloides experienced range contraction and expansion during glacial and interglacial cycles based on mismatch distribution analysis and ecological niche modelling. Based on several pieces of evidence, our study supports the idea that historical range shifts associated with Pleistocene climatic oscillations favoured multiple and regional introgression events in Schizocodon. © 2013 The Linnean Society of London, Botanical Journal of the Linnean Society, 2013, 173 , 46–63.