Atypical chemokine receptors (ACKRs) are a subclass of G protein-coupled receptors characterized by promiscuity of ligand binding and an obvious inability to signal after ligand binding. Although some discoveries regarding this family in Homo sapiens and other species have been reported in some studies, the evolution and function of multiple ACKR in mammals have not yet been clearly understood. We performed an evolutionary analysis of ACKR genes (ACKR1, ACKR2, ACKR3, and ACKR4) in mammals. Ninety-two full-length ACKR genes from 27 mammal species were retrieved from the Genbank and Ensemble databases. Phylogenetic analysis showed that there were four well-conserved subfamilies in mammals. Synteny analysis revealed that ACKR genes formed conserved linkage groups with their adjacent genes across mammalian species, facilitating the identification of ACKRs in as yet unannotated genome datasets. Analysis of the site-specific profiles established by posterior probability revealed the positive-selection sites to be distributed mainly in the ligand binding region of ACKR1. This study highlights the molecular evolution of the ACKR gene family in mammals and identifies the critical amino acid residues likely to be relevant to ligand binding. Further experimental verification of these findings may provide valuable information regarding the ACKR’s biochemical and physiological functions.
National Key Research and Development Program of China2016YFD05015032017YFD0501100Key Technology R&D Program of Gansu Province1604NKCA045-2National Pig Industrial SystemCARS-36-06B1. Introduction
The chemokine system plays an important role in mammalian immunity, which can guide immune effector cells to sites of infection or inflammation and coordinate interactions between immune cells. The chemokine family is classified into four main subfamilies (XC, CC, CXC, CX3C) based on the positioning of their initial cysteine residues (here C indicates cysteine and X/X3 indicates one or there non-cysteine amino acids) [1]. All of chemokines exert their functions by interacting with chemokine receptors that are selectively expressed on the surfaces of their target immune cells [2]. Chemokine receptors belong to the largest rhodopsin family of G protein-coupled receptors (GPCRs) and they consist of seven transmembrane domains and multiple extracellular and intracellular loops. These are involved in ligand binding and signaling [3]. Chemokine receptors are mainly divided into four subfamilies depending on the type of chemokine they bind, including CC chemokine receptors (CCRs), CXC chemokine receptors (CXCRs), XC chemokine receptors (XCRs), and CX3C chemokine receptors (CX3CRs) [4]. Apart from signaling receptors, chemokine receptors also include some atypical receptors. These are similar in structure to the conventional receptors, but lack an intracellular motif required for signaling [5]. The conventional chemokine receptors represent a larger subgroup (eighteen in humans) of G protein coupled leukocyte chemotactic receptors, and the atypical chemokine receptors represent a smaller subgroup (four in humans) of chemotactic receptors that do not transduce signals through G proteins and lack chemotactic activity [6]. Atypical chemokine receptors mainly include four types of receptors and named atypical chemokine receptor 1 (ACKR1), ACKR2, ACKR3, and ACKR4 on the basis of the new standardized nomenclature system [7]. Members of ACKR have been previously referred to by scientists in this filed as ‘chemokine-binding proteins’, ‘scavengers’, ‘interceptors’, and ‘decoys’, which has caused confusion and has delayed recognition of these molecules as a functionally related groups. For instance, ACKR1 is previously known as Duffy antigen, Fy antigen, and CD234. ACKR2 is previously known as D6, CCBP2, CCR9, and CMKBR9. ACKR3 is previously known as CXCR7 and RDC1. ACKR4 is previously known as CCRL1 and CCX-CKR.
Atypical chemokine receptors are characterized by promiscuity of ligand binding and an apparent inability to signal after ligand binding. This inability to signal is mainly because of alterations in the canonical DRY motif located in the second intracellular loop of the typical chemokine receptors [8]. This motif is responsible for G protein coupling to the receptors and its absence is an essential unifying theme of the atypical receptors [9, 10]. These nonsignaling receptors efficiently internalize their cognate chemokines and act as chemokine scavengers. This is achieved by continuous internalization and constant surface expression of the atopic receptors in a β-arrestin signaling dependent manner [11–13].
ACKR1 is mainly expressed in red blood cells and endothelial cells and has very little sequence similarity to other atypical receptors or chemokine receptors [14, 15]. ACKR genes exist in a specific chromosomal location, unlike those of the other receptors [15]. The ligand binding profile of ACKR1 includes several inflammatory chemokines, including CXCL1, CXCL5-9, CXCL11, and CXCL13, suggesting that ACKR1 plays an important role in maintaining chemokine receptors in the blood [16, 17]. ACKR2 is expressed in many different tissues including those of the skin, gut, and lung [18]. Within these classical barrier tissues, expression is largely confined to lymphatic endothelial cells, with no expression detected in blood endothelial cells [18]. The ligand binding profile includes CCL2-5, CCL7-8, CCL11, CCL13-14, CCL17, and CCL22 [19–24]. ACKR2, which serves as a constitutively internalizing and recycling receptor, is also capable of internalizing and degrading the chemokines of its binding profile [25]. ACKR3 is the best characterized receptor in the ACKR family and it is mainly expressed in hematopoietic cells, neurons, mesenchymal cells, endothelial cells, and cancer cells. The ligands of ACKR3 contain CXCL11 and CXCL12. Binding of ACKR3 to CXCL11, which is an IFN-γ-induced chemokine, is involved in the regulation of lymphocyte migration [26]. ACKR3 also forms a heterodimer with CXCR4 and acts as a scavenger for the ligand CXCL12 and plays key roles in organ development and tumor development and progression [27]. ACKR4 is an atypical receptor for homeostatic CC and CXC chemokines including CCL19, CCL21, CCL25, and CXCL13 [28]. Like ACKR3, ACKR4 is capable of internalizing its ligands and targeting them for intracellular degradation, and it does so in much the same way. ACKR4 is also able to antagonize CXCR3-induced chemokines through heterodimer formation with the CXCR3 [29].
Two other ACKRs, CCRL2 (ACKR5) and PITPNM3 (ACKR6), have been proposed, but functional confirmation is pending, so they have been provisionally assigned ACKR designations [30].
Compared with numerous studies on typical chemokine receptors, especially in fish and mammals, little is known about ACKRs in mammals [30]. The increasing wealth of sequence data available from sequenced genome databases has allowed researchers to perform evolutionary analyses of ACKRs in mammals. In the present study, we performed an evolutionary analysis of 92 full-length ACKR (ACKR1, ACKR2, ACKR3, and ACKR4) genes from 27 mammal species retrieved from the Genbank and Ensemble databases.
2. Materials and Methods2.1. Acquisition of ACKR Family Sequences and Synteny Analysis
For some well-annotated genomes, the amino acid sequences of ACKRs (ACKR1, ACKR2, ACKR3, and ACKR4) were retrieved directly from the Genbank and Ensemble databases. PSI-BLAST was performed to search these databases using Homo sapiens ACKRs (NP_001136269, NP_001008540, NP_001707, and NP_006555) as query sequences. The predicted coding sequences of the best hits were collected when the hits presented more than 70% in length and 50% in identity were aligned with the query sequence (with E values < e−10). These settings distinguished the potential ACKR members from different species but avoid involving other chemokine receptors effectively. After removal of redundant and incomplete sequences, the initial data set (S1 Table) for ACKR contained 92 protein sequences from 27 mammals. Because the functional confirmation of the two provisionally assigned ACKRs (ACKR5 and ACKR6) is pending, they were excluded from analysis in the present study.
Synteny analysis was conducted using the GENOMICUS v80.01 browser, which allows integration of the data available on the Ensemble database to provide a better visualization of conserved synteny blocks and to facilitate reconstruction of the organization of ancient genomes [31, 32]. Genes not annotated on the GENOMICUS browser were searched within the respective species by BLASTP and TBLASTN over the Genbank and Ensemble databases.
2.2. Sequence Alignment and Phylogenetic Analysis
A codon-based coding sequence alignment was constructed using MUSCLE with default parameters and manually adjusted using MEGA 6 [33] and viewed and edited in Jalview 2.0 [34]. The alignment was subsequently processed using Gblocks v0.91b [35] for phylogenetic reconstruction with default parameters. To access the selective pressures acting on the four mammals ACKR subfamilies, seven different alignments were produced: one for each paralog and a seventh with all sequences excluding outgroups. The substitution model that best fit the dataset was selected using Akaike Information Criterion (AIC) implemented in ProTest 3.2 [36], starting with 14 substitution matrices and using the fixed BIONJ tree for likelihood calculations. The phylogeny was estimated using the Maximum Likelihood (ML) methods. The ML phylogenetic tree was constructed in PhyML 3.0 [37], with 1000 bootstrap replicates and the NNI branch search algorithm. Finally, the phylogenetic trees were displayed using TreeView [38]. Besides, the neighbor joining (NJ), minimal evolution (ME) and maximum parsimony (MP) methods were used individually to reconstruct another three phylogenetic trees with MEGA v3.1 from the Gblocks alignment.
2.3. Codon-Based Analyses of Positive Selection
The selective pressures acting on coding region were evaluated across the phylogeny using a phylogenetic-based ML analysis. Accurate nucleotide sequences and related amino acid sequence alignments were retrieved with PAL2NAL [39], a program that constructs multiple codon alignments form matching protein sequences. The codeml program in PAML4.5 [40] was used to estimate the rates of synonymous (dS) and nonsynonymous substitution (dN) and the dN/dS ratio (omega, ω). ω>1 indicates positive selection, ω<1 indicates negative selection, and ω=1 is neutrality. Accurate nucleotide sequence alignments were constructed from matching related protein sequence with MUSCLE (MEGA 6). Then the resulting codon alignments and ML tree were used in the codeml program. The site-specific models were tested: Models M0 (one ratio), M1a (nearly neutral), M2a (positive selection), M3 (discrete), M7 (beta), and M8 (beta+ω) were all used in this analysis [41, 42]. Model M0 assumes one ratio for all sites. M1a presupposes a proportion p0 of conserved sites with ω<1 and p1=1-p0 of neutral sites with ω=1. M2a adds an additional class of sites with the frequency p2=1-p0-p1, and ω2 is estimated form the data. In the M3, the probabilities (p0, p1, and p2) of each site being submitted to purifying, neutral, and positive selection, respectively, and their corresponding ω ratios (ω0, ω1, and ω2) are inferred from the data. M7 and M8 assume a β-distribution for ω between 0 and 1, and M8 adds one extra class with the same ratio ω1. Subsequent likelihood rate comparisons of M0 and M3, M1a with M2a, and M7 with M8 were performed to determine which model fit the data best. The LRT was used to test positive selection of the two pairs of site model [43, 44]. Finally, the BEB approach was used to calculate the posterior probability that each site would belong to the site class of positive selection under each model [45].
3. Results3.1. Identification and Distribution of ACKR Genes across Mammals
The final data sets contained 92 ACKR gene sequences from 27 representative species of mammals, including two primates, two rodents, one monotreme, one hyracoidean, one edentate, one pilosa, two perissodactyls, three artiodactyls, two carnivorans, one cetacean, one proboscid, two erinaceidae, one insectivoran, two lagomorphs, one dasyuromorph, one scandentia, one diprotodont, and two chiropterans. The results of genomic database searches showed that the majority of mammals investigated in this study (from orders Primate, Rodentia, Edentata, Perissodactyla, Carnivora, Proboscidea, Erinaceidae, Insectivora, Lagomorpha, Dasyuromorphia, Chiroptera, Artiodactylas) possess 4 members of ACKR family. No ACKR genes were identified in the representative of Pilosa, Choloepus hoffmanni. The remaining mammals only possessed 2 ACKRs and all of these species lacked the ACKR1 gene (S1 Table).
3.2. Synteny Analysis of ACKR Genes in Mammal Genomes
As shown in the additional file 1, several ACKR genes could not be identified in some mammal genomes using the sequence collection method. Synteny analysis was performed to determine why some ACKR genes were missed. We observed that the ACKR1 gene formed a conserved linkage group with AIM2, CADM3, and FCER1A genes in the most mammal genomes (Figure 1(a)). In the genome databases of C. hoffmanni, only two genes of the conserved linkage group, the AIM2 and CADM3, were found. These were located in Scaffold_33705 and Scaffold_ 5395, respectively. CADM3 and FCER1A, which were in the conserved group, were identified in the genomes of Tupaia belangeri but AIM2 and ACKR1 were not. In the genome of Ornithorhynchus anatinus, no members of the conserved linkage group were found. ACKR2 formed a conserved linkage group with CCDC13, HIGD1A and CYP8B1 in mammalian genomes (Figure 1(b)). However, no ACKR2 genes were found in the genomes of C. hoffmanni or O. anatinus. Only CCDC13 and HIGD1A of the conserved gene group were identified. The conserved ACKR3-specific gene group consisted of ASB18, IQCA, ACKR3, COPS8, and COL6A3 in mammal genomes (Figure 1(c)). Nevertheless, neither ASB18 nor ACKR3 were found in the genome of C. hoffmanni or T. belangeri. In mammalian genomes, the ACKR4-specific conserved linkage group was composed of ACPP, DNAJC13, ACAD11, ACKR4, UBA5 and NPHP5 (Figure 1(d)). Among these genes, AKCR4, ACAD11, UBA5 and NAHP shared the same transcript: “ACAD11-NPHP5”. The ACAD11 and ACKR4 of this conserved gene group were absent from the genomes of C. hoffmanni and T. belangeri.
Gene synteny analysis of ACKR1-4 in representative species of mammals. The gene abbreviations were taken from the Ensemble database. The directions of arrows indicate the transcriptional orientation of genes.
3.3. Phylogenetic Analysis of ACKR Genes in Mammals
After the exclusion of partial and unfinished sequences, 92 sequences were retrieved from 27 mammal species. To determine the phylogenetic relationship of mammal ACKR genes, a rooted ML phylogenetic tree was constructed based on amino acids alignment under the best-fit model JTT+I+G+F. Here, the best-fit model (JTT+I+G+F) for amino acid substitution was selected by ProTest3.2 with discrete gamma distribution in four categories. All parameters (gamma shape = 1.687; proportion of invariants = 0.042) were estimated from the dataset. Tree topology was assessed using MEGA 6 with neighbor joining (NJ), minimum evolution (ME) and maximum parsimony (MP) methods, and it was found to be substantially similar to the ML tree (data not shown). Using the relaxin receptor from Ciona intestinalis as the outgroup of mammal ACKR genes [30], the ML trees showed the ACKRs of mammals to be grouped into four lineages: ACKR1 subfamily, ACKR2 subfamily, ACKR3 subfamily and ACKR4 subfamily (Figure 2). Our data suggested that two major duplications had occurred in mammal lineages. The first duplication led to the emergence of two lineages that evolved into ACKR1 and the ancestor of ACKR2, ACKR3 and ACKR4. The second duplication led to divergence of ACKR2, ACKR3 and ACKR4.
A maximum likelihood tree of mammal ACKRs. The tree was constructed using maximum likelihood method. The number indicates bootstrap values. See additional file 1 for details of gene accession numbers and species abbreviations.
3.4. Adaptive Evolution of ACKR Genes in Mammals
To detect signatures of adaptive evolution over the ACKR1, ACKR2, ACKR3, and ACKR4 codon sequences, four smaller phylogenetic trees were built for each group and the topology used for each site-specific model was implanted using the codeml program of PAML v4.0 package. Parameter estimates and log-likelihood values under model of variable ω ratios among sites were shown in Table 1. In all cases, the LRT did not differ significantly between M1a and M2a, but the LRTs did show significant differences between M0 and M3 and between M7 and M8 for all receptors except the ACKR4 lineage, indicating that M3 and M8 fit the data better. However, no selected sites were detectable in M3. In model M8, one site (154 G) from ACKR1 lineage was found to be a positively selected site, showing a P-value over 99%.
Likelihood ratio test statistics for models of variable selective pressure among sites.
Gene
Model
Parameters
Models compared
LRT(2Δl)
P-value
df
Positively selected sites (BEB)
ACKR1
M0
ω0=0.43574
M0 vs M3
534.02
0.0000
4
M3
ω0=0.03476 ω1= 0.35759ω2= 1.15972
p0= 0.23127 p1= 0.47276p2= 0.29597
M1a
ω0=0.16327 ω1=1.00000
M1a vs M2a
0
1.0000
2
p0=0.55005 p1=0.44995
M2a
ω0=0.16327 ω1= 1.00000ω2=1.00000
p0=0.55005 p1=0.28372p2=0.16623
M7
p =0.53074 q=0.60339
M7 vs M8
21.26
0.0000
2
154 G∗∗
M8
p0=0.85617 p=0.68793q=1.15451
p1=0.14383 ω=1.45019
ACKR2
M0
ω0=0.20556
M0 vs M3
652.48
0.0000
4
M3
ω0=0.01987 ω1=0.21352ω2= 0.77475
p0= 0.34607 p1=0.48100p2= 0.17293
M1a
ω0=0.12737 ω1=1.00000
M1a vs M2a
0
1.0000
2
p0=0.79265 p1=0.20735
M2a
ω0=0.12737 ω1= 1.00000ω2=1.00000
p0=0.79265 p1=0.18741p2=0.01994
M7
p =0.52179 q=1.61923
M7 vs M8
8.88
0.0117
2
M8
p0=0.97423 p=0.59529q=2.13409
p1=0.02577 ω=1.39256
ACKR3
M0
ω0=0.04691
M0 vs M3
569.64
0.0000
4
M3
ω0=0.00786 ω1=0.12501ω2= 0.49733
p0= 0.71583 p1=0.24584p2=0.03833
M1a
ω0=0.03806 ω1=1.00000
M1a vs M2a
0
1.0000
2
p0=0.95014 p1=0.04986
M2a
ω0=0.03806 ω1= 1.00000ω2=28.52785
p0=0.95014 p1=0.04986p2=0.00000
M7
p =0.24148 q=3.44384
M7 vs M8
7.78
0.0204
2
M8
p0=0.98766 p=0.27944q=5.15904
p1=0.01234 ω=1.00000
ACKR4
M0
ω0=0.12972
M0 vs M3
416.55
0.0000
4
M3
ω0=0.03081 ω1=0.32186ω2=1.16947
p0= 0.65355 p1=0.33146p2=0.01499
M1a
ω0=0.08280 ω1=1.00000
M1a vs M2a
0
1.0000
2
p0=0.84629 p1=0.15371
M2a
ω0=0.08280 ω1= 1.00000ω2=31.41467
p0=0.84629 p1=0.15371p2=0.00000
M7
p =0.41473 q=2.39665
M7 vs M8
5.54
0.0626
2
139 V∗
M8
p0=0.99543 p=0.43468q=2.65383
p1=0.00457 ω=1.77610
4. Discussion
Chemokines are important regulators of leukocyte migration and play key roles in diverse physiological and pathological immune and inflammatory contexts [28]. In addition to the typical signaling chemokine receptors, a recently discovered subclass of atypical chemokine receptors are characterized by promiscuity of ligand binding and an obvious inability to signal after ligand binding [46]. The inability to signal is largely a consequence of alterations in the canonical DRY motif in the second intracellular loop of the typical chemokine receptors [47]. The motif is responsible for G-protein coupling to the receptors and its absence is the key unifying theme of these atypical receptors [8]. The DRY motif of the ACKRs was identified using multiple sequence alignment. No DRY motif was found in the ACKR1 subgroup. The DKYLEIV motif, DRYLSVT motif, and DRYWAVT motif were identified in the mammal ACKR2, ACKR3, and ACKR4 subgroups, respectively. As in the DRY motif of typical chemokine receptors, the last three amino acids were essential to maintaining the function of signal transduction [8, 48] (Figure 3).
Sequence alignment of the DRY motif in the ACKR2, ACKR3, and ACKR4 proteins of mammals. Multiple alignments were performed using the full length protein sequences with the Jalview software. Identical amino acids are indicated by asterisks whereas those with high or low similarity are indicated by “:” and “.” respectively.
The number of the ACKR genome loci varied across several mammalian genomes. Synteny analysis was performed to determine the reason for the absence of some ACKR genes. ACKR genes formed conserved linkage groups with their adjacent genes across mammalian genomes. The genome sequence datasets of O. anatinus, P. capensis, and C. hoffmanni available in Genbank and Ensemble databases were limited and presented in scaffold form. These sequence data did not meet the requirements for assembly into chromosomes. Partial segments of the conserved ACKR-specific blocks of genes were here identified. In this way, the absence of ACKRs from some of the mammalian genomes investigated here may be attributed to the incomplete information available in genome databases rather than to gene loss during evolution.
The positively selective sites 154 G of ACKR1 is located within the extracellular domain between the fourth trans-membrane and the fifth trans-membrane (Figure 4). This region is responsible for the direct interaction between ACKR1 and the ligands [15]. ACKR1 is mainly expressed in red blood cell. ACKR1 serves as the chemokine buffer for the blood, and it can bind to many different chemokines. Increasing amounts of evidence have shown that ACKR1 possesses a larger ligand binding profile than the other ACKRs [49, 50]. The positively selective site within the binding region of ACKR1 may provide direct evidence for extended ligand binding profile.
The ACKR1 receptors do not signal in response to ligand binding but rather internalize the ligand and target it for intracellular degradation. The key to the lack of signaling is the absence of DRY motif in the second intracellular loop of ACKR1. The positively selective site was located in the second extracellular loop, which is responsible for ligand recognition and binding.
Data Availability
The sequences of ACKR family analyzed in this study are deposited in the NCBI and Ensemble databases. The accession number is listed in Supplementary file 1.
Conflicts of Interest
The authors declare that they have no conflicts of interest.
Authors’ Contributions
Li Pan and Yongguang Zhang designed and wrote the manuscript; Jianliang Lv and Zhongwang Zhang conducted the bioinformatic analyses. All authors read and approved the final manuscript.
Acknowledgments
This work was supported by grants from the National Key Research and Development Program of China (grant numbers 2016YFD0501503 and 2017YFD0501100), the Key Technology R&D Program of Gansu Province of China (grant number 1604NKCA045-2), and the National Pig Industrial System (CARS-36-06B).
Supplementary Materials
Table S1. List of species used for the Multiple Sequences Alignment (MSA).
ZlotnikA.YoshieO.NomiyamaH.The chemokine and chemokine receptor superfamilies and their molecular evolution2006712, article 24310.1186/gb-2006-7-12-2432-s2.0-34249707951Mélik-ParsadaniantzS.RostèneW.Chemokines and neuromodulation20081981-2626810.1016/j.jneuroim.2008.04.0222-s2.0-47249163193CharoI. F.RansohoffR. M.The many roles of chemokines and chemokine receptors in inflammation2006354661062110.1056/NEJMra0527232-s2.0-32144454172ZlotnikA.YoshieO.The chemokine superfamily revisited201236570571210.1016/j.immuni.2012.05.0082-s2.0-84861426163SokolC. L.LusterA. D.The chemokine system in innate immunity20157510.1101/cshperspect.a016303BachelerieF.GrahamG. J.LocatiM.MantovaniA.MurphyP. M.NibbsR.RotA.SozzaniS.ThelenM.New nomenclature for atypical chemokine receptors20141532072082-s2.0-8489445627310.1038/ni.2812GrahamG. J.LocatiM.MantovaniA.RotA.ThelenM.The biochemistry and biology of the atypical chemokine receptors20121451-230382-s2.0-8486218233810.1016/j.imlet.2012.04.004CancellieriC.VacchiniA.LocatiM.BonecchiR.BorroniE. M.Atypical chemokine receptors: From silence to sound20134112312362-s2.0-8487310497110.1042/BST2012024623356288KraftK.OlbrichH.MajoulI.MackM.ProudfootA.OppermannM.Characterization of sequence determinants within the carboxyl-terminal domain of chemokine receptor CCR5 that regulate signaling and receptor internalization200127637344083441810.1074/jbc.M1027822002-s2.0-0035860777DamajB. B.McCollS. R.NeoteK.SongqingN. A.OgbornK. T.HébertC. A.NaccacheP. H.Identification of G-protein binding sites of the human interleukin-8 receptors by functional mapping of the intracellular loops19961012142614342-s2.0-0029966772890351310.1096/fasebj.10.12.8903513GallieraE.JalaV. R.TrentJ. O.BonecchiR.SignorelliP.LefkowitzR. J.MantovaniA.LocatiM.HaribabuB.β-arrestin-dependent constitutive internalization of the human chemokine decoy receptor D620042792425590255972-s2.0-294256760610.1074/jbc.M40036320015084596WuF. Y.FanJ.TangL.ZhaoY. M.ZhouC. C.Atypical chemokine receptor D6 inhibits human non-small cell lung cancer growth by sequestration of chemokines20136191952-s2.0-8487833210210.3892/ol.2013.1358BorroniE. M.CancellieriC.VacchiniA.BenureauY.LaganeB.BachelerieF.Arenzana-SeisdedosF.MizunoK.MantovaniA.BonecchiR.LocatiM.β-arrestin-dependent activation of the cofilin pathway is required for the scavenging activity of the atypical chemokine receptor D620136273S31S3310.1126/scisignal.20048372-s2.0-84877591484RotA.Contribution of Duffy antigen to chemokine function20051666876942-s2.0-2714444611010.1016/j.cytogfr.2005.05.01116054417GrahamG. J.D6 and the atypical chemokine receptor family: Novel regulators of immune and inflammatory processes20093923423512-s2.0-6074912523110.1002/eji.20083885819130487HansellC. A. H.HursonC. E.NibbsR. J. B.DARC and D6: Silent partners in chemokine regulation?20118921972062-s2.0-7995192418310.1038/icb.2010.147GardnerL.PattersonA. M.AshtonB. A.StoneM. A.MiddletonJ.The human Duffy antigen binds selected inflammatory but not homeostatic chemokines200432123063122-s2.0-324275430310.1016/j.bbrc.2004.06.14615358176NibbsR. J. B.KriehuberE.PonathP. D.ParentD.QinS.CampbellJ. D. M.HendersonA.KerjaschkiD.MaurerD.GrahamG. J.RotA.The β-chemokine receptor D6 is expressed by lymphatic endothelium and a subset of vascular tumors200115838678772-s2.0-00351092021123803610.1016/S0002-9440(10)64035-7NibbsR. J. B.WylieS. M.PragnellI. B.GrahamG. J.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-5199727219124951250410.1074/jbc.272.19.124952-s2.0-0030924949NibbsR. J. B.WylieS. M.YangJ.LandauN. R.GrahamG. J.Cloning and characterization of a novel promiscuous human β-chemokine receptor D619972725132078320832-s2.0-003143480910.1074/jbc.272.51.320789405404BoniniJ. A.SteinerD. F.Molecular cloning and expression of a novel rat CC-chemokine receptor (rCCR10rR) that binds MCP-1 and MIP-1β with high affinity1997169102310302-s2.0-003141443210.1089/dna.1997.16.10239324304BonecchiR.LocatiM.GallieraE.VulcanoM.SironiM.FraA. M.GobbiM.VecchiA.SozzaniS.HaribabuB.Van DammeJ.MantovaniA.Differential recognition and scavenging of native and truncated macrophage-derived chemokine (macrophage-derived chemokine/CC Chemokine Ligand 22) by the D6 decoy receptor200417284972497610.4049/jimmunol.172.8.49722-s2.0-11144354216SavinoB.BorroniE. M.TorresN. M.ProostP.StruyfS.MortierA.MantovaniA.LocatiM.BonecchiR.Recognition versus adaptive up-regulation and degradation of CC chemokines by the chemokine decoy receptor D6 are determined by their N-terminal sequence20092843826207262152-s2.0-7035001533310.1074/jbc.M109.02924919632987HansellC. A. H.SchieringC.KinstrieR.FordL.BordonY.McInnesI. B.GoodyearC. S.NibbsR. J. B.Universal expression and dual function of the atypical chemokine receptor D6 on innate-like B cells in mice201111720541354242-s2.0-7995626575110.1182/blood-2010-11-31711521450903BachelerieF.Ben-BaruchA.BurkhardtA. M.CombadiereC.FarberJ. M.GrahamG. J.HorukR.Sparre-UlrichA. H.LocatiM.LusterA. D.MantovaniA.MatsushimaK.MurphyP. M.NibbsR.NomiyamaH.PowerC. A.ProudfootA. E. I.RosenkildeM. M.RotA.SozzaniS.ThelenM.YoshieO.ZlotnikA.International union of pharmacology. LXXXIX. Update on the extended family of chemokine receptors and introducing a new nomenclature for atypical chemokine receptors20146611792-s2.0-8488751297010.1124/pr.113.007724NaumannU.CameroniE.PruensterM.MahabaleshwarH.RazE.ZerwesH.-G.RotA.ThelenM.CXCR7 functions as a scavenger for CXCL12 and CXCL1120105210.1371/journal.pone.00091752-s2.0-77949482428LevoyeA.BalabanianK.BaleuxF.BachelerieF.LaganeB.CXCR7 heterodimerizes with CXCR4 and regulates CXCL12-mediated G protein signaling2009113246085609310.1182/blood-2008-12-1966182-s2.0-67650349084NibbsR. J. B.GrahamG. J.Immune regulation by atypical chemokine receptors2013131181582910.1038/nri35442-s2.0-84886768227VinetJ.Van ZwamM.DijkstraI. M.BrouwerN.Van WeeringH. R. J.WattsA.MeijerM.FokkensM. R.KannanV.VerzijlD.VischerH. F.SmitM. J.LeursR.BiberK.BoddekeH. W. G. M.Inhibition of CXCR3-mediated chemotaxis by the human chemokine receptor-like protein CCX-CKR20131686137513872-s2.0-8487444094710.1111/bph.1204223121557ZouJ.RedmondA. K.QiZ.DooleyH.SecombesC. J.The CXC chemokine receptors of fish: Insights into CXCR evolution in the vertebrates20152151171312-s2.0-8496429683810.1016/j.ygcen.2015.01.00425623148LouisA.NguyenN. T. T.MuffatoM.CrolliusH. R.Genomicus update 2015: KaryoView and MatrixView provide a genome-wide perspective to multispecies comparative genomics2015431D682D6892-s2.0-8494109244010.1093/nar/gku1112LouisA.MuffatoM.CrolliusH. R.Genomicus: Five genome browsers for comparative genomics in eukaryota2013411D700D7052-s2.0-8487655407610.1093/nar/gks1156TamuraK.StecherG.PetersonD.FilipskiA.KumarS.MEGA6: Molecular Evolutionary Genetics Analysis version 6.0201330122725272910.1093/molbev/mst197WaterhouseA. M.ProcterJ. B.MartinD. M. A.ClampM.BartonG. J.Jalview Version 2-A multiple sequence alignment editor and analysis workbench2009259118911912-s2.0-6544918823210.1093/bioinformatics/btp033TalaveraG.CastresanaJ.Improvement of phylogenies after removing divergent and ambiguously aligned blocks from protein sequence alignments20075645645772-s2.0-3454748908410.1080/10635150701472164DarribaD.TaboadaG. L.DoalloR.PosadaD.ProtTest 3: fast selection of best-fit models of protein evolution20112781164116510.1093/bioinformatics/btr0882-s2.0-79954529507GuindonS.DufayardJ.LefortV.AnisimovaM.HordijkW.GascuelO.New algorithms and methods to estimate maximum-likelihood phylogenies: assessing the performance of PhyML 3.0201059330732110.1093/sysbio/syq0102-s2.0-77950806408PageR. D.Visualizing phylogenetic trees using treeview20030016.2.16.2.1510.1002/0471250953.bi0602s01SuyamaM.TorrentsD.BorkP.PAL2NAL: robust conversion of protein sequence alignments into the corresponding codon alignments200634W609W61210.1093/nar/gkl3152-s2.0-33747847779XuB.YangZ.PAMLX: a graphical user interface for PAML201330122723272410.1093/molbev/mst1792-s2.0-84890381233AnisimovaM.BielawskiJ. P.YangZ.Accuracy and power of the likelihood ratio test in detecting adaptive molecular evolution2001188158515922-s2.0-003489770910.1093/oxfordjournals.molbev.a003945WongW. S. W.YangZ.GoldmanN.NielsenR.Accuracy and power of statistical methods for detecting adaptive evolution in protein coding sequences and for identifying positively selected sites200416821041105110.1534/genetics.104.0311532-s2.0-8544225780YangZ.NielsenR.GoldmanN.PedersenA. K.Codon-substitution models for heterogeneous selection pressure at amino acid sites200015514314492-s2.0-0034097381YangZ.Likelihood ratio tests for detecting positive selection and application to primate lysozyme evolution19981555685732-s2.0-003189796410.1093/oxfordjournals.molbev.a0259579580986PorterM. L.CroninT. W.McClellanD. A.CrandallK. A.Molecular characterization of crustacean visual pigments and the evolution of pancrustacean opsins20072412532682-s2.0-3384590661410.1093/molbev/msl152MahoA.BensimonA.VassartG.ParmentierM.Mapping of the CCXCR1, CX3CR1, CCBP2 and CCR9 genes to the CCR cluster within the 3p21.3 region of the human genome1999873-42652682-s2.0-00333682581070268910.1159/000015443NibbsR.GrahamG.RotA.Chemokines on the move: Control by the chemokine "interceptors" Duffy blood group antigen and D620031552872942-s2.0-074232199910.1016/j.smim.2003.08.00615001178NomiyamaH.YoshieO.Functional roles of evolutionary conserved motifs and residues in vertebrate chemokine receptors201597139472-s2.0-8492528706110.1189/jlb.2RU0614-290RNovitzky-BassoI.RotA.Duffy antigen receptor for chemokines and its involvement in patterning and control of inflammatory chemokines201232662-s2.0-84874217603BonecchiR.GrahamG. J.Atypical chemokine receptors and their roles in the resolution of the inflammatory response2016722410.3389/fimmu.2016.00224