The evolutionary trajectory of mitochondrial carrier family during metazoan evolution
BMC Evolutionary Biology volume 10, Article number: 282 (2010)
Exploring metabolic evolution is a way to understand metabolic complexity. The substrate transport of mitochondrial carrier family (MCF) influences direct metabolic activities, making it possible to understand indirectly metabolic evolution from the evolution of substrate transport of MCF. However, the evolutionary study of substrate transport of MCF does not mean that all the concrete structures of mitochondrial carriers (MCs) must first be gained.
Here we studied the alternation of MCF structure and potential correlated functions of MCF during metazoan evolution. The data analysis indicates that the types of substrates transported by MCF as a whole were maintained during metazoan evolution. However, the size of the substrates transported by members of MCs continuously diminished during the evolutionary process. We have found that the ratio of hydrophobic amino acids at specific helix-helix interfaces increases significantly during vertebrate evolution. Amino acid's spatial positioning and the calculating of packing values both indicate the increase in the number of hydrophobic amino acids would lead to a more "tight" structure of the TR domain, which is in agreement with the trend of diminishing size of substrates transported by MCs. In addition, there was a significant increase in the number of carriers of MCF during vertebrate evolution.
We propose that the more "tight" TR structure generated by the increase of the hydrophobic amino acids at specific helix-helix interfaces during vertebrate evolution enhances the substrate selectivity of MCF, reflecting the evolutionary trajectory of MCF during metazoan evolution.
Mitochondrial carriers (MCs), located at the inner membrane of mitochondria (IMM), are responsible for transporting a variety of metabolites between the cytosol and the mitochondria and influence the activities of the metabolic pathway direct [1, 2]. The substrates transported by MCs, such as nucleotides, amino acids, cofactors, and H+, are involved in oxidative phosphorylation, gluconeogenesis, the synthesis and degradations of amino and fatty acids, the macromolecular synthesis of proteins and nucleic acids, sterol metabolism, and thermogenesis . The mutations of the carriers are closely correlated with diseases, including Amish microcephaly, CAC deficiency, NICCD/CTLN2, neonatal myoclonic epilepsy, Sengers' syndrome, severe obesity, and type II diabetes [1, 3].
The substrates transported by MCs vary widely in their structure and size from the smallest, H+, to the largest and most highly charged species, such as ATP. By considering conservation of amino acids, distance and chemical constraints , and by modelling family members [5–12] on the known structure of the ADP/ATP carrier, recent research has identified a common substrate binding site located at approximately the midpoint of the membrane . However, the molecular mechanism of substrate selectivity of MCs remains incompletely understood. It is not possible to obtain the concrete structures of all mitochondrial carrier family (MCF) members through crystallization [1, 13], and thus uncovering the mechanism of substrate selectivity of MCF direct is impossible solely by linear superposition.
A powerful source of functional innovation is gene duplication and loss . Some reports have suggested that many of gene duplication events occur during the earlier period of vertebrate evolution [15, 16], indicating that if gene duplications exist in MCF [17, 18], the substrate selectivity changes may be a form of functional innovation closely associated with metabolic evolution. Here we investigate MCF, a gene family with common functional categories whose members do not have clearly defined functions, and subject it to The Big Experiment  to explore the evolutionary trajectory of MCF. We further investigate the synchronization between molecular evolution and macroevolution to search for evolutionary rules of sequences corresponding to structures of MCF in order to determine the characteristics of structural changes. The molecular mechanism of change of substrate selectivity was then analyzed to identify the transport functions of MCs.
Gene duplication of MCF is under strong purifying selection
According to the classification of the MCF secondary structure region (Additional file 1), TR has the most conservative sequence with an identity of 70.44%, while the average identities of LOOP135 (loop1, loop3, loop5) and NCLOOP24 (n, c, loop2, loop4) (Additional file 1) are only 42.24% and 19.25%, respectively. On the one hand, the high evolutionary conservation in TR reflects a stable structure maintaining the basic function of the transport channel. On the other hand, this evolutionary conservation indicates that the relatively small changes (29.56%) of TR sequences may have a large effect on the substrate transport of MCF.
Based on phylogenetic trees constructed from TR sequences, orthologous subfamilies were identified from gene families along with the relative gene duplication events. Fifteen gene duplication events were found in twelve gene families (Additional file 2). For the gene family tested, the alternative model (specifically testing for positive selection) was not a significantly better fit (2λ < 2.71, df = 1, p > 0.05) than the null model (Additional file 3). This result shows that there is no positive selection in the evolution of MCF. According to purifying selection (d N /d S < 1), the smaller the d N /d S ratio is, the greater the number of eliminated substitutions and the greater the selective constraint . Therefore, strong purifying selection is characteristic of a high number of eliminated substitutions and a strong selective constraint, indicating that its d N /d S values must be very small. The cumulative distribution of the d N /d S of MCF (Figure 1) and the mean value of d N /d S (0.023) both show that strong purifying selection occurs in the evolution of MCF in metazoan. The strong purifying selection indicates that MCF maintains a relatively stable substrate transport and its evolution is characteristic of continuity during metazoan evolution.
The evolution of MCF keeps pace with metazoan evolution
The distribution of gene duplication events was found mainly 750-430 Myr ago (Figure 2A), meeting approximately the earlier stage of vertebrate evolution along with a large scale of events of gene duplication [15, 16]. According to the trend of the cumulative distribution of age distribution (Figure 2B), an explosion of gene duplication events occurred approximately 600 Myr ago during the Ediacarian and Cambrian periods (approximately 630 - 490 Myr ago). This timing indicates that the functional evolution of MCF began at approximately the same time as metazoan evolution.
The correlation coefficient is 0.999 (P < 0.001) with well congruence by comparing the divergence time of species from the local clock estimation with the divergence time inferred from the fossil and molecular data (Figure 3B), indicating that the evolution of MCF keeps pace with metazoan evolution. The metabolic evolution of metazoans may be embodied in the functional evolution of MCF.
The significant rise of hydrophobic residues at helix-helix interfaces
The evolution of MCF was divided into three consecutive phases (see Methods) based on the fact that the evolution of MCF occurs in parallel to metazoan evolution. Statistics on the components of amino acids in different structural regions of MCF showed that the hydrophobic residues occupying a large proportion of TR had a significant rise from the invertebrate phase to the vertebrate phase (P < 0.001; Additional file 4A). There was a significant rise in the hydrophobic residues in TRDOWN (P < 0.001; Additional file 4B) too, while no significant changes in TRUP was observed. Thus, the region with a significant rise in hydrophobic residues was focused on TRDOWN135 (P < 0.001; Additional file 4C). There was no significant change in the proportion of charged residues and turn residues in TR. According to the method of deselecting species, three sequence pools were applied into the analysis of AACs. Non-ver was clearly separated from Ver-eco and Ver-endo according to the cumulative distribution of the proportion of the hydrophobic residues in TRDOWN135 (Figure 4A), while TRDOWN246 did not show a similar variation (Figure 4B). This result further demonstrated that the number of hydrophobic residues was significantly increased in TRDOWN135 during vertebrate evolution.
Interestingly, TRUP246 showed a similar distribution (Figure 4D) with TRDOWN135, indicating that the number of hydrophobic residues was increased in TRUP246 during vertebrate evolution too. Although TRDOWN135 and TRUP246 faced opposite directions (Additional file 1 and Figure 4E), the spatial distributions in the structure of ADP/ATP carrier clearly showed that they were close to each other and formed three helix-helix interfaces (Figure 4E). The conserved P and G (Figure 4F) in TR246 could serve as the hinges between TRUP246 and TRDOWN246 (Figure 4E) and made the helices in TRUP246 tilt to the TRDOWN135 sides (Figure 4E). TRDOWN246 participated in the formation of helix-helix interfaces with TRDOWN135 as well as TRUP246 (Additional file 5A). However, TRUP135 had no similar properties according to its space distribution (Additional file 5B). So, TRDOWN135 served as links between TRUP246 and TRDOWN246 and played a vital role in the formation of helix-helix interfaces, indicating the potential functional significance of the increased hydrophobic residues in TRDOWN135.
The "tight" structure of TR during vertebrate evolution
According to the orientation of the three dimensional structure of the ADP/ATP carrier, TRDOWN mainly consists of hydrophobic residues and some hydrophilic residues (Figure 5A and 5B). The changes of the type and number of hydrophobic residues in TRDOWN affected the connection among α-helices and further determined direct the degree of compactness of TRDOWN since hydrophobic residues in TRDOWN accounted for the connection among α-helices (Figure 5A and 5B). These hydrophobic residues were distributed in approximately the vertical planes along the central axis from the cytosol to the matrix (Figure 5B).
Furthermore, the logo  analysis based on all the sequences in the TRDOWN of MCF showed not only the enrichment of the hydrophobic residues at the same site or contiguous site in TRDOWN, but also the transformation of hydrophilic residues to hydrophobic residues in TRDOWN135 at sites 2, 3, 4, 6 and 7 during metazoan evolution (Figure 5C and 5E). The hydrophobic effect was lost once the helices were inserted into hydrophobic bilayers as a driving force for the helix association of integral membrane proteins  such as the mitochondrial inner membrane protein, MCF. The helix association then occurs through a combination of hydrogen-bonding, electrostatic, and van der Waals interactions . Due to the enrichment of the hydrophobic residues in position 3, 4, 6, 8 of TRDOWN246 (Figure 5D and 5F) and the enrichment of the hydrophobic residues in position 4, 5, 7, 8 of TRUP246 (Additional file 6), we postulated that the methyl groups of the side chain of the hydrophobic residues might interact with the methyl groups of the hydrophobic residues from the neighboring residues in TRDOWN135 by van der Waals interactions to maintain the helix-helix interactions. The experimental evidence showed that A134 of bovine mitochondrial oxoglutarate carrier potentially engaged in inter-helical interactions as well as other key residues in the odd-numbered transmembrane α-helices . Further, we found that A134 was located rightly at the TRDOWN135 and evolved from S at the same site of invertebrate (Additional file 7), which indirectly supported our hypothesis that the increase of hydrophobic residues would bring in more methyl groups, further strengthening the accumulated effect of the van der Waals interaction among them and generating a more "tight" TR structure.
According to the calculation of packing values of the residues in membrane proteins , TRDOWN135 in vertebrates had a significantly higher packing value (0.396) than that of invertebrates (0.383). This result thus further indicates that vertebrates have a more "tight" TR structure of MCF than invertebrates.
Enhancement of substrate selectivity
The salt network formed by the positive acids and negative acids (Additional file 8A) from three motif sequences (Px(D/E)xx(K/R)) constituting the conical pit of MCF (Additional file 8C) is important for transporting substrates [23–25]. Thus, the tandem concatenated sequences from three motifs were selected to serve as the indicators selecting the substrates. The Logo analysis on the mark sequences of MCF as a whole in metazoans showed a high similarity (Figure 6A), indicating that the types of substrates transported by MCF are not different. This is in agreement with the strong purifying selection (Additional file 3) occurring in the evolution of MCF in metazoan.
A phylogenetic tree was constructed using MEGA4.0  based on TR sequences from 44 carriers in SLC25 and their orthologous sequences in nine species (44 Χ 9) (Figure 6B). As can be seen, its branches were labelled according to previous reports related to substrates transported by MCs [1, 2]. The tree clearly showed three clusters that correspond to three categories of substrates, namely ADP/ATP, amino acid and cofactor, and H+. According to the E-value of each carrier obtained by BLAST in all the fungal genomes database, we calculated the average E-values of representative carriers for each cluster, showing that the ADP/ATP cluster was 10-50.56, the amino acid and cofactor cluster was 10-43, and the H+ cluster was 10-38.35. The results showed that the ADP/ATP cluster was an older cluster, while the H+ cluster was relatively younger cluster. Combined with the size of substrates transported by MCs , we found that the size of substrates transported by members of MCs continuously diminished throughout the evolution (Figure 6B).
Proline served as the hinge between TRDOWN135  and the label sequences of Px(D/E)xx(K/R) [24, 27] (Additional file 8). The conserved P and G (Figure 4F) in TR246 (Figure 4E and Figure 4F) could also serve as the hinges between TRUP246 and TRDOWN246. Furthermore, these conserved P and G in the transmembrane α-helices (Figure 4F) constituted the P-G level 1 located at approximately the midpoint of transmembrane region and the P-G level 2 located at approximately the bottom of transmembrane region . The two P-G levels  made the TRDOWN relatively independent and less affected by the structural changes due to helix tilt during the substrate transport process. So, the increased hydrophobic residues at the specific helix-helix interfaces could enhance the compactness degree in TRDOWN and strengthened the interaction between substrates and substrates contacting sites  of the MCF, thus allowing the passing of the substrates with the given size, in agreement with the trend of continuously diminishing size of substrates (Figure 6B).
MCF acts as the molecular fossil of metazoan
The rise of oxygen concentration appears to have precipitated a massive increase in metabolic complexity, culminating in the rise of metazoans around 600 million years ago . Feedback in the rise of oxygen concentrations has led to increasing metabolic complexity  and shows the existence of the intrinsic molecular mechanism of continuously adapting to the rise of oxygen concentration. We found that many gene duplication events of MCF occurred during the early period of vertebrate evolution, especially during the Ediacarian and Cambrian periods (Additional file 2). Further, the molecular evolution of MCF keeps pace with macroevolution (Figure 3). MCF located in mitochondria has a direct effect on metabolic activity [1, 2], and thus this feedback mechanism may be embodied in the functional evolution of MCF. These data make the MCF act as the molecular fossil for exploring metabolic complexity in metazoan.
Enhancement of substrates selectivity correlates to the rise of oxygen concentration
The enhancement of the substrate selectivity of MCF could partially explain the feedback mechanism evolved in metabolic evolution from the viewpoint of providing energy and reducing damage. From the view of providing energy, if mitochondria can be regarded as an energy-producing factory, its core department is the electron-transport chain. The enhancement of substrate selectivity of MCF allowed for more frequent communication between the mitochondrial matrix and the IMM. In addition, this increase of substrate selectivity provided more energy materials to the energy-producing factory by promoting the operation efficiency of this transport machine. Thus, the increase in the concentration of oxygen as the last electron receiver in the electron-transport chain [30, 31] would produce more abundant energy to satisfy physiological needs.
Oxygen participates in the provision of energy for the evolution of life but also brings by-products, such as superoxide and peroxides . The enhancement of substrate selectivity of MCF offers a possible way to reduce oxygen damage. It is the tendency towards strengthening of the substrate selectivity of MCF that makes the enhancement of selectivity of the H+ substrate. Uncoupling protein is generally considered as responsible for H+ conductance [32, 33]. No specific uncoupling protein is responsible for the proton conductance of yeast mitochondrion . The recent research represented that lamprey, generally considered as the most primitive vertebrate, had the uncoupling protein (UCP) transporting H+ , at least indicating the enhancement of selectivity of the H+ substrate. We also found the UCP subfamily was diverged approximately 637 million years ago (Additional file 2), in good agreement with around 600 million years ago when the rise of oxygen concentration culminates in the rise of metazoans. This good correlation further indicate that the rise of oxygen concentration may induces the enhancement of selectivity of the H+ substrate, which reduces oxygen damage by affecting the IMM potential, such as the role of the UCP located in the H+ cluster in reducing superoxide production and production of ROS [36–38].
In view of the important role of MCF as a molecular fossil in metabolism, we postulate that oxygen could be an environmental factor affecting the variation of substrate transport of MCF. Our studies explore the functional mechanism of oxygen on the evolution of MCF from the viewpoint of the enhancement of substrate selectivity of MCF.
The more "tight" TR structure generated by the increase of the hydrophobic amino acids at specific helix-helix interfaces enhances the substrate selectivity of MCF, which reflects the evolutionary trajectory of MCF. Because the goal of this study was to understand the variation in the structure of MCF from the level of sequence, our findings might provide the more specific mutation sites for laboratory experiments to improve experimental efficiency.
Partitioning evolutional processes and obtaining orthologous sequences
The orthologous sequences of human MCF (SLC25) were obtained from Swiss-Prot/TrEMBL and GenBank using the BLAST network service [39, 40]. The complexity of metabolism was reflected by the evolutional process of metazoan, and thus the orthologous sequences from metazoan were studied exclusively. In view of the facts that sequence similarity is close to 30% between members of SLC25 , the shorter process of metazoan evolution has a lower sequence variation and the orthologous gene itself has the a higher sequence similarity, the E-value of BLAST was set to 10-50 to exclude the interference from many nonhomologous genes. For cases when there was more than one sequence with a high similarity to human members, the sequence with the highest similarity was selected as the orthologous gene.
The emergences of vertebrates and endotherms indicate the two great leaps in metazoan evolution, and thus metazoan evolution was divided into three consecutive phases named the Non-ver (invertebrate) phase, the Ver-eco (vertebrate ectotherm) phase and the Ver-endo (vertebrate endotherm) phase. As current databases do not provide sequence data of MCF for more than three ectotherm species, we only chose three fully annotated species from three partitioned phases to meet the statistical requirements of ANOVA. Caenorhabditis elegans, Drosophila melanogaster and Strongylocentrotus purpuratus were used to represent the Non-ver phase. Tetraodon nigroviridis, Danio rerio and Xenopus tropicalis were used to represent the Ver-eco phase. Bos Taurus, Mus musculus, and Homo sapiens were used to represent the Ver-endo phase. The number of carriers from the nine species are as follows: human (44), mouse (44), bovine (40), frog (43), zebrafish (44), fugu (41), sea urchin (22), fly (26), and worm (25). A list of sequence accessions could be found in Additional file 9.
Classification of the whole structural region
In accordance with the generally accepted secondary structure of MCF [1, 2, 13, 23, 27, 41], we defined amino-ends, carboxyl-ends and two hydrophilic regions (n, c, loop2, loop4) facing the cytosol space as NCLOOP24; transmembrane regions consisting of six α-helices as TR; and three hydrophilic regions facing matrix space as LOOP135 (loop1, loop3, loop5) (Additional file 1). Six transmembrane α-helices were defined according to the classification of TR  and the c-terminal definition of three odd α-helices . The example of TR classification could be found in the CLUSTAL alignment among the oxoglutarate/malate sequences in metazoan species (Additional file 7). TR135 represents three odd α-helices. TR246 represents three even α-helices. TR was separated into two parts along the vertical direction from cytosol to matrix (Additional file 1). TRUP indicates the upper part facing cytosol. TRUP was further divided into two parts: TRUP135 was made up of three odd α-helices, while TRUP246 consisted of three even α-helices. TRDOWN represents the lower part of TR which faces the mitochondria matrix. TRDOWN was further divided into two parts: TRDOWN135 was made up of three odd α-helices, while TRDOWN246 consisted of three even α-helices.
Gene duplication and natural selection detection
Gene duplication events were detected through the construction of phylogenetic trees of gene families . The acceptance criterion was that the orthologous subfamily emerging from a duplication event had at least two species. The topology of the phylogenetic tree of the candidate family was constructed by the neighbor-joining method  (bootstrap = 500, JTT matrix).
The nonsynonymous substitution (those causing amino acid-alterations) to synonymous substitution (silent) rate ratio (ω = d N /d S ) provides a sensitive measure of natural selection at the protein level, with ω values of 1, > 1, and < 1, indicating neutral evolution, positive selection and purifying selection, respectively . We used branch-site model A [44, 45] to detect the type of natural selection of TR. This model has been reported to be more sensitive to the detection of positive selection than previous models including branch models  and site models , which may be useful in detecting positive selection after gene duplication in gene family evolution . We used a branch-site test of positive selection  to construct a likelihood ratio test between Model A with ω2≧1(alternative model,) and Model A1 with ω2 = 1 fixed (null model). The test was done by comparing the difference of likelihood values 2λ to a χ2 distribution of 1 degree of freedom. The critical values of 2λ at the 5% and 1% level are 2.71 and 5.41, respectively . Sites that could be under positive selection were identified with a Bayes Empirical Bayes approach . In addition, synonymous (d S ) and nonsynonymous (d N ) substitution rates were estimated using the methods of Yang and Nielson  as implemented in yn00 using PAML software .
Time scale and molecular clock estimation
The time scale estimation of gene duplication with multiple calibrations was implemented in the codeml program  of the PAML package (version 4) . We adopted several calibrations inferred from both fossil data and molecular data [52, 53], such as the primate-rodent (91 Myr ago), and vertebrate-fly (993 Myr ago) calibrations. The mtREV24 model  was applied as a model of amino acid substitution. We used both global clock and local clock methods to conduct the time scale estimation of gene duplication . Linear regression was applied to test the congruence between global clock and local clock.
We applied the above-mentioned method to the molecular clock estimation to infer whether the evolution of MCF keeps pace with macroevolution. Nine species could be found in the small number of orthologous sequences of human MCF, barely enough on behalf of MCF. Therefore, six species (worm, fruitfly, zebrafish, frog, mouse, and human) were selected to narrow the range of species. Then 33 orthologous sequences were found for each species. Further, the tandem concatenated sequence, which consisted of 33 orthologous TR sequences in each species, was applied in the molecular clock estimation.
Analysis of the feature of the sequence from the structure of MCF
The statistics of the amino acid component (AAC) were conducted according to the classification of the structure of MCF. Hydrophobic residues were based on the hydrophobicity table . Our compiled program using perl 5.8.7 was applied to the AAC statistic with two approaches of selecting and de-selecting species, respectively. Firstly, the statistics with selecting species were based on the above-mentioned nine species. A one-way ANOVA, followed by the Holm-Sidak test for multiple comparisons, was conducted using SigmaStat Windows Version 3.5 to evaluate differences in mean values among groups. The downloaded orthologous sequences according to 10-50 (not including sequence fragments and the sequences with obvious deletion in TR or conical pit region) were divided into the three sequences pools that were in the Non-ver pool (192), the Ver-eco pool (164), and the Ver-endo pool (208). A list of sequence accessions could be found in Additional file 10. The mean value of the percentage of the relative AAC was calculated from 132 (44 carriers × 3 species) sequences sampled randomly from three pools, respectively. This process was repeated for 132 times. The cumulative distribution function in the R language was applied to the analysis of the cumulative distribution of the mean values obtained.
Helix packing is important in the folding, stability, and association of membrane proteins [22, 57]. The packing value correlates positively with the degree of compactness of TR [22, 58]. Therefore, the quantitative statistics of packing values of the membrane protein can partially reflect the degree of compactness of TR. By studying the packing trait on the membrane protein helix from the obtained structure data of the membrane complex , it was able to obtain the packing values of the helix residues of TR. These data were applied to determine the average packing values of TR in MCF, indicating changes in the degree of compactness of TR during metazoan evolution.
The evolution of substrate transport of MCs
A phylogenetic tree was constructed using MEGA4.0  based on TR sequences from 44 carriers in SLC25 and their orthologous sequences in nine species (44 Χ 9). In view of the large number of sequences, the relatively low number of amino acids (the average length of the sequences used was 111 amino acids), and the relatively high divergence among the sequences, the neighbor-joining method was selected to generate the phylogenetic tree with the Poisson correction used to measure the pairwise sequence distances [59, 60]. Other tree-generating algorithms (minimal evolution, UPMGA and maximum likelihood) resulted in phylogenetic trees with a similar topology regarding the main branches generated by the neighbor-joining method. According to previous reports about the function of the members of MCF [1, 2], we could obtain the pattern of substrate transport of MCF.
As this family is exclusive to eukaryotes , we took all the fungal genomes in NCBI, the primitive species in eukaryotes, as a reference from which to estimate the ancient level of MCs. The sequences of SLC25 members were applied into the BLAST network service in all the fungal genomes in NCBI to obtain the corresponding list of E-values. The smallest 10 E-values on the list were used to determine the mean value statistics. The smaller the mean value, the higher the similarity between the sequences, and the older the carrier.
mitochondrial carrier family
solute carrier family 25
inner membrane of mitochondria
reactive oxygen species
amino acid component.
Palmieri F: The mitochondrial transporter family (SLC25): physiological and pathological implications. Pflugers Arch. 2004, 447: 689-709. 10.1007/s00424-003-1099-7.
Arco AD, Satrustegui J: New mitochondrial carriers: an overview. Cell Mol Life Sci. 2005, 62: 2204-2227. 10.1007/s00018-005-5197-x.
Palmieri F: Diseases caused by defects of mitochondrial carriers: a review. Biochim Biophys Acta. 2008, 1777: 564-578. 10.1016/j.bbabio.2008.03.008.
Robinson AJ, Kunji ER: Mitochondrial carriers in the cytoplasmic state have a common substrate binding site. Proc Natl Acad Sci USA. 2006, 103: 2617-2622. 10.1073/pnas.0509994103.
Walters DE, Kaplan RS: Homology-modeled structure of the yeast mitochondrial citrate transport protein. Biophys J. 2004, 87: 907-911. 10.1529/biophysj.104.042127.
Ma C, Kotaria R, Mayor JA, Remani S, Walters DE, Kaplan RS: The yeast mitochondrial citrate transport protein: characterization of transmembrane domain III residue involvement in substrate translocation. J Biol Chem. 2005, 280: 2331-2340. 10.1074/jbc.M411474200.
Morozzo Della Rocca B, Miniero DV, Tasco G, Dolce V, Falconi M, Ludovico A, Cappello AR, Sanchez P, Stipani I, Casadio R, et al: Substrate-induced conformational changes of the mitochondrial oxoglutarate carrier: a spectroscopic and molecular modelling study. Mol Membr Biol. 2005, 22: 443-452. 10.1080/09687860500269335.
Tonazzi A, Giangregorio N, Palmieri F, Indiveri C: Relationships of Cysteine and Lysine residues with the substrate binding site of the mitochondrial ornithine/citrulline carrier: an inhibition kinetic approach combined with the analysis of the homology structural model. Biochim Biophys Acta. 2005, 1718: 53-60. 10.1016/j.bbamem.2005.10.007.
Cappello AR, Curcio R, Valeria Miniero D, Stipani I, Robinson AJ, Kunji ER, Palmieri F: Functional and structural role of amino acid residues in the even-numbered transmembrane alpha-helices of the bovine mitochondrial oxoglutarate carrier. J Mol Biol. 2006, 363: 51-62. 10.1016/j.jmb.2006.08.041.
Cappello AR, Miniero DV, Curcio R, Ludovico A, Daddabbo L, Stipani I, Robinson AJ, Kunji ER, Palmieri F: Functional and structural role of amino acid residues in the odd-numbered transmembrane alpha-helices of the bovine mitochondrial oxoglutarate carrier. J Mol Biol. 2007, 369: 400-412. 10.1016/j.jmb.2007.03.048.
Marobbio CM, Giannuzzi G, Paradies E, Pierri CL, Palmieri F: alpha-Isopropylmalate, a leucine biosynthesis intermediate in yeast, is transported by the mitochondrial oxalacetate carrier. J Biol Chem. 2008, 283: 28445-28453. 10.1074/jbc.M804637200.
Tessa A, Fiermonte G, Dionisi-Vici C, Paradies E, Baumgartner MR, Chien YH, Loguercio C, de Baulny HO, Nassogne MC, Schiff M, et al: Identification of novel mutations in the SLC25A15 gene in hyperornithinemia-hyperammonemia-homocitrullinuria (HHH) syndrome: a clinical, molecular, and functional study. Hum Mutat. 2009, 30: 741-748. 10.1002/humu.20930.
Kunji ER: The role and structure of mitochondrial carriers. FEBS Lett. 2004, 564: 239-244. 10.1016/S0014-5793(04)00242-X.
Wapinski I, Pfeffer A, Friedman N, Regev A: Natural history and evolutionary principles of gene duplication in fungi. Nature. 2007, 449: 54-61. 10.1038/nature06107.
Gu X, Wang Y, Gu J: Age distribution of human gene families shows significant roles of both large- and small-scale duplications in vertebrate evolution. Nat Genet. 2002, 31: 205-209. 10.1038/ng902.
McLysaght A, Hokamp K, Wolfe KH: Extensive genomic duplication during early chordate evolution. Nat Genet. 2002, 31: 200-204. 10.1038/ng884.
Walker JE, Runswick MJ, Saraste M: Subunit equivalence in Escherichia coli and bovine heart mitochondrial F1F0 ATPases. FEBS Lett. 1982, 146: 393-396. 10.1016/0014-5793(82)80960-5.
Fiermonte G, Dolce V, Arrigoni R, Runswick MJ, Walker JE, Palmieri F: Organization and sequence of the gene for the human mitochondrial dicarboxylate carrier: evolution of the carrier family. Biochem J. 1999, 344 (Pt 3): 953-960. 10.1042/0264-6021:3440953.
Bartel DP: MicroRNAs: target recognition and regulatory functions. Cell. 2009, 136: 215-233. 10.1016/j.cell.2009.01.002.
Wagner A: Selection and gene duplication: a view from the genome. Genome Biol. 2002, 3: reviews1012.1011-1012.1013. 10.1186/gb-2002-3-5-reviews1012.
Crooks GE, Hon G, Chandonia JM, Brenner SE: WebLogo: a sequence logo generator. Genome Res. 2004, 14: 1188-1190. 10.1101/gr.849004.
Eilers M, Shekar SC, Shieh T, Smith SO, Fleming PJ: Internal packing of helical membrane proteins. Proc Natl Acad Sci USA. 2000, 97: 5796-5801. 10.1073/pnas.97.11.5796.
Pebay-Peyroula E, Dahout-Gonzalez C, Kahn R, Trezeguet V, Lauquin GJ, Brandolin G: Structure of mitochondrial ADP/ATP carrier in complex with carboxyatractyloside. Nature. 2003, 426: 39-44. 10.1038/nature02056.
Nelson DR, Felix CM, Swanson JM: Highly conserved charge-pair networks in the mitochondrial carrier family. J Mol Biol. 1998, 277: 285-308. 10.1006/jmbi.1997.1594.
Robinson AJ, Overy C, Kunji ER: The mechanism of transport by mitochondrial carriers based on analysis of symmetry. Proc Natl Acad Sci USA. 2008, 105: 17766-17771. 10.1073/pnas.0809580105.
Tamura K, Dudley J, Nei M, Kumar S: MEGA4: Molecular Evolutionary Genetics Analysis (MEGA) software version 4.0. Mol Biol Evol. 2007, 24: 1596-1599. 10.1093/molbev/msm092.
Walker JE, Runswick MJ: The mitochondrial transport protein superfamily. J Bioenerg Biomembr. 1993, 25: 435-446. 10.1007/BF01108401.
Palmieri F, Pierri CL: Structure and function of mitochondrial carriers - role of the transmembrane helix P and G residues in the gating and transport mechanism. FEBS Lett. 2010, 584: 1931-1939. 10.1016/j.febslet.2009.10.063.
Falkowski PG: Evolution. Tracing oxygen's imprint on earth's metabolic evolution. Science. 2006, 311: 1724-1725. 10.1126/science.1125937.
De Marais DJ: Evolution. When did photosynthesis emerge on Earth?. Science. 2000, 289: 1703-1705.
Dismukes GC, Klimov VV, Baranov SV, Kozlov YN, DasGupta J, Tyryshkin A: The origin of atmospheric oxygen on Earth: the innovation of oxygenic photosynthesis. Proc Natl Acad Sci USA. 2001, 98: 2170-2175. 10.1073/pnas.061514798.
Echtay KS, Winkler E, Frischmuth K, Klingenberg M: Uncoupling proteins 2 and 3 are highly active H(+) transporters and highly nucleotide sensitive when activated by coenzyme Q (ubiquinone). Proc Natl Acad Sci USA. 2001, 98: 1416-1421. 10.1073/pnas.98.4.1416.
Krauss S, Zhang CY, Lowell BB: The mitochondrial uncoupling-protein homologues. Nat Rev Mol Cell Biol. 2005, 6: 248-261. 10.1038/nrm1592.
Roussel D, Harding M, Runswick MJ, Walker JE, Brand MD: Does any yeast mitochondrial carrier have a native uncoupling protein function?. J Bioenerg Biomembr. 2002, 34: 165-176. 10.1023/A:1016027302232.
Wang C, Sun G, Chen K, Lv Z, Peng S, Jiang X, Xiang Y, Zhang C: Molecular cloning of lamprey uncoupling protein and assessment of its uncoupling activity using a yeast heterologous expression system. Mitochondrion. 2010, 10: 54-61. 10.1016/j.mito.2009.09.003.
Echtay KS, Roussel D, St-Pierre J, Jekabsons MB, Cadenas S, Stuart JA, Harper JA, Roebuck SJ, Morrison A, Pickering S, et al: Superoxide activates mitochondrial uncoupling proteins. Nature. 2002, 415: 96-99. 10.1038/415096a.
Brand MD, Affourtit C, Esteves TC, Green K, Lambert AJ, Miwa S, Pakay JL, Parker N: Mitochondrial superoxide: production, biological effects, and activation of uncoupling proteins. Free Radic Biol Med. 2004, 37: 755-767. 10.1016/j.freeradbiomed.2004.05.034.
Vidal-Puig AJ, Grujic D, Zhang CY, Hagen T, Boss O, Ido Y, Szczepanik A, Wade J, Mootha V, Cortright R, et al: Energy metabolism in uncoupling protein 3 gene knockout mice. J Biol Chem. 2000, 275: 16258-16266. 10.1074/jbc.M910179199.
Altschul SF, Madden TL, Schaffer AA, Zhang J, Zhang Z, Miller W, Lipman DJ: Gapped BLAST and PSI-BLAST: a new generation of protein database search programs. Nucleic Acids Res. 1997, 25: 3389-3402. 10.1093/nar/25.17.3389.
Schaffer AA, Aravind L, Madden TL, Shavirin S, Spouge JL, Wolf YI, Koonin EV, Altschul SF: Improving the accuracy of PSI-BLAST protein database searches with composition-based statistics and other refinements. Nucleic Acids Res. 2001, 29: 2994-3005. 10.1093/nar/29.14.2994.
Klingenberg M: Mechanism and evolution of the uncoupling protein of brown adipose tissue. Trends Biochem Sci. 1990, 15: 108-112. 10.1016/0968-0004(90)90194-G.
Trenker M, Malli R, Fertschai I, Levak-Frank S, Graier WF: Uncoupling proteins 2 and 3 are fundamental for mitochondrial Ca2+ uniport. Nat Cell Biol. 2007, 9: 445-452. 10.1038/ncb1556.
Saitou N, Nei M: The neighbor-joining method: a new method for reconstructing phylogenetic trees. Mol Biol Evol. 1987, 4: 406-425.
Yang Z, Nielsen R: Codon-substitution models for detecting molecular adaptation at individual sites along specific lineages. Mol Biol Evol. 2002, 19: 908-917.
Zhang J, Nielsen R, Yang Z: Evaluation of an improved branch-site likelihood method for detecting positive selection at the molecular level. Mol Biol Evol. 2005, 22: 2472-2479. 10.1093/molbev/msi237.
Yang Z: Likelihood ratio tests for detecting positive selection and application to primate lysozyme evolution. Mol Biol Evol. 1998, 15: 568-573.
Yang Z, Nielsen R, Goldman N, Pedersen AM: Codon-substitution models for heterogeneous selection pressure at amino acid sites. Genetics. 2000, 155: 431-449.
Yang Z: PAML 4: phylogenetic analysis by maximum likelihood. Mol Biol Evol. 2007, 24: 1586-1591. 10.1093/molbev/msm088.
Yang Z, Wong WS, Nielsen R: Bayes empirical bayes inference of amino acid sites under positive selection. Mol Biol Evol. 2005, 22: 1107-1118. 10.1093/molbev/msi097.
Yang Z, Nielsen R: Estimating synonymous and nonsynonymous substitution rates under realistic evolutionary models. Mol Biol Evol. 2000, 17: 32-43.
Goldman N, Yang Z: A codon-based model of nucleotide substitution for protein-coding DNA sequences. Mol Biol Evol. 1994, 11: 725-736.
Hedges SB: The origin and evolution of model organisms. Nat Rev Genet. 2002, 3: 838-849. 10.1038/nrg929.
Blair Hedges S, Kumar S: Genomic clocks and evolutionary timescales. Trends Genet. 2003, 19: 200-206. 10.1016/S0168-9525(03)00053-2.
Adachi J, Hasegawa M: Model of amino acid substitution in proteins encoded by mitochondrial DNA. J Mol Evol. 1996, 42: 459-468. 10.1007/BF02498640.
Yang Z, Yoder AD: Comparison of likelihood and Bayesian methods for estimating divergence times using multiple gene Loci and calibration points, with application to a radiation of cute-looking mouse lemur species. Syst Biol. 2003, 52: 705-716. 10.1080/10635150390235557.
Kyte J, Doolittle RF: A simple method for displaying the hydropathic character of a protein. J Mol Biol. 1982, 157: 105-132. 10.1016/0022-2836(82)90515-0.
Fleming KG, Ackerman AL, Engelman DM: The effect of point mutations on the free energy of transmembrane alpha-helix dimerization. J Mol Biol. 1997, 272: 266-275. 10.1006/jmbi.1997.1236.
Liu W, Eilers M, Patel AB, Smith SO: Helix packing moments reveal diversity and conservation in membrane protein structure. J Mol Biol. 2004, 337: 713-729. 10.1016/j.jmb.2004.02.001.
Nei M, Kumar S, Takahashi K: The optimization principle in phylogenetic analysis tends to give incorrect topologies when the number of nucleotides or amino acids used is small. Proc Natl Acad Sci USA. 1998, 95: 12390-12397. 10.1073/pnas.95.21.12390.
Takahashi K, Nei M: Efficiencies of fast algorithms of phylogenetic inference under the criteria of maximum parsimony, minimum evolution, and maximum likelihood when a large number of sequences are used. Mol Biol Evol. 2000, 17: 1251-1258.
Combet C, Jambon M, Deleage G, Geourjon C: Geno3D: automatic comparative molecular modelling of protein. Bioinformatics. 2002, 18: 213-214. 10.1093/bioinformatics/18.1.213.
Humphrey W, Dalke A, Schulten K: VMD: visual molecular dynamics. J Mol Graph. 1996, 14 (33-38): 27-38.
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-4680. 10.1093/nar/22.22.4673.
This work was supported by grants from the National Natural Science Foundation of China (no. 90813035, 30890044). We are very grateful to the three anonymous reviewers for valuable comments and suggestions.
MG and JL performed the analyses and drafted the manuscript. MW helped to draft the manuscript. MG, JW and KZ revised the manuscript. JW and KZ supervised the study. CZ initiated and supervised the project. All authors approved the final manuscript.
Ming Gong, Jie Li contributed equally to this work.
Electronic supplementary material
Additional file 1: Two-dimensional structure of MCF. The structure of MCF consists of the three segments that are NCLOOP24 facing the space of cytosol, TR located in the inner membrane, and LOOP135 facing the space of matrix, respectively. TR is made up of six α-helices represented by six bars. Dark dot represents individual amino acid. (TIFF 3 MB)
Additional file 2: Gene duplication events and the divergence time of MCs. Orthologous TR sequences from 44 members of SLC25 were applied into the construction of the phylogenetic trees to identify gene duplications events. 15 gene duplication events were detected in 12 independent trees. Arrow indicates a gene duplication event with the divergence time (Myr ago). The index in SLC25 (the appellation of MCF in human) was used to denote the corresponding orthologous sequences, such as A7, 8, 9 (representing the orthologous sequences of UCP1, 2, 3 in human, respectively). A47 represents the subfamily of Hepatocellular carcinoma-down-regulated mitochondrial carrier protein in human. (DOC 82 KB)
Additional file 3: Likelihood scores, model parameters, and likelihood-ratio test results for branch-site model analyses. A indicates the orthologous sequences of member of SLC25. The index in SLC25 (the appellation of MCF in human) was used to denote the corresponding orthologous sequences, such as A7, 8, 9 (representing the orthologous sequences of UCP1, 2, 3, respectively). We used a branch-site test of positive selection  to construct a likelihood ratio test between Model A with ω2≧1(alternative model,) and Model A1 with ω2 = 1 fixed (null model). The test was done by comparing the difference of likelihood values 2λ to a χ2 distribution of 1 degree of freedom. The critical values of 2λ at the 5% and 1% level are 2.71 and 5.41, respectively . (DOC 45 KB)
Additional file 4: Analysis of the AAC (amino acid component) in MCF. (A) The mean percentage of the hydrophobic residues in TR. (B) The comparison of the mean percentage of the hydrophobic residues in TRUP and TRDOWN. (C) The comparison of the mean percentage of the hydrophobic residues in TRDOWN135 and TRDOWN246. The groups were evaluated using a one way ANOVA followed by the Holm-Sidak test for multiple comparisons. n = 132. Bars represent the mean ± s.e.m, while an asterisk indicates P < 0.001 in comparison to the Non-ver phase. (TIFF 2 MB)
Additional file 5: The spatial distribution of transmembrane α-helices in the human ADP/ATP carrier. (A) The human ADP/ATP carrier viewed from the lateral side. Yellow represents TRDOWN135. Blue represents TR246. (B) The human ADP/ATP carrier viewed from the cytosol. Yellow represents TRUP135.(TIFF 4 MB)
Additional file 6: The logo chart of sequences in TR UP246 . (A) The logo chart of sequences in TRUP246 from Ver. (B) The logo chart of sequences in TRUP246 from Non-ver. Red indicates hydrophobic residues, blue indicates hydrophilic residues, and black indicates other residues. The logo chart describes the residue changes in the same vertical plane according to the space structure orientation of TRUP246 for metazoan evolution. WebLogo  was applied into the Logo analysis. Ver indicates vertebrates, while Non-ver indicates invertebrates. (TIFF 2 MB)
Additional file 7: The CLUSTAL alignment among the oxoglutarate/malate sequences in metazoan species. The sequences with fragment and obvious deletion in TR and conical pit region (CPR) were not included. CLUSTAL W  was applied to multiple sequence alignment. The marks "*", ".", and ":" stand for the identical amino acids, relative, and similar amino acids. CPR consists of three motif sequences (Px(D/E)xx(K/R)) that are located at the C-terminal of TR1, TR3 and TR5, respectively. The sequence accessions according to the species order above were as followed: Q02978, Q9CR62, P97700, NP_777096.1, NP_001090497.1, AAH71521.1, CAF90256.1, XP_002610854.1, XP_001639936.1, XP_001867726.1, NP_651703.1, NP_493694.2, XP_002571870.1 and XP_001893008.1. The detailed classifications of TR sequences and CPR in the oxoglutarate/malate carrier are applicable to other MCs too. (TIFF 11 MB)
Additional file 8: The structures of the CPR in the ADP/ATP carrier. (A) The ADP/ATP carrier in humans viewed from the lateral side. The surf region represents the CPR. Red indicates negative acids, blue indicates positive acids, white indicates non-polar residues, and green indicates polar residues. (B) The ADP/ATP carrier in human viewed from the cytosol. (C) TRDOWN and CPR viewed from the cytosol. Proline serves as the hinge between TRDOWN and CPR. (D) TRDOWN and CPR viewed from the matrix. (TIFF 6 MB)
Additional file 9: The accessions list from the selecting species. The details about selecting species are available in the method section. (TXT 6 KB)
Additional file 10: The accessions list from the deselecting species. The details about deselecting species are available in the method section. (TXT 12 KB)
Authors’ original submitted files for images
About this article
Cite this article
Gong, M., Li, J., Wang, M. et al. The evolutionary trajectory of mitochondrial carrier family during metazoan evolution. BMC Evol Biol 10, 282 (2010). https://doi.org/10.1186/1471-2148-10-282
- Hydrophobic Residue
- Gene Duplication Event
- Vertebrate Evolution
- Orthologous Sequence
- Substrate Transport