Pseudogenization of CCL14 in the Ochotonidae (pika) family

Abstract

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).

Keywords

Chemokine ligands CCL14 evolution Ochotonidae pseudogenization

Introduction

Chemokines are proteins that have several functions in the immune system and are found only in vertebrates.¹ These proteins are encoded by a multigene family that emerged by gene duplication events.¹ Gene duplication is considered a powerful force in the adaptive evolution of genes, in which one gene gives rise to new copies.^2–4 These new copies can evolve by pseudogenization, neofunctionalization or subfunctionalization.^3,5,6 Pseudogenization is the process in which a redundant functional copy passes through mutational decay and loses its ability to produce a functional protein.^6,7 In the remaining processes, both copies escape to mutation as either one copy emerges into new functions (neofunctionalization) or either divides the functions with the other copy (subfunctionalization).^3,5,8

Chemokines can only exert their functions through binding to specific receptors, having the ability to bind and activate different receptors; for example, CCL14 binds with high affinity to the chemokine receptors CCR1, CCR3 and CCR5.⁹ CCR5 is used as co-receptor by HIV type-1 (HIV-1) for infection and, despite controversy, it has been also associated with myxoma virus infection in rabbits.^10–17 In some mammalian species, CCR5 underwent gene conversion with CCR2.^18–23 Several studies have been conducted in the CCR5 ligands to determine the consequences of the CCR5–CCR2 gene conversion and its putative association with the myxoma virus infection. Indeed, the study of CCL8, a potential CCR5 ligand, revealed that in some leporid species, such as Oryctolagus and Bunolagus, this ligand is pseudogenized, whilst in Lepus and Sylvilagus is intact.¹² A study on CCL3, CCL4 and CCL5 showed that, unlike CCL8, these genes are all functional in leporids.²⁴ Although mouse CCL12 is the true ortholog to human CCL8,²⁵ in leporids, the gene identified as CCL8 is the true ortholog, as supported by their phylogenetic position within the mammalian CCL8 cluster.¹²

In leporids, other CCR5 ligands remain to be studied, including CCL14. CCL14, also known as hemofiltrate CC-chemokine 1 (HCC-1), is found in high concentrations in blood plasma and acts as an inflammatory chemokine after N-terminus cleavage.^15,26,27 In humans, this protein is located in the macrophage inflammatory protein (MIP) region of CC chemokine genes,^1,28 and is expressed in numerous tissues.^9,15,27,29 In humans, CCL14 encodes a protein with 93 aa with three predicted cleavage sites (Thr25, Glu31 and Ser33) that are necessary for the protein to become active.³⁰ There are two described variants: one encoded by three exons in which the active CCL14 consists of a 74 aa protein after the signal peptide cleavage, and another variant encoded by four exons, which has an insertion of 48 nucleotides between exons 1 and 2 encoding 16 aa.^1,15,17,27 These forms have no effect against HIV-1 infection; however, the cleaved form CCL14 (9–74) has the ability to block HIV-1 entry and replication.¹⁶ CCL14 (9–74) is an agonist of CCR1, CCR3 and CCR5, which promotes calcium influx and migration of T lymphocytes, monocytes and eosinophils; through the interaction with the second external loop of CCR5, it has also the ability to inhibit this receptor.^15–17,27 In CC chemokines there are two main interactions important for receptor binding and activation, both located in the N-terminus close to the CC motif. In CCL14, the first interaction corresponds to the high-affinity binding to the receptor and includes residues located after the CC motif, and the second interaction comprises the residues located before the CC motif and that are essential for the receptor activation.^15,31–34 The N-terminus is also important for CCL14 degradation.³⁵ Degradation occurs following the binding of D6 to the N-terminus of CCL14 and then through the interaction with the residue Pro53 of CCL14.^35,36

The order Lagomorpha includes two families, Ochotonidae (pikas) and Leporidae (hares and rabbits) that diverged ∼35 million years ago.³⁷ The Ochotonidae family is composed of only one genus, Ochotona, which is divided into four subgenera, Pika, Ochotona, Conothoa and Lagotona.³⁸ The Leporidae family is divided into 11 genera.³⁷ Rodents and lagomorphs form the superorder Glires that is considered a sister group of the superorder Euarchonta that includes primates.^39,40 In rodents, CCL14 was shown to be absent in mouse and rat, while present in squirrel and guinea pig;²⁸ in lagomorphs, information is only available for the European rabbit (Oryctolagus cuniculus), where a functional gene exists.

By comparing the O. princeps CCL14 coding sequence (CDS) available in GenBank/Ensembl with the remaining mammalian species it was possible to observe that the initiation codon typically present in mammals is mutated into a Thr (ACG) (Figure 1). Following this observation, and considering the apparent ‘random’ absence of CCL14 in rodents and its role as a ligand for CCR5 which in some lagomorphs underwent gene conversion with CCR2, we aimed at determining the presence and function conservation of CCL14 in lagomorphs. In addition, we studied the evolution of CCL14 in mammals by identifying codons that are evolving under positive selection.

Figure 1.

(A) Alignment of CCL14 for the different mammalian species (GenBank, Ensembl and Uniprot accession numbers are indicated for the retrieved sequences). Positively selected amino acids are boxed and the three possible cleavage sites are shaded in dark grey. (*) Represent stop codons, (–) represent deletions, *1 and *2 represent alleles. (B) Alignment of CCL14 for the different Ochotona species studied. The characteristic mammalian initiation codon is boxed. (C) Alignment of the CCL14 for O. princeps and O. pusilla with the insertion of the seven base pairs shaded in grey.

Material and methods

Mammalian CCL14 sequences were retrieved from public databases (accession numbers are given in Figure 1). Sequences were aligned using MUltiple Sequence Comparison by Log-Expectation (MUSCLE) available at http://www.ebi.ac.uk/.⁴¹ Similar results were obtained using Multiple Alignment using Fast Fourier Transform (MAFFT); however, MUSCLE alignment introduced less gaps, leading to a more conserved alignment.

In order to determine the impact of positive selection on CCL14 sequence evolution, we estimated ω, i.e, the ratio of nonsynonymous (dN) to synonymous (dS) differences in functional CCL14 orthologs (final dataset of 55 sequences). The codon-based maximum likelihood (ML) method (CODEML) implemented in PAML v4.4 was used.⁴² An un-rooted neighbor-joining tree was constructed using MEGA6,⁴³ with p-distance as substitution model and the pairwise deletion option for gaps/missing data. The topology of the phylogenetic tree obtained is in accordance with the accepted mammalian phylogeny. Two pairs of site-based models were compared: M1 (nearly neutral) vs. M2 (selection) and M7 (neutral, beta) vs. M8 (selection, beta and ω), where M1 and M7 correspond to the null hypothesis and M2 and M8 to the alternative hypothesis by allowing positive selection. A likelihood ratio test with two degrees of freedom determined if a selection model fitted better the data than a neutral model.^44,45 Codons under positive selection were identified by using a Bayes Empirical Bayes approach with probability >95%. Codon-based ML methods available in the HYPHY package implemented in the DataMonkey webserver were also used:^46,47 Single Likelihood Ancestor Counting (SLAC), Fixed-Effect Likelihood (FEL), Internal Branch Fel (iFEL), Mixed Effects Model of Evolution (MEME), Random Effect Likelihood (REL) and Fast Unconstrained Bayesian AppRoximation (FUBAR).^47–49 For the first four methods the P-value was set to ≤ 0.05; for FUBAR we used a P-value ≥ 0.95; and for REL we used a Bayes factor >95. The best fitting model for nucleotide substitution was determined by the automatic model selection available in the webserver. As done previously for other immunity genes,^50–55 only the codons detected by more than one method were considered as being under positive selection.

Samples from pikas (Ochotona alpina, O. dauurica, O. hyperborea, O. hoffmanni, O. mantchurica, O. pallasi, O. princeps, O. pusilla, O. rufescens and O. turuchanensis) were obtained from the Zoological Museum of Moscow University collection, and samples of O. cuniculus (subspecies O. c. cuniculus and O. c. algirus), Brachylagus idahoensis, Pentalagus furnessi, Sylvilagus bachmani, Sylvilagus floridanus, Lepus europaeus, Romerolagus diazi and Ochotona princeps were provided by the CIBIO Lagomorph tissue collection. For RNA extraction, a liver sample of O. princeps was kindly provided by Jay Storz from School of Biological Sciences, University of Nebraska, USA. Genomic DNA was extracted using the EasySpin Genomic DNA Minipreps Tissue Kit (Citomed, Torun, Poland) according to manufacturer’s instructions. Total RNA was extracted for one specimen of Oryctolagus cuniculus cuniculus, O. c. algirus, Lepus europaeus, S. floridanus and O. princeps using a RNeasy Mini Kit according to the manufacturer’s instructions (Qiagen, Hilden, Germany). cDNA was synthesized using oligo(dT) as primers and SuperScript III reverse transcriptase (Invitrogen, Carlsbad, CA, USA). CCL14 was PCR-amplified and sequenced; sequences were submitted to GenBank under the following accession numbers: KP407602-KP407627 and KP641178. The European rabbit (Oryctolagus cuniculus) and Pika (Ochotona princeps) CCL14 sequences available in GenBank were used as templates for primer design. For amplification of CCL14 from gDNA, amplification was partitioned and amplified with Multiplex PCR Kit (Qiagen). Primers and conditions are listed in Table 1. Sequencing was performed on an ABI PRISM 310 Genetic Analyser (PE Applied Biosystems, Foster City, CA, USA) and PCR products were sequenced in both directions.

Table 1.

Primers and conditions used in CCL14 PCR amplification.

Species amplified	Primers sequence (5’–3’)	Primer name	Exons amplified	PCR conditions	Fragment length
Leporids	CAGCTCTCCTGCAGCAGC	CCL14OrcuF1	Exon1	Initial denaturation: 95℃ (15 min) 40 cycles [95℃ (45 s), 58℃ (20 s), 72℃ (10 s)] Final extension: 60℃ (10 min)	∼ 182 bp
	GGAGTTGGATAGCTGAGAG	CCL14OrcuR1	Exon1		∼ 182 bp
	CAGAAGCAGCCAGCAATC	CCL14OrcuF2	Exon 2 + exon 3	Initial denaturation: 95℃ (15 min) 40 cycles [95℃ (45 s), 58℃ (20 s), 72℃ (45 s)] Final extension: 60℃ (10 min)	∼ 732 bp
	CAAAGGCACAGCTCAGAG	CCL14OrcuR2	Exon 2 + exon 3		∼ 732 bp
Ochotona	CTCAGCACATCCTTACCCAG	CCL14OcprF1	Exon1	Initial denaturation: 95℃ (15 min) 40 cycles [95℃ (45 s), 60℃ (20 s), 72℃ (10 s)] Final extension: 60℃ (10 min)	∼ 200 bp
	CTTCTAGCCTCGTCTCCAG	CCL14OcprR1	Exon1		∼ 200 bp
	CAGCAATCATGTCTCACCAG	CCL14OcprF2	Exon 2 + exon 3	Initial denaturation: 95℃ (15 min) 40 cycles [95℃ (45 s), 62℃ (20 s), 72℃ (45 s)] Final extension: 60℃ (10 min)	∼ 612 bp
	GAAGGCACAGTTCACAGACC	CCL14OcprR2			∼ 612 bp
	CCAGGGTAGCAGAAGCAGC	CCL14OcprF2a		Initial denaturation: 95℃ (15 min) 40 cycles [95℃ (45 s), 60℃ (20 s), 72℃ (45 s)] Final extension: 60℃ (10 min)	∼ 692 bp
	CTATGCCAGGAATTGCTCTG	CCL14OcprR2a			∼ 692 bp
	GTTCCTCTAACACCTCTCTGG	CCL14OcprF2b		Initial denaturation: 95℃ (15 min) 40 cycles [95℃ (45 s), 62℃ (20 s), 72℃ (1 min)] Final extension: 60℃ (10 min)	∼ 801 bp
	CTTGTGCACCTCACTGAG	CCL14OcprR2b			∼ 801 bp

The program PHASE, built into the DnaSP software,⁵⁶ was used to reconstruct the haplotype phases of the obtained sequences. These sequences were aligned and translated using BioEdit.⁵⁷

The Tajima’s Relative Rate Test was conducted in order to check the evolutionary rate of CCL14 between lagomorphs using MEGA6.^43,58

Results and discussion

In our study, six aa were identified as being under positive selection: Ala22, Thr25, Thr62, Lys67, Met74 and His96 (aa residues were numbered from the ﬁrst methionine residue in human CCL14. The signal peptide and indels, indicated as (–), were included in the numbering) (Figure 1). Ala22 and Thr25 are located in the pre-peptide, which is cleaved in order for the protein to become active;¹⁷ however, the modifications observed at residue 22 are not expected to alter the CCL14 protein. In contrast, Thr25 is located within the region where CCL14 cleavage occurs. As cleavage is required for protein activation, alterations in this site might negatively affect the CCL14 protein. Thr62 and Lys67 are in the region reported as to be essential to the binding of CCL14 to the several receptors;^15,31–34 thus, changes in this region may not only affect both the ligand–receptor binding, but also CCL14 degradation.^35,36 Met74 and His96 are located in the β1 and β3 sheets, respectively. Alterations in these sites may lead to alterations in the protein conformation.

The comparison of the O. princeps CCL14 CDS available in GenBank and Ensembl with the remaining mammalian species showed that in O. princeps the initiation codon, in contrast to the other mammalian sequences, is mutated into a Thr (ACG) (Figure 1A). Following this observation, we attempted the amplification of CCL14 from gDNA of O. princeps. Our results indicated that, as for the sequences publicly available, the amplified fragment also presents the ACG mutation. Next, we attempted to amplify CCL14 from cDNA of O. princeps. Despite having used different Taq polymerases and different PCR conditions (e.g. including the addition of DMSO in the reaction, increasing extension times) and different pairs of primers designed according to the O. princeps CDS, and having previously successfully amplified other CDS of other genes from the same sample,⁵⁹ we were unable to amplify the CCL14. This may suggest that CCL14 is a pseudogene in O. princeps. We amplified and sequenced the CCL14 gene for eight leporid species: O. cuniculus cuniculus and O. c. algirus, Brachylagus idahoensis, Pentalagus furnessi, Sylvilagus bachmani and S. floridanus, Lepus europaeus, Romerolagus diazi.

To support the putative pseudogenization of CCL14 in O. princeps we performed a Tajima’s Relative Rate Test in order to determine the rates of evolution of CCL14 in lagomorphs. This test determines if two genes (paralogs or orthologs) are evolving at similar rates, i.e. it determines statistically if the two sequences follow the molecular clock hypothesis of similar constant rates of evolution (null hypothesis) or if one sequence evolved at an accelerated or decelerated rate (alternative hypothesis).⁵⁸ By comparing the O. c. cuniculus CCL14 gene with the O. princeps orthologs, and using as an outgroup each of the different leporids studied, the P-values of the χ² were always statistically significant. This rejects the null hypothesis and favors the hypothesis of an acceleration of the mutation rate in Ochotona princeps. Thus, our results support a significant bias to neutrality between leporids and O. princeps (Table 2), compatible with a pseudogenization process in this species.

Table 2.

Results obtained in Tajima’s Relative Rate Test.

Outgroup	No. divergent sites in the three sequences^a	No. differences in O. c. cuniculus	Nr. differences in O. princeps	χ^2b	P-value^c
Sylvilagus bachmani	1	5	37	24.38	0.00000
Sylvilagus floridanus	2	7	34	17.78	0.00002
Brachylagus idahoensis	2	4	37	26.56	0.00000
Lepus europaeus	1	6	36	21.43	0.00000
Pentalagus furnessi	3	2	38	32.40	0.00000
Romerolagus diazi	1	9	33	13.71	0.00021

Nr means Number of.

Chi-square is a statistical test used to determine the substitution rates between species.

A P-value < 0.05 is used to reject the null hypothesis of equal rates between lineages.

We attempted to amplify the CCL14 gene from gDNA of several pika species that encompass the different Ochotona subgenera: O. alpina, O. dauurica, O. hyperborea, O. hoffmanni, O. mantchurica, O. pallasi, O. pusilla, O. rufescens and O. turuchanensis. Interestingly, we found different results within the Ochotonidae family (Figure 1B, C). Indeed, O. alpina, O. pallasi and O. turuchanensis, share the ATG > ACG mutation with O. princeps; in contrast, in O. dauurica, O. hoffmanni, O. mantchurica and O. rufescens, as in all the other mammals examined, the putative initiation codon is conserved. In addition, in O. hyperborea, we could define two alleles, one with the putative initiation codon and the other with ACG. In O. pusilla, despite presenting the putative start codon, CCL14 seems to be a pseudogene due to an insertion of seven base pairs (base pairs 43–49) that leads to a frameshift that disrupts the CDS (Figure 1C). Our results seem to indicate that in the Ochotonidae family subgenus Pika CCL14 is under a process of pseudogenization with inactivation of the gene in some, but not all, species; this pseudogenization, however, does not reflect the taxonomic relationships in the Ochotona genus (Figure 2). Indeed, the subgenus Pika where the mutation in the start codon emerged encompasses species from different lineages that present both the ATG and the ACG.^38,60 Ochotona pusilla is the only species studied where CCL14 is obviously a pseudogene (Figure 1C), and belongs to a different subgenus (Lagotona).^38,60–63

Figure 2.

Evolutionary topology showing the molecular phylogeny within the Ochotonidae family (adapted from Melo-Ferreira et al.⁶⁰). In the AYG found in O. hyperborea, Y is a pyrimidine indicative that this species has two alleles, one with ATG and other with ACG.

Taking into account the CCL14 alterations described in the Ochotona species, we also determined if those alterations could have had any influence in the receptors. For this we sequenced the CCR5 and the CCR1 genes (data not shown). However, these sequences did not reveal alterations that could explain the pseudogenization of CCL14.

In the human genome, pseudogenes are quite common, with high prevalence in multigene families such as chemokines,⁶⁴ but little is known for other mammals. For CCL14, it was previously shown that in Rodentia, CCL14 is absent in rat and mouse of the Muridae family,²⁸ but is a functional gene in the families Dipodidae (Jaculus jaculus), Spalacidae (Nannospalax galili), Sciuridae (Spermophilus tridecemlineatus), Bathyergidae (Heterocephalus glaber), Chinchillidae (Chinchilla laniger) and Caviidae (Cavia porcellus). In addition, in the Trichechus manatus latirostris (Florida manatee) that belongs to the superorder Afrotheria, the CCL14 predicted sequence available in GenBank (accession number XM_004385435.1) seems to be a pseudogene due to an early stop codon, while for the other species within this order for which the CCL14 sequences are available (Loxodonta africana and Echinops telfairi), the gene seems to encode a functional protein. This indicates that although CCL14 is present as a functional gene in the common ancestor of these species, it was later inactivated or deleted in some species, even within the same clade. This suggests that CCL14 might have functions that overlap with the function of other genes, making its presence unnecessary and leading to its pseudogenization.

Conclusions

In most of the mammals we observed that CCL14 encodes a functional gene that is evolving under positive selection, with positively selected codons located in regions important for ligand–receptor binding. In contrast, in two of the Ochotonidae family subgenera (Pika and Lagotona), CCL14 is a pseudogene due to distinct disabling mutations and the detection of an acceleration of the mutation rates. While in some mammals CCL14 retained its biological significance in the Ochotonidae family (subgenera Pika and Lagotona) it seems to have been lost.

Footnotes

Acknowledgements

We would like to thank Jeff Wilcox and Dr Michael Hamilton from Blue Oak Ranch Reserve, University of California, Berkeley, for providing brush rabbit (Sylvilagus bachmani) tissue samples. The authors are grateful to Dr Janet Rachlow, Dr Lisette Waits and Dr Caren Goldberg from Department of Fish and Wildlife Sciences, University of Idaho, USA, for providing Brachylagus idahoensis tissue samples.

Funding

This work was partially funded by FEDER (Fundo Europeu de Desenvolvimento Regional) funds through the Programa Operacional Factores de Competitividade (COMPETE program; FCOMP-01-0124-FEDER-028286) and Portuguese national funds through FCT (Fundação para a Ciência e a Tecnologia; research project PTDC/BIA-ANM/3963/2012) – Quadro de Referência Estratégico Nacional (QREN) funds from the European Social Fund and Portuguese Ministério da Educação e Ciência. FCT also supported the doctoral grants of Fabiana Neves (ref.:SFRH/BD/81916/2011) and the FCT Investigator grant of Joana Abrantes (ref.: IF/01396/2013). ‘Genomics Applied To Genetic Resources’ co-financed by North Portugal Regional Operational Programme 2007/2013 (ON.2 – O Novo Norte), under the National Strategic Reference Framework (NSRF), through the European Regional Development Fund (ERDF), also supported this work.

Conflict of interest

The authors do not have any potential conflicts of interest to declare.

References

Nomiyama

Osada

Yoshie

. The evolution of mammalian chemokine genes. Cytokine Growth Factor Rev 2010; 21: 253–262.

Ding

Zhou

Wang

. Origins of new genes and evolution of their novel functions. Annu Rev Ecol Evol Syst 2012; 43: 345–363.

Levasseur

Pontarotti

. The role of duplications in the evolution of genomes highlights the need for evolutionary-based approaches in comparative genomics. Biol Direct 2011; 6: 11–11.

Nomiyama

Osada

Yoshie

. Systematic classification of vertebrate chemokines based on conserved synteny and evolutionary history. Genes Cells 2013; 18: 1–16.

Kaessmann

. Origins, evolution, and phenotypic impact of new genes. Genome Res 2010; 20: 1313–1326.

Zhang

. Evolution by gene duplication: an update. Trends Ecol Evol 2003; 18: 292–298.

Clark

. Invasion and maintenance of a gene duplication. Proc Natl Acad Sci U S A 1994; 91: 2950–2954.

Lynch

Force

. The probability of duplicate gene preservation by subfunctionalization. Genetics 2000; 154: 459–473.

Zlotnik

Yoshie

. The chemokine superfamily revisited. Immunity 2012; 36: 705–716.

10.

Lalani

Masters

Zeng

. Use of chemokine receptors by poxviruses. Science 1999; 286: 1968–1971.

11.

Masters

Hinek

Uddin

. Poxvirus infection rapidly activates tyrosine kinase signal transduction. J Biol Chem 2001; 276: 48371–48375.

12.

van der Loo

Afonso

de Matos

. Pseudogenization of the MCP-2/CCL8 chemokine gene in European rabbit (genus Oryctolagus), but not in species of Cottontail rabbit (Sylvilagus) and Hare (Lepus). BMC Genet 2012; 13: 72–72.

13.

Abrantes

Esteves

Carmo

. Genetic characterization of the chemokine receptor CXCR4 gene in lagomorphs: comparison between the families Ochotonidae and Leporidae. Int J Immunogenet 2008; 35: 111–117.

14.

Lederman

Sieg

. CCR5 and its ligands: a new axis of evil? Nat Immunol 2007; 8: 1283–1285.

15.

Blain

Kwiatkowski

Zhao

. Structural and functional characterization of CC chemokine CCL14. Biochemistry 2007; 46: 10008–10015.

16.

Munch

Standker

Pohlmann

. Hemofiltrate CC chemokine 1[9–74] causes effective internalization of CCR5 and is a potent inhibitor of R5-tropic human immunodeficiency virus type 1 strains in primary T cells and macrophages. Antimicrob Agents Chemother 2002; 46: 982–990.

17.

Detheux

Standker

Vakili

. Natural proteolytic processing of hemofiltrate CC chemokine 1 generates a potent CC chemokine receptor (CCR)1 and CCR5 agonist with anti-HIV properties. J Exp Med 2000; 192: 1501–1508.

18.

Abrantes

Carmo

Matthee

. A shared unusual genetic change at the chemokine receptor type 5 between Oryctolagus, Bunolagus and Pentalagus. Conserv Genet 2011; 12: 325–330.

19.

Carmo

Esteves

Ferrand

van der Loo

. Genetic variation at chemokine receptor CCR5 in leporids: alteration at the 2nd extracellular domain by gene conversion with CCR2 in Oryctolagus, but not in Sylvilagus and Lepus species. Immunogenetics 2006; 58: 494–501.

20.

Esteves

Abrantes

van der Loo

. Extensive gene conversion between CCR2 and CCR5 in domestic cat (Felis catus). Int J Immunogenet 2007; 34: 321–324.

21.

Perelygin

Zharkikh

Astakhova

. Concerted evolution of vertebrate CCR2 and CCR5 genes and the origin of a recombinant equine CCR5/2 gene. J Hered 2008; 99: 500–511.

22.

Shields

. Gene conversion among chemokine receptors. Gene 2000; 246: 239–245.

23.

Vazquez-Salat

Yuhki

Beck

. Gene conversion between mammalian CCR2 and CCR5 chemokine receptor genes: a potential mechanism for receptor dimerization. Genomics 2007; 90: 213–224.

24.

de Matos

Lanning

Esteves

. Genetic characterization of CCL3, CCL4 and CCL5 in leporid genera Oryctolagus, Sylvilagus and Lepus. Int J Immunogenet 2014; 41: 154–158.

25.

Islam

Chang

Colvin

. Mouse CCL8, a CCR8 agonist, promotes atopic dermatitis by recruiting IL-5+ T(H)2 cells. Nat Immunol 2011; 12: 167–177.

26.

Schulz-Knappe

Magert

Dewald

. HCC-1, a novel chemokine from human plasma. J Exp Med 1996; 183: 295–299.

27.

Richter

Casarosa

Standker

. Significance of N-terminal proteolysis of CCL14a to activity on the chemokine receptors CCR1 and CCR5 and the human cytomegalovirus-encoded chemokine receptor US28. J Immunol 2009; 183: 1229–1237.

28.

Shibata

Nomiyama

Yoshie

Tanase

. Genome diversification mechanism of rodent and Lagomorpha chemokine genes. Biomed Res Int 2013; 2013: 856265–856265.

29.

D'Ambrosio

Panina-Bordignon

Sinigaglia

. Chemokine receptors in inflammation: an overview. J Immunol Methods 2003; 273: 3–13.

30.

Tsou

Gladue

Carroll

. Identification of C-C chemokine receptor 1 (CCR1) as the monocyte hemofiltrate C-C chemokine (HCC)-1 receptor. J Exp Med 1998; 188: 603–608.

31.

Lolis

Murphy

The structural biology of chemokines. In: Harrison

Lukacs

(eds). The receptors: the chemokine receptors, Totowa, NJ: Humana Press Inc, 2007, pp. 9–30.

32.

Pakianathan

Kuta

Artis

. Distinct but overlapping epitopes for the interaction of a CC-chemokine with CCR1, CCR3 and CCR5. Biochemistry 1997; 36: 9642–9648.

33.

Rajagopalan

Rajarathnam

. Structural basis of chemokine receptor function—a model for binding affinity and ligand selectivity. Biosci Rep 2006; 26: 325–339.

34.

Steen

Larsen

Thiele

Rosenkilde

. Biased and g protein-independent signaling of chemokine receptors. Front Immunol 2014; 5: 277–277.

35.

Savino

Borroni

Torres

. Recognition versus adaptive up-regulation and degradation of CC chemokines by the chemokine decoy receptor D6 are determined by their N-terminal sequence. J Biol Chem 2009; 284: 26207–26215.

36.

Nibbs

Wylie

Pragnell

Graham

. Cloning and characterization of a novel murine beta chemokine receptor, D6. Comparison to three other related macrophage inflammatory protein-1alpha receptors, CCR-1, CCR-3, and CCR-5. J Biol Chem 1997; 272: 12495–12504.

37.

Matthee

van Vuuren

Bell

Robinson

. A molecular supermatrix of the rabbits and hares (Leporidae) allows for the identification of five intercontinental exchanges during the Miocene. Syst Biol 2004; 53: 433–447.

38.

Lissovsky

. Taxonomic revision of pikas Ochotona (Lagomorpha, Mammalia) at the species level. Mammalia 2014; 78: 199–216.

39.

Horner

Lefkimmiatis

Reyes

. Phylogenetic analyses of complete mitochondrial genome sequences suggest a basal divergence of the enigmatic rodent Anomalurus. BMC Evol Biol 2007; 7: 16–16.

40.

Murphy

Eizirik

Johnson

. Molecular phylogenetics and the origins of placental mammals. Nature 2001; 409: 614–618.

41.

Edgar

. MUSCLE: multiple sequence alignment with high accuracy and high throughput. Nucleic Acids Res 2004; 32: 1792–1797.

42.

Yang

. PAML 4: phylogenetic analysis by maximum likelihood. Mol Biol Evol 2007; 24: 1586–1591.

43.

Tamura

Stecher

Peterson

. MEGA6: Molecular Evolutionary Genetics Analysis version 6.0. Mol Biol Evol 2013; 30: 2725–2729.

44.

Yang

. Inference of selection from multiple species alignments. Curr Opin Genet Dev 2002; 12: 688–964.

45.

Yang

Nielsen

Goldman

Pedersen

. Codon-substitution models for heterogeneous selection pressure at amino acid sites. Genetics 2000; 155: 431–449.

46.

Delport

Poon

Frost

Kosakovsky Pond

. Datamonkey 2010: a suite of phylogenetic analysis tools for evolutionary biology. Bioinformatics 2010; 26: 2455–2457.

47.

Poon

Frost

Pond

. Detecting signatures of selection from DNA sequences using Datamonkey. Methods Mol Biol 2009; 537: 163–183.

48.

Murrell

Moola

Mabona

. FUBAR: a fast, unconstrained bayesian approximation for inferring selection. Mol Biol Evol 2013; 30: 1196–1205.

49.

Murrell

Wertheim

Moola

. Detecting individual sites subject to episodic diversifying selection. PLoS Genet 2012; 8: e1002764–e1002764.

50.

Areal

Abrantes

Esteves

. Signatures of positive selection in Toll-like receptor (TLR) genes in mammals. BMC Evol Biol 2011; 11: 368–368.

51.

Neves

Abrantes

Steinke

Esteves

. Maximum-likelihood approaches reveal signatures of positive selection in IL genes in mammals. Innate Immun 2014; 20: 184–191.

52.

Lemos de Matos

McFadden

Esteves

. Positive evolutionary selection on the RIG-I-like receptor genes in mammals. PloS One 2013; 8: e81864–e81864.

53.

Lemos de Matos

McFadden

Esteves

. Evolution of viral sensing RIG-I-like receptor genes in Leporidae genera Oryctolagus, Sylvilagus, and Lepus. Immunogenetics 2014; 66: 43–52.

54.

Pinheiro

Woof

Abi-Rached

. Computational analyses of an evolutionary arms race between mammalian immunity mediated by immunoglobulin A and its subversion by bacterial pathogens. PloS One 2013; 8: e73934–e73934.

55.

Pinheiro

Woof

Almeida

. Leporid immunoglobulin G shows evidence of strong selective pressure on the hinge and CH3 domains. Open Biol 2014; 4: 140088–140088.

56.

Librado

Rozas

. DnaSP v5: a software for comprehensive analysis of DNA polymorphism data. Bioinformatics 2009; 25: 1451–1452.

57.

Hall

. BioEdit: a user-friendly biological sequence alignment editor and analysis program for Windows 95/98/NT. Nucl Acids Symp Ser 1999; 41: 95–98.

58.

Tajima

. Simple methods for testing the molecular evolutionary clock hypothesis. Genetics 1993; 135: 599–607.

59.

Neves

Abrantes

Pinheiro

. Convergent evolution of IL-6 in two leporids (Oryctolagus and Pentalagus) originated an extended protein. Immunogenetics 2014; 66: 589–595.

60.

Melo-Ferreira

de Matos

Areal

. The phylogeny of pikas (Ochotona) inferred from a multilocus coalescent approach. Mol Phylogenet Evol 2015; 84C: 240–244.

61.

Fostowicz-Frelik

Frelik

. The earliest occurrence of the steppe pika (Ochotona pusilla) in Europe near the Pliocene/Pleistocene boundary. Naturwissenschaften 2010; 97: 325–329.

62.

Fostowicz-Frelik

Frelik

Gasparik

. Morphological phylogeny of pikas (Lagomorpha: Ochotona), with a description of a new species from the Pliocene/Pleistocene transition of Hungary. Proc Natl Acad Nat Sci Phila 2010; 159: 97–118.

63.

Lanier

Olson

. Inferring divergence times within pikas (Ochotona spp.) using mtDNA and relaxed molecular dating techniques. Mol Phylogenet Evol 2009; 53: 1–12.

64.

Ondrej P and Jianzhi Z (2010) Pseudogenes and their evolution. In: eLS. John Wiley & Sons Ltd, Chichester. http://www.els.net. DOI: 10.1002/9780470015902.a0005118.pub2.