Molecular systematics of pikas (genus Ochotona) inferred from mitochondrial DNA sequences

Phylogeography of the Plateau Pika (Ochotona curzoniae) in Response to the Uplift of the Qinghai-Tibet Plateau

The evolution and current distribution of species on the Qinghai-Tibet Plateau have been significantly impacted by historical occurrences, including the uplift of the plateau and the Quaternary climate upheaval. As a remnant species, the plateau pika (Ochotona curzoniae) is a great model for researching historical events. In this study, 302 samples from 42 sample sites were utilized to analyze the impact of historical events on the evolution and distribution pattern of plateau pikas. The genetic diversity, patterns of differentiation, and historical dynamics of the plateau pika were investigated using molecular markers that included four mitochondrial genes (COI, D-loop, Cytb, and 12S rRNA) and three nuclear genes (GHR, IRBP, and RAG1). The results showed that: (1) The genetic diversity of the plateau pika was high in the Tibetan Plateau (Hd = 0.9997, π = 0.01205), and the plateau pika evolved into five lineages that occupied different geographical areas, with lineage 1 (Group 1) in the south of the Yarlung Zangbo River, lineage 2 (Group 2) in the hinterland of the plateau, lineage 3 (Group 3) in the northeastern part of the plateau, lineage 4 (Group 4) in the Hengduan Mountains, and lineage 5 (Group 5) in the eastern part of the plateau. (2) The gene flow among the five lineages was low, and the differentiation level was high (Nm < 0.25; Fst > 0.25), indicating that the geographical barriers between the five lineages, such as the Yarlung Zangbo River, the Qaidam-Ghuong-Guide Basin, and the Lancang River, effectively promoted the population differentiation of the plateau pika. (3) The plateau pika first spread from the Hengduan Mountains to the entire Qinghai-Tibet Plateau and then conducted small-scale migration and dispersal in several refuges across the plateau in response to climate changes during the glacial and interglacial periods. (4) Except for Group 1 and Group 4, all the other populations exhibited a rapid expansion between 0.06 and 0.01 Mya, but the expansion was considerably delayed or halted by the effects of climate change during the last glacial maximum (0.02 Mya). Overall, the plateau pika on the Qinghai-Tibet Plateau exhibits high genetic diversity, and topographic obstacles, including mountains, valleys, and basins, created by the uplift of the plateau and climatic changes since the Quaternary period have played an important role in the differentiation and historical dynamics of the plateau pika population.

A highland-adaptation mutation of the Epas1 protein increases its stability and disrupts the circadian clock in the plateau pika

The Qinghai-Tibet Plateau (QTP) harbors hundreds of species well adapted to its extreme conditions, including its low-oxygen (hypoxic) atmosphere. Here, we show that the plateau pika-a keystone mammal of the QTP-lacks robust circadian rhythms. The major form of the plateau pika Epas1 protein includes a 24-residue insert caused by a point mutation at the 5' juncture site of Intron14 and is more stable than other mammalian orthologs. Biochemical studies reveal that an Epas1-Bmal1 complex with lower trans-activation activity occupies the E1/E2 motifs at the promoter of the core-clock gene Per2, thus explaining how an Epas1 mutation-selected in the hypoxic conditions of the QTP-disrupts the molecular clockwork. Importantly, experiments with hypoxic chambers show that mice expressing the plateau pika Epas1 ortholog in their suprachiasmatic nucleus have dysregulated central clocks, and pika Epas1 knockin mice reared in hypoxic conditions exhibit dramatically reduced heart damage compared with wild-type animals.

Metagenomic Comparisons between Soft and Hard Feces of Plateau Pikas (Ochotona curzoniae)

Simple Summary

Plateau pika produces hard and soft feces with different morphology, component and microbial structure. Hard feces had more abundant Firmicutes, while soft feces had more abundant Akkermansia. The differences of microbial communities between hard and soft feces were mainly driven by core microbomes. Soft feces had a comprehensive advances in predict functional pathways compared to hard feces, these strengthened functional pathways were closely associated with metabolism of energy, vitamins, and amino acid. Our study preliminarily explored the differences in microbial structure and function between hard and soft feces, provided a foundation for future systematic explorations of the cecotrophy.

Abstract

The division of hard and soft feces is an effective digestion strategy in the order Lagomorpha. Although previous studies have reported that hard and soft feces differ in morphology and component, the discrepancy in the microbiome remains unclear. This study explored the microbiomes of hard and soft feces in plateau pikas by sequencing the V3 and V4 regions of 16S rDNA. We found that hard feces harbored higher Firmicutes, while soft feces harbored higher Akkermansia. Increased rare bacterial taxa were observed in hard feces compared with soft feces. Moreover, hard and soft feces displayed a greater difference in terms of core operational taxonomy units (OTUs) compared to the total OTUs. The soft feces showed enhancements in all predicted Kyoto Encyclopedia of Genes and Genomes (KEGG) functions, indicating an advancing microbial metabolism compared to hard feces. The significantly upregulated pathways in soft feces were mainly enriched in metabolism of energy and carbohydrate, glycan biosynthesis, cofactors and vitamins, and amino acids—all of which are associated with increased contents of microbial proteins, vitamins, and short-chain fatty acids. Our study reports, for the first time, the differential microbiomes between hard and soft feces of pikas and provides direction for the future studies on cecotrophy.

Sympatric Yaks and Plateau Pikas Promote Microbial Diversity and Similarity by the Mutual Utilization of Gut Microbiota

Interactions between species provide the basis for understanding coexisting mechanisms. The plateau pika (Ochotona curzoniae) and the yak (Bos grunniens) are considered competitors because they have shared habitats and consumed similar food on the Qinghai–Tibetan Plateau for more than 1 million years. Interestingly, the population density of plateau pikas increases with yak population expansion and subsequent overgrazing. To reveal the underlying mechanism, we sequenced the fecal microbial 16S rDNA from both sympatric and allopatric pikas and yaks. Our results indicated that sympatry increased both gut microbial diversity and similarity between pikas and yaks. The abundance of Firmicutes, Proteobacteria, Cyanobacteria, and Tenericutes decreased, while that of Verrucomicrobia increased in sympatric pikas. As for sympatric yaks, Firmicutes, Bacteroidetes, and Spirochaetes significantly increased, while Cyanobacteria, Euryarchaeota, and Verrucomicrobia significantly decreased. In sympatry, plateau pikas acquired 2692 OTUs from yaks, and yaks obtained 453 OTUs from pikas. The predominant horizontally transmitted bacteria were Firmicutes, Bacteroidetes, Verrucomicrobia, and Proteobacteria. These bacteria enhanced the enrichment of pathways related to prebiotics and immunity for pikas, such as heparin sulfate, heparin, chitin disaccharide, chondroitin-sulfate-ABC, and chondroitin-AC degradation pathways. In yaks, the horizontally transmitted bacteria enhanced pathways related to hepatoprotection, xenobiotic biodegradation, and detoxification. Our results suggest that horizontal transmission is a process of selection, and pikas and yaks tend to develop reciprocity through the horizontal transmission of gut microbiota.

Lagomorpha as a Model Morphological System

Due to their global distribution, invasive history, and unique characteristics, European rabbits are recognizable almost anywhere on our planet. Although they are members of a much larger group of living and extinct mammals [Mammalia, Lagomorpha (rabbits, hares, and pikas)], the group is often characterized by several well-known genera (e.g., Oryctolagus, Sylvilagus, Lepus, and Ochotona). This representation does not capture the extraordinary diversity of behavior and form found throughout the order. Model organisms are commonly used as exemplars for biological research, but there are a limited number of model clades or lineages that have been used to study evolutionary morphology in a more explicitly comparative way. We present this review paper to show that lagomorphs are a strong system in which to study macro- and micro-scale patterns of morphological change within a clade that offers underappreciated levels of diversity. To this end, we offer a summary of the status of relevant aspects of lagomorph biology.

Environment and host species identity shape gut microbiota diversity in sympatric herbivorous mammals

The previous studies have reported that the mammalian gut microbiota is a physiological consequence; nonetheless, the factors influencing its composition and function remain unclear. In this study, to evaluate the contributions of the host and environment to the gut microbiota, we conducted a sequencing analysis of 16S rDNA and shotgun metagenomic DNA from plateau pikas and yaks, two sympatric herbivorous mammals, and further compared the sequences in summer and winter. The results revealed that both pikas and yaks harboured considerably more distinct communities between summer and winter. We detected the over-representation of Verrucomicrobia and Proteobacteria in pikas, and Archaea and Bacteroidetes in yaks. Firmicutes and Actinobacteria, associated with energy-efficient acquisition, significantly enriched in winter. The diversity of the microbial community was determined by the interactive effects between the host and season. Metagenomic analysis revealed that methane-metabolism-related pathway of yaks was significantly enriched in summer, while some pathogenic pathways were more abundant in pikas. Both pikas and yaks had a higher capacity for lipid degradation in winter. Pika and yak shared more OTUs when food shortage occurred in winter, and this caused a convergence in gut microbial composition and function. From winter to summer, the network module number increased from one to five in pikas, which was different in yaks. Our study demonstrates that the host is a dominant factor in shaping the microbial communities and that seasonality promotes divergence or convergence based on dietary quality across host species identity.

Out of Tibet: Genomic Perspectives on the Evolutionary History of Extant Pikas

Pikas are widely distributed in the Northern Hemisphere and are highly adapted to cold and alpine environments. They are one of the most complex and problematic groups in mammalian systematics, and the origin and evolutionary history of extant pikas remain controversial. In this study, we sequenced the whole coding sequences of 105 pika samples (29 named species and 1 putative new species) and obtained DNA data for more than 10,000 genes. Our phylogenomic analyses recognized four subgenera of extant pikas: Alienauroa, Conothoa, Ochotona, and Pika. The interrelationships between the four subgenera were strongly resolved as (Conothoa, (Alienauroa, (Ochotona, Pika))), with the mountain group Conothoa being the sister group of all other pikas. Our divergence time and phylogeographic analyses indicated that the last common ancestor of extant pikas first occurred on in the middle Miocene, ∼14 Ma. The emergence of opportunities related to the climate, food supply, and spreading paths in concert promoted the dispersal of pikas from the Qinghai-Tibetan Plateau (QTP) to other parts of Eurasia and North America. We found that the genes that were positively selected in the early evolution of pikas were most concentrated in functional categories related to cold tolerance. These results suggest that the QTP may have served as a training ground for cold tolerance in early pikas, which gives pikas a great advantage when the climate continued to cool after the middle Miocene. Our study highlights the importance of the QTP as a center of origin for many cold-adapted animals.

Phylogeography of the Tibetan hamster Cricetulus kamensis in response to uplift and environmental change in the Qinghai-Tibet Plateau

AIM

The evolutionary process of an organism provides valuable data toward an understanding of the Earth evolution history. To investigate the relationship between the uplift of the Qinghai-Tibet Plateau (QTP) and mammalian evolution since the late Cenozoic, the geographic distribution of genetic variations in the Tibetan hamster Cricetulus kamensis was investigated using phylogeographical methods. In particular, population divergence, demographic history, genetic variation, and the prediction of species distribution area were investigated.

LOCATION

The Qinghai-Tibet Plateau.

METHODS

A total of 53 specimens, representing 13 geographic populations, were collected from the QTP. The phylogeographical pattern and demographic history of C. kamensis were analyzed, and the probable factors in the QTP uplift and the Quaternary glacial periods were inferred from one nuclear and four mitochondrial genes. Furthermore, the species distribution model (SDM) was used to predict changes in potentially suitable habitats since the last Interglacial.

RESULTS

Phylogenetic analysis demonstrated that two major genetic differentiations of the C. kamensis population occurred during the Early Pleistocene that were influenced by the Qing-Zang tectonic movement from the Middle Pliocene to the Early Pleistocene. Genetic distance between two major clades indicated low genetic divergence. Demographic history analysis showed that the C. kamensis population was affected by the Quaternary glacial period. SDM analysis indicated that C. kamensis was endemic to the QTP and the suitable habitat was affected by climate change, especially during the Last Glacial Maximum (LGM).

MAIN CONCLUSION

Our results indicated that the QTP uplift led to the population divergence of C. kamensis, and vicariance well accounted for the geographic distribution of genetic variation in C. kamensis as a result of genetic divergence and lack of gene flow. The genetic distance shows that C. alticola may be a subspecies of C. kamensis. Demographic history analysis suggests that the QTP was affected by the last glacial period. SDM analysis supports that almost the entire QTP is covered by a huge ice sheet during the LGM.

Independence between coping style and stress reactivity in plateau pika

The concept of coping style represents the way individual animals react to a stressful situation, both behaviourally and neurophysiologically. Over the last decades coping style has been linked to the development of research on animal personality. Based on this concept, we should find a proactive-reactive continuum in animal populations, with proactive individuals being fast explorer, bold, aggressive, and show high sympathetic reactivity (higher heart rate), as well as low hypothalamus-pituitary-adrenocortical (PHA) axis reactivity to external stressor (higher plasma glucocorticoid level). At the other extreme, shy, lowly aggressive, reactive individuals should be slow in their exploration, and show a low sympathetic reactivity and a high HPA axis reactivity. However, a recent two-tier model proposed that coping style and stress reactivity should be independent of each other. In this study, we tested the two-tier model in a wild plateau pika (Ochotona curzoniae) population on the Tibetan Plateau, by quantifying the associations between several behavioural and physiological traits at the among- and within-individual levels. We repeatedly measured exploration, docility, boldness, heart rate and plasma cortisol concentration in individuals between April and September of 2013. All traits tested were repeatable. At the among-individual level, all behavioural traits were correlated with each other and with heart rate, but were independent of both basal level and variation of plasma cortisol concentration. Most correlations were negligible at the within-individual level. In support of the two-tier model, these results suggest that coping style (i.e. behaviour and heart rate associations) is independent of stress reactivity (i.e. glucocorticoid reactivity) in that species.

The complete mitochondrial genome of the Caulerpa lentillifera (Ulvophyceae, Chlorophyta): Sequence, genome content, organization structure and phylogenetic consideration

The complete mitochondrial genome is greatly important for studies on genetic structure and phylogenetic relationship at various taxonomic levels. To obtain information about the evolutionary trends of mtDNA in the Ulvophyceae and also to gain insights into the phylogenetic relationships between ulvophytes and other chlorophytes, we determined the mtDNA sequence of Caulerpa lentillifera (sea grape) using de novo mitochondrial genome sequencing. The complete genomic DNA of C. lentillifera was circular and 209,034 bp in length, and it was the largest green-algal mitochondrial genome sequenced to date, with a low gene density of 65.2%, which is reminiscent of the "expanded" pattern of evolution exhibited by embryophyte mtDNAs. The C. lentillifera mtDNA consisted of a typical set of 17 protein-coding genes (PCGs), 20 transfer RNA (tRNA) genes, three ribosomal RNA (rRNA) genes, 42 putative open reading frames (ORFs) and 29 introns, which had homologs in green-algal mtDNAs displaying an "ancestral" or a "reduced-derived" pattern of evolution. The overall base composition of its mitochondrial genome was 24.19% for A, 24.94% for T, 25.80% for G, 25.07% for C and 50.87% for GC. The mitochondrial genome of C. lentillifera was characterized by numerous small intergenic regions and introns, which was clearly different from other green algae. With the exception of the NADH dehydrogenase subunit 6 (ND6), ND1, ATP and three tRNA genes (tRNA-His, tRNA-Thr and tRNA-Ala), all other mitochondrial genes were encoded on the heavy strand. All of the PCGs had ATG as their start codon and employed TAA, TGA or TAG as their termination codon. To gain insights into the evolutionary trends of mtDNA in the Ulvophyceae, we inferred the complete mtDNA sequence of C. lentillifera, an ulvophyte belonging to a distinct, early-diverging lineage. Taken together, our data offered useful information for the studies on phylogenetic hypotheses and phylogenetic relationships of C. lentillifera within the Chlorophyta.

The geography and timing of genetic divergence in the lizard Phrynocephalus theobaldi on the Qinghai-Tibetan plateau

The Qinghai-Tibetan Plateau (QTP) represents one of the earth’s most significant physical features and there is increasing interest in the historical generation of biodiversity within this region. We hypothesized that there should be clear geographically coherent genetic structuring within one of the world’s highest altitude lizards, Phrynocephalus theobaldi, due to considerable historical population fragmentation in this environment. This was tested using a major mitochondrial DNA (mtDNA) survey and sequencing of two nuclear markers (AME and RAG-1) from P. theobaldi, from across the southern QTP. A Bayesian method (BPEC) was used to detect four geographically structured mtDNA clusters. A Bayesian phylogenetic tree, together with associated dating analyses, supported four corresponding evolutionary lineages with a timing of 3.74–7.03 Ma for the most basal P. theobaldi split and Pliocene splits of 2.97–5.79 Ma and 2.40–5.39 Ma in the two daughter lineages. Himalayan uplift and changes in the Jilong basin may have contributed to these divergences, but uplift of the Gangdese mountains is rejected due to its timing. The nuclear markers appeared to be sorted between the four mtDNA groups, and species delimitation analyses supported the four phylogeographical groups as candidate species. The study contributes to our understanding of biodiversity on the QTP.

Multilocus approaches reveal underestimated species diversity and inter-specific gene flow in pikas (Ochotona) from southwestern China

The phylogeny of living pikas (Ochotonidae, Ochotona) remains obscure, and pika species diversity in southwestern China has never been well explored. In this study, 96 tissue samples from 11 valid species in three classified subgenera (Pika, Ochotona and Conothoa) from 23 locations were characterized using multilocus sequences of 7031bp. Two mitochondrial (CYT B and COI) and five nuclear gene segments (RAG1, RAG2, TTN, OXAIL and IL1RAPL1) were sequenced. We analysed evolutionary histories using maximum likelihood (RAxML) and Bayesian analyses (BEAST), and we also used molecular species delimitation analyses (BPP) to explore species diversity. Our study supported O. syrinx (O. huangensis) as a distinct clade from all named subgenera. Relationships among subgenera were not fully resolved, which may be due to a rapid diversification in the middle Miocene (∼13.90Ma). Conflicting gene trees implied mitochondrial introgression from O. cansus to O. curzoniae. We uncovered three cryptic species from Shaanxi, Sichuan and Yunnan with strong support, suggesting an underestimation of species diversity in the "sky-island" mountains of southwest China.

Opposite companion effect on flight initiation distance in sympatric species: plateau pika (Ochotona curzoniae) and White-rumped Snowfinch (Onychostruthus taczanowskii)

There are many factors influencing prey’s risk perception and escape decision during predator–prey encounters. The distance at which animals move away from perceived danger (often quantified as flight initiation distance or FID) has been used by behavioral ecologists to understand the economics of antipredator behavior. Using general linear models, we investigated escape decision-making processes in plateau pika (Ochotona curzoniae (Hodgson, 1858)) and White-rumped Snowfinch (Onychostruthus taczanowskii (Prjevalsky, 1876), formerly known as Montifringilla taczanowskii Przewalski, 1876) together and we found that (i) there are significant positive correlations between starting distance and FID in both species; (ii) pika escapes at a longer distance from an approaching intruder when it is far from its burrow; (iii) foraging animals tolerate closer distances than watching ones, both in pikas and in Snowfinches; (iv) conspecifics seem to have no effect on pikas’ escaping behavior, while the appearance of Snowfinches dramatically decreases the FID of pikas. On the contrary, conspecifics significantly decrease the FID of Snowfinches, while the presence of pikas has no effects. These findings provide new evidence, which are consistent with optimal escape theory. These novel results in multispecies interactions prompt us to pay attention to the potential relationship between pikas and Snowfinches, as well as the “bird–pika in one hole” phenomenon.

Comprehensive Transcriptome Analysis of Six Catfish Species from an Altitude Gradient Reveals Adaptive Evolution in Tibetan Fishes

Glyptosternoid fishes (Siluriformes), one of the three broad fish lineages (the two other are schizothoracines and Triplophysa), have a limited distribution in the rivers in the Tibetan Plateau and peripheral regions. To investigate the genetic mechanisms underlying adaptation to the Tibetan Plateau in several fish species from gradient altitudes, a total of 20,659,183–37,166,756 sequence reads from six species of catfish were generated by Illumina sequencing, resulting in six assemblies. Analysis of the 1,656 orthologs among the six assembled catfish unigene sets provided consistent evidence for genome-wide accelerated evolution in the three glyptosternoid lineages living at high altitudes. A large number of genes refer to functional categories related to hypoxia and energy metabolism exhibited rapid evolution in the glyptosternoid lineages relative to yellowhead catfish living in plains areas. Genes showing signatures of rapid evolution and positive selection in the glyptosternoid lineages were also enriched in functions associated with energy metabolism and hypoxia. Our analyses provide novel insights into highland adaptation in fishes and can serve as a foundation for future studies aiming to identify candidate genes underlying the genetic basis of adaptation in Tibetan fishes.

Biogeographic history and high-elevation adaptations inferred from the mitochondrial genome of Glyptosternoid fishes (Sisoridae, Siluriformes) from the southeastern Tibetan Plateau

Background

The distribution of the Chinese Glyptosternoid catfish is limited to the rivers of the Tibetan Plateau and peripheral regions, especially the drainage areas of southeastern Tibet. Therefore, Glyptosternoid fishes are ideal for reconstructing the geological history of the southeastern Tibet drainage patterns and mitochondrial genetic adaptions to high elevations.

Results

Our phylogenetic results support the monophyly of the Sisoridae and the Glyptosternoid fishes. The reconstructed ancestral geographical distribution suggests that the ancestral Glyptosternoids was widely distributed throughout the Brahmaputra drainage in the eastern Himalayas and Tibetan area during the Late Miocene (c. 5.5 Ma). We found that the Glyptosternoid fishes lineage had a higher ratio of nonsynonymous to synonymous substitutions than those found in non-Glyptosternoids. In addition, ω_pss was estimated to be 10.73, which is significantly higher than 1 (p-value 0.0002), in COX1, which indicates positive selection in the common ancestral branch of Glyptosternoid fishes in China. We also found other signatures of positive selection in the branch of specialized species. These results imply mitochondrial genetic adaptation to high elevations in the Glyptosternoids.

Conclusions

We reconstructed a possible scenario for the southeastern Tibetan drainage patterns based on the adaptive geographical distribution of the Chinese Glyptosternoids in this drainage. The Glyptosternoids may have experienced accelerated evolutionary rates in mitochondrial genes that were driven by positive selection to better adapt to the high-elevation environment of the Tibetan Plateau.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-015-0516-9) contains supplementary material, which is available to authorized users.

The phylogeny of pikas (Ochotona) inferred from a multilocus coalescent approach

The clarification of the systematics of pikas (genus Ochotona) has been hindered by largely overlapping morphological characters among species and the lack of a comprehensive molecular phylogeny. Here we estimate the first multilocus phylogeny of the genus to date, by analysing 12 nuclear DNA markers (total of 7.5Kb) in 11 species of pikas from the four classified subgenera (Pika, Ochotona, Lagotona and Conothoa) using a multispecies coalescent-based framework. The species-tree confirmed the subgeneric classification by retrieving as monophyletic the subgenera represented here by more than one species. Contrary to previous phylogenies based on mtDNA alone, Lagotona was found to be sister to Pika. Also, support for the monophyly of the alpina group was not strong, thus caution should be used in future analyses of this group. A relaxed molecular clock calibrated using the Ochotonidae-Leporidae divergence resulted in more recent estimates of divergence times relative to previous studies. Strong concordance with inferences based on fossil records was found, suggesting that the initial diversification of the genus took place by the end of late Miocene. Finally, this work sets up methodologies and gathers molecular markers that can be used to extend the understanding of the evolutionary history of the genus.

Multilocus phylogeny and cryptic diversity in Asian shrew-like moles (Uropsilus, Talpidae): implications for taxonomy and conservation

BACKGROUND

The genus Uropsilus comprises a group of terrestrial, montane mammals endemic to the Hengduan and adjacent mountains. These animals are the most primitive living talpids. The taxonomy has been primarily based on cursory morphological comparisons and the evolutionary affinities are little known. To provide insight into the systematics of this group, we estimated the first multi-locus phylogeny and conducted species delimitation, including taxon sampling throughout their distribution range.

RESULTS

We obtained two mitochondrial genes (~1, 985 bp) and eight nuclear genes (~4, 345 bp) from 56 specimens. Ten distinct evolutionary lineages were recovered from the three recognized species, eight of which were recognized as species/putative species. Five of these putative species were found to be masquerading as the gracile shrew mole. The divergence time estimation results indicated that climate change since the last Miocene and the uplift of the Himalayas may have resulted in the diversification and speciation of Uropsilus.

CONCLUSIONS

The cryptic diversity found in this study indicated that the number of species is strongly underestimated under the current taxonomy. Two synonyms of gracilis (atronates and nivatus) should be given full species status, and the taxonomic status of another three potential species should be evaluated using extensive taxon sampling, comprehensive morphological, and morphometric approaches. Consequently, the conservation status of Uropsilus spp. should also be re-evaluated, as most of the species/potential species have very limited distribution.

An intriguing shift occurs in the novel protein phosphatase 1 binding partner, TCTEX1D4: evidence of positive selection in a pika model

T-complex testis expressed protein 1 domain containing 4 (TCTEX1D4) contains the canonical phosphoprotein phosphatase 1 (PPP1) binding motif, composed by the amino acid sequence RVSF. We identified and validated the binding of TCTEX1D4 to PPP1 and demonstrated that indeed this protein is a novel PPP1 interacting protein. Analyses of twenty-one mammalian species available in public databases and seven Lagomorpha sequences obtained in this work showed that the PPP1 binding motif ₉₀RVSF₉₃ is present in all of them and is flanked by a palindromic sequence, PLGS, except in three species of pikas (Ochotona princeps, O. dauurica and O. pusilla). Furthermore, for the Ochotona species an extra glycosylation site, motif ₉₆NLS₉₈, and the loss of the palindromic sequence were observed. Comparison with other lagomorphs suggests that this event happened before the Ochotona radiation. The d_N/d_S for the sequence region comprising the PPP1 binding motif and the flanking palindrome highly supports the hypothesis that for Ochotona species this region has been evolving under positive selection. In addition, mutational screening shows that the ability of pikas TCTEX1D4 to bind to PPP1 is maintained, although the PPP1 binding motif is disrupted, and the N- and C-terminal surrounding residues are also abrogated. These observations suggest pika as an ideal model to study novel PPP1 complexes regulatory mechanisms.

Evolutionary history of lagomorphs in response to global environmental change

Although species within Lagomorpha are derived from a common ancestor, the distribution range and body size of its two extant groups, ochotonids and leporids, are quite differentiated. It is unclear what has driven their disparate evolutionary history. In this study, we compile and update all fossil records of Lagomorpha for the first time, to trace the evolutionary processes and infer their evolutionary history using mitochondrial genes, body length and distribution of extant species. We also compare the forage selection of extant species, which offers an insight into their future prospects. The earliest lagomorphs originated in Asia and later diversified in different continents. Within ochotonids, more than 20 genera occupied the period from the early Miocene to middle Miocene, whereas most of them became extinct during the transition from the Miocene to Pliocene. The peak diversity of the leporids occurred during the Miocene to Pliocene transition, while their diversity dramatically decreased in the late Quaternary. Mantel tests identified a positive correlation between body length and phylogenetic distance of lagomorphs. The body length of extant ochotonids shows a normal distribution, while the body length of extant leporids displays a non-normal pattern. We also find that the forage selection of extant pikas features a strong preference for C₃ plants, while for the diet of leporids, more than 16% of plant species are identified as C₄ (31% species are from Poaceae). The ability of several leporid species to consume C₄ plants is likely to result in their size increase and range expansion, most notably in Lepus. Expansion of C₄ plants in the late Miocene, the so-called ‘nature’s green revolution’, induced by global environmental change, is suggested to be one of the major ‘ecological opportunities’, which probably drove large-scale extinction and range contraction of ochotonids, but inversely promoted diversification and range expansion of leporids.

The chromosomes of Tsing-Ling pika, Ochotona huangensis Matschie, 1908 (Lagomorpha, Ochotonidae)

The karyotype of the Tsing-Ling (Huanghe) pika, Ochotona huangensis Matschie, 1908 from the forest habitats of the Qinling Mountains (Shaanxi Province, China) was described for the first time. The chromosome set contains 42 chromosomes (NFa=80). The autosomes are 15 meta-submetacentric pairs and 5 subtelocentric pairs. The X chromosome is a medium sized submetacentric; the Y chromosome is a small sized acrocentric. C-banding revealed a localization of heterochromatin in the pericentromeric regions of all autosomes.

Low, complex and probably reticulated chromosome evolution of Sciuromorpha (Rodentia) and Lagomorpha

Lagomorpha (rabbits and pikas) and Sciuromorpha (squirrels) are grouped in the Glires superorder. Their chromosome diversification, since their separation from the eutherian mammalian common ancestor, was characterized by a low rate of chromosome rearrangements. Consequently, the structure of some chromosomes was either conserved or only slightly modified, making their comparison easy at the genus, family and even order level. Interspecific in situ hybridization (Zoo-FISH) largely corroborates classical cytogenetic data but provides much more reliability in comparisons, especially for distant species. We reconstructed common ancestral karyotypes for Glires, Lagomorpha, Sciuromorpha, and Sciuridae species, and then, determined the chromosome changes separating these ancestors from their common eutherian ancestor. We propose that reticulated evolution occurred during the diversification of Glires, which implies that several pericentric inversions and Robertsonian translocations were conserved in the heterozygous status for an extensive period. Finally, among Lagomorpha and Sciuromorpha, we focused on Leporidae and Sciuridae chromosome evolution. In the various attempts to establish dichotomic evolutionary schemes, it was necessary to admit that multiple homoplasies (convergent and reverse rearrangements) occurred in Sciuridae and in a lesser degree, in Leporidae. In Leporidae, additional rearrangements were sufficient to propose a resolved phylogeny. However, a resolved phylogeny was not possible for Sciuridae because most of the rearrangements occurred in terminal branches. We conclude that a reticulated evolution took place early during the evolution of both families and lasted longer in Sciuridae than in Leporidae. In Sciuridae, most chromosome rearrangements were pericentric inversions involving short fragments. Such rearrangements have only mild meiotic consequences, which may explain the long persistence of the heterozygous status characterizing reticulated evolution.

Northern host-parasite assemblages: history and biogeography on the borderlands of episodic climate and environmental transition

Diversity among assemblages of mammalian hosts and parasites in northern terrestrial ecosystems was structured by a deep history of biotic and abiotic change that overlies a complex geographic arena. Since the Pliocene, Holarctic ecosystems assembled in response to shifting climates (glacial and interglacial stages). Cycles of episodic dispersal/isolation and diversification defined northern diversity on landscape to regional scales. Episodes of geographic expansion and colonisation linked Eurasia and North America across Beringia and drove macroevolutionary structure of host and parasite associations. Asynchronous dispersal from centres of origin in Eurasia into the Nearctic resulted in gradients in parasite diversity in the carnivoran, lagomorph, rodent and artiodactyl assemblages we reviewed. Recurrent faunal interchange and isolation in conjunction with episodes of host colonisation have produced a mosaic structure for parasite faunas and considerable cryptic diversity among nematodes and cestodes. Mechanisms of invasion and geographic colonisation leading to the establishment of complex faunal assemblages are equivalent in evolutionary and ecological time, as demonstrated by various explorations of diversity in these high-latitude systems. Our ability to determine historical responses to episodic shifts in global climate may provide a framework for predicting the cascading effects of contemporary environmental change.

When cold is better: climate-driven elevation shifts yield complex patterns of diversification and demography in an alpine specialist (American pika, Ochotona princeps)

The genetic consequences of climate-driven range fluctuation during the Pleistocene have been well studied for temperate species, but cold-adapted (e.g., alpine, arctic) species that may have responded uniquely to past climatic events have received less attention. In particular, we have no a priori expectation for long-term evolutionary consequences of elevation shifts into and out of sky islands by species adapted to alpine habitats. Here, we examined the influence of elevation shifts on genetic differentiation and historical demography in an alpine specialist, the American pika (Ochotona princeps). Pika populations are divided into five genetic lineages that evolved in association with separate mountain systems, rather than lineages that reflect individual sky islands. This suggests a role for glacial-period elevation shifts in promoting gene flow among high-elevation populations and maintaining regional cohesion of genetic lineages. We detected a signature of recent demographic decline in all lineages, consistent with the expectation that Holocene climate warming has driven range retraction in southern lineages, but unexpected for northern populations that presumably represent postglacial expansion. An ecological niche model of past and future pika distributions highlights the influence of climate on species range and indicates that the distribution of genetic diversity may change dramatically with continued climate warming.

Identification of a RELIK orthologue in the European hare (Lepus europaeus) reveals a minimum age of 12 million years for the lagomorph lentiviruses

The retroviral genus Lentivirus comprises retroviruses characterised from five mammalian orders. Lentiviruses typically undergo rapid rates of evolution, a feature that has allowed recent evolutionary relationships to be elucidated, but has also obscured their distant evolutionary past. However, the slowdown in the rate of evolution associated with genome invasion, as has occurred in the European rabbit, enables longer-term lentiviral evolutionary history to be inferred. Here we report the identification of orthologous RELIK proviruses in the European hare, demonstrating a minimum age of 12 million years for the lagomorph lentiviruses. This finding indicates an association between lentiviruses and their hosts covering much of the evolutionary history of the lagomorphs, and taking place within species with a worldwide distribution.

Natural selection and adaptive evolution of leptin in the ochotona family driven by the cold environmental stress

BACKGROUND

Environmental stress can accelerate the evolutionary rate of specific stress-response proteins and create new functions specialized for different environments, enhancing an organism's fitness to stressful environments. Pikas (order Lagomorpha), endemic, non-hibernating mammals in the modern Holarctic Region, live in cold regions at either high altitudes or high latitudes and have a maximum distribution of species diversification confined to the Qinghai-Tibet Plateau. Variations in energy metabolism are remarkable for them living in cold environments. Leptin, an adipocyte-derived hormone, plays important roles in energy homeostasis.

METHODOLOGY/PRINCIPAL FINDINGS

To examine the extent of leptin variations within the Ochotona family, we cloned the entire coding sequence of pika leptin from 6 species in two regions (Qinghai-Tibet Plateau and Inner Mongolia steppe in China) and the leptin sequences of plateau pikas (O. curzonia) from different altitudes on Qinghai-Tibet Plateau. We carried out both DNA and amino acid sequence analyses in molecular evolution and compared modeled spatial structures. Our results show that positive selection (PS) acts on pika leptin, while nine PS sites located within the functionally significant segment 85-119 of leptin and one unique motif appeared only in pika lineages-the ATP synthase alpha and beta subunit signature site. To reveal the environmental factors affecting sequence evolution of pika leptin, relative rate test was performed in pikas from different altitudes. Stepwise multiple regression shows that temperature is significantly and negatively correlated with the rates of non-synonymous substitution (Ka) and amino acid substitution (Aa), whereas altitude does not significantly affect synonymous substitution (Ks), Ka and Aa.

CONCLUSIONS/SIGNIFICANCE

Our findings support the viewpoint that adaptive evolution may occur in pika leptin, which may play important roles in pikas' ecological adaptation to extreme environmental stress. We speculate that cold, and probably not hypoxia, may be the primary environmental factor for driving adaptive evolution of pika leptin.

Genetic characterization of the chemokine receptor CXCR4 gene in lagomorphs: comparison between the families Ochotonidae and Leporidae

Chemokines receptors are transmembrane proteins that bind chemokines. Chemokines and their receptors are known to play a crucial role in the immune system and in pathogen entry. There is evidence that myxoma virus, the causative agent of myxomatosis, can use the chemokine receptor CXCR4 to infect cells. This virus causes a benign disease in its natural host, Sylvilagus, but in the European rabbit (Oryctolagus cuniculus) it causes a highly fatal and infectious disease known as myxomatosis. We have characterized the chemokine receptor CXCR4 gene in five genera of the order Lagomorpha, Ochotona (Ochotonidae), and Oryctolagus, Lepus, Bunolagus and Sylvilagus (Leporidae). In lagomorphs, the CXCR4 is highly conserved, with most of the protein diversity found at surface regions. Five amino acid replacements were observed, two in the intracellular loops, one in the transmembrane domain and two in the extracellular loops. Oryctolagus features unique amino acid changes at the intracellular domains, putting this genus apart of all other lagomorphs. Furthermore, in the 37 European rabbits analysed, which included healthy rabbits and rabbits with clinical symptoms of myxomatosis, 14 nucleotide substitutions were obtained but no amino acid differences were observed.