Skip to main content

Molecular evolution and expression profile of the chemerine encoding gene RARRES2 in baboon and chimpanzee

Abstract

Background

Chemerin, encoded by the retinoic acid receptor responder 2 (RARRES2) gene is an adipocytesecreted protein with autocrine/paracrine functions in adipose tissue, metabolism and inflammation with a recently described function in vascular tone regulation, liver, steatosis, etc. This molecule is believed to represent a critical endocrine signal linking obesity to diabetes. There are no data available regarding evolution of RARRES2 in non-human primates and great apes. Expression profile and orthology in RARRES2 genes are unknown aspects in the biology of this multigene family in primates. Thus; we attempt to describe expression profile and phylogenetic relationship as complementary knowledge in the function of this gene in primates. To do that, we performed A RT-PCR from different tissues obtained during necropsies. Also we tested the hypotheses of positive evolution, purifying selection, and neutrality. And finally a phylogenetic analysis was made between primates RARRES2 protein.

Results

RARRES2 transcripts were present in liver, lung, adipose tissue, ovary, pancreas, heart, hypothalamus and pituitary tissues. Expression in kidney and leukocytes were not detectable in either species. It was determined that the studied genes are orthologous.

Conclusions

RARRES2 evolution fits the hypothesis of purifying selection. Expression profiles of the RARRES2 gene are similar in baboons and chimpanzees and are also phylogenetically related.

Background

Chemerine is a protein that initiates chemotaxis via the ChemR23 - G protein-coupled seven-trans-membrane domain receptor ligand, which has been classified as an adipokine due to its role in adipocyte differentiation and glucose uptake [1]. It also plays a potential role in controlling immune responses at sites of tissue injury and inflammation [2], including chronic inflammation of adipose tissue in obesity [1, 35]. Chemerine has also been suggested as an essential endocrine signal, linking obesity to insulin resistance [3, 68], therefore it is an independent biomarker of metabolic syndrome [912]. In addition to adipose tissue, chemerine plays an important role in metabolic regulation in the liver and skeletal muscle [6, 13]. Recently a novel role for chemerine as a stimulator of angiogenesis was identified [9].

The retinoic acid receptor responder protein 2 (RARRES2) gene (also named RAR-responsive protein TIG2, chemerin and tazarotene-induced gene 2 protein), which encodes chemerine, is located in chromosome 7 at 7q36.1 in humans. RARRES2 mRNA is highly expressed in white adipose tissue, liver and lungs, while the mRNA for chemerine receptor is predominantly expressed in immune cells and adipose tissue [1420].

The study of the evolution of Old World primates (OWM) and great apes has been an excellent approach to understand human pathology such as metabolic syndrome. Currently, the baboon (Papio spp) has been proven to be an ideal model to study metabolism disturbances. Previous studies have found a substantial variation in weight and body composition in adult baboons sharing the same diet and living conditions. Baboons spontaneously develop obesity [21], type 2 diabetes mellitus (T2DM) [22] and a metabolic syndrome-like phenotype has been described in this species [23]. While the RARRES2 gene sequence has been described in humans, there is no information available regarding baboon and chimpanzee. The present study analyzed the expression profile and phylogenetic relationship of the RARRES2 gene from baboon and chimpanzee.

Results

Expression profile

In the PCR from genomic DNA, a single band of 3,257 bp was visualized for both primates (data not shown). In a similar manner in the PCR from cDNA, only one product of 641 bp was observed for both primates. PCR product sizes and nucleotide sequence are not sufficient to determine absence of paralogous genes, but with the evidence given at the moment could be that there is only one RARRES2 gene for each of the two species studied, baboon and chimpanzee. Both genes are organized in six exons and five introns of variable size: exon 1,110 pb; intron A, 948 pb; exon 2, 194 pb; intron B, 231 pb; exon 3, 105 pb; intron C, 1018 pb; exon 4, 96 pb; intron D, 244 pb; exon 5, 127 pb; intron E, 175 pb; and exon 6, 121 pb.

For each species, RARRES2 transcripts were amplified from liver, lung, adipose tissue, ovary, pancreas, heart, hypothalamus, and pituitary gland. Expression in kidney and leucocytes was not detected. Negative and positive controls gave the expected results (Fig. 1). All PCR products were cloned and sequenced. For every primate, at least three independent clones were sequenced for each gene and transcript. Baboon and chimpanzee mRNA sequences resulted identical for all tissues.

Fig. 1
figure 1

RT-PCR profile expression screening of RARRES2 gene. Top panel shows RT-PCR assays from baboon’s tissues. Button panel shows RT-PCR assays from chimpanzee’s tissues. In both panels, the lanes 1 to 10 represent: liver, lung, adipose tissue, ovary, pancreas, heart, hypothalamus, pituitary gland, kidney and leucocytes. Lanes (−) indicate the negative RT-PCR controls. Positives controls (324 bp) were co amplified with samples (641 bp)

Phylogenetic analysis

The phylogenetic tree (Fig. 2) shows four clades in a lineage specific manner. These clades correspond to apes, OWM, NWM and lemur (out-group). It confirms the orthology between primate RARRES2 genes. Bootstrap values are shown on the tree’s branches. Similar results were obtained using Maximum Likelihood (ML), Neighbor-Joining (NJ) and UPGMA phylogenetic methods. We confirmed that RARRES2 evolution fits the hypothesis of purifying selection (P < 0.05) for OWM (except for green monkey) and apes. Since for RARRES2 from NWM the tests did not yield statistically significant differences (P > 0.05), at this moment it is unclear which forces govern their evolution. Neutrality is observed for orangutan, rhesus monkey and crab-eating macaque RARRES2 genes (Table 1).

Fig. 2
figure 2

Phylogenetic tree of RARRES2 proteins from various primates. The tree was built using MEGA version 6.06 by the ML, NJ and UPGMA methods and further bootstrap analysis of 1000 replicas. Numbers in the branches indicate Bootstrap value (bold) and below the branches length, both estimated by three methods (ML/NJ/UPGMA). Clades are in linage specific manner, apes, OWM, NWM and lemur (out-group)

Table 1 Test to determine the evolutionary forces that underlie the RARRES2 primate genes using the Li-Wu-Luo method (Kimura 2-parameters)

Discussions

mRNA tissue expressions

This is the first report showing the expression profile of the novel RARRES2 mRNA in different tissues from baboon and chimpanzee. The novel finding of this study is that in baboon and chimpanzee, there is expression in the pituitary gland, while in humans RARRES2 was not detected in this tissue. In humans, expression in kidney has been reported (ETS profile at the NCBI, UniGene), while in baboon and chimpanzee in our study it was not detected. To address this difference studies such as; promoters alignment, bioinformatics analyses to predict transcptional factors bindings sites, searching CpG islands susceptible to be methylated, as others; are needed.

Similarity in the expression pattern suggests that mechanisms of regulation of gene expression are similar in the three species. It also suggests that the RARRES2 gene might have similar physiological effects. However more studies in functional differences are needed, such as studying the expression profile under different conditions such as insulin resistance, obesity, fasting or others associated with metabolic syndrome.

Gene structure

We found that the RARRES2 cDNA sequence in baboon and chimpanzee have a high similarity to others apes. The conservation of exon–intron boundaries and absence of obvious mutations suggest that the baboon and chimpanzee RARRES2 genes are functional.

Phylogenetic analysis

Positive selection (dN > dS) implies that the substitutions mostly non-synonymous are functional and benefit the organism, conferring some evolutionary advantage. While purifying selection (dN < dS) indicates that evolutionary pressure has been relaxed. The dN and dS rates show that the evolutionary force actually acting on these is the purification of the selection (P < 0.05). Fit hypothesis of purifying of selection is a clue that these genes are functional in the studied species, because there are not functional genes that do not fit this hypothesis. Although we did not do an analysis between the expression profiles of primate RARRES2 genes and their putative binding sites for known transcriptional factors, we speculate that this evolutionary force has been relaxed perhaps because in apes and OWM the functions and/or regulatory mechanisms are not necessarily the same as those governing NWM or lemurs. However, more research is needed to investigate the role of this molecule in physiology and regulatory expression profiles.

The phylogenetic relationship between NWM, OWM and apes RARRRES proteins was determined to evaluate their evolution in primates. The phylogenetic tree (Fig. 2) showed three clades in a linage-specific manner. These clades corresponded to NWM, OWM and apes, finally galago (out-group). Tree’s topology, length’s branches and Bootstrap values are similar using either phylogenetic method (ML/NJ/UPGMA). This confirms orthology within the RARRES gene. Bootstrap values are shown bold on the branches of the tree; below of Bootstrap values are shown branches length estimated by the same methods (ML/NJ/UPGMA).

Conclusions

We found that there is only one RARRES2 gene for baboon and chimpanzee. RARRES2 transcripts were present in liver, lung, adipose tissue, ovary, pancreas, heart, hypothalamus and pituitary tissues from both primates. The dN and dS rates conclusively show that these genes have no positive selection. RARRES2 evolution fits the hypothesis of purifying selection. The studied sequences are clearly orthologous.

Methods

Animal housing

The baboons (Papio anubis) and chimpanzees (Pan troglodytes) were housed at the Southwest National Primate Research Center in San Antonio, TX; an AAALAC-accredited facility at the Texas Biomedical Research Institute. These animals are housed with indoor-outdoor access in standard stainless steel cages with covered shelters equipped with external containers for food and ad libitum access to water. The protocol was approved by the IACUC of the Texas Biomedical Research Institute for all procedures.

Tissue harvesting and processing

The liver, lung, adipose tissue, ovary, pancreas, heart, hypothalamus, pituitary gland, kidney and total leukocytes were collected at routine necropsy and stored at −80 °C until further evaluation. Baboon tissues came from three different animals, while chimpanzee tissues came from a single female that underwent humane euthanasia due to a complication of diabetes. Total RNA and DNA of the different tissues and leukocytes were extracted with Trizol reagent according to the manufacturer’s instructions (Invitrogen, Carlsbad, CA). RNA was treated with RQ1 DNase (Promega, Madison, WI) for 15 min at 37 °C to remove traces of genomic DNA. The purity and integrity of RNA and DNA were assessed using standard spectrophotometry methods with NanoDrop equipment (Thermo Scientific, Wilmington, DE) and agarose gel electrophoresis.

Reverse transcription (RT), polymerase chain reaction (PCR), molecular cloning and nucleotide sequencing

Total RNA from frozen tissues (~1.0 μg) was retro-transcribed according to the manufacturer’s instructions with a High Capacity cDNA Reverse Transcription (RT) kit (Applied Biosystems, Foster City, CA) into a 50 μL final volume. A primer set to amplify baboon and chimpanzee’s RARRES2 transcripts was designed using RARRES2 sequences from rhesus monkey (Macaca mulatta) and humans (Homo sapiens) as templates. Sequence templates are deposited in the NCBI (Table 2) [http://www.ncbi.nlm.nih.gov/]. Such design used an online tool [24]. Sense primer (5′-GGACAGCGAGGCCAAGAT-3′) hybridizes in exon one at nucleotide 68 prior to the start codon (AUG). Antisense primer (5′-CTGGGGTCTTCCACTGGTTAC-3′) hybridizes in exon six at nucleotide after the stop triplet (UAA) (positions referent to the RNA molecule). For PCR we used 5 μL of RT reaction or 200 ng of genomic DNA as templates, 10 μM of each primer and the 2X PCR master mix kit (Qiagen, Valencia, CA). The reactions were performed in a 25 μL final volume carried out in a Veriti thermal cycler (Applied Biosystems, Foster City, CA), with the following amplification programs: first for RT samples, an initial denaturalization step for 3 min at 94 °C; then 30 cycles of 30 s at 94 °C, 30 s at 58 °C, and 30 s at 72 °C; finally an elongation step of 6 min at 72 °C. Second for genomic DNA samples, an initial denaturalization step for 4 min at 94 °C; then 30 cycles of 1 min at 94 °C, 30 s at 58 °C, 2.5 min at 72 °C; finally an elongation step of 15 min at 72 °C. Primers for ribosomal RNA fraction 18S (Applied Biosystems, Foster City, CA), were used as RT positive controls (324 bp in size). The amplification reactions were confirmed in 1 % agarose electrophoresis gel stained with ethidium bromide and visualized with UV light. PCR products were cloned using a 3.5 kb-XL-TOPO-vector kit and transformed in the electrocompetent bacterium E. coli strain Top 10. These procedures were performed according to the manufacturer’s specifications (Invitrogen, Carlsbad, CA). The positive clones were sequenced using a Big Dye terminator kit and either specific or universal (M13) primers. The reactions were analyzed with an ABI PRISM 3100 Genetic Analyzer and in its software (Applied Biosystems, Foster City, CA).

Table 2 Primate RARRES2 sequences from NCBI GenBank used in this study

Phylogenetic analysis

The information obtained from the sequencing assays was subjected to a BLAST test to determine identity. The structures of intron-exon boundaries of each one of the baboon and chimpanzee RARRES2 genes were determined by directly retrieving information given at the NCBI server. Alignments were performed using the CLUSTAL W program [25]. GenBank accession numbers of the sequences used in this study are provided in Table 2. From amino acidic sequence, a phylogenetic tree was built with MEGA 6.06 software [26] using the Maximum Likelihood (ML), Neighbor-Joining (NJ) and UPGMA methods, thena bootstrap test was done with 1,000 replicates [27].

Seeking to identify the evolutionary forces that underlie the process of divergence in the RARRES2 primate genes, we tested the hypothesis of positive or adaptive evolution (dN > dS), purifying selection (dN < dS) and neutrality (dN = dS). For this purpose; first, we calculated the non-synonymous (causes an amino acid change) and synonymous (does not cause an amino acid change) dN, dS distances, respectively, by the the Li-Wu-Luo method (Kimura 2-parameters) [28] from RARRES2’ coding sequences from apes, OWM and NWM with their lemur’s counterpart. Second, we tested whether dN is significantly greater, lower or equal, respectively, than dS using codon-based Z test of selection as implemented in MEGA 6.06 software [26]. Differences were considered statistically significant at a P <0.05.

Abbreviations

RARRES2:

Retinoic acid receptor responder protein 2

OWM:

Old World primates

T2DM:

Type 2 diabetes mellitus

References

  1. Goralski KB, McCarthy TC, Hanniman EA, Zabel BA, Butcher EC, Parlee SD, et al. Chemerin, a novel adipokine that regulates adipogenesis and adipocyte metabolism. J BiolChem. 2007;282:28175–88.

    CAS  Google Scholar 

  2. Allen SJ, Zabel BA, Kirkpatrick J, Butcher EC, Nietlispach D, Handel TM. NMR assignment of human chemerin, a novel chemoattractant. Biomol NMR Assign. 2007;1:171–3.

    Article  CAS  PubMed  Google Scholar 

  3. Ernst MC, Sinal CJ. Chemerin: at the crossroads of inflammation and obesity. Trends EndocrinolMetab. 2010;21:660–7.

    Article  CAS  Google Scholar 

  4. Hart R, Greaves DR. Chemerin contributes to inflammation by promoting macrophage adhesion to VCAM-1 and fibronectin through clustering of VLA-4 and VLA-5. J Immunol. 2010;185:3728–39.

    Article  CAS  PubMed  Google Scholar 

  5. Parolini S, Santoro A, Marcenaro E, Luini W, Massardi L, Facchetti F, et al. The role of chemerin in the colocalization of NK and dendritic cell subsets into inflamed tissues. Blood. 2007;109:3625–32.

    Article  CAS  PubMed  Google Scholar 

  6. Becker M, Rabe K, Lebherz C, Zugwurst J, Goke B, Parhofer KG, et al. Expression of human chemerin induces insulin resistance in the skeletal muscle but does not affect weight, lipid levels, and atherosclerosis in LDL receptor knockout mice on high-fat diet. Diabetes. 2010;59:2898–903.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  7. Parlee SD, Ernst MC, Muruganandan S, Sinal CJ, Goralski KB. Serum chemerin levels vary with time of day and are modified by obesity and tumor necrosis factor-{alpha}. Endocrinology. 2010;151:2590–602.

    Article  CAS  PubMed  Google Scholar 

  8. Pfau D, Stepan H, Kratzsch J, Verlohren M, Verlohren HJ, Drynda K, et al. Circulating levels of the adipokinechemerin in gestational diabetes mellitus. Horm Res Paediatr. 2010;74:56–61.

    Article  CAS  PubMed  Google Scholar 

  9. Bozaoglu K, Segal D, Shields KA, Cummings N, Curran JE, Comuzzie AG, et al. Chemerin is associated with metabolic syndrome phenotypes in a Mexican-American population. J ClinEndocrinolMetab. 2009;94:3085–8.

    Article  CAS  Google Scholar 

  10. Lehrke M, Becker A, Greif M, Stark R, Laubender RP, von Ziegler F, et al. Chemerin is associated with markers of inflammation and components of the metabolic syndrome but does not predict coronary atherosclerosis. Eur J Endocrinol. 2009;161:339–44.

    Article  CAS  PubMed  Google Scholar 

  11. Lijnen HR. Angiogenesis and obesity. Cardiovasc Res. 2008;78:286–93.

    Article  CAS  PubMed  Google Scholar 

  12. Stejskal D, Karpisek M, Hanulova Z, Svestak M. Chemerin is an independent marker of the metabolic syndrome in a Caucasian population--a pilot study. Biomed Pap Med FacUnivPalacky Olomouc Czech Repub. 2008;152:217–21.

    Article  CAS  Google Scholar 

  13. Ernst MC, Issa M, Goralski KB, Sinal CJ. Chemerin exacerbates glucose intolerance in mouse models of obesity and diabetes. Endocrinology. 2010;151:1998–2007.

    Article  CAS  PubMed  Google Scholar 

  14. Lin W, Chen YL, Jiang L, Chen JK. Reduced expression of chemerin is associated with a poor prognosis and a lowed infiltration of both dendritic cells and natural killer cells in human hepatocellular carcinoma. Clin Lab. 2011;57:879–85.

    CAS  PubMed  Google Scholar 

  15. Bozaoglu K, Curran JE, Stocker CJ, Zaibi MS, Segal D, Konstantopoulos N, et al. Chemerin, a novel adipokine in the regulation of angiogenesis. J ClinEndocrinolMetab. 2010;95:2476–85.

    Article  CAS  Google Scholar 

  16. Du XY, Leung LL. Proteolytic regulatory mechanism of chemerin bioactivity. Acta BiochimBiophys Sin (Shanghai). 2009;41:973–9.

    Article  CAS  Google Scholar 

  17. Kaur J, Adya R, Tan BK, Chen J, Randeva HS, Kaur J, et al. Identification of chemerin receptor (ChemR23) in human endothelial cells: chemerin-induced endothelial angiogenesis. BiochemBiophys Res Commun. 2010;391:1762–8.

    Article  CAS  Google Scholar 

  18. Knight A. The beginning of the end for chimpanzee experiments? Philos Ethics Humanit Med. 2008;3:16.

    Article  PubMed Central  PubMed  Google Scholar 

  19. Muruganandan S, Roman AA, Sinal CJ. Role of chemerin/CMKLR1 signaling in adipogenesis and osteoblastogenesis of bone marrow stem cells. J Bone Miner Res. 2010;25:222–34.

    Article  CAS  PubMed  Google Scholar 

  20. Yoshimura T, Oppenheim JJ. Chemerin reveals its chimeric nature. J Exp Med. 2008;205:2187–90.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  21. Comuzzie AG, Cole SA, Martin L, Carey KD, Mahaney MC, Blangero J, et al. The baboon as a nonhuman primate model for the study of the genetics of obesity. Obes Res. 2003;11:75–80.

    Article  PubMed  Google Scholar 

  22. Banks WA, Altmann J, Sapolsky RM, Phillips-Conroy JE, Morley JE. Serum leptin levels as a marker for a syndrome X-like condition in wild baboons. J ClinEndocrinolMetab. 2003;88:1234–40.

    Article  CAS  Google Scholar 

  23. Cole SA, Martin LJ, Peebles KW, Leland MM, Rice K, VandeBerg JL, et al. Genetics of leptin expression in baboons. Int J ObesRelatMetabDisord. 2003;27:778–83.

    CAS  Google Scholar 

  24. Boutros R, Stokes N, Bekaert M, Teeling EC. UniPrime2: a web service providing easier Universal Primer design. Nucleic Acids Res. 2009;37:W209–213.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  25. Thompson JD, Higgins DG, Gibson TJ. CLUSTAL W: improving the sensitivity of progressive multiple sequence alignment through sequence weighting, position-specific gap penalties and weight matrix choice. Nucleic Acids Res. 1994;22:4673–80.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  26. Tamura K, Stecher G, Peterson D, Filipski A, Kumar S. MEGA6: Molecular Evolutionary Genetics Analysis version 6.0. MolBiolEvol. 2013;30:2725–9.

    CAS  Google Scholar 

  27. Saitou N, Nei M. The neighbor-joining method: a new method for reconstructing phylogenetic trees. MolBiolEvol. 1987;4:406–25.

    CAS  Google Scholar 

  28. Zhang J, Rosenberg HF, Nei M. Positive Darwinian selection after gene duplication in primate ribonuclease genes. ProcNatlAcadSci U S A. 1998;95:3708–13.

    Article  CAS  Google Scholar 

Download references

Acknowledgements

The present work was supported by grants from the Mexican Council of Science and Technology, CONACyT (167697 and 157965). The authors gratefully acknowledge the critical reading of the manuscript by Sergio Lozano-Rodríguez, MD.

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Irám Pablo Rodríguez-Sánchez.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors’ contributions

RGA development the alignment and phylogenetic analysis; MLGR, IDG, VMTA, LEMDV, RCDC, ALT and MET carried out primer designed, nucleic acid extraction, performed RT-PCR assays, molecular cloning; RGG and OBQ sequenced the respective recombinant plasmids; NESL, MDRRP, SAC and DRP carried out the design of the study;MMA, AGC, HABS and IPRS participated in the design, coordination and writing of the manuscript. All authors read and approved the final manuscript.

Rights and permissions

Open Access  This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made.

The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder.

To view a copy of this licence, visit https://creativecommons.org/licenses/by/4.0/.

The Creative Commons Public Domain Dedication waiver (https://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

González-Alvarez, R., Garza-Rodríguez, M.d.L., Delgado-Enciso, I. et al. Molecular evolution and expression profile of the chemerine encoding gene RARRES2 in baboon and chimpanzee. Biol Res 48, 31 (2015). https://doi.org/10.1186/s40659-015-0020-0

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s40659-015-0020-0

Keywords