Phylogeographic studies of schizothoracine fishes on the central Qinghai-Tibet Plateau reveal the highest known glacial microrefugia

Across the highest mountain on earth: discordant phylogeographic patterns and recent dispersal of Tibetan stone loaches (Triplophysa) in the Himalayas

Phylogeographic congruence among co-distributed taxa is regarded as an inherent inference to vicariance events. Nonetheless, incongruent patterns of contemporary lineage divergence among taxa indicated that species differ in their response to common past events. To investigate the role of past events, ecological traits and lineage diversification time in shaping the contemporary phylogeographic patterns, comparative analyses were conducted for Tibetan stone loaches in the Himalayas using three gene markers and two ecological traits (depth of caudal peduncle in their length and presence/absence of posterior chamber of the air bladder). By a thorough sampling in two flanks of the Himalayas, the authors detected that phylogenetic breaks were spatially discordant and divergences of populations were also temporally asynchronous in co-distributed loaches. Estimated divergence time using fossil-calibrated node dating indicated that the Tibetan stone loaches colonised into the south flank of the Himalayas until the Pleistocene. The demographic expansions were also disconcerted between populations in north and south flanks, or east and west Himalayas. Ongoing gene flows between populations in north and south sides implied that the Himalayas do not strictly impede dispersal of cold-adapted species. The results highlight that the quaternary climatic oscillation, in conjunction with ecological traits and lineage diversification time, shaped contemporary phylogenetic patterns of stone loaches in the Himalayas and provide new insights into the biodiversity and composition of species in the Himalayas and surrounding region.

Pleistocene Landscape Dynamics Drives Lineage Divergence of a Temperate Freshwater Fish Gobio rivuloides in Coastal Drainages of Northern China

Understanding historical processes underlying lineage distribution patterns is a primary goal of phylogeography. We selected Gobio rivuloides (Cypriniformes: Gobionidae) as a model to improve our knowledge about how intraspecific genetic divergence of freshwater fishes arises in coastal drainages of northern China via statistical analysis using cytochrome b gene. The time-calibrated phylogeny of G. rivuloides showed the divergence of two major lineages (I and II) at ~0.98 Ma (million years ago). Lineage I can be divided into two sub-lineages (I-A and I-B) with a divergence time of ~0.83 Ma. Sub-lineage I-A inhabits the Amur River, and sub-lineage I-B lives in the Luan River and Liao River. Lineage II is distributed in the Yellow River and Hai River, with close genetic relationships between the two drainages, and can be split into two sub-lineages (II-C and II-D) with a divergence time of ~0.60 Ma. Our findings indicate that the splitting of lineages and sub-lineages could be attributed to geographic isolation caused by the formation of the Bohai Sea, river capture, and the episodic hydrologic closing of a paleolake during the late Lower-Middle Pleistocene. It is also the first report we know of displaying a clear phylogeographic break for freshwater fishes across coastal drainages in northern China.

Microsatellites as Molecular Markers with Applications in Exploitation and Conservation of Aquatic Animal Populations

A large number of species and taxa has been studied for genetic polymorphism. Microsatellites have been known as hypervariable neutral molecular markers with the highest resolution power in comparison with any other markers. However, the discovery of a new type of molecular marker—single nucleotide polymorphism (SNP) has put the existing applications of microsatellites to the test. To ensure good resolution power in studies of populations and individuals, a number of microsatellite loci from 14 to 20 was often used, which corresponds to about 200 independent alleles. Recently, these numbers have tended to be increased by the application of genomic sequencing of expressed sequence tags (ESTs), and the choice of the most informative loci for genotyping depends on the aims of research. Examples of successful applications of microsatellite molecular markers in aquaculture, fisheries, and conservation genetics in comparison with SNPs have been summarized in this review. Microsatellites can be considered superior markers in such topics as kinship and parentage analysis in cultured and natural populations, the assessment of gynogenesis, androgenesis and ploidization. Microsatellites can be coupled with SNPs for mapping QTL. Microsatellites will continue to be used in research on genetic diversity in cultured stocks, and also in natural populations as an economically advantageous genotyping technique.

Phylogenetic divergence associated with climate oscillations and topology illustrates the dispersal history of Formosan sambar deer (Rusa unicolor swinhoii) in Taiwan

The island of Taiwan represents an ideal context for studying the effects of climatic oscillations and topographic variation on large herbivores due to its varied tropical to sub-tropical climate zones at different elevational ranges. We explored the phylogenetics of Formosan sambar deer (Rusa unicolor swinhoii) using the control region of the mitochondrial genome. We detected 18 haplotypes among 454 sequences across the island and grouped them into six regions based on SAMOVA, with 68.78% variance among regions. A Bayesian phylogenetic dendrogram revealed two spatially segregated genetic clades. Neutrality tests and Bayesian skyline plots uncovered different demographic expansion histories for the two clades. We further tested divergence times and chronology to propose potential phylogenetic scenarios, which were examined using approximate Bayesian computation. Finally, we present a credible hypothesis for a glacial refugium in the northern part of the Central Mountain Range. Subsequent secondary contact between the two clades during interglacial periods has led to the extant genetic structure of Formosan sambar deer.

Changes in the Suitable Habitats of Three Endemic Fishes to Climate Change in Tibet

As one of the most sensitive regions to global climate change, Tibet is subject to remarkable changes in biota over the past decades, including endemic fish species. However, no study has attempted to predict the changes in the distribution of Tibetan fishes, leaving a great blank for aquatic conservation in Tibet. Based on the Maximum Entropy model (MaxEnt), this study predicted the changes in the suitable habitats of three endemic fish species, including two species mainly inhabiting the rivers (Glyptosternon maculatum, Oxygymnocypris stewartii) and one species mainly inhabiting lakes (Gymnocypris selincuoensis) in Tibet under two representative concentration pathways (RCP2.6 and RCP8.5) under two future scenarios (2050 and 2090), and explored the impact of the barrier effects of hydropower projects on the suitable habitats of fish. The results showed that under the four scenarios, the net change in the suitable habitats of the G. maculatum was negative (-2.0--18.8%), while the suitable habitats of the O. stewartii and G. selincuoensis would be expanded, with the net change of 60.0-238.3% and 46.4-56.9%, respectively. Under different scenarios, the suitable habitats of the three species had a tendency to migrate to a higher elevation, and the largest expansion in the range of migration was projected to occur under the 2090-RCP8.5 scenario. In addition, due to the impact of the hydropower projects, the ability of G. maculatum to obtain new suitable habitats from climate change would be reduced by 2.0-8.1%, which was less than the loss induced by climate change (5.5-25.1%), while the suitable habitats of O. stewartii would be reduced by 3.0-9.7%, which was more than the impact of climate change (about 1%). The results of this study have guiding significance for the conservation and management of fish resources diversity in the Yarlung Tsangpo River basin and Siling Co basin of Tibet, and also provide a reference for the coordination and scientific planning of hydropower projects in Tibet.

Population genetics analysis of Tolai hares (Lepus tolai) in Xinjiang, China using genome-wide SNPs from SLAF-seq and mitochondrial markers

The main topic of population genetics and evolutionary biology is the influence of the ecological environment, geographical isolation, and climatic factors on population structure and history. Here, we estimated the genetic diversity, genetic structure, and population history of two subspecies of Tolai hares (Lepus tolai Pallas, 1778), L. t. lehmanni inhabiting Northern and Northwest Xinjiang and L. t. centrasiaticus inhabiting Central and Eastern Xinjiang using SNP of specific-length amplified fragment sequencing (SLAF-seq) and four mitochondrial DNA (mtDNA). Our results showed a relatively high degree of genetic diversity for Tolai hares, and the diversity of L. t. lehmanni was slightly higher than that of L. t. centrasiaticus, likely due to the more favorable ecological environment, such as woodlands and plains. Phylogenetic analysis from SNP and mtDNA indicated a rough phylogeographical distribution pattern among Tolai hares. Strong differentiation was found between the two subspecies and the two geographical groups in L. t. centrasiaticus, possibly due to the geographical isolation of mountains, basins, and deserts. However, gene flow was also detected between the two subspecies, which might be attributed to the Tianshan Corridor and the strong migration ability of hares. Tolai hare population differentiation occurred at approximately 1.2377 MYA. Population history analysis based on SNP and mtDNA showed that the Tolai hare population has a complex history and L. t. lehmanni was less affected by the glacial event, possibly because its geographic location and terrain conditions weaken the drastic climate fluctuations. In conclusion, our results indicated that the joint effect of ecological environment, geographic events, and climatic factors might play important roles in the evolutionary process of L. t. lehmanni and L. t. centrasiaticus, thus resulting in differentiation, gene exchange, and different population history.

Geological and climatic influences on population differentiation of the Phrynocephalus vlangalii species complex (Sauria: Agamidae) in the northern Qinghai-Tibet Plateau

Extremely heterogeneous topography and complex paleoclimate history of the Qinghai-Tibet Plateau (QTP) have a key role in promoting genetic divergence and lineage/species formation. Here, we sequenced one nuclear and three mitochondrial markers of 532 individuals from the entire range of the Phrynocephalus vlangalii species complex including two species, P. putjatai and P. vlangalii, endemic to the northern QTP. We integrated multilocus phylogeny, demographic analysis and geographic barrier detection to evaluate the population structure and dynamics. We found a new mitochondrial clade (PV-I) in the Gonghe County population of P. vlangalii, partial mitochondrial DNA replacement within P. vlangalii and complete mitochondrial DNA replacement between P. putjatai and P. vlangalii. Neutrality test, mismatch distribution analysis and Extended Bayesian Skyline Plot (EBSP) analysis all supported a significant expansion of the Qaidam Basin population of P. vlangalii (PV-II-2) from 0.091 to 0.026 Ma after Penultimate Glaciation. The uplift of the Arjin and Anyemanqen Mountains during the Kunhuang Movement (∼1.2 Ma) split populations of P. vlangalii in Akesai, Qaidam Basin and source of the Yellow River. The uplift of the Elashan Mountains during the second phase of the Qingzang Movement (∼2.5 Ma) contributed to the divergence of the Gonghe County population of P. vlangalii from other conspecific populations. The third phase of the Qingzang Movement (∼1.7 Ma) contributed to the divergence of the Xinghai population of P. vlangalii from P. putjatai and to the divergence of the northern populations of P. putjatai from the southern conspecific populations. Our data support the idea that the geological and climatic changes following the orogeny of the QTP may have promoted population differentiation and shaped the current population patterns of the P. vlangalii species complex in the northeastern QTP.

Characterization of the complete mitogenome of Gymnocypris dobula (Günther, 1868) (Cypriniformes: Cyprinidae)

Gymnocypris dobula, classified into the highly specialized Schizothoracine fish, is endemic to Tibet, China. The complete mitochondrial DNA sequence of G. dobula was 16,728 base pairs in length and comprised 22 transfer RNA genes, 13 protein-coding genes, two ribosomal RNA genes as well as one control region as in a typical vertebrate mitochondrial DNA gene. The ML and BI trees showed that G. dobula was most closely related to Gymnocypris scleracanthus within the highly specialized group. This mitogenome provides new molecular data for further taxonomic and phylogenetic studies of Schizothoracine fish.

Melatonin: A smart molecule in the DNA repair system

DNA repair is an important pathway for the protection of DNA molecules from destruction. DNA damage can be produced by oxidative reactive nitrogen or oxygen species, irritation, alkylating agents, depurination and depyrimidination; in this regard, DNA repair pathways can neutralize the negative effects of these factors. Melatonin is a hormone secreted from the pineal gland with an antioxidant effect by binding to oxidative factors. In addition, the effect of melatonin on DNA repair pathways has been proven by the literature. DNA repair is carried out by several mechanisms, of which homologous recombination repair (HRR) and non-homologous end-joining (NHEJ) are of great importance. Because of the importance of DNA repair in DNA integrity and the anticancer effect of this pathway, we presented the effect of melatonin on DNA repair factors regarding previous studies conducted in this area.

Phylogeography of Triplophysa stenura (Nemacheilidae): Responded to the Mid-Pleistocene Climate Transition in the Qinghai-Tibetan Plateau

Pleistocene climatic oscillations have greatly influenced the evolutionary histories and distribution patterns of most extant species. In this study, the phylogeographic patterns and evolutionary history of Triplophysa stenura were investigated. A total of 359 individuals from 19 locations covering the species' distribution range were collected, and two mitochondrial genes (COI and Cyt b) and the recombination activating protein 1 gene (Rag1) were analyzed. Two deeply divergent lineages, corresponding geographically to the northern and southern Tanggula Mountain, were observed, and shallow divergences were found within the southern and central Qinghai-Tibetan plateau (QTP). The estimated divergence time between the northern and southern Tanggula lineages was around 1.07 Mya. Within the southern Tanggula Mountain lineages, the Nu and Subansiri rivers populations were split about 0.74 Mya, and the southwestern and central QTP populations were divided with the southeastern QTP ones approximately 0.67 Mya. The divergence times of the lineages were matched with two major glaciations in QTP (the Xixiabangma Glaciation, 1.17-0.80 Mya and Kunlun Glaciation, 0.72-0.62 Mya). Together with demographic history analyses, our results highlighted that major glacial cycles during the mid-Pleistocene played a major role in sculpting the distribution pattern of T. stenura, and led to the gene homogenization crossing the drainage populations in the southwestern and central QTP.

Processes that drive the population structuring of Jenynsia lineata (Cyprinidontiformes, Anablepidae) in the La Plata Basin

The distribution of genetic diversity across a species distribution range is rarely homogeneous, as the genetic structure among populations is related to the degree of isolation among them, such as isolation by distance, isolation by barrier, and isolation by environment. Jenynsia lineata is a small viviparous fish that inhabits a wide range of habitats in South America. To decipher the isolation processes that drive population structuring in J. lineata, we analyzed 221 sequences of the mitochondrial cytochrome c oxidase I gene (COI), from 19 localities. Then, we examined the influence of the three most common types of isolation in order to explain the genetic variation found in this species.Our results revealed a marked structuration, with three groups: (a) La Plata/Desaguadero Rivers (sampling sites across Argentina, Uruguay, and Southern Brazil), (b) Central Argentina, and (c) Northern Argentina. A distance-based redundancy analysis, including the explanatory variables geographical distances, altitude, latitude, and basin, was able to explain up to 65% of the genetic structure. A variance partitioning analysis showed that the two most important variables underlying the structuration in J. lineata were altitude (isolation by environment) and type of basin (isolation by barrier).Our results show that in this species, the processes of population diversification are complex and are not limited to a single mechanism. The processes that play a prominent role in this study could explain the high rate of diversity that characterizes freshwater fish species. And these processes in turn are the basis for possible speciation events.

When rice gets the chills: comparative transcriptome profiling at germination shows WRKY transcription factor responses

Rice is vital for food security. Due to its tropical origin, rice suffers from cold temperatures that affect its entire life cycle. Key genes have been identified involved in cold tolerance. WRKYs are generally downstream of the MAPK cascade and can act together with VQ proteins to regulate stress-responsive genes. Chilling treatment was applied at germination to two rice genotypes (tolerant and sensitive). Shoots at S3 stage were collected for RNA-seq to identify OsWRKY, OsMAPKs and OsVQs expression. Relationships among MAPKs, WRKYs and VQs were predicted through correlation analysis. OsWRKY transcriptional regulation was predicted by in silico analysis of cis-regulatory elements. A total of 39 OsWRKYs were differentially expressed. OsWRKY21, OsWRK24 and OsWRKY69 are potential positive regulators, while OsWRKY10, OsWRK47, OsWRKY62, OsWRKY72 and OsWRKY77 are potential negative regulators, of chilling tolerance. 12 OsMAPKs were differentially expressed. OsMAPKs were downregulated and negatively correlated with the upregulated OsWRKYs in the tolerant genotype. 19 OsVQs were differentially expressed, three and six OsVQs were positively correlated with OsWRKYs in the tolerant and sensitive genotypes, respectively. Seven differentially expressed OsWRKYs have cold-responsive elements in their promoters and five upregulated OsWRKYs in the tolerant genotype contained the W-box motif. Chilling causes changes in OsWRKY, OsMAPK and OsVQ gene expression at germination. OsWRKYs may not act downstream of the MAPK cascade to coordinate chilling tolerance, but OsWRKYs may act with VQs to regulate chilling tolerance. Candidate OsWRKYs are correlated and have a W-box in the promoter, suggesting an auto-regulation mechanism.

Comparative transcriptome analysis of scaled and scaleless skins in Gymnocypris eckloni provides insights into the molecular mechanism of scale degeneration

Background

The scale degeneration is thought to be related to the adaptation to the extreme environment with cold climate and high-altitude in schizothoracine fishes. Gymnocypris eckloni, a schizothoracine fish living in plateau waters with the elevation above 2500 m, is nearly esquamate and only covered with shoulder scales and anal scales, making it a good model species to study the molecular mechanism of scale degeneration.

Results

The transcriptomes of shoulder scaled skins (SSS), anal scaled skins (ASS) and scaleless skins (NSS) were sequenced and analyzed in G. eckloni at the age of 1 year. Histological examination showed that shoulder scale had completed its differentiation and anal scale just initiated the differentiation. A total of 578,046 unigenes were obtained from the transcriptomes, with 407,799 unigenes annotated in public databases. A total of 428 and 142 differentially expressed unigenes (DEUs) were identified between SSS and NSS, and between ASS and NSS, respectively, with 45 DEUs that were overlapped. Annotation analysis indicated that these DEUs were mainly enriched in Gene Ontology (GO) terms and KEGG pathways associated with bone and muscle formation, such as myofibril, contractile fiber, cytoskeletal protein binding, muscle structure development, cardiac muscle contraction, hypertrophic cardiomyopathy (HCM) and calcium signaling pathway.

Conclusions

Our results would provide insights into the molecular mechanisms of scale degeneration in G. eckloni and other congeneric fishes. In addition, the transcriptome data provides candidate genes and markers for future studies.

Changing climate reallocates the carbon debt of frequent-fire forests

Ongoing climate change will alter the carbon carrying capacity of forests as they adjust to climatic extremes and changing disturbance regimes. In frequent-fire forests, increasing drought frequency and severity are already causing widespread tree mortality events, which can exacerbate the carbon debt that has developed as a result of fire exclusion. Forest management techniques that reduce tree density and surface fuels decrease the risk of high-severity wildfire and may also limit drought-induced mortality by reducing competition. We used a long-term thinning and burning experiment in a mixed-conifer forest to investigate the effects of the 2012-2015 California drought on forest carbon dynamics in each treatment, including the carbon emissions from a second-entry prescribed fire that followed the drought. We assessed differences in carbon stability and tree survival across treatments, expecting that both carbon stability and survival probability would increase with increasing treatment intensity (decreasing basal area). Additionally, we analyzed the effects of drought- mortality on second-entry burn emissions and compared emissions for the first- and second-entry burns. We found a non-linear relationship between treatment intensity and carbon stability, which was in part driven by varying relationships between individual tree growing space and survival across treatments. Drought mortality increased dead tree and surface fuel carbon in all treatments, which contributed to higher second-entry burn emissions for two of the three burn treatments when compared to the first burn. Our findings suggest that restoration treatments will not serve as a panacea for ongoing climate change and that the carbon debt of these forests will become increasingly unstable as the carbon carrying capacity adjusts to severe drought events. Managing the carbon debt with prescribed fire will help reduce the risk of additional mortality from wildfire, but at an increasing carbon cost for forest management.

Phylogeography of Schizopygopsis malacanthus Herzenstein (Cypriniformes, Cyprinidae) in relation to the tectonic events and Quaternary climatic oscillations in the Shaluli Mountains Region

The tectonic events and Quaternary climatic oscillations in the Shaluli Mountains (the margin of the southeastern Tibetan Plateau) had an extensive effect on the genetic structure and distribution patterns of this region's terrestrial fauna and flora. It is not yet clear whether similar mechanisms influence this region's fish fauna. Schizopygopsis malacanthus is limited to high-elevation rivers and lakes, and exhibits distinct adaptations to the mountains of, and near to, the Qinghai-Tibet Plateau (QTP). Therefore, this species is a good candidate for investigating patterns of genetic variation resulting from palaeoenvironmental fluctuations in the Shaluli Mountains (China). Here, we used microsatellite and mitochondrial DNA control region sequences to analyze six populations of S. malacanthus collected from the Jinsha and Yalong River drainages. Genealogical analyses identified four maternal lineages and perhaps even four putative species, of which the Ouqu River lineage played a pivotal role during the course of the species' evolution. Two lineages from the Yalong River drainage did not cluster together, whereas those from different drainages grouped together, suggesting tectonic event impacts that possibly altered regional river drainages have been highly influential in population connectivity and gene flow. Population genetic analysis indicated that the geographic barriers and this species preference for higher elevations both played key roles in the divergence of S. malacanthus populations. Demographic tests suggested large-scale spatial synchrony in population fluctuations of S. malacanthus, accompanying dramatic Pleistocene climatic oscillations. It appears that palaeoenvironmental changes in the Shaluli Mountains influenced the distribution and evolution of studied S. malacanthus populations, which provide important information about the factors that influenced the phylogeographic history of this region's fish fauna. Additionally, our study also has implications for ongoing conservation of this vulnerable fish.

Prevalence and genetic diversity of Echinorhynchus gymnocyprii (Acanthocephala: Echinorhynchidae) in schizothoracine fishes (Cyprinidae: Schizothoracinae) in Qinghai-Tibetan Plateau, China

Background

The schizothoracine fishes, an excellent model for several studies, is a dominant fish group of the Qinghai-Tibet Plateau (QTP). However, species populations have rapidly declined due to various factors, and infection with Echinorhynchus gymnocyprii is cited as a possible factor. In the present study, the molecular characteristics of E. gymnocyprii in four species of schizothoracine fishes from the QTP were explored.

Methods

We investigated the infection status of E. gymnocyprii in 156 schizothoracine fishes from the upper Yangtze River, upper Yellow River, and Qinghai Lake in Qinghai Province, China. The complete internal transcribed spacer (ITS) of the ribosomal RNA (rRNA) gene and part of the mitochondrial cytochrome c oxidase subunit 1 (cox1) gene of 35 E. gymnocyprii isolates from these fishes were sequenced and their characteristics analyzed. In addition, we inferred phylogenetic relationships of the E. gymnocyprii populations based on the rRNA-ITS and cox1 sequences.

Results

The total prevalence of E. gymnocyprii in schizothoracine fishes was 57.69% (90/156). However, the prevalence among different species as well as that across the geographical locations of the schizothoracine fishes was significantly different. The results of sequence analysis showed that the four E. gymnocyprii populations from different hosts and regions of Qinghai Province were conspecific, exhibiting rich genetic diversity. Phylogenetic analysis based on rRNA-ITS and cox1 sequences supported the coalescence of branches within E. gymnocyprii; the cox1 gene of E. gymnocyprii populations inferred some geographical associations with water systems. In addition, three species of schizothoracine fishes were recorded as new definitive hosts for E. gymnocyprii.

Conclusions

To the best of our knowledge, this is the first molecular description of E. gymnocyprii populations in schizothoracine fishes from the Qinghai-Tibet Plateau that provides basic data for epidemiological surveillance and control of acanthocephaliasis to protect endemic fish stocks.

Pelagic-benthic resource polymorphism in Schizopygopsis thermalis Herzenstein 1891 (Pisces, Cyprinidae) in a headwater lake in the Salween River system on the Tibetan Plateau

Resource polymorphism is a ubiquitous phenomenon in vertebrates and may represent a critical intermediate stage in speciation. Freshwater lakes in high-altitude areas represent a natural system for understanding resource polymorphism in fishes benefiting from diverse lacustrine environments and species-poor fish assemblages. We report resource polymorphism in a cyprinid fish, Schizopygopsis thermalis, in Lake Amdo Tsonak Co, a headwater lake in the upper Salween River system. Two discrete intraspecific morphs, planktivorous and benthivorous, were identified according to geometric morphometrics and traditional univariate linear measures. The planktivorous morph exhibits a longer head, longer upper and lower jaw, larger asymptotic standard length (L∞ ), lower growth rate (k), and higher growth performance index (φ) than the benthivorous morph. With respect to descriptive traits, the planktivorous morph possesses a large, terminal mouth and obvious mucus pores on the cheek and chin, while the benthivorous morph is characterized by a more inferior mouth with a sharpen horny edge on the lower jaw and unconspicuous mucus pores. The discrete pelagic-benthic resources and low interspecific competition in the lake system might drive the initial differentiation of the two morphs, and partial spatial reproductive isolation in breeding further maintains and reinforces the differences between them.

Quaternary climate change and habitat preference shaped the genetic differentiation and phylogeography of Rhodiola sect. Prainia in the southern Qinghai-Tibetan Plateau

There are two long-standing biogeographic hypotheses regarding the glacial survival of plant species in the Qinghai-Tibetan Plateau (QTP): the in situ survival hypothesis and the tabula rasa hypothesis. We tested these two hypotheses in a phylogeographic study of Rhodiola sect. Prainia, a monophyletic section with ecologically divergent lineages. Molecular data from the nuclear internal transcribed spacer, six plastid markers and 13 nuclear microsatellite loci were analyzed for 240 individuals from 19 populations of this section. Environmental data were used to analyze the niches of major phylogenetic lineages within this section and to model changes in their distributions since the Last Glacial Maximum (LGM). We found that Rhodiola sect. Prainia consists of three evolutionary lineages: all populations of R. stapfii, R. prainii populations at the southern edge of the QTP, and R. prainii populations in the interior part of the QTP. During the LGM, the survival of R. prainii in the interior part of the QTP corresponded with the in situ survival hypothesis, while R. stapfii most probably survived the LGM in a manner corresponding with the tabula rasa hypothesis. The evolutionary history of different lineages of this section was shaped by topography, climate change, and lineage-specific habitat preferences.

DNA barcodes and their characteristic diagnostic sites analysis of Schizothoracinae fishes in Qinghai province

The Qinghai-Tibetan Plateau (QTP), the source and upper reaches of many Asian rivers, are crisscrossed by rivers and dotted with lakes. Schizothoracinae fishes, species native to the QTP, are distributed widely through these rivers and lakes. Over the past decades, ecological protection has become increasingly intense. The rapid acquisition of the genetic information and accurate gene sequence database are assumed to play an important role in the conservation of species diversity and biodiversity. In this study, 153 COI sequences (648bp in length) covering 13 species in 8 genera of Schizothoracinae fishes in Qinghai Province were used to determine whether barcode could identify Schizothoracinae species accurately. The average Kimura two parameter (K2P) genetic distances within and among species were 0.35% and 8.83%, respectively. The maximum K2P distance within species was observed in Gymnocypris eckloni (1.36%) while minimum K2P distance among species was observed between Chuanchia labiosa and Schizopygopsis pylzovi (0.23%). Overlaps existed in K2P distance intra- and inter- species based on both the genes. Eleven groups with 9 single-species groups and 2 multi-species groups were identified through Automatic Barcode Gap Discovery System, which were consistent with the overlaps of K2P distance. 96.7% as the accurate ratio for COI barcode was calculated and high solution was observed in the phylogenetic trees based on COI gene and Cyt b gene. Except for the similar results based on two genes above, COI barcode was more economical than Cyt b gene. The SOM model successfully predicted characteristic-diagnostic sites at species level: 36 characteristic-diagnostic sites from eight species, in which 12 from Gmnodiptychus pachycgeilus, 2 from Platypharodon extremus, 7 from Ptychobarbus kaznakovi, 2 from Schizopygopsis anteroventris, 2 from Schizopygopsis malacanthus, 3 from Schizopygopsis malacanthus chengi, 3 from Schizothorax dolichonema and 5 from Schizothorax lantsangensis. Our results show that Schizothoracinae fishes can be identified validly by using COI DNA barcode. Thirty-six characteristic-diagnostic sites were proposed to be applied into works of species identification for the Schizothoracinae fishes in Qinghai Province.

Genetic diversity in frogs linked to past and future climate changes on the roof of the world

Mountains, representing storehouses of biodiversity, endemism and threatened species, are biodiversity hotspots of great conservation importance. However, increasing evidence indicates that mountain species throughout the world are responding to climate change, past or contemporary, by shifting their geographic distributions and patterns of genetic diversity, potentially affecting their adaptive capacity and increasing risk of extinction. Using the iconic high-elevation frog Nanorana parkeri as indicator, we showed how spatial analyses of climatic stability combined with genetic data allow unravelling amphibian responses to past and future climate changes on 'the roof of the world'-the Qinghai-Tibetan Plateau. We found that areas along the Yarlung Tsangpo Valley were climatically more stable relative to other regions, apparently serving as a large climatic refugium during Quaternary glaciations, but that these areas will likely be affected by future climate change. As populations closer to Quaternary refugia usually had higher genetic diversity, current genetic diversity can be explained in the largest part by distance to historically stable areas, outweighing other historical and contemporary factors. Along with the dynamics of suitable range, a fluctuating habitat fragmentation supported the pattern of historical changes in genetic diversity (N_e ) over time. Our results emphasize strong relationships between amphibian genetic diversity, past range dynamics and where to preserve suitable habitats in the face of future climate changes. More generally, our findings highlighted a central role of refugia during Quaternary climatic fluctuations, and how isolation from refugia may have modulated amphibian genetic diversity across the Qinghai-Tibetan Plateau.