Evolution of proton pumping ATPases: Rooting the tree of life

Proton pumping ATPases are found in all groups of present day organisms. The F-ATPases of eubacteria, mitochondria and chloroplasts also function as ATP synthases, i.e., they catalyze the final step that transforms the energy available from reduction/oxidation reactions (e.g., in photosynthesis) into ATP, the usual energy currency of modern cells. The primary structure of these ATPases/ATP synthases was found to be much more conserved between different groups of bacteria than other parts of the photosynthetic machinery, e.g., reaction center proteins and redox carrier complexes.These F-ATPases and the vacuolar type ATPase, which is found on many of the endomembranes of eukaryotic cells, were shown to be homologous to each other; i.e., these two groups of ATPases evolved from the same enzyme present in the common ancestor. (The term eubacteria is used here to denote the phylogenetic group containing all bacteria except the archaebacteria.) Sequences obtained for the plasmamembrane ATPase of various archaebacteria revealed that this ATPase is much more similar to the eukaryotic than to the eubacterial counterpart. The eukaryotic cell of higher organisms evolved from a symbiosis between eubacteria (that evolved into mitochondria and chloroplasts) and a host organism. Using the vacuolar type ATPase as a molecular marker for the cytoplasmic component of the eukaryotic cell reveals that this host organism was a close relative of the archaebacteria.A unique feature of the evolution of the ATPases is the presence of a non-catalytic subunit that is paralogous to the catalytic subunit, i.e., the two types of subunits evolved from a common ancestral gene. Since the gene duplication that gave rise to these two types of subunits had already occurred in the last common ancestor of all living organisms, this non-catalytic subunit can be used to root the tree of life by means of an outgroup; that is, the location of the last common ancestor of the major domains of living organisms (archaebacteria, eubacteria and eukaryotes) can be located in the tree of life without assuming constant or equal rates of change in the different branches.A correlation between structure and function of ATPases has been established for present day organisms. Implications resulting from this correlation for biochemical pathways, especially photosynthesis, that were operative in the last common ancestor and preceding life forms are discussed.     

The Intriguing Evolution of the “b” and “G” Subunits in F-type and V-type ATPases: Isolation of the vma-10 Gene from Neurospora crassa

We have characterized the vma-10 gene which encodes the G subunit of the vacuolar ATPase in Neurospora crassa. The gene is somewhat unusual in filamentous fungi because it contains five introns, comprising 71% of the region between the translation start and stop codons. The 5 untranslated region of the gene contains several elements that have been identified in other genes that encode subunits of the vacuolar ATPase in N. crassa. A comparison of G subunits from N. crassa, S. cerevisiae, and animal cells showed that the N-terminal half of the polypeptide shows the highest degree of sequence conservation. Most striking is the observation that this region could form an alpha helix in which all of the conserved residues are clustered on one face. Subunit G appears to be homologous to the b subunit found in F-type ATPases. The major difference between the b and G subunits is the lack of a membrane-spanning region in the G subunit. We have also identified homologous subunits in the operons which encode V-type ATPases in a eubacterium, Enterrococcus hirae, and an archaebacterium, Methanococcus jannaschii. As in eukaryotic vacuolar ATPases the G subunit homologs lack a membrane-spanning region. Although the b and G subunits appear to be derived from a common ancestor, significant changes have evolved. In F-type and V-type ATPases these subunits can have zero, one, or two membrane-spanning regions and can also differ significantly in the number of copies per enzyme.     

Molecular Evolution of Nitrogen Fixation: The Evolutionary History of the nifD, nifK, nifE, and nifN Genes

The pairs of nitrogen fixation genes nifDK and nifEN encode for the α and β subunits of nitrogenase and for the two subunits of the NifNE protein complex, involved in the biosynthesis of the FeMo cofactor, respectively. Comparative analysis of the amino acid sequences of the four NifD, NifK, NifE, and NifN in several archaeal and bacterial diazotrophs showed extensive sequence similarity between them, suggesting that their encoding genes constitute a novel paralogous gene family. We propose a two-step model to reconstruct the possible evolutionary history of the four genes. Accordingly, an ancestor gene gave rise, by an in-tandem paralogous duplication event followed by divergence, to an ancestral bicistronic operon; the latter, in turn, underwent a paralogous operon duplication event followed by evolutionary divergence leading to the ancestors of the present-day nifDK and nifEN operons. Both these paralogous duplication events very likely predated the appearance of the last universal common ancestor. The possible role of the ancestral gene and operon in nitrogen fixation is also discussed. Received: 21 June 1999 / Accepted: 1 March 2000     

Gene duplication,genome duplication,and the functional diversification of vertebrate globins

The functional diversification of the vertebrate globin gene superfamily provides an especially vivid illustration of the role of gene duplication and whole-genome duplication in promoting evolutionary innovation. For example, key globin proteins that evolved specialized functions in various aspects of oxidative metabolism and oxygen signaling pathways (hemoglobin [Hb], myoglobin [Mb], and cytoglobin [Cygb]) trace their origins to two whole-genome duplication events in the stem lineage of vertebrates. The retention of the proto-Hb and Mb genes in the ancestor of jawed vertebrates permitted a physiological division of labor between the oxygen-carrier function of Hb and the oxygen-storage function of Mb. In the Hb gene lineage, a subsequent tandem gene duplication gave rise to the proto α- and β-globin genes, which permitted the formation of multimeric Hbs composed of unlike subunits (α₂β₂). The evolution of this heteromeric quaternary structure was central to the emergence of Hb as a specialized oxygen-transport protein because it provided a mechanism for cooperative oxygen-binding and allosteric regulatory control. Subsequent rounds of duplication and divergence have produced diverse repertoires of α- and β-like globin genes that are ontogenetically regulated such that functionally distinct Hb isoforms are expressed during different stages of prenatal development and postnatal life. In the ancestor of jawless fishes, the proto Mb and Hb genes appear to have been secondarily lost, and the Cygb homolog evolved a specialized respiratory function in blood-oxygen transport. Phylogenetic and comparative genomic analyses of the vertebrate globin gene superfamily have revealed numerous instances in which paralogous globins have convergently evolved similar expression patterns and/or similar functional specializations in different organismal lineages.     

Conserved intron positions in FGFR genes reflect the modular structure of FGFR and reveal stepwise addition of domains to an already complex ancestral FGFR

We have analyzed the evolution of fibroblast growth factor receptor (FGFR) tyrosine kinase genes throughout a wide range of animal phyla. No evidence for an FGFR gene was found in Porifera, but we tentatively identified an FGFR gene in the placozoan Trichoplax adhaerens. The gene encodes a protein with three immunoglobulin-like domains, a single-pass transmembrane, and a split tyrosine kinase domain. By superimposing intron positions of 20 FGFR genes from Placozoa, Cnidaria, Protostomia, and Deuterostomia over the respective protein domain structure, we identified ten ancestral introns and three conserved intron groups. Our analysis shows (1) that the position of ancestral introns correlates to the modular structure of FGFRs, (2) that the acidic domain very likely evolved in the last common ancestor of triploblasts, (3) that splicing of IgIII was enabled by a triploblast-specific insertion, and (4) that IgI is subject to substantial loss or duplication particularly in quickly evolving genomes. Moreover, intron positions in the catalytic domain of FGFRs map to the borders of protein subdomains highly conserved in other serine/threonine kinases. Nevertheless, these introns were introduced in metazoan receptor tyrosine kinases exclusively. Our data support the view that protein evolution dating back to the Cambrian explosion took place in such a short time window that only subtle changes in the domain structure are detectable in extant representatives of animal phyla. We propose that the first multidomain FGFR originated in the last common ancestor of Placozoa, Cnidaria, and Bilateria. Additional domains were introduced mainly in the ancestor of triploblasts and in the Ecdysozoa.     

Structures of P-type transporting ATPases and chromosomal locations of their genes

P-type ATPases (E1E2-ATPases) are primary active transporters which form phospho-intermediates during their catalytic cycle. They are classified into P1 to P4 based on the primary structure and potential transmembrane segments. Although the classic P-type ATPases are cation transporters, two new members have recently been found; one is a flippase catalyzing the flip-flop movement of aminophospholipids, but the substrate and function of the other one remain unknown. It would be interesting to determine whether the cations and aminophospholipids are transported by similar or different mechanisms. P-type ATPases are believed to have been derived from a common ancestor, and their genes are found to be distributed in various chromosomal loci. However, gene duplication events can be traced from the tandem arrangement of genes and their linkage map. Na+/K+- and H+/K+-ATPases have not only closely related a subunits but also similar beta subunits. Renal Na+/K+-ATPase has an additional subunit gamma. Similar small polypeptides (phospholemman, Mat-8 and CHIF), which induce Cl- and K+ currents, have been found. The idea of their functional and structural coupling with P-type ATPases, especially with H+/K+-ATPase, is intriguing. Each P-type ATPase must have specific domains or sequences for its intracellular trafficking (sorting, retention and recycling). Identification of such regions and studies on the molecules playing role in their recognition may facilitate the unveiling of various cellular processes regulated by P-type ATPases.     

The Prokaryote-to-Eukaryote Transition Reflected in the Evolution of the V/F/A-ATPase Catalytic and Proteolipid Subunits

Changes in the primary and quarternary structure of vacuolar and archaeal type ATPases that accompany the prokaryote-to-eukaryote transition are analyzed. The gene encoding the vacuolar-type proteolipid of the V-ATPase from Giardia lamblia is reported. Giardia has a typical vacuolar ATPase as observed from the common motifs shared between its proteolipid subunit and other eukaryotic vacuolar ATPases, suggesting that the former enzyme works as a hydrolase in this primitive eukaryote. The phylogenetic analyses of the V-ATPase catalytic subunit and the front and back halves of the proteolipid subunit placed Giardia as the deepest branch within the eukaryotes. Our phylogenetic analysis indicated that at least two independent duplication and fusion events gave rise to the larger proteolipid type found in eukaryotes and in Methanococcus. The spatial distribution of the conserved residues among the vacuolar-type proteolipids suggest a zipper-type interaction among the transmembrane helices and surrounding subunits of the V-ATPase complex. Important residues involved in the function of the F-ATP synthase proteolipid have been replaced during evolution in the V-proteolipid, but in some cases retained in the archaeal A-ATPase. Their possible implication in the evolution of V/F/A-ATPases is discussed. Received: 27 August 1997 / Accepted: 14 January 1998     

Ancestral paralogs and pseudoparalogs and their role in the emergence of the eukaryotic cell

Gene duplication is a crucial mechanism of evolutionary innovation. A substantial fraction of eukaryotic genomes consists of paralogous gene families. We assess the extent of ancestral paralogy, which dates back to the last common ancestor of all eukaryotes, and examine the origins of the ancestral paralogs and their potential roles in the emergence of the eukaryotic cell complexity. A parsimonious reconstruction of ancestral gene repertoires shows that 4137 orthologous gene sets in the last eukaryotic common ancestor (LECA) map back to 2150 orthologous sets in the hypothetical first eukaryotic common ancestor (FECA) [paralogy quotient (PQ) of 1.92]. Analogous reconstructions show significantly lower levels of paralogy in prokaryotes, 1.19 for archaea and 1.25 for bacteria. The only functional class of eukaryotic proteins with a significant excess of paralogous clusters over the mean includes molecular chaperones and proteins with related functions. Almost all genes in this category underwent multiple duplications during early eukaryotic evolution. In structural terms, the most prominent sets of paralogs are superstructure-forming proteins with repetitive domains, such as WD-40 and TPR. In addition to the true ancestral paralogs which evolved via duplication at the onset of eukaryotic evolution, numerous pseudoparalogs were detected, i.e. homologous genes that apparently were acquired by early eukaryotes via different routes, including horizontal gene transfer (HGT) from diverse bacteria. The results of this study demonstrate a major increase in the level of gene paralogy as a hallmark of the early evolution of eukaryotes.     

Defective assembly of a hybrid vacuolar H(+)-ATPase containing the mouse testis-specific E1 isoform and yeast subunits

Mammalian vacuolar-type proton pumping ATPases (V-ATPases) are diverse multi-subunit proton pumps. They are formed from membrane V(o) and catalytic V(1) sectors, whose subunits have cell-specific or ubiquitous isoforms. Biochemical study of a unique V-ATPase is difficult because ones with different isoforms are present in the same cell. However, the properties of mouse isoforms can be studied using hybrid V-ATPases formed from the isoforms and other yeast subunits. As shown previously, mouse subunit E isoform E1 (testis-specific) or E2 (ubiquitous) can form active V-ATPases with other subunits of yeast, but E1/yeast hybrid V-ATPase is defective in proton transport at 37 degrees C (Sun-Wada, G.-H., Imai-Senga, Y., Yamamoto, A., Murata, Y., Hirata, T., Wada, Y., and Futai, M., 2002, J. Biol. Chem. 277, 18098-18105). In this study, we have analyzed the properties of E1/yeast hybrid V-ATPase to understand the role of the E subunit. The proton transport by the defective hybrid ATPase was reversibly recovered when incubation temperature of vacuoles or cells was shifted to 30 degrees C. Corresponding to the reversible defect of the hybrid V-ATPase, the V(o) subunit a epitope was exposed to the corresponding antibody at 37 degrees C, but became inaccessible at 30 degrees C. However, the V(1) sector was still associated with V(o) at 37 degrees C, as shown immunochemically. The control yeast V-ATPase was active at 37 degrees C, and its epitope was not accessible to the antibody. Glucose depletion, known to dissociate V(1) from V(o) in yeast, had only a slight effect on the hybrid at acidic pH. The domain between Lys26 and Val83 of E1, which contains eight residues not conserved between E1 and E2, was responsible for the unique properties of the hybrid. These results suggest that subunit E, especially its amino-terminal domain, plays a pertinent role in the assembly of V-ATPase subunits in vacuolar membranes.     

Characterization of a membrane-associated ATPase from Methanococcus voltae, a methanogenic member of the Archaea.

A membrane-associated ATPase with an M(r) of approximately 510,000 and containing subunits with M(r)s of 80,000 (alpha), 55,000 (beta), and 25,000 (gamma) was isolated from the methanogen Methanococcus voltae. Enzymatic activity was not affected by vanadate or azide, inhibitors of P- and F1-ATPase, respectively, but was inhibited by nitrate and bafilomycin A1, inhibitors of V1-type ATPases. Since dicyclohexylcarbodiimide inhibited the enzyme when it was present in membranes but not after the ATPase was solubilized, we suggest the presence of membrane-associated component analogous to the F0 and V0 components of both F-type and V-type ATPases. N-terminal amino acid sequence analysis of the alpha subunit showed a higher similarity to ATPases of the V-type family than to those of the F-type family.     

Gene duplication and protein evolution. Case of aminoacyl-t-RNA synthetases]

Aminoacyl-tRNA synthetases are a functionally homologous group of enzymes which catalyse the first step of protein synthesis. The accumulation of data on the oligomeric structures of these enzymes has revealed a wide diversity in the sizes and organization of the protomers. This is not consistent with the idea of a family with homologous primary and tertiary structures. Howewer recent studies have shown that all these enzymes may have evolved from a common ancestor through gene duplication and fusion which has led to extensive repeating sequences in the corresponding subunits.     

Evolution of the soluble diiron monooxygenases

Based on structural, biochemical, and genetic data, the soluble diiron monooxygenases can be divided into four groups: the soluble methane monooxygenases, the Amo alkene monooxygenase of Rhodococcus corallinus B-276, the phenol hydroxylases, and the four-component alkene/aromatic monooxygenases. The limited phylogenetic distribution of these enzymes among bacteria, together with available genetic evidence, indicates that they have been spread largely through horizontal gene transfer. Phylogenetic analyses reveal that the alpha- and beta-oxygenase subunits are paralogous proteins and were derived from an ancient gene duplication of a carboxylate-bridged diiron protein, with subsequent divergence yielding a catalytic alpha-oxygenase subunit and a structural beta-oxygenase subunit. The oxidoreductase and ferredoxin components of these enzymes are likely to have been acquired by horizontal transfer from ancestors common to unrelated diiron and Rieske center oxygenases and other enzymes. The cumulative results of phylogenetic reconstructions suggest that the alkene/aromatic monooxygenases diverged first from the last common ancestor for these enzymes, followed by the phenol hydroxylases, Amo alkene monooxygenase, and methane monooxygenases.     

Evolution of catalytic proteins

Summary It is believed that all present-day organisms descended from a common cellular ancestor. Such a cell must have evolved from more primitive and simpler precursors, but neither their organization nor the route such evolution took are accessible to the molecular techniques available today. We propose a mechanism, based on functional properties of enzymes and the kinetics of growth, which allows us to reconstruct the general course of early enzyme evolution. A precursor cell containing very few multifunctional enzymes with low catalytic activities is shown to lead inevitably to descendants with a large number of differentiated monofunctional enzymes with high turnover numbers. Mutation and natural selection for faster growth are shown to be the only conditions necessary for such a change to have occurred.     

Isolation of genes encoding the Neurospora vacuolar ATPase. Analysis of vma-2 encoding the 57-kDa polypeptide and comparison to vma-1

In partially purified preparations of the vacuolar ATPase from Neurospora crassa, the two most prominent components are polypeptides of Mr = 70,000 and 60,000. We previously reported the isolation of the gene vma-1, which encodes the Mr = 70,000 polypeptide, and presented evidence that the polypeptide contains the site of ATP hydrolysis (Bowman, E. J., Tenney, K., and Bowman, B. J. (1988) J. Biol. Chem. 263, 13994-14001). We now report the isolation of a gene (designated vma-2), that encodes the Mr = 60,000 polypeptide. Analysis of the DNA sequence shows that the polypeptide has 513 amino acids and a molecular mass of 56,808 daltons (and will thus be referred to as the 57-kDa polypeptide). It is fairly rich in polar amino acids and has no apparent membrane-spanning domains. The vma-2 gene contains five short introns (55-71 bases), all clustered in the 5' end of the coding region. The gene maps to the right arm of linkage group II, near 5 S RNA gene 3. Thus, it is unlinked to vma-1 and to other known ATPase genes in N. crassa. The 57-kDa polypeptide shows 25% amino acid sequence identity with the vma-1 gene product. It shows essentially the same degree of similarity (25-28%) to both the alpha and beta subunits of F0F1 ATPases. Analysis of specific regions of the 57-kDa polypeptide, however, suggests it may have a function like that of the alpha subunit in F0F1 ATPases. The data indicate that all four types of ATPase polypeptides have evolved from a common ancestor and that the vacuolar-type ATPases have a structure surprisingly similar to that of the F0F1 ATPases.     

Malate decarboxylases: evolution and roles of NAD(P)-ME isoforms in species performing C(4) and C(3) photosynthesis

In the C(4) pathway of photosynthesis two types of malate decarboxylases release CO(2) in bundle sheath cells, NADP- and NAD-dependent malic enzyme (NADP-ME and NAD-ME), located in the chloroplasts and the mitochondria of these cells, respectively. The C(4) decarboxylases involved in C(4) photosynthesis did not evolve de novo; they were recruited from existing housekeeping isoforms. NADP-ME housekeeping isoforms would function in the control of malate levels during hypoxia, pathogen defence responses, and microspore separation, while NAD-ME participates in the respiration of malate in the tricarboxylic acid cycle. Recently, the existence of three enzymatic NAD-ME entities in Arabidopsis, occurring by alternative association of two subunits, was described as a novel mechanism to regulate NAD-ME activity under changing metabolic environments. The C(4) NADP-ME is thought to have evolved from a C(3) chloroplastic ancestor, which in turn would have evolved from an ancient cytosolic enzyme. In this way, the C(4) NADP-ME would have emerged through gene duplication, acquisition of a new promoter, and neo-functionalization. In contrast, there would exist a unique NAD-ME in C(4) plants, which would have been adapted to perform a dual function through changes in the kinetic and regulatory properties of the C(3) ancestors. In addition to this, for the evolution of C(4) NAD-ME, insertion of promoters or enhancers into the single-copy genes of the C(3) ancestors would have changed the expression without gene duplication.     

Birth and rapid subcellular adaptation of a hominoid-specific CDC14 protein

Gene duplication was prevalent during hominoid evolution, yet little is known about the functional fate of new ape gene copies. We characterized the CDC14B cell cycle gene and the functional evolution of its hominoid-specific daughter gene, CDC14Bretro. We found that CDC14B encodes four different splice isoforms that show different subcellular localizations (nucleus or microtubule-associated) and functional properties. A microtubular CDC14B variant spawned CDC14Bretro through retroposition in the hominoid ancestor 18–25 million years ago (Mya). CDC14Bretro evolved brain-/testis-specific expression after the duplication event and experienced a short period of intense positive selection in the African ape ancestor 7–12 Mya. Using resurrected ancestral protein variants, we demonstrate that by virtue of amino acid substitutions in distinct protein regions during this time, the subcellular localization of CDC14Bretro progressively shifted from the association with microtubules (stabilizing them) to an association with the endoplasmic reticulum. CDC14Bretro evolution represents a paradigm example of rapid, selectively driven subcellular relocalization, thus revealing a novel mode for the emergence of new gene function.     

Birth and Rapid Subcellular Adaptation of a Hominoid-Specific CDC14 Protein

Gene duplication was prevalent during hominoid evolution, yet little is known about the functional fate of new ape gene copies. We characterized the CDC14B cell cycle gene and the functional evolution of its hominoid-specific daughter gene, CDC14Bretro. We found that CDC14B encodes four different splice isoforms that show different subcellular localizations (nucleus or microtubule-associated) and functional properties. A microtubular CDC14B variant spawned CDC14Bretro through retroposition in the hominoid ancestor 18–25 million years ago (Mya). CDC14Bretro evolved brain-/testis-specific expression after the duplication event and experienced a short period of intense positive selection in the African ape ancestor 7–12 Mya. Using resurrected ancestral protein variants, we demonstrate that by virtue of amino acid substitutions in distinct protein regions during this time, the subcellular localization of CDC14Bretro progressively shifted from the association with microtubules (stabilizing them) to an association with the endoplasmic reticulum. CDC14Bretro evolution represents a paradigm example of rapid, selectively driven subcellular relocalization, thus revealing a novel mode for the emergence of new gene function.     

The phylogeny of human globin genes investigated by the maximum parsimony method

Gene phylogenetic trees were constructed by the maximum parsimony method for various sets of ninety six globin chain amino acid sequences spanning plant and animal kingdoms. The method, executed by several computer programs, constructed ancestor and descendant globin messengers on tree topologies which required the least number of nucleotide replacements to account for the evolution of the globins. The human myoglobin-hemoglobin divergence was traced to a gene duplication which occurred either in the first vertebrates or earlier yet in the common ancestor of chordates and annelids, the alpha-beta divergence to a gene duplication in the common ancestor of teleosts and tetrapods, the gamma divergence from typical beta chains to a gene duplication in basal therian mammals, and the delta separation from beta to a duplication in the basal catarrhine primates. Evidence was provided by the globin phylogenies for the hominoid affinities of the gibbon and the close phyletic relationship of the African apes to man. Over the period of teleos-tetrapod divergence the globin messengers evolved at an average rate of 18.5 nucleotide replacements per 100 codons per 10⁸ years, a faster rate than most previous estimates. Very fast and very slow rates were encountered in different globin lineages and at different stages of descent, reducing the effectiveness of globins as molecular clocks. Rates increased with gene duplication and decreased after selection discovered useful specializations in the products of genes which had previously been freer to accept mutations. The early eutherian radiation was characterized by rapid rates of globin evolution, but the later hominoid radiation by extremely slow rates. This pattern was related to more complicated grades of internal organization evolving in human ancestors. The types of nucleotide replacements in the globin messengers over the long course of globin evolution did not seem indicative of any special mutational mechanisms.