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

Evolutionary divergence on the Qinghai-Tibet Plateau: How life-history traits shape the diversity of plateau zokor and pika populations

Understanding how species diverge and adapt is fundamental to unraveling biodiversity. While environmental impacts on species evolution are well-documented, the roles of intrinsic life-history traits remain underexplored. The Qinghai-Tibet Plateau, with its harsh conditions and unique biodiversity, offers a natural laboratory for such investigations. Here, we examined two sympatric small mammals-the solitary, low-dispersal plateau zokor (Eospalax baileyi) and the social, high-dispersal plateau pika (Ochotona curzoniae)-to elucidate how life-history traits shape population structures and adaptive strategies. Through whole-genome sequencing and cardiac-blood phenotype analyses, we reveal striking differences in their evolutionary trajectories. Despite enduring similar environmental pressures, plateau zokor populations exhibit pronounced genetic subdivisions, high inbreeding, and distinct local adaptations. In contrast, plateau pika populations display genetic panmixia, widespread diversity, and adaptive uniformity. Demographic inference highlights plateau zokors experienced severe population bottlenecks and restricted gene flow during glacial periods, underscoring the impact of dispersal capacity on evolutionary outcomes. Our findings demonstrate that intrinsic biological traits, particularly dispersal ability, fundamentally influence genetic architecture, population connectivity, and local adaptation. This study not only provides empirical evidence of how life-history traits shape evolutionary dynamics but offers a framework for integrating intrinsic and extrinsic factors in understanding biodiversity formation.

Genomic Vulnerability of a Sentinel Mammal Under Climate Change

Climate change poses a significant threat to biodiversity, particularly in alpine ecosystems where species have already undergone elevational range shifts. Genomics can be used to estimate the adaptive potential of species, as well as the shift in adaptive genomic composition necessary for populations to adjust to climate change (e.g., genomic offset). Here, we investigated patterns of climate-mediated adaptive genetic variation and predicted the degree of genomic offset under multiple climate change scenarios for a sentinel alpine mammal, the American pika (Ochotona princeps). We collected genome-wide data (29,709 SNPs) from 363 individuals spanning the entire range in western North America and employed genotype-environment association analyses to identify 924 robust outlier SNPs, several of which were linked to genes previously associated with high elevation and hypoxia responses in various pika species (Ochotonidae). Adaptive genomic variation was most strongly influenced by mean warmest month temperature, followed by precipitation of the coldest quarter. Spatial patterns of genomic offset were heterogeneous, significantly predicted by levels of adaptive genetic variation, elevation and latitude. Sites within the Northern Rocky Mountains exhibited the highest genomic offset under projected climate change despite possessing high levels of adaptive genetic variation. As such, while our study provides an example of how genomic data can be used to explore the potential consequences of climate change, it further highlights the need for careful consideration of genomic offset values within their proper ecological context.

Molecular time estimates for the Lagomorpha diversification

Despite their importance as members of the Glires group, lagomorph diversification processes have seldom been studied using molecular data. Notably, only a few phylogenetic studies have included most of the examined lagomorph lineages. Previous studies that included a larger sample of taxa and markers used nonconservative tests to support the branches of their proposed phylogeny. The objective of this study was to test the monophyly of families and genera of lagomorphs and to evaluate the group diversification process. To that end, this work expanded the sampling of markers and taxa in addition to implementing the bootstrap, a more rigorous statistical test to measure branch support; hence, a more robust phylogeny was recovered. Our supermatrix included five mitochondrial genes and 14 nuclear genes for eighty-eight taxa, including three rodent outgroups. Our maximum likelihood tree showed that all tested genera and both families, Leporidae and Ochotonidae, were recovered as monophyletic. In the Ochotona genus, the subgenera Conothoa and Pika, but not Ochotona, were recovered as monophyletic. Six calibration points based on fossils were used to construct a time tree. A calibration test was performed (via jackknife) by removing one calibration at a time and estimating divergence times for each set. The diversification of the main groups of lagomorphs indicated that the origin of the order's crown group was dated from the beginning of the Palaeogene. Our diversification time estimates for Lagomorpha were compared with those for the largest mammalian order, i.e., rodent lineages in Muroidea. According to our time-resolved phylogenetic tree, the leporids underwent major radiation by evolving a completely new morphospace-larger bodies and an efficient locomotor system-that enabled them to cover wide foraging areas and outrun predators more easily than rodents and pikas.

Unique microbial communities of parasitic fleas on wild animals from the Qinghai-Tibet Plateau

Fleas, one of the most significant ectoparasites, play a crucial role as vectors in spreading zoonotic diseases globally. The Qinghai Province, as part of the Qinghai-Tibet Plateau, is one of the provinces in China with the largest number of flea species. In this study, we characterized the microbial communities of eighty-five adult fleas, belonging to nineteen species within four families (Ceratophyllidae, Ctenophthalmidae, Leptopsyllidae, and Pulicidae). We identified a total of 1162 unique operational taxonomic units at the genus level, with flea-borne pathogens such as Wolbachia, Bartonella, Rickettsia being the members of top abundant taxa. Except for comparison between Ctenophthalmidae and Leptopsyllidae families, the analyses of both alpha- and beta- diversity indicators suggested that bacterial diversity varied among flea families. This could be attributed to flea phylogeny, which also influenced by their geographical sites and animal hosts. Results of Linear discriminant analysis effect size (LEfSe) indicated that 29 genera in Ceratophylloidea, 11 genera in Ctenophthalmidae, 15 genera in Leptopsyllidae, and 22 genera in Pulicidae were significantly responsible for explaining the differences among the four flea families (linear discriminant analysis score > 2, P < 0.05). Phylogenetic Investigation of Communities by Reconstruction of Unobserved States (PICRUSt2) analyses showed that the functional pathways varied significantly across flea families, which was supported by the significant correlation between the functional pathways and the microbial communities.

Phylogenomics of American pika (Ochotona princeps) lineage diversification

Quaternary climate oscillations have profoundly influenced current species distributions. For many montane species, these fluctuations were a prominent driver in species range shifts, often resulting in intraspecific diversification, as has been the case for American pikas (Ochotona princeps). Range shifts and population declines in this thermally-sensitive lagomorph have been linked to historical and contemporary environmental changes across its western North American range, with previous research reconstructing five mitochondrial DNA lineages. Here, we paired genome-wide data (25,244 SNPs) with range-wide sampling to re-examine the number and distribution of intra-specific lineages, and investigate patterns of within- and among-lineage divergence and diversity. Our results provide genomic evidence of O. princeps monophyly, reconstructing six distinct lineages that underwent multiple rounds of divergence (0.809-2.81 mya), including a new Central Rocky Mountain lineage. We further found evidence for population differentiation across multiple spatial scales, and reconstructed levels of standing variation comparable to those found in other small mammals. Overall, our findings demonstrate the influence of past glacial cycles on O. princeps lineage diversification, suggest that current subspecific taxonomy may need to be revisited, and provide an important framework for investigations of American pika adaptive potential in the face of anthropogenic climate change.

GeneMiner: A tool for extracting phylogenetic markers from next-generation sequencing data

The advancement of next-generation sequencing (NGS) technologies has been revolutionary for the field of evolutionary biology. This technology has led to an abundance of available genomes and transcriptomes for researchers to mine. Specifically, researchers can mine for various types of molecular markers that are vital for phylogenetic, evolutionary and ecological studies. Numerous tools have been developed to extract these molecular markers from NGS data. However, due to an insufficient number of well-annotated reference genomes for non-model organisms, it remains challenging to obtain these markers accurately and efficiently. Here, we present GeneMiner, an improved and expanded version of our previous tool, Easy353. GeneMiner combines the reference-guided de Bruijn graph assembly with seed self-discovery and greedy extension. Additionally, it includes a verification step using a parameter-bootstrap method to reduce the pitfalls associated with using a relatively distant reference. Our results, using both experimental and simulation data, showed GeneMiner can accurately acquire phylogenetic molecular markers for plants using transcriptomic, genomic and other NGS data. GeneMiner is designed to be user-friendly, fast and memory-efficient. Further, it is compatible with Linux, Windows and macOS. All source codes are publicly available on GitHub (https://github.com/sculab/GeneMiner) and Gitee (https://gitee.com/sculab/GeneMiner) for easy accessibility and transparency.

The first complete mitochondrial genome data of the Afghan pika Ochotona rufescens (Lagomorpha, Ochotonidae), near the type locality

The Afghan pika Ochotona rufescens (Gray, 1842) is widely distributed across the mountains of Afghanistan, Iran, Pakistan, and southwestern Turkmenistan, most often at elevations between 2,000 and 3,000 m. Here we present, for the first time, the complete mitochondrial genomes of two specimens of the nominotypical subspecies Ochotona rufescens rufescens, de novo assembled from Illumina short reads of fragmented probe-enriched DNA. The lengths of the circular mitogenomes are 16,408 bp and 16,407 bp, respectively. Both mitogenomes contain 13 protein-coding genes (PCGs), two ribosomal RNAs (16S rRNA and 12S rRNA), 22 transfer RNA genes, and a control region. The gene NAD6 and the tRNA (Gln), tRNA (Ala), tRNA (Asn), tRNA (Cys), tRNA (Tyr), tRNA (Ser), tRNA (Glu), and tRNA (Pro) are encoded on the light strand while the rest are encoded on the heavy strand. The overall nucleotide composition was ∼30% for A, 25% for T, 15% for G, and 29% for C. The mitogenome data are available in the GenBank under the accession numbers ON859136 and ON859137.

Dysbiosis of Gut Microbiome Aggravated Male Infertility in Captivity of Plateau Pika

Captivity is an important and efficient technique for rescuing endangered species. However, it induces infertility, and the underlying mechanism remains obscure. This study used the plateau pika (Ochotona curzoniae) as a model to integrate physiological, metagenomic, metabolomic, and transcriptome analyses and explore whether dysbiosis of the gut microbiota induced by artificial food exacerbates infertility in captive wild animals. Results revealed that captivity significantly decreased testosterone levels and the testicle weight/body weight ratio. RNA sequencing revealed abnormal gene expression profiles in the testicles of captive animals. The microbial α-diversity and Firmicutes/Bacteroidetes ratio were drastically decreased in the captivity group. Bacteroidetes and Muribaculaceae abundance notably increased in captive pikas. Metagenomic analysis revealed that the alteration of flora increased the capacity for carbohydrate degradation in captivity. The levels of microbe metabolites' short-chain fatty acids (SCFAs) were significantly high in the captive group. Increasing SCFAs influenced the immune response of captivity plateau pikas; pro-inflammatory cytokines were upregulated in captivity. The inflammation ultimately contributed to male infertility. In addition, a positive correlation was observed between Gastranaerophilales family abundance and testosterone concentration. Our results provide evidence for the interactions between artificial food, the gut microbiota, and male infertility in pikas and benefit the application of gut microbiota interference in threatened and endangered species.

Novel trypanosomatid species detected in Mongolian pikas (Ochotona pallasi) and their fleas in northwestern China

BACKGROUND

In the family Trypanosomatidae, the genus Trypanosoma contains protozoan parasites that infect a diverse range of hosts, including humans, domestic animals, and wildlife. Wild rodents, as natural reservoir hosts of various pathogens, play an important role in the evolution and emergence of Trypanosomatidae. To date, no reports are available on the trypanosomatid infection of pikas (Lagomorpha: Ochotonidae).

METHODS

In this study, Mongolian pikas and their fleas were sampled at the China-Mongolia border, northwestern China. The samples were analyzed with polymerase chain reaction (PCR) and sequencing for the presence of Trypanosomatidae on the basis of both the 18S ribosomal RNA (18S rRNA) gene and the glyceraldehyde-3-phosphate dehydrogenase (gGAPDH) gene. The morphology of trypomastigotes was also observed in peripheral blood smears by microscopy.

RESULTS

Molecular and phylogenetic analyses revealed a new genotype of the Trypanosoma lewisi clade that was found both in pika blood and flea samples. This genotype, which probably represents a new species, was provisionally designated as "Trypanosoma sp. pika". In addition, a novel genotype belonging to the genus Blechomonas of Trypanosomatidae was detected in fleas. On the basis of its molecular and phylogenetic properties, this genotype was named Blechomonas luni-like, because it was shown to be the closest related to B. luni compared with other flea-associated trypanosomatids.

CONCLUSIONS

To the best of our knowledge, this is the first study to report any trypanosomatid species in Mongolian pikas and their fleas. Further studies are needed to investigate the epidemiology of these protozoan parasites, as well as to evaluate their pathogenicity for humans or domestic animals.

Dispersal from the Qinghai-Tibet plateau by a high-altitude butterfly is associated with rapid expansion and reorganization of its genome

Parnassius glacialis is a typical "Out of the QTP" alpine butterfly that originated on the Qinghai-Tibet Plateau (QTP) and dispersed into relatively low-altitude mountainous. Here we assemble a chromosome-level genome of P. glacialis and resequence 9 populations in order to explore the genome evolution and local adaptation of this species. These results indicated that the rapid accumulation and slow unequal recombination of transposable elements (TEs) contributed to the formation of its large genome. Several ribosomal gene families showed extensive expansion and selective evolution through transposon-mediated processed pseudogenes. Additionally, massive structural variations (SVs) of TEs affected the genetic differentiation of low-altitude populations. These low-altitude populations might have experienced a genetic bottleneck in the past and harbor genes with selective signatures which may be responsible for the potential adaptation to low-altitude environments. These results provide a foundation for understanding genome evolution and local adaptation for "Out of the QTP" of P. glacialis.

Establishment of a secondary infection laboratory model of Echinococcus shiquicus metacestode using BALB/c mice and Mongolian jirds (Meriones unguiculatus)

Echinococcus shiquicus is peculiar to the Qinghai–Tibet plateau of China. Research on this parasite has mainly focused on epidemiological surveys and life cycle studies. So far, limited laboratory studies have been reported. Here, experimental infection of E. shiquicus metacestode in BALB/c mice and Mongolian jirds (Meriones unguiculatus) was carried out to establish alternative laboratory animal models. Intraperitoneal inoculation of metacestode material containing protoscoleces (PSCs) obtained from infected plateau pikas were conducted on BALB/c mice. Furthermore, metacestode material without PSCs deriving from infected BALB/c mice was intraperitoneally inoculated to Mongolian jirds. Experimental animals were dissected for macroscopic and histopathological examination. The growth of cysts in BALB/c mice was infiltrative, and they invaded the murine entire body. Most of the metacestode cysts were multicystic, but a few were unilocular. The cysts contained sterile vesicles, which had no PSCs. The metacestode materials were able to successfully infect new mice. In the jirds model, E. shiquicus cysts were typically formed freely in the peritoneal cavity; the majority of these cysts were free while a small portion adhered loosely to nearby organs. The proportion of fertile cysts was high, and contained many PSCs. The PSCs produced in Mongolian jirds also successfully infected new ones, which confirms that jirds can serve as an alternative experimental intermediate host. In conclusion, a laboratory animal infection was successfully established for E. shiquicus using BALB/c mice and Mongolian jirds. These results provide new models for the in-depth study of Echinococcus metacestode survival strategy, host interactions and immune escape mechanism.

Gene flow drives genomic diversity in Asian Pikas distributed along the core and range-edge habitats in the Himalayas

Studying the genetic variation among different species distributed across their core and range-edge habitats can provide valuable insights into how genetic variation changes across the species' distribution range. This information can be important for understanding local adaptation, as well as for conservation and management efforts. In this study, we have carried out genomic characterization of six species of Asian Pikas distributed along their core and range-edge habitats in the Himalayas. We utilized a population genomics approach using ~28,000 genome-wide SNP markers obtained from restriction-site associated DNA sequencing. We identified low nucleotide diversity and high inbreeding coefficients in all six species across their core and range-edge habitats. We also identified evidence of gene flow among genetically diverse species. Our results provide evidence of reduced genetic diversity in Asian pikas distributed across the Himalayas and the neighboring regions and indicate that recurrent gene flow is possibly a key mechanism for maintaining genetic diversity and adaptive potential in these pikas. However, full-scale genomics studies that utilize whole-genome sequencing approaches will be needed to quantify the direction and timing of gene flow and functional changes associated with introgressed regions in the genome. Our results represent an important step toward understanding the patterns and consequences of gene flow in species, sampled at the least studied, yet climatically vulnerable part of their habitat that can be further used to inform conservation strategies that promote connectivity and gene flow between populations.

The pikas of China: a review of current research priorities and challenges for conservation

Of the approximately 30 extant species of pika (Ochotona; Ochotonidae; Lagomorpha), at least 23 occur in China. Of these, 8 are endemic, and 3 are classified as Endangered by the IUCN. Research into most Chinese species is scarce, and there is much to learn about basic aspects of their ecology. We reviewed the literature on the 23 Chinese Ochotona species, with the aim of highlighting knowledge gaps and biases in research allocation. Specifically, we identify (1) which pika species receive the most attention from researchers, (2) which topics are the most frequently studied, and (3) how well research satisfies conservation priorities. We covered a total of 879 articles, most of which were written in English or Chinese. Around 75% of all publications focused on the plateau pika Ochotona curzoniae. Seven species were represented by 10 or fewer publications each. Endangered and endemic species were particularly poorly represented. 204 papers treated pikas as a pest, while 13 studied conservation issues. Nine species were considered possible targets for pest control, including some poorly known endemics. We make some recommendations on how research into Chinese Ochotona could be improved, including: (1) recognizing that the current species-level taxonomy may be an ineffective way to prioritize conservation research in Ochotona, (2) directing more research toward threatened and endemic species, subspecies, and populations, (3) researching the ecosystem engineering effects of pikas, (4) collecting basic data on natural history of the many understudied species, and (5) researching threats, including climate change and pest control campaigns.

Out of the Qinghai-Tibet plateau: Genomic biogeography of the alpine monospecific genus Megadenia (Biscutelleae, Brassicaceae)

Numerous high-elevation alpine plants of the Qinghai-Tibet Plateau (QTP) also have disjunct distribution in adjacent low-altitude mountains. The out-of-QTP versus into-the-QTP hypothesis of alpine plants provide strong evidence for the highly disputed assumption of the massive ice sheet developed in the central plateau during the Last Glacial Maximum (LGM). In this study, we sequenced the genomes of most known populations of Megadenia, a monospecific alpine genus of Brassicaceae distributed primarily in the QTP, though rarely found in adjacent low-elevation mountains of north China and Russia (NC-R). All sequenced samples clustered into four geographic genetic groups: one pair was in the QTP and another was in NC-R. The latter pair is nested within the former, and these findings support the out-of-QTP hypothesis. Dating the four genetic groups and niche distribution suggested that Megadenia migrated out of the QTP to adjacent regions during the LGM. The NC-R group showed a decrease in the effective population sizes. In addition, the genes with high genetic divergences in the QTP group were mainly involved in habitat adaptations during low-altitude colonization. These findings reject the hypothesis of development massive ice sheets, and support glacial survival of alpine plants within, as well as further migration out of, the QTP.

Evolution and diversification of Mountain voles (Rodentia: Cricetidae)

The systematics of the Cricetid genus Neodon have long been fraught with uncertainty due to sampling issues and a lack of comprehensive datasets. To gain better insights into the phylogeny and evolution of Neodon, we systematically sampled Neodon across the Hengduan and Himalayan Mountains, which cover most of its range in China. Analyses of skulls, teeth, and bacular structures revealed 15 distinct patterns corresponding to 15 species of Neodon. In addition to morphological analyses, we generated a high-quality reference genome for the mountain vole and generated whole-genome sequencing data for 47 samples. Phylogenomic analyses supported the recognition of six new species, revealing a long-term underestimation of Neodon diversity. We further identified positively selected genes potentially related to high-elevation adaptation. Together, our results illuminate how climate change caused the plateau to become the centre of Neodon origin and diversification and how mountain voles have adapted to the hypoxic high-altitude plateau environment.

Genomic Consequences of and Demographic Response to Pervasive Hybridization Over Time in Climate-Sensitive Pikas

Rare and geographically restricted species may be vulnerable to genetic effects from inbreeding depression in small populations or from genetic swamping through hybridization with common species, but a third possibility is that selective gene flow can restore fitness (genetic rescue). Climate-sensitive pikas (Ochotona spp.) of the Qinghai-Tibetan Plateau (QHTP) and its vicinity have been reduced to residual populations through the movement of climatic zones during the Pleistocene and recent anthropogenic disturbance, whereas the plateau pika (O. curzoniae) remains common. Population-level whole-genome sequencing (n = 142) of six closely related species in the subgenus Ochotona revealed several phases of ancient introgression, lineage replacement, and bidirectional introgression. The strength of gene flow was the greatest from the dominant O. curzoniae to ecologically distinct species in areas peripheral to the QHTP. Genetic analyses were consistent with environmental reconstructions of past population movements. Recurrent periods of introgression throughout the Pleistocene revealed an increase in genetic variation at first but subsequent loss of genetic variation in later phases. Enhanced dispersion of introgressed genomic regions apparently contributed to demographic recovery in three peripheral species that underwent range shifts following climate oscillations on the QHTP, although it failed to drive recovery of northeastern O. dauurica and geographically isolated O. sikimaria. Our findings highlight differences in timescale and environmental background to determine the consequence of hybridization and the unique role of the QHTP in conserving key evolutionary processes of sky island species.

Transcriptome analysis of pika heart tissue reveals mechanisms underlying the adaptation of a keystone species on the roof of the world

High-altitude environments impose intense stresses on living organisms and drive striking phenotypic and genetic adaptations, such as hypoxia resistance, cold tolerance, and increases in metabolic capacity and body mass. As one of the most successful and dominant mammals on the Qinghai-Tibetan Plateau (QHTP), the plateau pika (Ochotona curzoniae) has adapted to the extreme environments of the highest altitudes of this region and exhibits tolerance to cold and hypoxia, in contrast to closely related species that inhabit the peripheral alpine bush or forests. To explore the potential genetic mechanisms underlying the adaptation of O. curzoniae to a high-altitude environment, we sequenced the heart tissue transcriptomes of adult plateau pikas (comparing specimens from sites at two different altitudes) and Gansu pikas (O. cansus). Differential expression analysis and weighted gene co-expression network analysis (WGCNA) were used to identify differentially expressed genes (DEGs) and their primary functions. Key genes and pathways related to high-altitude adaptation were identified. In addition to the biological processes of signal transduction, energy metabolism and material transport, the identified plateau pika genes were mainly enriched in biological pathways such as the negative regulation of smooth muscle cell proliferation, the apoptosis signalling pathway, the cellular response to DNA damage stimulus, and ossification involved in bone maturation and heart development. Our results showed that the plateau pika has adapted to the extreme environments of the QHTP via protection against cardiomyopathy, tissue structure alterations and improvements in the blood circulation system and energy metabolism. These adaptations shed light on how pikas thrive on the roof of the world.

Phylogenomics reveals the evolution, biogeography, and diversification history of voles in the Hengduan Mountains

The Hengduan Mountains (HDM) of China are a biodiversity hotspot whose temperate flora and fauna are among the world's richest. However, the origin and evolution of biodiversity in the HDM remain poorly understood, especially in mammals. Given that the HDM shows the highest richness of vole species in the world, we used whole-exome capture sequencing data from the currently most comprehensive sampling of HDM voles to investigate their evolutionary history and diversification patterns. We reconstructed a robust phylogeny and re-estimated divergence times of the HDM voles. We found that all HDM voles could be divided into a western lineage (Volemys, Proedromys, and Neodon) and an eastern lineage (Caryomys and Eothenomys), and the two lineages originated from two migration events from North Eurasia to the HDM approximately 9 Mya. Both vole lineages underwent a significant acceleration of net diversification from 8-5 Mya, which was temporally congruent with the orogeny of the HDM region. We also identified strong intertribal gene flow among the HDM voles and hypothesized that frequent gene flow might have facilitated the speciation burst of the HDM voles. Our study highlights the importance of both environmental and biotic factors in shaping the biodiversity of mammals in mountain ecosystems.

Comparative genomics reveals putative evidence for high-elevation adaptation in the American pika (Ochotona princeps)

High-elevation environments have lower atmospheric oxygen content, reduced temperatures, and higher levels of UV radiation than found at lower elevations. As such, species living at high elevations must overcome these challenges to survive, grow, and reproduce. American pikas (Ochotona princeps) are alpine lagomorphs that are habitat specialists typically found at elevations >2,000 m. Previous research has shown putative evidence for high-elevation adaptation; however, investigations to date have been limited to a fraction of the genome. Here, we took a comparative genomics approach to identify putative regions under selection using a chromosomal reference genome assembly for the American pika relative to 8 other mammalian species targeted based on phylogenetic relatedness and (dis)similarity in ecology. We first identified orthologous gene groups across species and then extracted groups containing only American pika genes as well as unclustered pika genes to inform functional enrichment analyses; among these, we found 141 enriched terms with many related to hypoxia, metabolism, mitochondrial function/development, and DNA repair. We identified 15 significantly expanded gene families within the American pika across all orthologous gene groups that displayed functionally enriched terms associated with hypoxia adaptation. We further detected 196 positively selected genes, 41 of which have been associated with putative adaptation to hypoxia, cold tolerance, and response to UV following a literature review. In particular, OXNAD1, NRDC, and those genes critical in DNA repair represent important targets for future research to examine their functional implications in the American pika, especially as they may relate to adaptation to rapidly changing environments.

Expansion of Cyclophyllidea Biodiversity in Rodents of Qinghai-Tibet Plateau and the "Out of Qinghai-Tibet Plateau" Hypothesis of Cyclophyllideans

The Cyclophyllidea comprises the most species-rich order of tapeworms (Platyhelminthes, Cestoda) and includes species with some of the most severe health impact on wildlife, livestock, and humans. We collected seven Cyclophyllidea specimens from rodents in Qinghai-Tibet Plateau (QTP) and its surrounding mountain systems, of which four specimens in QTP were unsequenced, representing “putative new species.” Their complete mitochondrial (mt) genomes were sequenced and annotated. Phylogenetic reconstruction of partial 28S rDNA, cox1 and nad1 datasets provided high bootstrap frequency support for the categorization of three “putative new species,” assigning each, respectively, to the genera Mesocestoides, Paranoplocephala, and Mosgovoyia, and revealing that some species and families in these three datasets, which contain 291 species from nine families, may require taxonomic revision. The partial 18S rDNA phylogeny of 29 species from Taeniidae provided high bootstrap frequency support for the categorization of the “putative new species” in the genus Hydatigera. Combined with the current investigation, the other three known Taeniidae species found in this study were Taenia caixuepengi, T. crassiceps, and Versteria mustelae and may be widely distributed in western China. Estimates of divergence time based on cox1 + nad1 fragment and mt protein-coding genes (PCGs) showed that the differentiation rate of Cyclophyllidea species was strongly associated with the rate of change in the biogeographic scenarios, likely caused by the uplift of the QTP; i.e., species differentiation of Cyclophyllidea might be driven by host-parasite co-evolution caused by the uplift of QTP. We propose an “out of QTP” hypothesis for the radiation of these cyclophyllidean tapeworms.

Conservation genetics as a management tool: The five best-supported paradigms to assist the management of threatened species

About 50 y ago, Crow and Kimura [An Introduction to Population Genetics Theory (1970)] and Ohta and Kimura [Genet. Res. 22, 201-204 (1973)] laid the foundations of conservation genetics by predicting the relationship between population size and genetic marker diversity. This work sparked an enormous research effort investigating the importance of population dynamics, in particular small population size, for population mean performance, population viability, and evolutionary potential. In light of a recent perspective [J. C. Teixeira, C. D. Huber, Proc. Natl. Acad. Sci. U.S.A. 118, 10 (2021)] that challenges some fundamental assumptions in conservation genetics, it is timely to summarize what the field has achieved, what robust patterns have emerged, and worthwhile future research directions. We consider theory and methodological breakthroughs that have helped management, and we outline some fundamental and applied challenges for conservation genetics.

Chromosome-Level Reference Genome Assembly for the American Pika (Ochotona princeps)

The American pika (Ochotona princeps) is an alpine lagomorph found throughout western North America. Primarily inhabiting talus slopes at higher elevations (>2000 m), American pikas are well adapted to cold, montane environments. Warming climates on both historical and contemporary scales have contributed to population declines in American pikas, positioning them as a focal mammalian species for investigating the ecological effects of climate change. To support and expand ongoing research efforts, here, we present a highly contiguous and annotated reference genome assembly for the American pika (OchPri4.0). This assembly was produced using Dovetail de novo proximity ligation methods and annotated through the NCBI Eukaryotic Genome Annotation pipeline. The resulting assembly was chromosome- scale, with a total length of 2.23 Gb across 9350 scaffolds and a scaffold N50 of 75.8 Mb. The vast majority (>97%) of the total assembly length was found within 36 large scaffolds; 33 of these scaffolds correlated to whole autosomes, while the X chromosome was covered by 3 large scaffolds. Additionally, we identified 17 enriched gene ontology terms among American pika-specific genes putatively related to adaptation to high-elevation environments. This high-quality genome assembly will serve as a springboard for exploring the evolutionary underpinnings of behavioral, ecological, and taxonomic diversification in pikas as well as broader-scale eco-evolutionary questions pertaining to cold-adapted species in general.

Demographic History and Genomic Response to Environmental Changes in a Rapid Radiation of Wild Rats

For organisms to survive and prosper in a harsh environment, particularly under rapid climate change, poses tremendous challenges. Recent studies have highlighted the continued loss of megafauna in terrestrial ecosystems and the subsequent surge of small mammals, such as rodents, bats, lagomorphs, and insectivores. However, the ecological partitioning of these animals will likely lead to large variation in their responses to environmental change. In the present study, we investigated the evolutionary history and genetic adaptations of white-bellied rats (Niviventer Marshall, 1976), which are widespread in the natural terrestrial ecosystems in Asia but also known as important zoonotic pathogen vectors and transmitters. The southeastern Qinghai-Tibet Plateau was inferred as the origin center of this genus, with parallel diversification in temperate and tropical niches. Demographic history analyses from mitochondrial and nuclear sequences of Niviventer demonstrated population size increases and range expansion for species in Southeast Asia, and habitat generalists elsewhere. Unexpectedly, population increases were seen in N. eha, which inhabits the highest elevation among Niviventer species. Genome scans of nuclear exons revealed that among the congeneric species, N. eha has the largest number of positively selected genes. Protein functions of these genes are mainly related to olfaction, taste, and tumor suppression. Extensive genetic modification presents a major strategy in response to global changes in these alpine species.

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.

Genetic Evolution and Implications of the Mitochondrial Genomes of Two Newly Identified Taenia spp. in Rodents From Qinghai-Tibet Plateau

The larva of Taeniidae species can infect a wide range of mammals, causing major public health and food safety hazards worldwide. The Qinghai-Tibet Plateau (QTP), a biodiversity hotspot, is home to many species of rodents, which act as the critical intermediate hosts of many Taeniidae species. In this study, we identified two new larvae of Taenia spp., named T. caixuepengi and T. tianguangfui, collected from the plateau pika (Ochotona curzoniae) and the Qinghai vole (Neodon fuscus), respectively, in QTP, and their mitochondrial genomes were sequenced and annotated. Phylogenetic trees based on the mitochondrial genome showed that T. caixuepengi has the closest genetic relationship with T. pisiformis, while T. tianguangfui was contained in a monophyletic group with T. crassiceps, T. twitchelli, and T. martis. Biogeographic scenarios analysis based on split time speculated that the speciation of T. caixuepengi (∼5.49 Mya) is due to host switching caused by the evolution of its intermediate host. Although the reason for T. tianguangfui (∼13.11 Mya) speciation is not clear, the analysis suggests that it should be infective to a variety of other rodents following the evolutionary divergence time of its intermediate host and the range of intermediate hosts of its genetically close species. This study confirms the species diversity of Taeniidae in the QTP, and speculates that the uplift of the QTP has not only a profound impact on the biodiversity of plants and animals, but also that of parasites.

When parasites persist: tapeworms survive host extinction and reveal waves of dispersal across Beringia

Investigations of intercontinental dispersal between Asia and North America reveal complex patterns of geographic expansion, retraction and isolation, yet historical reconstructions are largely limited by the depth of the record that is retained in patterns of extant diversity. Parasites offer a tool for recovering deep historical insights about the biosphere, improving the resolution of past community-level interactions. We explored biogeographic hypotheses regarding the history of dispersal across Beringia, the region intermittently linking Asia and North America, through large-scale multi-locus phylogenetic analyses of the genus Schizorchis, an assemblage of host-specific cestodes in pikas (Lagomorpha: Ochotonidae). Our genetic data support palaeontological evidence for two separate geographic expansions into North America by Ochotona in the late Tertiary, a history that genomic evidence from extant pikas does not record. Pikas descending from the first colonization of Miocene age persisted into the Pliocene, subsequently coming into contact with a second wave of Nearctic colonists from Eurasia before going extinct. Spatial and temporal overlap of historically independent pika populations provided a window for host colonization, allowing persistence of an early parasite lineage in the contemporary fauna following the extinction of its ancestral hosts. Empirical evidence for ancient 'ghost assemblages' of hosts and parasites demonstrates how complex mosaic faunas are assembled in the biosphere through episodes of faunal mixing encompassing parasite lineages across deep and shallow time.