Europe PMC

This website requires cookies, and the limited processing of your personal data in order to function. By using the site you are agreeing to this as outlined in our privacy notice and cookie policy.

Phosphoenolpyruvate carboxykinase 1 gene (Pck1) displays parallel evolution between Old World and New World fruit bats.

Author information

Affiliations

  • 1. Institute of Molecular Ecology and Evolution, SKLEC & IECR & IAIR, East China Normal University, Shanghai, China.
      Authors
    • Zhu L 1
    • Yin Q 1
    • Zhang S 1
    (3 authors)
  • 2. Department of Laboratory Medicine and Pathobiology, University of Toronto, Toronto, Canada.
      Authors
    • Irwin DM 2
    (1 author)

ORCIDs linked to this article

Plos one, 25 Mar 2015, 10(3):e0118666
https://doi.org/10.1371/journal.pone.0118666 PMID: 25807515 PMCID: PMC4373879

This article is in the Europe PMC Open access subset. Refer to the copyright information in the article for licensing details.
Free full text in Europe PMC

Abstract 

Bats are an ideal mammalian group for exploring adaptations to fasting due to their large variety of diets and because fasting is a regular part of their life cycle. Mammals fed on a carbohydrate-rich diet experience a rapid decrease in blood glucose levels during a fast, thus, the development of mechanisms to resist the consequences of regular fasts, experienced on a daily basis, must have been crucial in the evolution of frugivorous bats. Phosphoenolpyruvate carboxykinase 1 (PEPCK1, encoded by the Pck1 gene) is the rate-limiting enzyme in gluconeogenesis and is largely responsible for the maintenance of glucose homeostasis during fasting in fruit-eating bats. To test whether Pck1 has experienced adaptive evolution in frugivorous bats, we obtained Pck1 coding sequence from 20 species of bats, including five Old World fruit bats (OWFBs) (Pteropodidae) and two New World fruit bats (NWFBs) (Phyllostomidae). Our molecular evolutionary analyses of these sequences revealed that Pck1 was under purifying selection in both Old World and New World fruit bats with no evidence of positive selection detected in either ancestral branch leading to fruit bats. Interestingly, however, six specific amino acid substitutions were detected on the ancestral lineage of OWFBs. In addition, we found considerable evidence for parallel evolution, at the amino acid level, between the PEPCK1 sequences of Old World fruit bats and New World fruit bats. Test for parallel evolution showed that four parallel substitutions (Q276R, R503H, I558V and Q593R) were driven by natural selection. Our study provides evidence that Pck1 underwent parallel evolution between Old World and New World fruit bats, two lineages of mammals that feed on a carbohydrate-rich diet and experience regular periods of fasting as part of their life cycle.

Free full text 

Logo of plosoneLink to Publisher's site
PLoS One. 2015; 10(3): e0118666.
Published online 2015 Mar 25. https://doi.org/10.1371/journal.pone.0118666
PMCID: PMC4373879
PMID: 25807515

Phosphoenolpyruvate Carboxykinase 1 Gene (Pck1) Displays Parallel Evolution between Old World and New World Fruit Bats

Lei Zhu,# 1 Qiuyuan Yin,# 1 David M. Irwin, 2 and Shuyi Zhang 1 ,*
Marc Robinson-Rechavi, Academic Editor
Author information Article notes Copyright and License information Disclaimer

Associated Data

Supplementary Materials
Data Availability Statement

Abstract

Bats are an ideal mammalian group for exploring adaptations to fasting due to their large variety of diets and because fasting is a regular part of their life cycle. Mammals fed on a carbohydrate-rich diet experience a rapid decrease in blood glucose levels during a fast, thus, the development of mechanisms to resist the consequences of regular fasts, experienced on a daily basis, must have been crucial in the evolution of frugivorous bats. Phosphoenolpyruvate carboxykinase 1 (PEPCK1, encoded by the Pck1 gene) is the rate-limiting enzyme in gluconeogenesis and is largely responsible for the maintenance of glucose homeostasis during fasting in fruit-eating bats. To test whether Pck1 has experienced adaptive evolution in frugivorous bats, we obtained Pck1 coding sequence from 20 species of bats, including five Old World fruit bats (OWFBs) (Pteropodidae) and two New World fruit bats (NWFBs) (Phyllostomidae). Our molecular evolutionary analyses of these sequences revealed that Pck1 was under purifying selection in both Old World and New World fruit bats with no evidence of positive selection detected in either ancestral branch leading to fruit bats. Interestingly, however, six specific amino acid substitutions were detected on the ancestral lineage of OWFBs. In addition, we found considerable evidence for parallel evolution, at the amino acid level, between the PEPCK1 sequences of Old World fruit bats and New World fruit bats. Test for parallel evolution showed that four parallel substitutions (Q276R, R503H, I558V and Q593R) were driven by natural selection. Our study provides evidence that Pck1 underwent parallel evolution between Old World and New World fruit bats, two lineages of mammals that feed on a carbohydrate-rich diet and experience regular periods of fasting as part of their life cycle.

Introduction

Food limitation is a common physiological challenge faced by animals and leads to a risk of death due to starvation. Prolonged fasting, and starvation, causes a reduction in blood glucose levels, which can result in hypoglycemia and threaten life [1,2]. In this study, we focused on the molecular evolution of a gene involved in the resistance to the consequences of fasting. Maintaining glucose homeostasis under fasting conditions is essential [3], and many tissues are known to be involved in the regulation of blood glucose levels [4]. The liver plays a major role, as hepatic gluconeogenesis is the primary metabolic mechanism in the fasting state, in generating fuel to maintain the basic functions of other tissues, including skeletal muscle, red blood cells, and the brain [5]. In addition to the liver, renal and intestinal gluconeogeneses also contribute to the maintenance of glucose homeostasis during fasting [6,7]. The gluconeogenic pathway contains 11 steps that are catalyzed in turn by 11 enzymes [8] (S1 Fig.). The enzyme phosphoenolpyruvate carboxykinase (PEPCK), which catalyzes an irreversible reaction, is the rate-limiting step in gluconeogenesis [9]. Mammals have two isozymes of PEPCK, which are encoded by separate genes: phosphoenolpyruvate carboxykinase1 (PEPCK1, encoded by the Pck1 gene), which is located in the cytosol, and phosphoenolpyruvate carboxykinase 2 (PEPCK2, encoded by the Pck2 gene), which is restricted to the mitochondria [10]. PEPCK1 accounts for about 95% of the PEPCK enzymatic activity in rats and mice, which are popular model organisms [11,12], and its single copy gene consisting of 10 exons [13]. Studies in Pck1-knockout mice confirmed the importance of Pck1 in hepatic gluconeogenesis and glucose homeostasis, as Pck1 -/- mice display severe hypoglycemia [14]. The level of transcription of PEPCK1 determines the rate of gluconeogenesis [15].

Mammals consume a variety of different diets, which has led to the development of differing approaches to store and utilize body energy reserves, with different responses to fasting evolving in different species [16]. Previous studies have shown that the plasma glucose levels of animals fed a carbohydrate-rich diets fall severely during fasting, which contrasts to the steady blood glucose concentrations observed in vertebrates fed a protein-rich diet [17–19], suggesting that animals with a carbohydrate-rich diet should be less capable of resisting the consequences of a fast.

Among mammalian orders, bats have the widest range of diets, with diets based on a variety of specific food sources including insects, large arthropods, small vertebrates, blood, nectar, fruits, or leaves [20,21]. Old World fruit bats (OWFB) and New World fruit bats (NWFB), belonging to the families Pteropodidae and Phyllostomidae, respectively, have independently acquired carbohydrate-rich diets that are comprised mostly of fruit and/or nectar [22,23]. Bats dedicate only a very small portion (up to 4 h) of their daily activity cycle to foraging, and spend the remainder of their day in sheltered roosts, where they neither eat nor drink [24]. Therefore, bats experience both short and long term fasts as regular occurrences [25]. The evolution of flight in bats has also brought high energetic costs, as bats in flight have metabolic rates about 15 times greater than those at rest [26–29]. Since bats in flight have extremely high metabolic rates, one might expect that bats that feed on carbohydrate rich diets should show a greater consequence to fasting than other bats. Glucose homeostasis during a fast can be maintained by two primary mechanisms: glycogenolysis and gluconeogenesis [8]. Recent studies have demonstrated that in frugivorous bats that gluconeogenesis, rather than glycogenolysis, is the main mechanism for regulating glucose homeostasis during a fast [30,31]. These observations make it reasonable to hypothesize that key genes involved in gluconeogenesis, through adaptations for fasting, may have undergone adaptive evolution in the fruit eating bats. The enzymatic activity levels of all key enzymes in gluconeogenesis, including Glucose-6-phosphatase, fructose-1,6-bisphosphatase, PEPCK1 and PEPCK2, were measured in two species of fruit bats (Artibeus lituratus and Artibeus jamaicensis) before and after a 24 hour fasts by Pinheiro et al [32], who concluded that a significant change in the activity of only PEPCK1 could be demonstrated. These results indicated that PEPCK1 likely has a crucial role in the adaptation of gluconeogenesis for fasting in fruit bats.

We speculated that Pck1 experienced adaptive evolution in OWFBs and NWFBs in association with their acquired high carbohydrate diets. To test our hypothesis, we obtained the coding region of Pck1 from 20 species of bats (including five OWFBs species and two NWFBs species), and assessed whether these sequences display molecular evolutionary patterns consistent with adaptive evolution. In addition, we also examined whether the Pck1 gene has experienced convergent and/or parallel evolution between the OWFBs and NWFBs at the amino acid level.

Materials and Methods

Ethics statement

Protocols, including bat species collection and tissue sampling in China, were approved by the Animal Ethics Committee of East China Normal University (ID No: AR2012/03001). We captured Cynopterus sphinx from Zhongshan park (21°28′N, 109°7′E) in Guangxi province in Nov. 2009, and Rousettus leschenaultii were captured during Nov. 2009 from Jinlun cave (23°33′N, 108°15′E) located in Guangxi province, China. Eonycteris spelaea was captured in Mar. 2009 from Mengla park (21°55′N, 101°15′E) located in Yunnan province, China. Rhinolophus pusillus and Rhinolophus ferrumequinum were captured in Apr. 2013 from Lianhua cave (36°50′N, 113°54′E) located in Hebei province, China. Hipposideros armiger was captured in Apr. 2013 from Shihua cave (30°5′N, 119°5′E) located in Zhejiang province, China. Hipposideros pratti was captured in Apr. 2013 from Tianzidi cave (29°46′N, 119°21′E) located in Zhejiang province, China. Miniopterus schreibersi was captured in Apr. 2013 from Yunhua cave (33°18′N, 111°25′E) located in Henan province, China. Myotis ricketti was captured in Nov. 2012 from a cave (28°4′N, 116°58′E) located Jiangxi province, China. Pipistrellus abramus was captured in Jul. 2011 from a park (23°26′N, 111°30′E) located in Jiangxi province, China. Ia io was captured in Aug. 2013 near Xilin road (24°23′N, 106°17′E) located in Yunnan province, China. Scotophilus heathi was captured in Jul. 2011 from a building (23°1′N, 114°56′E) located in Guangdong province, China. Artibeus lituratus, Leptonycteris yerbabuenae and Pteronotus parnellii were captured in Apr. 2010 in Mexico. Bats captured in these locations were permitted by the local Protection and Research Center. No bat species used for this study are considered endangered. Bats were sacrificed by decapitation, quickly after capture, and all efforts were made to minimize potential pain and suffering.

Taxonomic coverage

In this study, we sequenced the coding region (97%) of the Pck1 gene from 15 species of bats, which were then combined with sequences from five other bat species obtained from NCBI and Ensembl. Our dataset included sequences from the following species. From the suborder Yinpterochiroptera, we included five OWFBs from the family pteropodidae (Pteropusvampyrus, Pteropus alecto, Cynopterus sphinx, Eonycteris spelaea, and Rousettus leschenaultii) and four insectivorous species from sister families to OWFBs (Rhinolophus pusillus, Rhinolophus ferrumequinum, Hipposideros armiger and Hipposideros pratti). From the suborder Yangochiroptera, we sampled two NWFBs from the family Phyllostomidae (Artibeus lituratus and Leptonycteris yerbabuenae) and nine insectivorous bat species (Pteronotus parnellii, Miniopterus schreibersi, Scotophilus heathi, Ia io, Pipistrellus abramus, Myotis ricketti, Myotis davidii, Myotis brandtii, and Myotis lucifugus). All new Pck1 sequences were deposited into GenBank with the following accession numbers: KJ957749-KJ957763.

To expand our dataset across mammals we incorporated published Pck1 sequences from nine other mammalian species from NCBI and Ensembl: Homo sapiens (NM_002591), Mus musculus (NM_011044), Rattus norvegicus (NM_198780), Equus caballus (XM_001489771), Bos taurus (NM_174737), Sus scrofa (NM_001123158), Canis familiaris (NM_001197143), Loxodonta Africana (XM_003419943), and Monodelphis domestica (XM_001377770). Details on all species and their corresponding accession numbers are listed in S1 Table.

Isolation, Amplification and Sequencing of Pck1 Coding Sequences

For the 15 bat species used to isolate new sequences, total RNA was isolated from liver tissue (stored at -80°C) using the Trizol reagent. Following the standard protocol, 5 μg of total RNA was reverse transcribed into cDNA using the Reverse Transcriptase kit. Primers (F: ATGCCTCCTCAGCTGCAAAACG and R: TACATCTGGCTGACTCTCTGCCTC) were designed to amplify the coding sequences of Pck1. All PCR products were separated using 1% agarose gels, with the correct sized fragments purified with Gel Extraction Kits, ligated into pGEM-T easy vector, cloned, and sequenced using Terminator kits on an ABI 3730 DNA sequencer. To avoid artifacts, multiple clones (>5) were sequenced for each species.

Phylogenetic Reconstruction

Pck1 nucleotide sequences from 29 mammalian species were aligned using ClustalX [33] and checked for accuracy by eye, with the coding sequences translated into amino acids using MEGA5 [34]. A Bayesian phylogenetic tree was reconstructed based on the aligned nucleotide sequences using MrBayes 3.1.2 [35] with the GTR+I+Γ nucleotide substitution model selected as the best-fitting model by jModelTest0.1 [36]. For the Bayesian analysis, we performed 10,000,000 generations of MCMC and sampled every 100 generations, with the first 2,000,000 generations discarded as burn-in. All other options and priors were the default settings of the MrBayes 3.1.2 software. The standard deviations of split frequencies were stably below 0.01 after 2,000,000 generations of MCMC performances.

Molecular evolution analyses

Before beginning the molecular evolution analyses, we first tested the sequences for evidence of recombination, since recombination may lead to the incorrect results in tests for positive selection. We used GARD [37] from the HyPhy package [38] to detect evidence for recombination breakpoints and tested their statistical significance. For tests of selection we used a topology for the phylogenetic tree that was based on the accepted species relationships [39–42]. For positive selection in Pck1, we estimated the rates of nonsynonymous (dN) and synonymous (dS) substitutions using CODEML from the PAML package [43]. We first used the two-ratio model, where the dN/dS ratio (termed omega or ω) was allowed to be different on the foreground and background branches. To detect differences in the selective pressure acting on the ancestral branches leading to OWFBs and NWFBs, each of these two groups were separately set as the foreground branch. The one ratio-model, where the omega values were equal among branches, served as the null hypothesis for these comparisons. We also employed the free-ratio model, where the value of dN/dS varies among all branches, to explore the overall selection pressure acting on Pck1 in the 29 mammalian species used in this study [44]. Similarly, the one-ratio model served as the null expectation for this comparison. “PairwiseRelativeRate.bf” from the HyPhy package [38] was also used to detect whether there were significant differences in the rates of synonymous and nonsynonymous substitutions between the frugivorous and insectivorous bats (see S1 Text for details).

The branch-site model A (test1 and test2) was then applied to detect positively selected changes along particular branches [45]. In this model, the branch of interest is set as the foreground, with the remaining branches set as the background. We conducted both test1 and test2 of the branch-site models [46] on the ancestral branches for the OWFBs and NWFBs. In both tests, the branch-site model A is the alternative hypothesis. In test1, the null hypothesis is the M1a (Nearly Neutral) model. For test2, the null hypothesis is the branch-site null model. Results from the alternative and null hypothesis were then compared in likelihood ratio tests (LRTs).

Ancestral sequences were reconstructed using CODEML from the PAML package [47], and the amino acid substitutions on each branch were then inferred. The identity and direction of the substitutions was also assessed using the maximum parsimony method with Mesquite 2.74 [48]. To further explore the heterogeneous selection pressure acting across codons of the Pck1 gene in frugivorous bats, a sliding window analysis of the ω values was performed using the program SWAAP 1.0.2 [49] with the Nei and Gojobori method [50]. Window and step sizes for this analysis were set as 90 and 36 nucleotides, respectively.

To detect evidence of parallel or convergent amino acid substitutions between OWFBs and NWFBs, we used a method described by Castoe et al. implemented in the codeMLancestral package [51] (see S1 Text). To test whether the convergence detected between pairs of focal branches was significant, we used a method described by Rossiter et al [52] (see S1 Text). We also used CONVERGE2 [53] to test whether the number of observed parallel or convergent amino acid substitutions shared by the OWFBs and NWFBs branches were significantly different from those expected by random chance. Finally, we examined whether the locations of the positively selected, convergently evolved and other specific sites in the protein sequences were associated with particular functional domains. The protein domains were defined from previous research [54,55].

Results

Our Pck1 sequence dataset contained 29 taxa, including five Old World fruit bats (Pteropodidae), two New World fruit bats (Phyllostomidae) and 13 insectivorous bats. The alignment of the Pck1 coding sequence spanned 1818 nucleotides (representing ~97% of the complete coding sequence) and equates to 606 amino acids. No insertions, deletions or stop codons were identified in the isolated coding sequences from the 20 bat species, and no evidence for recombination breakpoints were detected in our dataset.

Bayesian phylogenetic reconstruction of the nucleotides sequences of Pck1 recovered a tree with the main groupings of species in agreement with the accepted species topology [39–41] (Fig. 1B). The five OWFB species (P. vampyrus, P. alecto, C. sphinx, E. spelaea, and R. leschenaultii) grouped with the species from the family Rhinolophidae (R. pusillus and R. ferrumequinum) and Hipposideridae (H. armiger and H. pratti) to comprise the suborder Yinpterochiroptera [Bayesian posterior probability (BPP) = 0.56] (Fig. 1A). The other bat species, including the NWFBs (A. lituratus and L. yerbabuenae), grouped together to comprise the suborder Yangochiroptera (BPP = 1.00) (Fig. 1A). The two species of NWFBs grouped together by a high nodal support (BPP = 1.00) (Fig. 1A).

An external file that holds a picture, illustration, etc. Object name is pone.0118666.g001.jpg
Unconstrained Pck1 Bayesian phylogenetic tree and the species tree.

(A) Unconstrained Bayesian phylogenetic tree based on the Pck1 coding sequences, under the GTR+I+Γ model. Values on the nodes are posterior probabilities. (B) Species tree of the 29 mammals used in this study based on their accepted relationships. Red and green thick lines labeled ‘OWFBs’ and ‘NWFBs’ represent the ancestral branches for the Old World fruit bats and New World fruit bats, respectively. Nonsynonymous amino acid substitutions were mapped onto select branches. 13 sites on the ancestral branch of Old World fruit bats having omega values >1 are shown in blue. Convergent substitutions between OWFBs and NWFBs are shown in orange.

To detect evidence for changes in the selection pressure acting on Pck1 in OWFBs and NWFBs, we set either the ancestral branch to OWFBs or NWFBs as the foreground branch and performed two-ratio model tests on these branches. Our results showed that the two-ratio model with the ancestral branch to OWFBs as the foreground was a better fit than the one-ratio model (null hypothesis) [likelihood ratio test (LRT) statistic (2Δ[ell]) = 3.99, df = 1, P = 0.046] (Table 1), and had an omega value (dN/dS) on the ancestral branch leading to OWFBs that was higher than on the background (0.1090 vs. 0.0676, respectively). These results indicate that a change in selection pressure on the Pck1 gene occurred on the OWFB lineage. However, the two-ratio model test with the ancestral branch for NWFBs as the foreground was not significantly different from the one-ratio model (2Δ[ell] = 3.09, df = 1, P = 0.079) (Table 1), which implied that no change in selection pressure occurred on this NWFB lineage.

Table 1

Result of the branch model tests of selection pressure in the Pck1 gene of bats.
Model np [ell] ω 0 a ω OW a ω NW a Model Compare 2Δ[ell] P
A. One ratio:ω0 = ωOW = ωNW 57-13778.820.069= ω0 = ω0
B. Two ratio:ω0 = ωNW,ωOW 58-13776.830.0680.109= ω0 B vs. A3.990.046
C. Two ratio:ω0 = ωOW,ωNW 58-13777.280.068= ω0 0.163C vs. A3.090.079
D. Free ratio111-13718.94---D vs. A119.77<0.01

aωOW,ωNW,ω0, are the ω ratios for the OWFB, NWFB, and other branches, respectively.

Pairwise relative rate tests were also conducted to determine whether there were significant differences in the rates of synonymous and nonsynonymous substitutions between frugivorous and insectivorous bats. Our results showed that a significant difference in the rate of nonsynonymous substitutions occurs between OWFBs and insectivorous bats; however, the same result was not found for synonymous substitutions. These results from the pairwise relative rate tests indicated that a change in selection pressure on the Pck1 gene occurred in OWFBs, and that this change mainly effected the rate of nonsynonymous substitutions (Table 2).

Table 2

Results of the pairwise relative rate tests.
Group 1 a Group 2 b Syn rate test c Nonsyn rate test c
(Old World fruit bats)(insectivorous bats)LRT d P-valueLRT d p-value
Pteropus vampyrus Rhinolophus ferrumequinum 1.2360.26612.352<0.001
Hipposideros armiger 1.2010.27312.733<0.001
Miniopterus schreibersi 0.9110.3396.769<0.01
Pteronotus parnellii 2.9270.0879.368<0.001
Myotis lucifugus 0.4420.5066.5860.010
Pteropus alecto Rhinolophus ferrumequinum 0.3440.55812.628<0.001
Hipposideros armiger 0.3160.57413.430<0.001
Miniopterus schreibersi 0.1920.6616.954<0.01
Pteronotus parnellii 1.5040.2209.7090.002
Myotis lucifugus 0.0280.8676.929<0.01
Cynopterus sphinx Rhinolophus ferrumequinum 9.0550.00316.261<0.001
Hipposideros armiger 609690.00816.124<0.001
Miniopterus schreibersi 6.6030.0109.4090.002
Pteronotus parnellii 10.3820.00112.042<0.001
Myotis lucifugus 6.0170.0149.4490.002
Eonycteris spelaea Rhinolophus ferrumequinum 1.9720.1608.9130.003
Hipposideros armiger 1.1870.2769.1480.002
Miniopterus schreibersi 1.6600.1984.1560.041
Pteronotus parnellii 3.1190.0776.879<0.01
Myotis lucifugus 1.0030.3174.2600.039
Rousettus leschenaultii Rhinolophus ferrumequinum 5.2650.02210.5580.001
Hipposideros armiger 4.3160.03810.847<0.001
Miniopterus schreibersi 4.3210.0385.2900.021
Pteronotus parnellii 7.2640.0078.4020.004
Myotis lucifugus 3.5270.0605.3180.021

aFive Old World fruit bats are used as ingroup 1 to separately compare with the insectivorous bats.

bFive typical Insectivorous bats are used separately as ingroup 2 for the relative rate tests to compare with the Old World fruit bats.

cSyn rate test, Synonymous rate tests. Nonsyn rate test, Nonsynonymous rate tests.

dLRT, Likelihood Ratio Test.

Branch-site model A test was performed to detect evidence for positive selection on the ancestral branches of OWFBs and NWFBs. The result of branch-site model A revealed 13 sites on the ancestral branch of OWFBs with elevated omega values (10G, 13G, 62A, 66L, 95V, 188D, 270V, 349E, 352R, 375E, 534V, 578L and 591V, amino acid positions refer to the human PEPCK1 sequence), and the fit of branch-site model A was better than the M1a model (2Δ[ell] = 7.95, df = 2, P = 0.0188; test1) (Table 3). However, the in-depth test for positive selection (test2) did not reveal a similar result, as the branch-site model A of OWFBs was not significantly better than the null branch-site null model (2Δ[ell] = 0.0, df = 1, P = 1.0; test2) (Table 3). Since test1 cannot distinguish relaxed constraints from positive selection [46], we were unable to rule out the possibility of relaxed evolution. No evidence for a change in selection pressure was found in the ancestral branch of NWFBs with both test1 and test2, as the alternative hypothesis of both the M1a model and branch-site null model were not significantly better than the null models.

Table 3

Result of the branch-site model A test for the detection of positively selected sites on the ancestral branches of Old World fruit bats and New World fruit bats.
Branch-site modelnp a ParametersLRT b [ell] P-valueSite with elevated omega values c
M1a(nearly Neutral)58 P 0 = 0.906,P 1 = 0.094, ω0 = 0.038, ω1 = 1.00-13500.24
Model A (alternative hypothesis) for Old World fruit bats60 P 0 = 0.862, P 1 = 0.090, P 2a = 0.044, P 2b = 0.005 Background: ω0 = 0.037, ω1 = 1.00, ω2a = 0.037, ω2b = 1.00 Foreground: ω0 = 0.037, ω1 = 1.00, ω2a = 1.00,ω2b = 1.00Test 1-13496.270.019 10G, 13G, 62A, 66L, 95V, 188D, 270V, 349E, 352R, 375E, 534V, 578L, 591V
Model A (null hypothesis) for Old World fruit bats59 P 0 = 0.862, P 1 = 0.090, P 2a = 0.044,P 2b = 0.005 Background: ω0 = 0.037, ω1 = 1.00, ω2a = 0.037, ω2b = 1.00 Foreground: ω0 = 0.037, ω1 = 1.00, ω2a = 1.00,ω2b = 1.00Test 2-13496.271Not allowed
Model A (alternative hypothesis) for New World fruit bats60 P 0 = 0.837, P 1 = 0.087, P 2a = 0.068,P 2b = 0.007 Background: ω0 = 0.038, ω1 = 1.000, ω2a = 0.038, ω2b = 1.00 Foreground: ω0 = 0.038, ω1 = 1.00, ω2a = 1.00,ω2b = 1.00Test 1-13498.580.191Not allowed
Model A (null hypothesis) for New World fruit bats59 P 0 = 0.837, P 1 = 0.087, P 2a = 0.068,P 2b = 0.007 Background: ω0 = 0.038, ω1 = 1.00, ω2a = 0.038, ω2b = 1.00 Foreground: ω0 = 0.038, ω1 = 1.00, ω2a = 1.00,ω2b = 1.00Test 2-13498.581Not allowed

anp, number of parameters.

bLRT, likelihood ratio test.

clocation of sites with elevated omega values detected by branch-site model A test refer to Homo sapiens sequence.

The sites with BEB posterior probabilities > 0.95 were highlight by underline.

To identify amino acid changes on particular branches, we use the software Mesquite 2.74 [48] to reconstruct the ancestral amino acid sequences. These reconstructions showed that there were six specific amino acid replacements (A13G, S252G, K349E, F578L, D591V, and I608V) on the ancestral branch of OWFBs (Fig. 2), and four (A13G, K349E, F578L, and D591V) of them were detected by branch-site model A as having elevated omega values (Table 3). To further explore the selection pressure changes among Pck1 codons, we calculated the values of ω among OWFBs, NWFBs, and insectivorous bats using SWAAP 1.0.2. The results of this analysis showed that the estimated ω values for OWFBs are much larger than those of insectivorous bats, with a maximal ω value of 1.9 observed on the nucleotide binding domain 1 (Fig. 3). The ω values of the NWFBs were also larger than those of the insect-eating bats, but not as large as those of the OWFBs.

An external file that holds a picture, illustration, etc. Object name is pone.0118666.g002.jpg
The species tree of 29 mammals with Old World fruit bat specific amino acid replacements identified by ancestral sequence reconstruction using the maximum parsimony method.

(A) A13G, (B) S252G, (C) K349E, (D) F578L, (E) D591V and (F) I608V. Branch lengths are not drawn to scale.

An external file that holds a picture, illustration, etc. Object name is pone.0118666.g003.jpg
Sliding window analysis of the variation in omega values along Pck1 genes in OWFBs, NWFBs, and insectivorous bats.

Window size and the step size were set to 90bp and 36bp, respectively. Schematic of the PEPCK1 protein domain structure is shown beneath the plot.

To examine parallel and/or convergent amino acid substitutions between the OWFB and NWFB ancestral branches, we reconstructed ancestral sequences using PAML [43], with the amino acid replacements on the major branches of the fruit-eating bats inferred. Five parallel amino acid substitutions were shared between OWFBs and NWFBs (N127S, Q276R, R503H, I558V, and Q593R). The total posterior probabilities of divergence versus convergence for focus branch-pairs were shown in Fig. 4. We used two methods to test the probability of acquiring the observed parallel substitutions in the frugivorous branches by chance. Both methods suggest that the probabilities of the substitutions Q276R and I558V, shared between the ancestral branches of OWFBs and A. lituratus, the Q593R substitution, shared between the ancestral branch of NWFBs and C. sphinx, and the R503H substitution, shared between the ancestral branch of the genera Pteropus and L. yerbabuenae, all deviated significantly from the random expectation (Table 4 and Table 5). The N127S substitution, shared between the ancestral branches of OWFBs and L. yerbabuenae, however, was not significant (P > 0.05), suggesting that this amino acid substitution likely occurred in parallel by chance (Table 4 and Table 5).

An external file that holds a picture, illustration, etc. Object name is pone.0118666.g004.jpg
Plot of the total posterior probabilities of divergence versus convergence for all pairs of branches in the tree.

Pairwise comparison for the ancestral branch of Old World fruit bats versus Artibeus lituratus, the ancestral branch of NWFBs versus Cynopterus sphinx, ancestral branch of Pteropus versus Leptonycteris yerbabuenae, and the ancestral branch of OWFBs versus Leptonycteris yerbabuenae are highlighted by green, blue, purple and orange, respectively.

Table 4

Total probabilities of convergence for parallel amino acid substitutions in the Pck1 gene for key pairs of branches of frugivorous bats.
Branch Pair a Total Convergence ProbabilityParallel Substitutions with BPP Values b P-value c
1 vs. 52.834Q276R (0.873), I558V (0.996)0.007
2 vs. 41.003Q593R (0.977)0.025
3 vs. 60.983R503H (0.983)0.033
1 vs. 60.892N127S (0.887)0.229

aSee Fig. 1(B) for branch labels.

bBPP, Bayesian posterior probability.

cSignificance levels of the branch-wise convergences are comparisons of the observed probabilities against the null distribution based on 1000 simulated data sets.

Table 5

Statistical tests for parallel evolution between Old World fruit bats and New World fruit bats.
Branch pair a Parallel substitutionObserved numberExpected number P-value b
1 vs. 5Q276R, I558V20.0307<0.01
2 vs. 4Q593R10.00750.030
3 vs. 6R503H10.00960.038
1 vs. 6N127S10.02570.103

aSelected branches are labeled with Arabic numerals in Fig. 1(B)

bThe p value of each test was multiplied by four (branch pairs) to correct for multiple testing.

Discussion

In this study, we examined the molecular evolution of the Pck1 gene, which plays an important role in the maintenance of glucose homeostasis during fasting in frugivirous bats, which eat carbohydrate-rich diets. As fasting is a regular part of the life cycle of fruit eating bats, they require an enhanced ability to resist hypoglycemia [24,25]. By comparing Pck1 sequences among frugivorous and insectivorous bats, and other mammals, we provide the first data on the molecular evolution of this gene in frugivorous bats. We found that the omega value (dN/dS) on the ancestral branch of OWFBs is higher than for the remaining branches of the species tree, and that the rate of nonsynonymous substitutions in the Pck1 gene of OWFBs is significantly higher than that of insectivorous bats. Although there was a change in the selection pressure acting on the Pck1 gene on the ancestral branch of NWFBs, the omega value is still quite low indicating that Pck1 is still under purifying selection.

Although the results of test1 of the branch-site model A test revealed 13 amino acid replacement with elevated omega value, we are unable to rule out the possibility that these were due to a relaxation of the selection acting on the pck1 gene, as the result of the explicit test of positive selection (test2) was negative. Test1, by itself, is unable to distinguish between relaxation of constraints and positive selection [46]. However, the likelihood that Pck1 experienced relaxed selection in Old World fruit bats is probably low as PEPCK1 has a very important function across mammals [56]. Even between bacteria and mammals, a number of residues at the active site are conserved and the reaction mechanism for this enzyme in both bacteria and mammals is thought to be same [57]. From our sequenced Pck1 genes we found that the crucial residues in the active site loops, including the R-loop, P-loop and Ω-loop, are completely conserved. These results indicated that PEPCK1 retains an important role in glucose homeostasis in frugivore bats [58]. Previous studies have shown that knockout of Pck1 in mice leads to hypoglycemia and death [14], supporting the necessity of PEPCK for life in mammals, and making it unlikely that this gene in OWFBs underwent relaxed evolution. All these results indicate that Pck1 has experienced strong purifying selection not only in fruit bats, but also in all other mammalian species examined. Although the results are congruent with the fact that the PEPCK1protein is very critical for hepatic glycometabolism, there are two possible explanations that may explain the failure to detect positive selection. First, OWFBs diverged from the other Yinpterochiroptera lineages nearly 60 million years ago [40], and previous studies have shown that signals of positive selection can be swamped by a history of purifying selection [46]. Second, the low level of sequence divergence among the sequences limits the power of statistical methods to detect signals of positive selection [59–61]. In addition, we only had two NWFB species, which further limits our ability to detect positive selection in this group.

Results from the ancestral reconstruction of the amino acid sequences using the maximum parsimony method with Mesquite 2.74 [48] identified six specific amino acid replacement that occur on the OWFBs ancestral lineage. A change from the basic amino acid lysine to the acidic amino acid glutamic acid (K349E) is located in the PEPCK-specific active site domain (residues 326–425) [54]. The A13G and S252G substitutions occur in the N-terminal domain (residues 1–259) associated with degradation and acetylation of PEPCK1 [62]. F578L, D591V and I608V substitutions lie near the nucleotide-binding domain 2 (residues 426–622) (Fig. 5) responsible for GTP binding, which is necessary for enzymatic activity [55]. Sliding window analysis also found a significantly elevated omega value for the sequences encoding the nucleotide-binding domain 1 (residues 260–325). The elevated omega value is likely due to the I270V substitution, which was detected by the branch-site model A. Replacement of I270V is a OWFBs-specific change among species of the order Chiroptera, and is only shared with the pig (data not shown).

An external file that holds a picture, illustration, etc. Object name is pone.0118666.g005.jpg
Distribution of the six OWFBs-specific amino acid substitutions and four parallel amino acid substitutions occurring between OWFBs and NWFBs in the secondary structure of PEPCK1 protein.

OWFB substitutions are labeled in black, while parallel substitutions are labeled in red within a model of the domain structure of PEPCK1.

In contrast to OWFBs, we did not find any specific amino acid replacements in the NWFBs. It is interesting to note that other recent molecular evolution studies focusing on diet-related genes in bats also found similar patterns [63–65]. For example, GLUT4 which play an important role in glucose homeostasis had two specific amino acid substitutions only in the OWFBs but not in NWFBs [64]. This evolutionary discrepancy in patterns of sequence evolution between OWFBs and NWFBs for the Pck1 gene may have two explanations. First, OWFBs evolved their carbohydrate-rich diet earlier than the NWFBs (28 mya and 20 mya, respectively) [40], thus have had a greater amount of time to accumulated accelerated changes. Second, the specific sequence changes may happen in the regulatory regions in NWFBs, instead of the coding region of these genes, as the promoter can regulate the expression of PEPCK1 [66], change may only occur at the protein abundance and not the protein sequence level.

Although no evidence for positive selection was found in the Pck1 gene of OWFBs and NWFBs, this does not mean that this gene did not undergo adaptive evolution. Parallel evolution between OWFBs and NWFBs is also a strong indicator of adaptive evolution [53]. Previous molecular evolution studies in bats have also found significant evidence for parallel or convergent evolution with no evidence for positive selection [67,68]. A possible explanation for the parallel evolution at the amino acid level between Pck1 genes in OWFBs and NWFBs is molecular adaption to the regular fasting by these fruit bats. It is widely accepted that the hepatic gluconeogenesis plays a critical role in the maintenance of glucose homeostasis during fasting [5,8], and it has been shown that hepatic gluconeogenesis is the primary metabolic mechanism for fasting resistance in the Old World Fruit Bats [31]. PEPCK1, the key enzyme of gluconeogenesis, has been shown in a previous study to have a crucial role in fasting in New World fruit bats [32]. Thus, it is reasonable to assume that the parallel evolution of Pck1 in OWFBs and NWFBs might relate to the evolution of adaptive mechanisms for fasting resistance in fruit bats. Here we found four credible parallel evolution sites (Q276R, R503H, I558V and Q593R) in PEPCK1. All of these sites are located near the nucleotide binding domain, which might influence the binding of GTP to PEPCK1 and thus enzymatic activity [54] (Fig. 5). It is interesting to note that, besides the four convergently evolving sites, the Old World fruit bats also have six other specific amino acid substitutions, which may reflect the fact that the evolutionary history of OWFBs is longer than that of the NWFBs [40].

In conclusion, our results show that the Pck1 gene, which plays an important role in the maintenance of glucose homeostasis during fasting, has undergone parallel evolution between the OWFBs and NWFBs, which provides evidence for adaptive evolution in both lineages, although our evidence for NWFBs is limited to only two species of Phyllostomid bats. Greater sampling of NWFBs is needed to confirm this finding. Apart from gluconeogenesis, PEPCK1 is also the rate-limiting enzyme in glyceroneogenesis, which is responsible for recycling free fatty acid back to triglyceride [69]. Chronic release of free fatty acid by adipose tissue is a critical factor in the development of Type 2 diabetes [70], therefore, PEPCK1 could prevent diabetes by controlling the rate of glyceroneogenesis [70]. Here we found evidence for parallel evolution of the Pck1 gene between OWFBs and NWFBs, which could possibly be explained by its function in glyceroneogenesis. However, additional in-depth studies are needed to address this issue. We speculate that additional deeper studies will reveal that numerous other genes related to carbohydrate metabolism will be found to have experienced adaptive evolution in OWFBs and NWFBs.

Supporting Information

S1 Fig

Enzymatic pathway for gluconeogenesis.

Four unique enzymes required in gluconeogenesis are shown in the rectangles. All of the remaining steps are catalyzed by glycolytic enzymes.

(TIF)

S1 Table

Information on species examined in the study.

(DOC)

S1 Text

Details on some methods used in this study.

(DOC)

Acknowledgments

We thank Bin Shen and Fanxing Meng for helpful comments on the molecular evolutionary analyses, and Lei Ming for help with the bioinformatic analyses.

Funding Statement

This work was supported by the Key Construction Program of the National ‘985 Project’ and ‘211 Project’ and Shanghai Excellent Academic Leaders Project (Grant No. 11XD1402000) to SZ. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.

Data Availability

All new pck1 sequences files are available from the NCBI database (accession number: KJ957749-KJ957763).

References

1. Boyle PJ, Shah SD, Cryer PE (1989) Insulin, glucagon, and catecholamines in prevention of hypoglycemia during fasting. Am J Physiol 256: E651–661. [Abstract] [Google Scholar]
2. Houten SM, Herrema H, Te Brinke H, Denis S, Ruiter JP, et al. (2013) Impaired amino acid metabolism contributes to fasting-induced hypoglycemia in fatty acid oxidation defects. Hum Mol Genet 22: 5249–5261. 10.1093/hmg/ddt382 [Abstract] [CrossRef] [Google Scholar]
3. Shrayyef MZ, Gerich JE (2010) Normal glucose homeostasis Principles of Diabetes Mellitus: Springer; pp. 19–35. [Google Scholar]
4. Brunner Y, Schvartz D, Priego-Capote F, Couté Y, Sanchez J-C (2009) Glucotoxicity and pancreatic proteomics. Journal of proteomics 71: 576–591. 10.1016/j.jprot.2008.10.002 [Abstract] [CrossRef] [Google Scholar]
5. Nordlie RC, Foster JD, Lange AJ (1999) Regulation of glucose production by the liver. Annu Rev Nutr 19: 379–406. [Abstract] [Google Scholar]
6. Segura J, Ruilope L (2013) Contribution of the kidney to glucose homeostasis. Medicina Clinica 141: 26–30. 10.1016/S0025-7753(13)70060-5 [Abstract] [CrossRef] [Google Scholar]
7. Mithieux G, Andreelli F, Magnan C (2009) Intestinal gluconeogenesis: key signal of central control of energy and glucose homeostasis. Current Opinion in Clinical Nutrition & Metabolic Care 12: 419–423. [Abstract] [Google Scholar]
8. Ferrier D, Harvey R (2011) Lippincotts illustrated review biochemistry Walters Kluwer: Lippincott Williams and Wilkins. [Google Scholar]
9. Hanson RW, Patel YM (1994) Phosphoenolpyruvate carboxykinase (GTP): the gene and the enzyme. Advances in Enzymology and Related Areas of Molecular Biology, Volume 69: 203–281. [Abstract] [Google Scholar]
10. Beale EG, Harvey BJ, Forest C (2007) PCK1 and PCK2 as candidate diabetes and obesity genes. Cell Biochemistry and Biophysics 48: 89–95. [Abstract] [Google Scholar]
11. Hanson RW, Reshef L (1997) Regulation of phosphoenolpyruvate carboxykinase (GTP) gene expression. Annual Review of Biochemistry 66: 581–611. [Abstract] [Google Scholar]
12. Hanson RW, Mehlman MA (1976) Gluconeogenesis: its regulation in mammalian species: John Wiley & Sons. [Google Scholar]
13. Yu H, Thun R, Chandrasekharappa SC, Trent JM, Zhang J, et al. (1993) Human PCK1 encoding phosphoenolpyruvate carboxykinase is located on chromosome 20q13. 2. Genomics 15, 219–221 [Abstract] [Google Scholar]
14. She P, Shiota M, Shelton KD, Chalkley R, Postic C, et al. (2000) Phosphoenolpyruvate carboxykinase is necessary for the integration of hepatic energy metabolism. Molecular and Cellular Biology 20: 6508–6517. [Europe PMC free article] [Abstract] [Google Scholar]
15. Goto M, Yoshioka T, Battelino T, Ravindranath T, Zeller WP (2001) TNF&agr; Decreases Gluconeogenesis in Hepatocytes Isolated from 10-Day-Old Rats. Pediatric Research 49: 552–557. [Abstract] [Google Scholar]
16. McCue MD (2012) Comparative physiology of fasting, starvation, and food limitation: Springer. [Google Scholar]
17. Kettelhut IC, Foss MC, Migliorini R (1980) Glucose homeostasis in a carnivorous animal (cat) and in rats fed a high-protein diet. American Journal of Physiology-Regulatory, Integrative and Comparative Physiology 239: R437–R444. [Abstract] [Google Scholar]
18. Veiga J, Roselino E, Migliorini R (1982) Effects of fasting and adrenalectomy on the kinetics of glucose metabolism in granivorous and carnivorous birds. Braz. J. Med. Biol. Res., 15 (1982), pp. 175–180. [Abstract] [Google Scholar]
19. da Silva RS, Migliortni RH (1990) Effects of starvation and refeeding on energy-linked metabolic processes in the turtle (iTphrynops hilarii). Comparative Biochemistry and Physiology Part A: Physiology 96: 415–419. [Google Scholar]
20. Neuweiler G (2000) The biology of bats: Oxford University Press. [Google Scholar]
21. Kunz TH, Fenton MB (2006) Bat ecology: University of Chicago Press. [Google Scholar]
22. Altringham JD, McOwat T, Hammond L (1996) Bats: Biology and Behaviour: Oxford University Press. [Google Scholar]
23. Voigt CC, Zubaid A, Kunz TH, Kingston T (2011) Sources of assimilated proteins in Old and New World phytophagous bats. Biotropica 43: 108–113. 10.1093/jncimonographs/lgr038 [Abstract] [CrossRef] [Google Scholar]
24. Bassett JE (2004) Role of urea in the postprandial urine concentration cycle of the insectivorous bat(Antrozous pallidus). Comparative Biochemistry and Physiology Part A: Molecular & Integrative Physiology 137: 271–284. [Abstract] [Google Scholar]
25. Ben-Hamo M, Muñoz-Garcia A, Pinshow B (2012) Physiological responses to fasting in bats Comparative Physiology of Fasting, Starvation, and Food Limitation: Springer; pp. 257–275. [Google Scholar]
26. Austad SN, Fischer KE (1991) Mammalian aging, metabolism, and ecology: evidence from the bats and marsupials. Journal of Gerontology 46: B47–B53. [Abstract] [Google Scholar]
27. Harrison JF, Roberts SP (2000) Flight respiration and energetics. Annual Review of Physiology 62: 179–205. [Abstract] [Google Scholar]
28. Munshi-South J, Wilkinson GS (2010) Bats and birds: exceptional longevity despite high metabolic rates. Ageing Research Reviews 9: 12–19. 10.1016/j.arr.2009.07.006 [Abstract] [CrossRef] [Google Scholar]
29. Szewczak JM (1997) Matching gas exchange in the bat from flight to torpor. American zoologist 37: 92–100. [Google Scholar]
30. Keegan D (1977) Aspects of the assimilation of sugars by(Rousettus aegyptiacus). Comparative Biochemistry and Physiology Part A: Physiology 58: 349–352. [Google Scholar]
31. Korine C, Zinder O, Arad Z (1999) Diurnal and seasonal changes in blood composition of the free-living Egyptian fruit bat (Rousettus aegyptiacus). Journal of Comparative Physiology B 169: 280–286. [Abstract] [Google Scholar]
32. Pinheiro EC, Taddei VA, Migliorini RH, Kettelhut IC (2006) Effect of fasting on carbohydrate metabolism in frugivorous bats (Artibeus lituratus and Artibeus jamaicensis). Comparative Biochemistry and Physiology Part B: Biochemistry and Molecular Biology 143: 279–284. [Abstract] [Google Scholar]
33. Thompson JD, Gibson TJ, Plewniak F, Jeanmougin F, Higgins DG (1997) The CLUSTAL_X windows interface: flexible strategies for multiple sequence alignment aided by quality analysis tools. Nucleic Acids Research 25: 4876–4882. [Europe PMC free article] [Abstract] [Google Scholar]
34. Tamura K, Peterson D, Peterson N, Stecher G, Nei M, et al. (2011) MEGA5: molecular evolutionary genetics analysis using maximum likelihood, evolutionary distance, and maximum parsimony methods. Molecular Biology and Evolution 28: 2731–2739. 10.1093/molbev/msr121 [Europe PMC free article] [Abstract] [CrossRef] [Google Scholar]
35. Ronquist F, Huelsenbeck JP (2003) MrBayes 3: Bayesian phylogenetic inference under mixed models. Bioinformatics 19: 1572–1574. [Abstract] [Google Scholar]
36. Posada D (2008) jModelTest: phylogenetic model averaging. Molecular biology and evolution 25: 1253–1256. 10.1093/molbev/msn083 [Abstract] [CrossRef] [Google Scholar]
37. Pond SLK, Posada D, Gravenor MB, Woelk CH, Frost SD (2006) GARD: a genetic algorithm for recombination detection. Bioinformatics 22: 3096–3098. [Abstract] [Google Scholar]
38. Pond SLK, Muse SV (2005) HyPhy: hypothesis testing using phylogenies Statistical methods in molecular evolution: Springer; pp. 125–181. [Google Scholar]
39. Datzmann T, von Helversen O, Mayer F (2010) Evolution of nectarivory in phyllostomid bats (Phyllostomidae Gray, 1825, Chiroptera: Mammalia). BMC Evolutionary Biology 10: 165 10.1186/1471-2148-10-165 [Europe PMC free article] [Abstract] [CrossRef] [Google Scholar]
40. Teeling EC, Springer MS, Madsen O, Bates P, O'brien SJ, et al. (2005) A molecular phylogeny for bats illuminates biogeography and the fossil record. Science 307: 580–584. [Abstract] [Google Scholar]
41. Zhou X, Xu S, Xu J, Chen B, Zhou K, et al. (2012) Phylogenomic analysis resolves the interordinal relationships and rapid diversification of the Laurasiatherian mammals. Systematic Biology 61: 150–164. 10.1093/sysbio/syr089 [Europe PMC free article] [Abstract] [CrossRef] [Google Scholar]
42. Bickham JW, Patton JC, Schlitter DA, Rautenbach IL, Honeycutt RL (2004) Molecular phylogenetics, karyotypic diversity, and partition of the genus Myotis (Chiroptera: Vespertilionidae). Molecular Phylogenetics and Evolution 33: 333–338. [Abstract] [Google Scholar]
43. Yang Z (2007) PAML 4: phylogenetic analysis by maximum likelihood. Molecular Biology and Evolution 24: 1586–1591. [Abstract] [Google Scholar]
44. Yang Z (1998) Likelihood ratio tests for detecting positive selection and application to primate lysozyme evolution. Molecular Biology and Evolution 15: 568–573. [Abstract] [Google Scholar]
45. Yang Z, Wong WS, Nielsen R (2005) Bayes empirical Bayes inference of amino acid sites under positive selection. Molecular Biology and Evolution 22: 1107–1118. [Abstract] [Google Scholar]
46. Zhang J, Nielsen R, Yang Z (2005) Evaluation of an improved branch-site likelihood method for detecting positive selection at the molecular level. Molecular Biology and Evolution 22: 2472–2479. [Abstract] [Google Scholar]
47. Yang Z, Kumar S, Nei M (1995) A new method of inference of ancestral nucleotide and amino acid sequences. Genetics 141: 1641–1650. [Europe PMC free article] [Abstract] [Google Scholar]
48. Maddison W, Maddison D (2010) Mesquite: a modular system for evolutionary analysis. 2010. Version 2.74. Available at: mesquiteproject org/mesquite/download/download html.
49. Pride D (2000) SWAAP: a tool for analyzing substitutions and similarity in multiple alignments. Version 1.0. 2. http://www.bacteriamuseum.org/SWAAP/SwaapPage.htm.
50. Nei M, Gojobori T (1986) Simple methods for estimating the numbers of synonymous and nonsynonymous nucleotide substitutions. Molecular Biology and Evolution 3: 418–426. [Abstract] [Google Scholar]
51. Castoe TA, de Koning AJ, Kim H-M, Gu W, Noonan BP, et al. (2009) Evidence for an ancient adaptive episode of convergent molecular evolution. Proceedings of the National Academy of Sciences 106: 8986–8991. 10.1073/pnas.0900233106 [Europe PMC free article] [Abstract] [CrossRef] [Google Scholar]
52. Davies K, Cotton J, Kirwan J, Teeling E, Rossiter S (2012) Parallel signatures of sequence evolution among hearing genes in echolocating mammals: an emerging model of genetic convergence. Heredity 108: 480–489. 10.1038/hdy.2011.119 [Europe PMC free article] [Abstract] [CrossRef] [Google Scholar]
53. Zhang J, Kumar S (1997) Detection of convergent and parallel evolution at the amino acid sequence level. Molecular Biology and Evolution 14: 527–536. [Abstract] [Google Scholar]
54. Johnson TA, Holyoak T (2012) The Ω-Loop Lid Domain of Phosphoenolpyruvate Carboxykinase Is Essential for Catalytic Function. Biochemistry 51: 9547–9559. 10.1021/bi301278t [Europe PMC free article] [Abstract] [CrossRef] [Google Scholar]
55. Carlson GM, Holyoak T (2009) Structural insights into the mechanism of phosphoenolpyruvate carboxykinase catalysis. Journal of Biological Chemistry 284: 27037–27041. 10.1074/jbc.R109.040568 [Europe PMC free article] [Abstract] [CrossRef] [Google Scholar]
56. Utter MF, Kolenbrander HM (1972) 4 Formation of Oxalacetate by CO2 Fixation on Phosphoenolpyruvate. The Enzymes 6: 117–168. [Google Scholar]
57. Dunten P, Belunis C, Crowther R, Hollfelder K, Kammlott U, et al. (2002) Crystal structure of human cytosolic phosphoenolpyruvate carboxykinase reveals a new GTP-binding site. Journal of Molecular Biology 316: 257–264. [Abstract] [Google Scholar]
58. Klover PJ, Mooney RA (2004) Hepatocytes: critical for glucose homeostasis. The International Journal of Biochemistry & Cell Biology 36: 753–758. [Abstract] [Google Scholar]
59. Bielawski JP, Yang Z (2003) Maximum likelihood methods for detecting adaptive evolution after gene duplication Genome Evolution: Springer; pp. 201–212. [Abstract] [Google Scholar]
60. Anisimova M, Bielawski JP, Yang Z (2002) Accuracy and power of Bayes prediction of amino acid sites under positive selection. Molecular Biology and Evolution 19: 950–958. [Abstract] [Google Scholar]
61. Yang Z, Bielawski JP (2000) Statistical methods for detecting molecular adaptation. Trends in Ecology & Evolution 15: 496–503. [Europe PMC free article] [Abstract] [Google Scholar]
62. Jiang W, Wang S, Xiao M, Lin Y, Zhou L, et al. (2011) Acetylation regulates gluconeogenesis by promoting PEPCK1 degradation via recruiting the UBR5 ubiquitin ligase. Molecular Cell 43: 33–44. 10.1016/j.molcel.2011.04.028 [Europe PMC free article] [Abstract] [CrossRef] [Google Scholar]
63. Shen B, Han X, Jones G, Rossiter SJ, Zhang S (2013) Adaptive evolution of the Myo6 gene in Old World fruit bats (Family: Pteropodidae). PloS One 8: e62307 10.1371/journal.pone.0062307 [Europe PMC free article] [Abstract] [CrossRef] [Google Scholar]
64. Shen B, Han X, Zhang J, Rossiter SJ, Zhang S (2012) Adaptive evolution in the glucose transporter 4 gene Slc2a4 in Old World fruit bats (Family: Pteropodidae). PloS One 7: e33197 10.1371/journal.pone.0033197 [Europe PMC free article] [Abstract] [CrossRef] [Google Scholar]
65. Liu Y, Xu H, Yuan X, Rossiter SJ, Zhang S (2012) Multiple adaptive losses of alanine-glyoxylate aminotransferase mitochondrial targeting in fruit-eating bats. Molecular Biology and Evolution 29: 1507–1511. 10.1093/molbev/mss013 [Abstract] [CrossRef] [Google Scholar]
66. Zhang Q, Koser S, Donkin S (2014) Propionate is a dominant inducer of bovine cytosolic phosphoenolpyruvate carboxykinase (PCK1) expression (818.11). The FASEB Journal 28: 818.811. [Google Scholar]
67. Liu Y, Han N, Franchini LF, Xu H, Pisciottano F, et al. (2012) The voltage-gated potassium channel subfamily KQT member 4 (KCNQ4) displays parallel evolution in echolocating bats. Molecular Biology and Evolution 29: 1441–1450. 10.1093/molbev/msr310 [Europe PMC free article] [Abstract] [CrossRef] [Google Scholar]
68. Shen Y-Y, Liu J, Irwin DM, Zhang Y-P (2010) Parallel and convergent evolution of the dim-light vision gene RH1 in bats (Order: Chiroptera). PLoS One 5: e8838 10.1371/journal.pone.0008838 [Europe PMC free article] [Abstract] [CrossRef] [Google Scholar]
69. Cadoudal T, Fouque F, ccedil, oise, Benelli C, et al. (2008) Glyceroneogenesis and PEPCK-C: pharmacological targets in type 2 diabetes. Medecine Sciences: 24(4): 407–413 p. [Abstract] [Google Scholar]
70. Hanson RW, Reshef L (2003) Glyceroneogenesis revisited. Biochimie 85: 1199–1205. [Abstract] [Google Scholar]
Articles from PLOS ONE are provided here courtesy of PLOS

Full text links 

Read article at publisher's site: https://doi.org/10.1371/journal.pone.0118666

Read article for free, from open access legal sources, via Unpaywall: https://journals.plos.org/plosone/article/file?id=10.1371/journal.pone.0118666&type=printableUnpaywall

Citations & impact 

Impact metrics

3
Citations
Jump to Citations
20
Data citations
Jump to Data

Citations of article over time

Alternative metrics

Altmetric item for https://www.altmetric.com/details/3838130
Altmetric
Discover the attention surrounding your research
https://www.altmetric.com/details/3838130

Smart citations by scite.ai
Smart citations by scite.ai include citation statements extracted from the full text of the citing article. The number of the statements may be higher than the number of citations provided by EuropePMC if one paper cites another multiple times or lower if scite has not yet processed some of the citing articles.
Explore citation contexts and check if this article has been supported or disputed.
https://scite.ai/reports/10.1371/journal.pone.0118666

Supporting
Mentioning
Contrasting
0
8
0

Article citations

Data 

Data behind the article

This data has been text mined from the article, or deposited into data resources.

BioStudies: supplemental material and supporting data

Nucleotide Sequences (2)

RefSeq - NCBI Reference Sequence Database (Showing 9 of 9)

Similar Articles 

To arrive at the top five similar articles we use a word-weighted algorithm to compare words from the Title and Abstract of each citation.