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

Evolutionary and functional characterization of lagomorph guanylate-binding proteins: a story of gain and loss and shedding light on expression, localization and innate immunity-related functions

Guanylate binding proteins (GBPs) are an evolutionarily ancient family of proteins that are widely distributed among eukaryotes. They belong to the dynamin superfamily of GTPases, and their expression can be partially induced by interferons (IFNs). GBPs are involved in the cell-autonomous innate immune response against bacterial, parasitic and viral infections. Evolutionary studies have shown that GBPs exhibit a pattern of gene gain and loss events, indicative for the birth-and-death model of evolution. Most species harbor large GBP gene clusters that encode multiple paralogs. Previous functional and in-depth evolutionary studies have mainly focused on murine and human GBPs. Since rabbits are another important model system for studying human diseases, we focus here on lagomorphs to broaden our understanding of the multifunctional GBP protein family by conducting evolutionary analyses and performing a molecular and functional characterization of rabbit GBPs. We observed that lagomorphs lack GBP3, 6 and 7. Furthermore, Leporidae experienced a loss of GBP2, a unique duplication of GBP5 and a massive expansion of GBP4. Gene expression analysis by reverse transcriptase quantitative polymerase chain reaction (RT-qPCR) and transcriptome data revealed that leporid GBP expression varied across tissues. Overexpressed rabbit GBPs localized either uniformly and/or discretely to the cytoplasm and/or to the nucleus. Oryctolagus cuniculus (oc)GBP5L1 and rarely ocGBP5L2 were an exception, colocalizing with the trans-Golgi network (TGN). In addition, four ocGBPs were IFN-inducible and only ocGBP5L2 inhibited furin activity. In conclusion, from an evolutionary perspective, lagomorph GBPs experienced multiple gain and loss events, and the molecular and functional characteristics of ocGBP suggest a role in innate immunity.

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.

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.

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.

TLR7 and TLR8 evolution in lagomorphs: different patterns in the different lineages

Toll-like receptors (TLRs) are one of the most ancient and widely studied innate immune receptors responsible for host defense against invading pathogens. Among the known TLRs, TLR7 and TLR8 sense and recognize single-stranded (ss) RNAs with a dynamic evolutionary history. While TLR8 was lost in birds and duplicated in turtles and crocodiles, TLR7 is duplicated in some birds, but in other tetrapods, there is only one copy. In mammals, with the exception of lagomorphs, TLR7 and TLR8 are highly conserved. Here, we aim to study the evolution of TLR7 and TLR8 in mammals, with a special focus in the order Lagomorpha. By searching public sequence databases, conducting evolutionary analysis, and evaluating gene expression, we were able to confirm that TLR8 is absent in hares but widely expressed in the European rabbit. In contrast, TLR7 is absent in the European rabbit and quite divergent in hares. Our results suggest that, in lagomorphs, more in particular in leporids, TLR7 and TLR8 genes have evolved faster than in any other mammalian group. The long history of interaction with viruses and their location in highly dynamic telomeric regions might explain the pattern observed.

Convergent Loss of the Necroptosis Pathway in Disparate Mammalian Lineages Shapes Viruses Countermeasures

Programmed cell death is a vital process in the life cycle of organisms. Necroptosis, an evolutionary form of programmed necrosis, contributes to the innate immune response by killing pathogen-infected cells. This virus-host interaction pathway is organized around two components: the receptor-interacting protein kinase 3 (RIPK3), which recruits and phosphorylates the mixed lineage kinase-like protein (MLKL), inducing cellular plasma membrane rupture and cell death. Critically, the presence of necroptotic inhibitors in viral genomes validates necroptosis as an important host defense mechanism. Here, we show, counterintuitively, that in different mammalian lineages, central components of necroptosis, such as RIPK3 and MLKL, are deleted or display inactivating mutations. Frameshifts or premature stop codons are observed in all the studied species of cetaceans and leporids. In carnivores’ genomes, the MLKL gene is deleted, while in a small number of species from afrotheria and rodentia premature stop codons are observed in RIPK3 and/or MLKL. Interestingly, we also found a strong correlation between the disruption of necroptosis in leporids and cetaceans and the absence of the N-terminal domain of E3-like homologs (responsible for necroptosis inhibition) in their naturally infecting poxviruses. Overall, our study provides the first comprehensive picture of the molecular evolution of necroptosis in mammals. The loss of necroptosis multiple times during mammalian evolution highlights the importance of gene/pathway loss for species adaptation and suggests that necroptosis is not required for normal mammalian development. Moreover, this study highlights a co-evolutionary relationship between poxviruses and their hosts, emphasizing the role of host adaptation in shaping virus evolution.

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.

Not so unique to Primates: The independent adaptive evolution of TRIM5 in Lagomorpha lineage

The plethora of restriction factors with the ability to inhibit the replication of retroviruses have been widely studied and genetic hallmarks of evolutionary selective pressures in Primates have been well documented. One example is the tripartite motif-containing protein 5 alpha (TRIM5α), a cytoplasmic factor that restricts retroviral infection in a species-specific fashion. In Lagomorphs, similarly to what has been observed in Primates, the specificity of TRIM5 restriction has been assigned to the PRYSPRY domain. In this study, we present the first insight of an intra-genus variability within the Lagomorpha TRIM5 PRYSPRY domain. Remarkably, and considering just the 32 residue-long v1 region of this domain, the deduced amino acid sequences of Daurian pika (Ochotona dauurica) and steppe pika (O. pusilla) evidenced a high divergence when compared to the remaining Ochotona species, presenting values of 44% and 66% of amino acid differences, respectively. The same evolutionary pattern was also observed when comparing the v1 region of two Sylvilagus species members (47% divergence). However, and unexpectedly, the PRYSPRY domain of Lepus species exhibited a great conservation. Our results show a high level of variation in the PRYSPRY domain of Lagomorpha species that belong to the same genus. This suggests that, throughout evolution, the Lagomorpha TRIM5 should have been influenced by constant selective pressures, likely as a result of multiple different retroviral infections.

Strong selection of the TLR2 coding region among the Lagomorpha suggests an evolutionary history that differs from other mammals

Toll-like receptors (TLRs) are one of the first lines of defense against pathogens and are crucial for triggering an appropriate immune response. Among TLRs, TLR2 is functional in all vertebrates and has high ability in detecting bacterial and viral pathogen ligands. The mammals' phylogenetic tree of TLR2 showed longer branches for the Lagomorpha clade, raising the hypothesis that lagomorphs experienced an acceleration of the mutation rate. This hypothesis was confirmed by (i) Tajima's test of neutrality that revealed different evolutionary rates between lagomorphs and the remaining mammals with lagomorphs presenting higher nucleotide diversity; (ii) genetic distances were similar among lagomorphs and between lagomorphs and other mammals; and (iii) branch models reinforced the existence of an acceleration of the mutation rate in lagomorphs. These results suggest that the lagomorph TLR2 has been strongly involved in pathogen recognition, which probably caused a host-pathogen arms race that led to the observed acceleration of the mutation rate.

Evolution of CCL16 in Glires (Rodentia and Lagomorpha) shows an unusual random pseudogenization pattern

Background

The C-C motif chemokine ligand 16 (CCL16) is a potent pro-inflammatory chemokine and a chemoattractant for monocytes and lymphocytes. In normal plasma, it is present at high concentrations and elicits its effects on cells by interacting with cell surface chemokine receptors. In the European rabbit and in rodents such as mouse, rat and guinea pig, CCL16 was identified as a pseudogene, while in the thirteen-lined ground squirrel it appears to be potentially functional. To gain insight into the evolution of this gene in the superorder Glires (rodents and lagomorphs), we amplified the CCL16 gene from eleven Leporidae and seven Ochotonidae species.

Results

We compared our sequences with CCL16 sequences of twelve rodent species retrieved from public databases. The data show that for all leporid species studied CCL16 is a pseudogene. This is primarily due to mutations at the canonical Cys Cys motif, creating either premature stop codons, or disrupting amino acid replacements. In the Mexican cottontail, CCL16 is pseudogenized due to a frameshift deletion. Additionally, in the exon 1 (signal peptide), there are frameshift deletions present in all leporids studied. In contrast, in Ochotona species, CCL16 is potentially functional, except for an allele in Hoffmann’s pika. In rodents, CCL16 is functional in a number of species, but patterns of pseudogenization similar to those observed in lagomorphs also exist.

Conclusions

Our results suggest that while functional in the Glires ancestor, CCL16 underwent pseudogenization in some species. This process occurred stochastically or in specific lineages at different moments in the evolution of Glires. These observations suggest that the CCL16 had different evolutionary constrains in the Glires group that could be associated with the CCL16 biological function.

Electronic supplementary material

The online version of this article (10.1186/s12862-019-1390-7) contains supplementary material, which is available to authorized users.

Evolution for extreme living: variation in mitochondrial cytochrome c oxidase genes correlated with elevation in pikas (genus Ochotona)

The genus Ochotona (pikas) is a clade of cold-tolerant lagomorphs that includes many high-elevation species. Pikas offer a unique opportunity to study adaptations and potential limitations of an ecologically important mammal to high-elevation hypoxia. We analyzed the evolution of 3 mitochondrial genes encoding the catalytic core of cytochrome c oxidase (COX) in 10 pika species occupying elevations from sea level to 5000 m. COX is an enzyme highly reliant on oxygen and essential for cell function. One amino acid property, the equilibrium constant (ionization of COOH), was found to be under selection in the overall protein complex. We observed a strong relationship between the net value change in this property and the elevation each species occupies, with higher-elevation species having potentially more efficient proteins. We also found evidence of selection in low-elevation species for potentially less efficient COX, perhaps trading efficiency for heat production in the absence of hypoxia. Our results suggest that different pika species may have evolved elevation-specific COX proteins, specialization that may indicate limitations in their ability to shift their elevational ranges in response to future climate change.

Evolutionary studies on the betaretrovirus RERV-H in the Leporidae family reveal an endogenization in the ancestor of Oryctolagus, Bunolagus and Pentalagus at 9 million years ago

RERV-H was first identified in human tissues and mistaken for a human exogenous retrovirus. However, the integration sites carried by this virus showed that it was instead a European rabbit (Oryctolagus cuniculus) endogenous retrovirus. The first clones retrieved from European rabbit samples represented defective proviruses, although estimation of proviral copy numbers found in the European rabbit genome ranged from hundreds to thousands. Screening for the presence of RERV-H showed the absence of the virus in two other lagomorphs, pika (Ochotona) and hares (Lepus), which diverged from rabbits about 35 and 12 million years ago, respectively. Using a PCR-based approach, samples of seven different Lagomorph genera were tested for the presence of RERV-H. It was possible to amplify a proviral fragment corresponding to RNaseH from Oryctolagus, Bunolagus and Pentalagus genomic samples. The amplification of proviral DNA in species other than Oryctolagus revealed that this virus was endogenized in their common ancestor, roughly 9 million years ago. Using the European rabbit genome sequence OryCun2.0, it was possible to find multiple copies spread throughout the genome and several complete proviral genomes were retrieved. Some copies contained full open reading frames for all viral components. The lack of a complete genome in the other Lagomorph species did not allow further analyses of the provirus, although more deleterious mutations were found in Bunolagus and Pentalagus than in Oryctolagus RNaseH-amplified sequences. To what extent RERV-H and other endogenous viruses might have had an impact on the rabbit genome and its immune system remains elusive.

From rabbit antibody repertoires to rabbit monoclonal antibodies

In this review, we explain why and how rabbit monoclonal antibodies have become outstanding reagents for laboratory research and increasingly for diagnostic and therapeutic applications. Starting with the unique ontogeny of rabbit B cells that affords highly distinctive antibody repertoires rich in in vivo pruned binders of high diversity, affinity and specificity, we describe the generation of rabbit monoclonal antibodies by hybridoma technology, phage display and alternative methods, along with an account of successful humanization strategies.

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.

Genetics, morphology and ecology reveal a cryptic pika lineage in the Sikkim Himalaya

Asian pika species are morphologically ∼similar and have overlapping ranges. This leads to uncertainty and species misidentification in the field. Phylogenetic analyses of such misidentified samples leads to taxonomic ambiguity. The ecology of many pika species remains understudied, particularly in the Himalaya, where sympatric species could be separated by elevation and/or substrate. We sampled, measured, and acquired genetic data from pikas in the Sikkim Himalaya. Our analyses revealed a cryptic lineage, Ochotona sikimaria, previously reported as a subspecies of O. thibetana. The results support the elevation of this lineage to the species level, as it is genetically divergent from O. thibetana, as well as sister species, O. cansus (endemic to central China) and O. curzoniae (endemic to the Tibetan plateau). The Sikkim lineage diverged from its sister species' about 1.7-0.8myrago, coincident with uplift events in the Himalaya. Our results add to the recent spate of cryptic diversity identified from the eastern Himalaya and highlight the need for further study within the Ochotonidae.

LaGomiCs-Lagomorph Genomics Consortium: An International Collaborative Effort for Sequencing the Genomes of an Entire Mammalian Order

The order Lagomorpha comprises about 90 living species, divided in 2 families: the pikas (Family Ochotonidae), and the rabbits, hares, and jackrabbits (Family Leporidae). Lagomorphs are important economically and scientifically as major human food resources, valued game species, pests of agricultural significance, model laboratory animals, and key elements in food webs. A quarter of the lagomorph species are listed as threatened. They are native to all continents except Antarctica, and occur up to 5000 m above sea level, from the equator to the Arctic, spanning a wide range of environmental conditions. The order has notable taxonomic problems presenting significant difficulties for defining a species due to broad phenotypic variation, overlap of morphological characteristics, and relatively recent speciation events. At present, only the genomes of 2 species, the European rabbit (Oryctolagus cuniculus) and American pika (Ochotona princeps) have been sequenced and assembled. Starting from a paucity of genome information, the main scientific aim of the Lagomorph Genomics Consortium (LaGomiCs), born from a cooperative initiative of the European COST Action "A Collaborative European Network on Rabbit Genome Biology-RGB-Net" and the World Lagomorph Society (WLS), is to provide an international framework for the sequencing of the genome of all extant and selected extinct lagomorphs. Sequencing the genomes of an entire order will provide a large amount of information to address biological problems not only related to lagomorphs but also to all mammals. We present current and planned sequencing programs and outline the final objective of LaGomiCs possible through broad international collaboration.

An overview of the lagomorph immune system and its genetic diversity

Our knowledge of the lagomorph immune system remains largely based upon studies of the European rabbit (Oryctolagus cuniculus), a major model for studies of immunology. Two important and devastating viral diseases, rabbit hemorrhagic disease and myxomatosis, are affecting European rabbit populations. In this context, we discuss the genetic diversity of the European rabbit immune system and extend to available information about other lagomorphs. Regarding innate immunity, we review the most recent advances in identifying interleukins, chemokines and chemokine receptors, Toll-like receptors, antiviral proteins (RIG-I and Trim5), and the genes encoding fucosyltransferases that are utilized by rabbit hemorrhagic disease virus as a portal for invading host respiratory and gut epithelial cells. Evolutionary studies showed that several genes of innate immunity are evolving by strong natural selection. Studies of the leporid CCR5 gene revealed a very dramatic change unique in mammals at the second extracellular loop of CCR5 resulting from a gene conversion event with the paralogous CCR2. For the adaptive immune system, we review genetic diversity at the loci encoding antibody variable and constant regions, the major histocompatibility complex (RLA) and T cells. Studies of IGHV and IGKC genes expressed in leporids are two of the few examples of trans-species polymorphism observed outside of the major histocompatibility complex. In addition, we review some endogenous viruses of lagomorph genomes, the importance of the European rabbit as a model for human disease studies, and the anticipated role of next-generation sequencing in extending knowledge of lagomorph immune systems and their evolution.

Endogenization of mouse mammary tumor virus (MMTV)-like elements in genomes of pikas (Ochotona sp.)

Despite the finding in European rabbit and other leporid genomes of the first ever described endogenous lentivirus and of a European rabbit exclusive endogenous gammaretrovirus, until now no exogenous retroviruses have been isolated in Lagomorpha species. Nevertheless, looking for the presence of endogenous retroviruses (ERVs) in the species genomes could lead to the discovery of retroviral lineages yet to be found in Lagomorpha. Different mammalian genomes harbor endogenous viral sequences phylogenetically close to the betaretrovirus mouse mammary tumor virus (MMTV), propelling us to look for such retroviral "fossil" in American pika (Ochotona princeps) and European rabbit (Oryctolagus cuniculus) genomes. By performing genomic mining using MMTV gag and LTR as query sequences, we found that such viral elements were absent from the European rabbit genome. Oppositely, significant matches were found in American pika, and more importantly, a nearly complete MMTV-like virus (Pika-BERV) was identified. Using Pika-BERV gag and LTR as templates, we found similar sequences endogenized in different pika (Ochotona sp.) species. The orthology of the LTR flanking region between some pika species supported shared ancestry of specific endogenous betaretroviruses, while in other pika species similar sequences, but not orthologous, should have resulted from independent insertions. Our study supports the possible existence of infecting exogenous betaretroviruses for a long term, after the divergence of Ochotonidae from Leporidae, but yet to be identified.

Pseudogenization of CCL14 in the Ochotonidae (pika) family

The interaction between chemokines and their receptors is crucial for inflammatory cell trafficking. CCL14 binds with high affinity to CCR5. In leporids, CCR5 underwent gene conversion with CCR2. The study of CCR5 ligands in leporid species showed that CCL8 is pseudogenized, while CCL3, CCL4 and CCL5 are functional. Here, we study the evolution of CCL14 in mammals with emphasis in the order Lagomorpha. By employing maximum likelihood methods we detected six sites under positive selection. Some of these sites are located in regions crucial for CCL14 activation and binding to receptors. Sequencing of CCL14 in Ochotona species showed that O. princeps, O. pallasi, O. alpina and O. turuchanensis have a mutation at the start codon (Met > Thr), while O. hoffmanni, O. mantchurica, O. dauurica and O. rufescens present the mammalian conserved Met. Ochotona hyperborea has the two alleles. In O. pusilla, CCL14 is a pseudogene due to a seven base pair insertion. Like CCL3, CCL4 and CCL5, CCL14 is functional in all leporids but in the Ochotonidae family it underwent a pseudogenization process. This suggests that CCL14 has an important biological role in other mammals by evolving under positive selection that has been lost in Ochotonidae (subgenera Pika and Lagotona).