A genome phylogeny for mitochondria among alpha-proteobacteria and a predominantly eubacterial ancestry of yeast nuclear genes

Analyses of 55 individual and 31 concatenated protein data sets encoded in Reclinomonas americana and Marchantia polymorpha mitochondrial genomes revealed that current methods for constructing phylogenetic trees are insufficiently sensitive (or artifact-insensitive) to ascertain the sister of mitochondria among the current sample of eight alpha-proteobacterial genomes using mitochondrially-encoded proteins. However, Rhodospirillum rubrum came as close to mitochondria as any alpha-proteobacterium investigated. This prompted a search for methods to directly compare eukaryotic genomes to their prokaryotic counterparts to investigate the origin of the mitochondrion and its host from the standpoint of nuclear genes. We examined pairwise amino acid sequence identity in comparisons of 6,214 nuclear protein-coding genes from Saccharomyces cerevisiae to 177,117 proteins encoded in sequenced genomes from 45 eubacteria and 15 archaebacteria. The results reveal that approximately 75% of yeast genes having homologues among the present prokaryotic sample share greater amino acid sequence identity to eubacterial than to archaebacterial homologues. At high stringency comparisons, only the eubacterial component of the yeast genome is detectable. Our findings indicate that at the levels of overall amino acid sequence identity and gene content, yeast shares a sister-group relationship with eubacteria, not with archaebacteria, in contrast to the current phylogenetic paradigm based on ribosomal RNA. Among eubacteria and archaebacteria, proteobacterial and methanogen genomes, respectively, shared more similarity with the yeast genome than other prokaryotic genomes surveyed.     

A plant locus essential for phylloquinone (vitamin K1) biosynthesis originated from a fusion of four eubacterial genes

Phylloquinone is a compound present in all photosynthetic plants serving as cofactor for Photosystem I-mediated electron transport. Newly identified seedling-lethal Arabidopsis thaliana mutants impaired in the biosynthesis of phylloquinone possess reduced Photosystem I activity. The affected gene, called PHYLLO, consists of a fusion of four previously individual eubacterial genes, menF, menD, menC, and menH, required for the biosynthesis of phylloquinone in photosynthetic cyanobacteria and the respiratory menaquinone in eubacteria. The fact that homologous men genes reside as polycistronic units in eubacterial chromosomes and in plastomes of red algae strongly suggests that PHYLLO derived from a plastid operon during endosymbiosis. The principle architecture of the fused PHYLLO locus is conserved in the nuclear genomes of plants, green algae, and the diatom alga Thalassiosira pseudonana. The latter arose from secondary endosymbiosis of a red algae and a eukaryotic host indicating selective driving forces for maintenance and/or independent generation of the composite gene cluster within the nuclear genomes. Besides, individual menF genes, encoding active isochorismate synthases (ICS), have been established followed by splitting of the essential 3' region of the menF module of PHYLLO only in genomes of higher plants. This resulted in inactivation of the ICS activity encoded by PHYLLO and enabled a metabolic branch from the phylloquinone biosynthetic route to independently regulate the synthesis of salicylic acid required for plant defense. Therefore, gene fusion, duplication, and fission events adapted a eubacterial multienzymatic system to the metabolic requirements of plants.     

Evolution and Horizontal Transfer of dUTPase-Encoding Genes in Viruses and Their Hosts

dUTPase is a ubiquitous and essential enzyme responsible for regulating cellular levels of dUTP. The dut gene exists as single, tandemly duplicated, and tandemly triplicated copies. Crystallized single-copy dUTPases have been shown to assemble as homotrimers. dUTPase is encoded as an auxiliary gene in a number of virus genomes. The origin of viral dut genes has remained unresolved since their initial discovery. A comprehensive analysis of dUTPase amino acid sequence relationships was performed to explore the evolutionary dynamics of dut in viruses and their hosts. Our data set, comprised of 24 host and 51 viral sequences, includes representative sequences from available eukaryotes, archaea, eubacteria cells, and viruses, including herpesviruses. These amino acid sequences were aligned by using a hidden Markov model approach developed to align divergent data. Known secondary structures from single-copy crystals were mapped onto the aligned duplicate and triplicate sequences. We show how duplicated dUTPases might fold into a monomer, and we hypothesize that triplicated dUTPases also assemble as monomers. Phylogenetic analysis revealed at least five viral dUTPase sequence lineages in well-supported monophyletic clusters with eukaryotic, eubacterial, and archaeal hosts. We have identified all five as strong examples of horizontal transfer as well as additional potential transfer of dut genes among eubacteria, between eubacteria and viruses, and between retroviruses. The evidence for horizontal transfers is particularly interesting since eukaryotic dut genes have introns, while DNA virus dut genes do not. This implies that an intermediary retroid agent facilitated the horizontal transfer process between host mRNA and DNA viruses.     

Identification and Characterization of a Gene Cluster for Synthesis of the Polyketide Antibiotic 2,4-Diacetylphloroglucinol from Pseudomonas fluorescens Q2-87

The polyketide metabolite 2,4-diacetylphloroglucinol (2,4-DAPG) is produced by many strains of fluorescent Pseudomonas spp. with biocontrol activity against soilborne fungal plant pathogens. Genes required for 2,4-DAPG synthesis by P. fluorescens Q2-87 are encoded by a 6.5-kb fragment of genomic DNA that can transfer production of 2,4-DAPG to 2,4-DAPG-nonproducing recipient Pseudomonas strains. In this study the nucleotide sequence was determined for the 6.5-kb fragment and flanking regions of genomic DNA from strain Q2-87. Six open reading frames were identified, four of which (phlACBD) comprise an operon that includes a set of three genes (phlACB) conserved between eubacteria and archaebacteria and a gene (phlD) encoding a polyketide synthase with homology to chalcone and stilbene synthases from plants. The biosynthetic operon is flanked on either side by phlE and phlF, which code respectively for putative efflux and regulatory (repressor) proteins. Expression in Escherichia coli of phlA, phlC, phlB, and phlD, individually or in combination, identified a novel polyketide biosynthetic pathway in which PhlD is responsible for the production of monoacetylphloroglucinol (MAPG). PhlA, PhlC, and PhlB are necessary to convert MAPG to 2,4-DAPG, and they also may function in the synthesis of MAPG.     

The Evolution of Fungal Metabolic Pathways

Fungi contain a remarkable range of metabolic pathways, sometimes encoded by gene clusters, enabling them to digest most organic matter and synthesize an array of potent small molecules. Although metabolism is fundamental to the fungal lifestyle, we still know little about how major evolutionary processes, such as gene duplication (GD) and horizontal gene transfer (HGT), have interacted with clustered and non-clustered fungal metabolic pathways to give rise to this metabolic versatility. We examined the synteny and evolutionary history of 247,202 fungal genes encoding enzymes that catalyze 875 distinct metabolic reactions from 130 pathways in 208 diverse genomes. We found that gene clustering varied greatly with respect to metabolic category and lineage; for example, clustered genes in Saccharomycotina yeasts were overrepresented in nucleotide metabolism, whereas clustered genes in Pezizomycotina were more common in lipid and amino acid metabolism. The effects of both GD and HGT were more pronounced in clustered genes than in their non-clustered counterparts and were differentially distributed across fungal lineages; specifically, GD, which was an order of magnitude more abundant than HGT, was most frequently observed in Agaricomycetes, whereas HGT was much more prevalent in Pezizomycotina. The effect of HGT in some Pezizomycotina was particularly strong; for example, we identified 111 HGT events associated with the 15 Aspergillus genomes, which sharply contrasts with the 60 HGT events detected for the 48 genomes from the entire Saccharomycotina subphylum. Finally, the impact of GD within a metabolic category was typically consistent across all fungal lineages, whereas the impact of HGT was variable. These results indicate that GD is the dominant process underlying fungal metabolic diversity, whereas HGT is episodic and acts in a category- or lineage-specific manner. Both processes have a greater impact on clustered genes, suggesting that metabolic gene clusters represent hotspots for the generation of fungal metabolic diversity.     

Carotenoid biosynthetic pathway: molecular phylogenies and evolutionary behavior of crt genes in eubacteria

Phylogenetic analysis of carotenoid biosynthetic pathway genes and their evolutionary rate variations were studied among eubacterial taxa. The gene sequences for the enzymes involved in this pathway were obtained for major phylogenetic groups of eubacteria (green sulfur bacteria, green nonsulphur bacteria, Gram-positive bacteria, proteobacteria, flavobacteria, cyanobacteria) and archeabacteria. These gene datasets were distributed under five major steps of carotenoid biosynthesis in eubacteria; isoprenoid precursor biosynthesis, phytoene synthesis, dehydrogenation of phytoene, lycopene cyclization, formation of acyclic xanthophylls, formation of cyclic xanthophylls and carotenoid biosynthesis regulation. The NJ algorithm was used on protein coding DNA sequences to deduce the evolutionary relationship for the respective crt genes among different eubacterial lineages. The rate of nonsynonymous nucleotide substitutions per nonsynonymous site (d(N)) and synonymous nucleotide substitutions per synonymous site (d(S)) were calculated for different clades of the respective phylogenetic tree for specific crt genes. The phylogenetic analysis suggests that evolutionary pattern of crt genes in eubacteria is characterized by lateral gene transfer and gene duplication events. The d(N) values indicate that carotenoid biosynthetic genes are more conserved in proteobacteria than in any other eubacterial phyla. Furthermore, of the genes involved in carotenoid biosynthesis pathway, structural genes evolve slowly than the regulatory genes in eubacteria.     

Biosynthesis of riboflavin: an unusual riboflavin synthase of Methanobacterium thermoautotrophicum.

Riboflavin synthase was purified by a factor of about 1,500 from cell extract of Methanobacterium thermoautotrophicum. The enzyme had a specific activity of about 2,700 nmol mg(-1) h(-1) at 65 degrees C, which is relatively low compared to those of riboflavin synthases of eubacteria and yeast. Amino acid sequences obtained after proteolytic cleavage had no similarity with known riboflavin synthases. The gene coding for riboflavin synthase (designated ribC) was subsequently cloned by marker rescue with a ribC mutant of Escherichia coli. The ribC gene of M. thermoautotrophicum specifies a protein of 153 amino acid residues. The predicted amino acid sequence agrees with the information gleaned from Edman degradation of the isolated protein and shows 67% identity with the sequence predicted for the unannotated reading frame MJ1184 of Methanococcus jannaschii. The ribC gene is adjacent to a cluster of four genes with similarity to the genes cbiMNQO of Salmonella typhimurium, which form part of the cob operon (this operon contains most of the genes involved in the biosynthesis of vitamin B12). The amino acid sequence predicted by the ribC gene of M. thermoautotrophicum shows no similarity whatsoever to the sequences of riboflavin synthases of eubacteria and yeast. Most notably, the M. thermoautotrophicum protein does not show the internal sequence homology characteristic of eubacterial and yeast riboflavin synthases. The protein of M. thermoautotrophicum can be expressed efficiently in a recombinant E. coli strain. The specific activity of the purified, recombinant protein is 1,900 nmol mg(-1) h(-1) at 65 degrees C. In contrast to riboflavin synthases from eubacteria and fungi, the methanobacterial enzyme has an absolute requirement for magnesium ions. The 5' phosphate of 6,7-dimethyl-8-ribityllumazine does not act as a substrate. The findings suggest that riboflavin synthase has evolved independently in eubacteria and methanobacteria.     

Large-Scale Phylogenetic Classification of Fungal Chitin Synthases and Identification of a Putative Cell-Wall Metabolism Gene Cluster in Aspergillus Genomes

The cell wall is a protective and versatile structure distributed in all fungi. The component responsible for its rigidity is chitin, a product of chitin synthase (Chsp) enzymes. There are seven classes of chitin synthase genes (CHS) and the amount and type encoded in fungal genomes varies considerably from one species to another. Previous Chsp sequence analyses focused on their study as individual units, regardless of genomic context. The identification of blocks of conserved genes between genomes can provide important clues about the interactions and localization of chitin synthases. On the present study, we carried out an in silico search of all putative Chsp encoded in 54 full fungal genomes, encompassing 21 orders from five phyla. Phylogenetic studies of these Chsp were able to confidently classify 347 out of the 369 Chsp identified (94%). Patterns in the distribution of Chsp related to taxonomy were identified, the most prominent being related to the type of fungal growth. More importantly, a synteny analysis for genomic blocks centered on class IV Chsp (the most abundant and widely distributed Chsp class) identified a putative cell wall metabolism gene cluster in members of the genus Aspergillus, the first such association reported for any fungal genome.     

On the evolutionary descent of organisms and organelles: a global phylogeny based on a highly conserved structural core in small subunit ribosomal RNA

To probe the earliest evolutionary events attending the origin of the five known genome types (archaebacterial, eubacterial, nuclear, mitochondrial and plastid), we have analyzed sequences corresponding to a ubiquitous, highly conserved core of secondary structure in small subunit rRNA. Our results support (i) the existence of three primary lineages (archaebacterial, eubacterial, and nuclear), (ii) a specific eubacterial ancestry for plastids and mitochondria (plant, animal, fungal), and (iii) an endosymbiotic, evolutionary origin of the two types of organelle from within distinct groups of eubacteria (blue-green algae (cyanobacteria) in the case of plastids, nonphotosynthetic aerobic bacteria in the case of mitochondria). In addition, our analysis suggests (iv) a biphyletic origin of mitochondria, with animal and fungal mitochondria branching together but separately from plant mitochondria, and (v) a monophyletic origin of plastids. The method described here provides a powerful and generally applicable molecular taxonomic approach towards a global phylogeny encompassing all organisms and organelles.     

Cloning, expression, and substrate specificity of a fungal chymotrypsin. Evidence for lateral gene transfer from an actinomycete bacterium

Unlike trypsins, chymotrypsins have not until now been found in fungi. Expressed sequence tag analysis of the deuteromycete Metarhizium anisopliae identified two trypsins (family S1) and a novel chymotrypsin (CHY1). CHY1 resembles actinomycete (bacterial) chymotrypsins (family S2) rather than other eukaryote enzymes (family S1) in being synthesized as a precursor species (374 amino acids, pI/MW: 5.07/38,279) containing a large N-terminal fragment (186 amino acids). Chy1 was expressed in Pichia pastoris yielding an enzyme with a chymotrypsin specificity for branched aliphatic and aromatic C-terminal amino acids. This is predictable as key catalytic residues determining the specificity of Streptomyces griseus chymotrypsins are conserved with CHY1. Mature (secreted) CHY1 (pI/MW: 8.29/18,499) shows closest overall amino acid identity to S. griseus protease C (55%) and clustered with other secreted bacterial S2 chymotrypsins that diverged widely from animal and endocellular bacterial enzymes in phylogenetic trees of the chymotrypsin superfamily. Conversely, actinomycete chymotrypsins are much more closely related to fungal proteases than to other eubacterial sequences. Complete genomes of yeast, gram eubacteria, archaebacteria, and mitochondria do not contain paralogous genes. Expressed sequence tag data bases from other fungi also lack chymotrypsin homologs. In light of this patchy distribution, we conclude that chy1 probably arose by lateral gene transfer from an actinomycete bacterium.     

Evolutionary origins of Hsp90 chaperones and a deep paralogy in their bacterial ancestors

The 82-90 kD family of molecular chaperone proteins has homologs in eukaryotes (Hsp90) and many eubacteria (HtpG) but not in Archaebacteria. We used representatives of all four different eukaryotic paralogs (cytosolic, endoplasmic reticulum (ER), chloroplast, mitochondrial) together with numerous eubacterial HtpG proteins for phylogenetic analyses to investigate their evolutionary origins. Our trees confirm that none of the organellar Hsp90s derives from the endosymbionts of early eukaryotes. Contrary to previous suggestions of distant origins through lateral gene transfer (LGT) all eukaryote Hsp90s are related to Gram-positive eubacterial HtpG proteins. The nucleocytosolic, ER and chloroplast Hsp90 paralogs are clearly mutually related. The origin of mitochondrial Hsp90 is more obscure, as these sequences are deeply nested within eubacteria. Our trees also reveal a deep split within eubacteria into a group of mainly long-branching sequences (including the eukaryote mitochondrial Hsp90s) and another group comprising exclusively short-branching HtpG proteins, from which the cytosolic/ER versions probably arose. Both versions are present in several eubacterial phyla, suggesting gene duplication very early in eubacterial evolution and multiple independent losses thereafter. We identified one probable case of LGT within eubacteria. However, multiple losses can simply explain the evolutionary pattern of the eubacterial HtpG paralogs and predominate over LGT. We suggest that the actinobacterial ancestor of eukaryotes harbored genes for both eubacterial HtpG paralogs, as the actinobacterium Streptomyces coelicolor still does; one could have given rise to the mitochondrial Hsp90 and the other, following another duplication event in the ancestral eukaryote, to the cytosolic and ER Hsp90 homologs.     

Widespread occurrence of secondary lipid biosynthesis potential in microbial lineages

Bacterial production of long-chain omega-3 polyunsaturated fatty acids (PUFAs), such as eicosapentaenoic acid (EPA, 20:5n-3) and docosahexaenoic acid (DHA, 22:6n-3), is constrained to a narrow subset of marine γ-proteobacteria. The genes responsible for de novo bacterial PUFA biosynthesis, designated pfaEABCD, encode large, multi-domain protein complexes akin to type I iterative fatty acid and polyketide synthases, herein referred to as "Pfa synthases". In addition to the archetypal Pfa synthase gene products from marine bacteria, we have identified homologous type I FAS/PKS gene clusters in diverse microbial lineages spanning 45 genera representing 10 phyla, presumed to be involved in long-chain fatty acid biosynthesis. In total, 20 distinct types of gene clusters were identified. Collectively, we propose the designation of "secondary lipids" to describe these biosynthetic pathways and products, a proposition consistent with the "secondary metabolite" vernacular. Phylogenomic analysis reveals a high degree of functional conservation within distinct biosynthetic pathways. Incongruence between secondary lipid synthase functional clades and taxonomic group membership combined with the lack of orthologous gene clusters in closely related strains suggests horizontal gene transfer has contributed to the dissemination of specialized lipid biosynthetic activities across disparate microbial lineages.     

Evolutionary implications of bacterial polyketide synthases

Polyketide synthases (PKS) perform a stepwise biosynthesis of diverse carbon skeletons from simple activated carboxylic acid units. The products of the complex pathways possess a wide range of pharmaceutical properties, including antibiotic, antitumor, antifungal, and immunosuppressive activities. We have performed a comprehensive phylogenetic analysis of multimodular and iterative PKS of bacteria and fungi and of the distinct types of fatty acid synthases (FAS) from different groups of organisms based on the highly conserved ketoacyl synthase (KS) domains. Apart from enzymes that meet the classification standards we have included enzymes involved in the biosynthesis of mycolic acids, polyunsaturated fatty acids (PUFA), and glycolipids in bacteria. This study has revealed that PKS and FAS have passed through a long joint evolution process, in which modular PKS have a central position. They appear to have derived from bacterial FAS and primary iterative PKS and, in addition, share a common ancestor with animal FAS and secondary iterative PKS. Furthermore, we have carried out a phylogenomic analysis of all modular PKS that are encoded by the complete eubacterial genomes currently available in the database. The phylogenetic distribution of acyltransferase and KS domain sequences revealed that multiple gene duplications, gene losses, as well as horizontal gene transfer (HGT) have contributed to the evolution of PKS I in bacteria. The impact of these factors seems to vary considerably between the bacterial groups. Whereas in actinobacteria and cyanobacteria the majority of PKS I genes may have evolved from a common ancestor, several lines of evidence indicate that HGT has strongly contributed to the evolution of PKS I in proteobacteria. Discovery of new evolutionary links between PKS and FAS and between the different PKS pathways in bacteria may help us in understanding the selective advantage that has led to the evolution of multiple secondary metabolite biosyntheses within individual bacteria.