Next Article in Journal
Presence of Recombinant Bat Coronavirus GCCDC1 in Cambodian Bats
Next Article in Special Issue
Different but Not Unique: Deciphering the Immunity of the Jamaican Fruit Bat by Studying Its Viriome
Previous Article in Journal
Targeting the Virus Capsid as a Tool to Fight RNA Viruses
Previous Article in Special Issue
Exploring the Role of Innate Lymphocytes in the Immune System of Bats and Virus-Host Interactions
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Brief Report

Evolution of CCR5 and CCR2 Genes in Bats Showed Multiple Independent Gene Conversion Events

by
Alexandre P. Fernandes
1,2,3,
Ana Águeda-Pinto
1,2,3,
Ana Pinheiro
1,2,3,
Hugo Rebelo
1,3,4 and
Pedro J. Esteves
1,2,3,5,*
1
CIBIO, Centro de Investigação em Biodiversidade e Recursos Genéticos, InBIO Laboratório Associado, Campus de Vairão, Universidade do Porto, 4485-661 Vairao, Portugal
2
Departamento de Biologia, Faculdade de Ciências, Universidade do Porto, 4099-002 Porto, Portugal
3
BIOPOLIS Program in Genomics, Biodiversity and Land Planning, CIBIO, Campus de Vairão, 4485-661 Vairao, Portugal
4
CIBIO/InBIO, Universidade de Lisboa, Tapada da Ajuda, 1349-017 Lisboa, Portugal
5
CITS-Centro de Investigação em Tecnologias da Saúde, Instituto Politécnico de Saúde do Norte (IPSN), Cooperativa de Ensino Superior Politécnico e Universitário (CESPU), 4585-116 Gandra, Portugal
*
Author to whom correspondence should be addressed.
Viruses 2022, 14(2), 169; https://doi.org/10.3390/v14020169
Submission received: 22 December 2021 / Revised: 12 January 2022 / Accepted: 15 January 2022 / Published: 18 January 2022
(This article belongs to the Special Issue Host Diversity and Responses to Bat-Associated Viruses)

Abstract

:
Chemokine receptors are an important determinant for the infectiousness of different pathogens, which are able to target the host cells by binding to the extracellular domains of these proteins. This is the mechanism of infection of HIV-1, among other concerning human diseases. Over the past years, it has been shown that two chemokine receptors, CCR2 and CCR5, have been shaped by events of gene conversion in different mammalian lineages, which has been linked to a possible selective advantage against pathogens. Here, by taking advantage of available bat genomes, we present the first insight of CCR2 and CCR5 evolution within the Chiroptera order. In total, four independent events of recombination between CCR2 and CCR5 were detected: two in a single species, Miniopterus natalensis; one in two species from the Rhinolophoidea superfamily; and one in four species from the Pteropodidae family. The regions affected by the gene conversions were generally extensive and always encompassed extracellular domains. Overall, we demonstrate that CCR2 and CCR5 have been subject to extensive gene conversion in multiple species of bats. Considering that bats are known to be large reservoirs of virus in nature, these results might indicate that chimeric CCR2-CCR5 genes might grant some bat species a selective advantage against viruses that rely in the extracellular portions of either CCR2 or CCR5 as gateways into the cell.

1. Introduction

Chemokine receptors are a family of seven-transmembrane-spanning G proteins important to immunological pathway signaling due to their binding ability with chemokines. They are composed of an extracellular N-terminus, an intracellular C-terminus and seven transmembrane helices (TM1-7), which are intercalated by three intracellular (ICL1-3) and three extracellular (ECL1-3) loops [1]. Many of these proteins are known to form homodimers or heterodimers, enabling new signaling pathways or enhancing already existing ones [2,3]. Chemokine receptors are expressed in lymphocytes, neutrophils, dendritic cells, and many other cell types, and play a part in both the innate and adaptive immune responses, being mainly involved in the trafficking of leukocyte populations to a site of injury and/or infection [1]. Additionally, some of these receptors have been deemed as determinants of infectiousness for different pathogens [4,5,6,7,8,9,10,11]. Most notably, CCR5 and CXCR4 have both been identified as vital co-receptors that mediate human immunodeficiency virus (HIV) entry into human cells [10,11,12,13].
Over the past years, it was shown that gene conversion events have shaped two chemokine receptors, CCR2 and CCR5, in several clades of mammals [14,15,16,17,18,19]. These genes are particularly prone to events of recombination due to their high degree of similarity (~76% identity for the human CCR2 and CCR5) and because of their proximity in the genome [1]. Two hypotheses have been put forward as possible explanations for the fixation of chimeric CCR2-CCR5 genes in different orders of mammals: (i) they facilitate the formation of CCR2-CCR5 heterodimers, which have been shown to naturally occur in human cells [14]; or (ii) they provide a selective advantage against viruses which may have used either CCR5 or CCR2 as gateways to enter the cell (in a similar fashion to HIV-1 in humans) [15,16,17]. So far, such CCR2-CCR5 recombinants have been reported in rodent, feline, leporid, guinea pig, horse, pig, marsupial, and primate lineages [14,15,16,17,18,19]. However, no evolutionary studies regarding these genes have been performed in bat species.
Bats species (Chiroptera order) have been frequently suggested as the source of outbreaks of multiple human zoonotic diseases, making them useful animal models for the fields of immunology and virology. Some examples of such diseases include rabies, Marburg virus, Ebola virus, Middle East respiratory syndrome coronavirus, and, more recently, severe acute respiratory syndrome coronavirus-2, all of which are believed to have recent ancestry in bat species and to have later been transmitted to humans by virus spillover [20,21,22,23]. Because bat colonies represent one of the largest natural reservoirs of virus in the animal kingdom [24,25], it is of the utmost importance to understand how their immune system interacts and co-evolves with viruses. In this study, we investigate possible events of gene conversion happening in CCR2 and CCR5 genes from different bat species, which may pave the way to further clarify their immune system function on these mammals.

2. Materials and Methods

2.1. Data Retrieval

The coding DNA sequences (CDSs) of CCR2 and CCR5 genes from 18 bat species were manually extracted from the NCBI database (https://www.ncbi.nlm.nih.gov/, accessed on 22 December 2021). In one species, Myotis lucifugus, the CCR2 gene was found as a pseudogene and, therefore, it was not included in our dataset. For the remaining species, a single copy of both CCR2 and CCR5 was found. Additionally, Homo sapiens, Macaca mulatta, Canis lupus familiaris and Choloepus didactylus CCR2 and CCR5 genes were included in the dataset as outgroups, resulting in a final dataset of 43 sequences. The sequences were then aligned using the ClustalW2 algorithm [26] with a gap opening penalty of 10 and a gap extending penalty of 0.2, followed by manual corrections where necessary.

2.2. Detection of Recombination and Phylogenetic Analysis

The obtained nucleotide alignment was used to create a Maximum Likelihood (ML) phylogenetic tree with the MEGA software [27,28] and 1000 bootstrap repetitions. The Tamura 3-parameter model with five gamma categories was selected as the best model of substitution for the tree by the BIC (Bayesian Information Criterion) implemented in MEGA.
For the detection of recombination events, the sequence alignment dataset was analyzed under six different algorithms (3Seq, SiScan, Chimaera, MaxChi, GENECONV and RPD) implemented in the Recombination Detection Program (RDP) [29]. To keep a conservative approach, only events of recombination detected by at least five of the six algorithms implemented were considered significant (Table S1).

3. Results and Discussion

In this study, our main goal was to determine whether the CCR2 and CCR5 genes of bats were affected by gene conversion during their evolutionary history. For that, CCR2 and CCR5 CDSs from 18 different bat species were retrieved from public databases (Figure 1). Bats included in this study represent a wide variety of species, coming from seven different major Chiroptera families: Rhinolophidae (one species), Hipposideridae (one species), Pteropodidae (four species), Phyllostomidae (four species), Molossidae (one species), Miniopteridae (one species) and Vespertilionidae (six species). As can be observed from Figure 1, the obtained ML tree could not solve the phylogeny of the sequences at its root, causing it to be split in three different clusters supported by high bootstraps: one containing the CCR5 and CCR2 of Miniopterus natalensis (91 bootstrap value), one containing the remaining CCR5 sequences, and one containing the remaining CCR2 CDSs sequences. The fact that M. natalensis CCR2 and CCR5 genes were clustered together in a separate group indicate that these sequences are more similar to each other than to their respective orthologs, which prompted us to investigate whether gene conversion could explain the phylogeny of the obtained tree.
To search for signals of gene conversion in our dataset, we used the RDP software [29] to analyze the sequences under different algorithms (see Section 2). This method identified four different events of gene conversion that have independently occurred in ancestors of bats (Figure 2 and Table S1). Individual phylogenetic trees were then created for each event of gene conversion using only the regions of the alignment identified to be involved in recombination (Figure 3). These trees were used along with statistics provided by RDP to infer the direction of gene conversion (whether CCR5 was converted to CCR2 or the opposite) for each of the recombination events.
One of the events of gene conversion identified by RDP has affected the CCR2 of all the studied species from the Pteropodidae family (Pteropus alecto, P. vampyrus, P. giganteus, and Rousettus aegyptiacus). The region affected by this recombination included residues 284 to 344 of the CCR2 protein in these species, spanning over part of TM6, ECL3, TM7, and part of the C-terminus (Figure 2A). The CCR2 and CCR5 sequences of each Pteropodidae species are not grouped in pairs, but rather follow loosely the expected phylogenetic relations of this family [30], as shown by the phylogenetic tree generated for this alignment region (Figure 3A). This suggests that a single event of recombination has shaped the CCR2 gene in Pteropodidae before the radiation of this family (~25 mya) [30]. Another event of gene conversion was detected in two species from the Rhinolophoidea superfamily (Rhinolophus ferrumequinum and Hipposideros armiger). This time, the recombinant gene was identified to be CCR5, instead of CCR2. The region involved in this recombination was also considerably more extensive, ranging from residues 66 to 191 and spanning over part of TM1, ICL1, TM2, ECL1, TM3, ICL2, and part of TM4 (Figure 2B). In the ML tree generated for this event of gene conversion (Figure 3B), the CCR5 and CCR2 sequences of both R. ferrumequinum and H. armiger can be seen grouped together within their respective species, contrary to what was observed for the recombination in Pteropodidae. This suggests that concerted evolution is currently acting on CCR2 and CCR5 in Rhinolophoidea, homogenizing the region affected by recombination in CCR5 (Figure 4). Lastly, two events of gene conversion were detected in the CCR5 of M. natalensis, explaining the observed proximity between the CCR2 and CCR5 sequences of this species in our initial ML tree (Figure 1). The first event includes residues 157 to 214, which span over ICL2, TM4 and half of the ECL2 regions (Figure 2C). The second event was more extensive, affecting positions 248 to 384 and including part of TM5, ICL3, TM6, ECL3, and TM7 regions, as well as the C-terminus (Figure 2C). M. natalensis’ CCR2 and CCR5 are identical in the regions affected by the events of recombination (Figure 4), suggesting that concerted evolution is also ongoing there. Overall, our results support that gene conversion was significant in the evolution of CCR2 and CCR5 genes in Chiroptera, having affected large portions of these genes in different clades of bats. The regions of CCR2 and CCR5 affected by gene conversion in this group varied significantly between Chiroptera lineages, including both intracellular and extracellular domains of the two proteins.
In a previous study of recombination between the CCR2 and CCR5 in cats (Felis catus) and other lineages of mammals, Vàzquez-Salat et al. [14] have suggested that facilitation of CCR2-CCR5 heterodimerization could provide a selective advantage for the fixation of CCR2-CCR5 chimeric forms. This work only identified recombination to have occurred in structural domains (TMs and ICLs) of both CCR2 and CCR5, which supported their hypothesis since these regions are known to be the most relevant for CCR2-CCR5 heterodimerization [31,32]. However, posterior studies did not find such correlation, and have rather reported that gene conversion events spanned over both structural and functional domains, which is also in line with our results for bats (Figure 2). Furthermore, a more detailed analysis of the F. catus CCR2 and CCR5 has shown that gene conversion in this species is also not restricted to TMs and ICLs [17], suggesting that the specificity of recombination between these genes might have been initially overestimated.
Alternatively, host–pathogen interactions have also been proposed as a factor driving the selection of CCR2-CCR5 recombinants [15,16,17]. This hypothesis relies heavily on the existence of pathogens that are able to hijack chemokine receptors and use them to enter the cell. Some examples of such pathogens have indeed been identified, such as HIV-1, dengue virus (DENV), Epstein–Barr virus (EBV), and Plasmodium vivax [4,5,6,7,8,9,10,11]. Although DENV and EBV interactions with chemokine receptors remain somewhat elusive, it is well established that both HIV-1 and Plasmodium vivax interact directly with chemokine receptors in order to infect the cell [4,5,7,10,11]. Site-directed mutagenesis has shown that HIV infection mainly depends on residues located in extracellular domains of CCR5, such as ECL1-3 and the N-terminus [11,13]. Similarly, Plasmodium vivax infection has been shown to rely on the extracellular portion of another chemokine receptor, DARC [7]. The importance of the extracellular domains of chemokine receptors in ligand binding is further supported by evidence of positive selection restricted to these domains [33]. Yet, it must be noted that all events of gene conversion detected in our dataset comprised portions of extracellular domains of either CCR2 or CCR5, although the N-terminus was never part of a recombination. Therefore, it is possible that affected extracellular domains of these chemokine receptors have been selected in bats in order to restrict infection of different pathogens.
Because bats are known to be reservoirs of viruses in nature, it is tempting to speculate that the gene conversion events between CCR2 and CCR5 observed in this group is due to the existence of viruses that have used these proteins as gateways for infection in the past, driving a selective force that has resulted in the fixation of CCR2-CCR5 chimeric forms. However, no virus with these characteristics has yet been identified in bats. Additionally, it has been shown that gene conversion among clustered, closely related genes is not uncommon in nature, suggesting that the fixation of CCR recombinants might be due to random chance [19]. Therefore, whether these recombinant events represent an actual adaptive advantage or not is still dependent on more research. Particularly, by discovering more viruses that use chemokine receptors as a pathway to infection would be of great interest in studying the effects of gene conversion among these proteins and in the adaptive immune system as a whole.

Supplementary Materials

The following are available online at https://www.mdpi.com/article/10.3390/v14020169/s1, Table S1: RDP [29] statistical support (average p-value) for each detected recombination.

Author Contributions

Conceptualization, P.J.E.; methodology, A.P.F. and A.Á.-P.; formal analysis, A.P.F., A.P. and A.Á.-P.; investigation, A.Á.-P., A.P.F. and P.J.E.; resources, H.R. and P.J.E.; data curation, A.P.F. and A.Á.-P.; writing—original draft preparation, A.P.F., A.Á.-P. and A.P.; writing—review and editing, A.P.F., A.P., A.Á.-P., H.R. and P.J.E.; supervision, H.R. and P.J.E.; funding acquisition, P.J.E. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by FCT—Foundation for Science and Technology, which supported the doctoral fellowship of AAP (ref. SFRH/BD/128752/2017) and the investigator grant of PJE (CEECIND/CP1601/CT0005). This work is supported by National Funds through FCT-Fundação para a Ciência e a Tecnologia in the scope of the project UIDB/50027/2020 (Base).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study.

References

  1. Olson, T.S.; Ley, K. Chemokines and chemokine receptors in leukocyte trafficking. Am. J. Physiol. Regul. Integr. Comp. Physiol. 2002, 283, R7–R28. [Google Scholar] [CrossRef] [Green Version]
  2. Stephens, B.; Handel, T.M. Chemokine receptor oligomerization and allostery. Prog. Mol. Biol. Transl. Sci. 2013, 115, 375–420. [Google Scholar]
  3. Mellado, M. Chemokine receptor homo- or heterodimerization activates distinct signaling pathways. EMBO J. 2001, 20, 2497–2507. [Google Scholar] [CrossRef] [Green Version]
  4. Hoffman, T.L.; Doms, R.W. Chemokines and coreceptors in HIV/SIV-host interactions. AIDS 1998, 12, S17–S26. [Google Scholar]
  5. Kaur, H.; Sehgal, R.; Rani, S. Duffy antigen receptor for chemokines (DARC) and susceptibility to plasmodium vivax malaria. Parasitol. Int. 2019, 71, 73–75. [Google Scholar] [CrossRef]
  6. Marques, R.E.; Guabiraba, R.; Del Sarto, J.L.; Rocha, R.F.; Queiroz, A.L.; Cisalpino, D.; Marques, P.E.; Pacca, C.C.; Fagundes, C.T.; Menezes, G.B.; et al. Dengue virus requires the CC-chemokine receptor CCR5 for replication and infection development. Immunology 2015, 145, 583–596. [Google Scholar] [CrossRef]
  7. Ntumngia, F.B.; Thomson-Luque, R.; Pires, C.V.; Adams, J.H. The role of the human duffy antigen receptor for chemokines in Malaria susceptibility: Current opinions and future treatment prospects. J. Recept. Ligand Channel Res. 2016, 9, 1–11. [Google Scholar]
  8. Vischer, H.F.; Nijmeijer, S.; Smit, M.J.; Leurs, R. Viral Hijacking of human receptors through heterodimerization. Biochem. Biophys. Res. Commun. 2008, 377, 93–97. [Google Scholar] [CrossRef]
  9. Kozireva, S.; Rudevica, Z.; Baryshev, M.; Leonciks, A.; Kashuba, E.; Kholodnyuk, I. Upregulation of the chemokine receptor CCR2B in Epstein-Barr virus-positive Burkitt lymphoma cell lines with the latency III program. Viruses 2018, 10, 239. [Google Scholar] [CrossRef] [Green Version]
  10. Huang, C.-C.; Lam, S.N.; Acharya, P.; Tang, M.; Xiang, S.-H.; Hussan, S.S.-U.; Stanfield, R.L.; Robinson, J.; Sodroski, J.; Wilson, I.A.; et al. Structures of the CCR5 N terminus and of a tyrosine-sulfated antibody with HIV-1 gp120 and CD4. Science 2007, 317, 1930–1934. [Google Scholar] [CrossRef] [Green Version]
  11. Abayev, M.; Moseri, A.; Tchaicheeyan, O.; Kessler, N.; Arshava, B.; Naider, F.; Scherf, T.; Anglister, J. An extended CCR5 ECL2 peptide forms a helix that binds HIV-1 gp120 through non-specific hydrophobic interactions. FEBS J. 2015, 282, 1906–1921. [Google Scholar] [CrossRef]
  12. Dragic, T. An overview of the determinants of CCR5 and CXCR4 Co-receptor function. J. Gen. Virol. 2001, 82, 1807–1814. [Google Scholar] [CrossRef] [Green Version]
  13. Garcia-Perez, J.; Rueda, P.; Alcami, J.; Rognan, D.; Arenzana-Seisdedos, F.; Lagane, B.; Kellenberger, E. Allosteric model of maraviroc binding to CC chemokine receptor 5 (CCR5). J. Biol. Chem. 2011, 286, 33409–33421. [Google Scholar] [CrossRef] [Green Version]
  14. Vàzquez-Salat, N.; Yuhki, N.; Beck, T.; O’Brien, S.J.; Murphy, W.J. Gene conversion between mammalian CCR2 and CCR5 chemokine receptor genes: A potential mechanism for receptor dimerization. Genomics 2007, 90, 213–224. [Google Scholar] [CrossRef] [Green Version]
  15. Carmo, C.R.; Esteves, P.J.; Ferrand, N.; van der Loo, W. 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. [Google Scholar] [CrossRef]
  16. Abrantes, J.; Carmo, C.R.; Matthee, C.A.; Yamada, F.; van der Loo, W.; Esteves, P.J. A shared unusual genetic change at the chemokine receptor type 5 between oryctolagus, bunolagus and pentalagus. Conserv. Genet. 2011, 12, 325–330. [Google Scholar] [CrossRef]
  17. Esteves, P.J.; Abrantes, J.; van der Loo, W. Extensive gene conversion between CCR2 and CCR5 in domestic cat (Felis Catus). Int. J. Immunogenet. 2007, 34, 321–324. [Google Scholar] [CrossRef]
  18. Perelygin, A.A.; Zharkikh, A.A.; Astakhova, N.M.; Lear, T.L.; Brinton, M.A. Concerted evolution of vertebrate CCR2 and CCR5 genes and the origin of a recombinant equine CCR5/2 gene. J. Hered. 2008, 99, 500–511. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  19. Shields, D.C. Gene conversion among chemokine receptors. Gene 2000, 246, 239–245. [Google Scholar] [CrossRef]
  20. Chomel, B.B.; Stuckey, M.J.; Boulouis, H.-J.; Setién, A.A. Bat-related zoonoses. In Zoonoses-Infections Affecting Humans and Animals; Springer: Dordrecht, The Netherlands, 2015; pp. 697–714. [Google Scholar]
  21. Wang, Q.; Qi, J.; Yuan, Y.; Xuan, Y.; Han, P.; Wan, Y.; Ji, W.; Li, Y.; Wu, Y.; Wang, J.; et al. Bat origins of MERS-CoV supported by bat Coronavirus HKU4 usage of human receptor CD26. Cell Host Microbe 2014, 16, 328–337. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  22. Andersen, K.G.; Rambaut, A.; Lipkin, W.I.; Holmes, E.C.; Garry, R.F. The proximal origin of SARS-CoV-2. Nat. Med. 2020, 26, 450–452. [Google Scholar] [CrossRef] [Green Version]
  23. Seyran, M.; Pizzol, D.; Adadi, P.; El-Aziz, T.M.A.; Hassan, S.S.; Soares, A.; Kandimalla, R.; Lundstrom, K.; Tambuwala, M.; Aljabali, A.A.A.; et al. Questions concerning the proximal origin of SARS-CoV-2. J. Med. Virol. 2021, 93, 1204–1206. [Google Scholar] [CrossRef]
  24. Calisher, C.H.; Childs, J.E.; Field, H.E.; Holmes, K.V.; Schountz, T. Bats: Important reservoir hosts of emerging viruses. Clin. Microbiol. Rev. 2006, 19, 531–545. [Google Scholar] [CrossRef] [Green Version]
  25. Wong, S.; Lau, S.; Woo, P.; Yuen, K.-Y. Bats as a continuing source of emerging infections in humans. Rev. Med. Virol. 2007, 17, 67–91. [Google Scholar] [CrossRef] [PubMed]
  26. Larkin, M.A.; Blackshields, G.; Brown, N.P.; Chenna, R.; McGettigan, P.A.; McWilliam, H.; Valentin, F.; Wallace, I.M.; Wilm, A.; Lopez, R.; et al. Clustal W and Clustal X version 2.0. Bioinformatics 2007, 23, 2947–2948. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  27. Kumar, S.; Stecher, G.; Li, M.; Knyaz, C.; Tamura, K. MEGA X: Molecular evolutionary genetics analysis across computing platforms. Mol. Biol. Evol. 2018, 35, 1547–1549. [Google Scholar] [CrossRef]
  28. Cummings, M.P. BAMBE, DnaSP, ENCprime/SeqCount, LAMARC, MacClade, MEGA, Modeltest, MrBayes, PAML, PAUP*, PHYLIP, r8s, readseq, Seq-Gen, Sites, TreeView. In Dictionary of Bioinformatics and Computational Biology; Hancock, J.M., Zvelebil, M.J., Eds.; Wiley-Liss: Hoboken, NJ, USA, 2004; p. 318. ISBN 978-0-471-43622-5. [Google Scholar]
  29. Martin, D.P.; Murrell, B.; Golden, M.; Khoosal, A.; Muhire, B. RDP4: Detection and analysis of recombination patterns in virus genomes. Virus Evol. 2015, 1. [Google Scholar] [CrossRef] [Green Version]
  30. Almeida, F.C.; Simmons, N.B.; Giannini, N.P. A Species-level phylogeny of old world fruit bats with a new higher-level classification of the family pteropodidae. Am. Mus. Novit. 2020, 2020, 1. [Google Scholar] [CrossRef]
  31. Gahbauer, S.; Pluhackova, K.; Böckmann, R.A. Closely related, yet unique: Distinct homo-and heterodimerization patterns of G protein coupled chemokine receptors and their fine-tuning by cholesterol. PLoS Comput. Biol. 2018, 14, e1006062. [Google Scholar] [CrossRef]
  32. Hernanz-Falcón, P.; Rodríguez-Frade, J.M.; Serrano, A.; Juan, D.; del Sol, A.; Soriano, S.F.; Roncal, F.; Gómez, L.; Valencia, A.; Martínez-A, C.; et al. Identification of amino acid residues crucial for chemokine receptor dimerization. Nat. Immunol. 2004, 5, 216–223. [Google Scholar] [CrossRef]
  33. Metzger, K.J.; Thomas, M.A. Evidence of positive selection at codon sites localized in extracellular domains of mammalian CC motif chemokine receptor proteins. BMC Evol. Biol. 2010, 10, 139. [Google Scholar] [CrossRef] [PubMed] [Green Version]
Figure 1. Phylogenetic analysis of CCR5 and CCR2 nucleotide sequences from different bat species. The Maximum Likelihood tree of bat CCR2 and CCR5 nucleotide sequences was generated with MEGA [27,28] using 1000 bootstrap repetitions. Only bootstrap values >75% are shown.
Figure 1. Phylogenetic analysis of CCR5 and CCR2 nucleotide sequences from different bat species. The Maximum Likelihood tree of bat CCR2 and CCR5 nucleotide sequences was generated with MEGA [27,28] using 1000 bootstrap repetitions. Only bootstrap values >75% are shown.
Viruses 14 00169 g001
Figure 2. Diagrams of CCR2 and CCR5 genes detected to have undergone gene conversion in bats. RDP [29] analyses performed with the CDSs of CCR2 and CCR5 genes from Chiroptera species have identified four independent events of gene conversion to have occurred in this order. Sections of the proteins which have been acquired through gene conversion are shown in the diagrams as dislocated blocks with their starting and ending positions marked. Protein domains are colored according to the legend. (A) CCR2 protein of Pteropodidae family members. (B) CCR5 protein of the Rhinolophoidea superfamily members. (C) CCR5 protein of M. natalensis.
Figure 2. Diagrams of CCR2 and CCR5 genes detected to have undergone gene conversion in bats. RDP [29] analyses performed with the CDSs of CCR2 and CCR5 genes from Chiroptera species have identified four independent events of gene conversion to have occurred in this order. Sections of the proteins which have been acquired through gene conversion are shown in the diagrams as dislocated blocks with their starting and ending positions marked. Protein domains are colored according to the legend. (A) CCR2 protein of Pteropodidae family members. (B) CCR5 protein of the Rhinolophoidea superfamily members. (C) CCR5 protein of M. natalensis.
Viruses 14 00169 g002
Figure 3. ML trees generated using only regions affected by each event of gene conversion detected in bat species. ML trees were generated with 1000 bootstrap repetitions using MEGA [27,28]. Only bootstraps values above 75% are shown below each branch. Recombinant sequences are marked in red in the ML trees. (A) ML tree of the subsection of the alignment where recombination has shaped the CCR2 of the Pteropodidae family (positions 284 to 344). (B) ML tree of the subsection of the alignment where recombination has shaped the CCR5 of the Rhinolophoidea superfamily (positions 66 to 191). (C) ML tree of the first subsection of the alignment where recombination has shaped the CCR5 of M. natalensis (positions 157 to 214). (D) ML tree of the second subsection of the alignment where recombination has shaped the CCR5 of M. natalensis (positions 248 to 384).
Figure 3. ML trees generated using only regions affected by each event of gene conversion detected in bat species. ML trees were generated with 1000 bootstrap repetitions using MEGA [27,28]. Only bootstraps values above 75% are shown below each branch. Recombinant sequences are marked in red in the ML trees. (A) ML tree of the subsection of the alignment where recombination has shaped the CCR2 of the Pteropodidae family (positions 284 to 344). (B) ML tree of the subsection of the alignment where recombination has shaped the CCR5 of the Rhinolophoidea superfamily (positions 66 to 191). (C) ML tree of the first subsection of the alignment where recombination has shaped the CCR5 of M. natalensis (positions 157 to 214). (D) ML tree of the second subsection of the alignment where recombination has shaped the CCR5 of M. natalensis (positions 248 to 384).
Viruses 14 00169 g003
Figure 4. Amino acid alignment of CCR2 and CCR5 from bat species showing different regions where recombination was detected. The CCR5 sequence of Phyllostomus discolor was chosen as a reference and identity to this sequence is marked by dots (.). Colored boxes identify the different recombination events detected by RDP [29].
Figure 4. Amino acid alignment of CCR2 and CCR5 from bat species showing different regions where recombination was detected. The CCR5 sequence of Phyllostomus discolor was chosen as a reference and identity to this sequence is marked by dots (.). Colored boxes identify the different recombination events detected by RDP [29].
Viruses 14 00169 g004
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Fernandes, A.P.; Águeda-Pinto, A.; Pinheiro, A.; Rebelo, H.; Esteves, P.J. Evolution of CCR5 and CCR2 Genes in Bats Showed Multiple Independent Gene Conversion Events. Viruses 2022, 14, 169. https://doi.org/10.3390/v14020169

AMA Style

Fernandes AP, Águeda-Pinto A, Pinheiro A, Rebelo H, Esteves PJ. Evolution of CCR5 and CCR2 Genes in Bats Showed Multiple Independent Gene Conversion Events. Viruses. 2022; 14(2):169. https://doi.org/10.3390/v14020169

Chicago/Turabian Style

Fernandes, Alexandre P., Ana Águeda-Pinto, Ana Pinheiro, Hugo Rebelo, and Pedro J. Esteves. 2022. "Evolution of CCR5 and CCR2 Genes in Bats Showed Multiple Independent Gene Conversion Events" Viruses 14, no. 2: 169. https://doi.org/10.3390/v14020169

APA Style

Fernandes, A. P., Águeda-Pinto, A., Pinheiro, A., Rebelo, H., & Esteves, P. J. (2022). Evolution of CCR5 and CCR2 Genes in Bats Showed Multiple Independent Gene Conversion Events. Viruses, 14(2), 169. https://doi.org/10.3390/v14020169

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop