Evolutionary relationships within the fungi: analyses of nuclear small subunit rRNA sequences.

Nucleotide sequences of the small subunit ribosomal RNA (18S) gene were used to investigate evolutionary relationships within the Fungi. The inferred tree topologies are in general agreement with traditional classifications in the following ways: (1) the Chytridiomycota and Zygomycota appear to be basal groups within the Fungi. (2) The Ascomycota and Basidiomycota are a derived monophyletic group. (3) Relationships within the Ascomycota are concordant with traditional orders and divide the hemi- and euascomycetes into distinct lineages. (4) The Basidiomycota is divided between the holobasidiomycetes and phragmobasidiomycetes. Conflicts with traditional classification were limited to weakly supported branches of the tree. Strongly supported relationships were robust to minor changes in alignment, method of analysis, and various weighting schemes. Weighting, either of transversions or by site, did not convincingly improve the status of poorly supported portions of the tree. The rate of variation at particular sites does not appear to be independent of lineage, suggesting that covariation of sites may be an important phenomenon in these genes.     

Evolutionary conservation of a testes-specific proteasome subunit gene in Drosophila

Proteasomes are large multisubunit particles that act as the proteolytic machinery for the ubiquitin-dependent proteolytic pathway. The core of this complex, the 20S proteasome, is made up of seven α-type and seven β-type subunits, arranged in an (α1–α7)(β1–β7)(β1–β7)(α1–α7) configuration. Previous work had shown that there exist alternative isoforms of the Drosophila melanogaster α4-type subunit, encoded by two distinct genes, α4t1_dm and α4t2_dm, and that these are expressed exclusively in the germline of the testes. We sought to investigate the evolutionary conservation of this phenomenon by screening for orthologs of the α4-type gene family in the distantly related Drosophila species, D. virilis. We isolated the D. virilis orthologs of the somatically expressed gene, α4_dm, and the testes-specific gene, α4t2_dm. We failed to find an ortholog of the other testes-specific gene, α4t1_dm. The α4_dv gene maps to the X chromosome at 12A-C, its product shares 90% amino acid identity with α4_dm, and it is expressed at high levels in both males and females. The other gene, α4t_dv, encodes a protein most similar to the testes-specific α4t2_dm proteasome subunit (59% a.a. identity), and it maps to position 27 on chomosome 2. The expression of the α4t_dv gene is testes-specific, like that of α4t2_dm. The existence of testes-specific α4-type subunits in two widely diverged subgenera of Drosophila suggests that these subunit isoforms have important functional roles in spermatogenesis.     

Evolutionary origins of the eukaryotic shikimate pathway: gene fusions, horizontal gene transfer, and endosymbiotic replacements

Currently the shikimate pathway is reported as a metabolic feature of prokaryotes, ascomycete fungi, apicomplexans, and plants. The plant shikimate pathway enzymes have similarities to prokaryote homologues and are largely active in chloroplasts, suggesting ancestry from the plastid progenitor genome. Toxoplasma gondii, which also possesses an alga-derived plastid organelle, encodes a shikimate pathway with similarities to ascomycete genes, including a five-enzyme pentafunctional arom. These data suggests that the shikimate pathway and the pentafunctional arom either had an ancient origin in the eukaryotes or was conveyed by eukaryote-to-eukaryote horizontal gene transfer (HGT). We expand sampling and analyses of the shikimate pathway genes to include the oomycetes, ciliates, diatoms, basidiomycetes, zygomycetes, and the green and red algae. Sequencing of cDNA from Tetrahymena thermophila confirmed the presence of a pentafused arom, as in fungi and T. gondii. Phylogenies and taxon distribution suggest that the arom gene fusion event may be an ancient eukaryotic innovation. Conversely, the Plantae lineage (represented here by both Viridaeplantae and the red algae) acquired different prokaryotic genes for all seven steps of the shikimate pathway. Two of the phylogenies suggest a derivation of the Plantae genes from the cyanobacterial plastid progenitor genome, but if the full Plantae pathway was originally of cyanobacterial origin, then the five other shikimate pathway genes were obtained from a minimum of two other eubacterial genomes. Thus, the phylogenies demonstrate both separate HGTs and shared derived HGTs within the Plantae clade either by primary HGT transfer or secondarily via the plastid progenitor genome. The shared derived characters support the holophyly of the Plantae lineage and a single ancestral primary plastid endosymbiosis. Our analyses also pinpoints a minimum of 50 gene/domain loss events, demonstrating that loss and replacement events have been an important process in eukaryote genome evolution.     

Cloning of a yeast gene coding for the glutamate synthase small subunit (GUS2) by complementation of Saccharomyces cerevisiae and Escherichia coli glutamate auxotrophs

A Saccharomyces cerevisiae glutamate auxotroph, lacking NADP-glutamate dehydrogenase (NADP-GDH) and glutamate synthase (GOGAT) activities, was complemented with a yeast genomic library. Clones were obtained which still lacked NADP-GDH but showed GOGAT activity. Northern analysis revealed that the DNA fragment present in the complementing plasmids coded for a 1.5kb mRNA. Since the only GOGAT enzyme so far purified from S. cerevisiae is made up of a small and a large subunit, the size of the mRNA suggested that the cloned DNA fragment could code for the GOGAT small subunit. Plasmids were purified and used to transform Escherichia coli glutamate auxotrophs. Transformants were only recovered when the recipient strain was an E. coli GDH-less mutant lacking the small GOGAT subunit. These data show that we have cloned the structural gene coding for the yeast small subunit (GUS2). Evidence is also presented indicating that the GOGAT enzyme which is synthesized in the E. coli transformants is a hybrid comprising the large E. coli subunit and the small S. cerevisiae subunit.     

Simultaneous identification of duplications and lateral gene transfers

The incongruency between a gene tree and a corresponding species tree can be attributed to evolutionary events such as gene duplication and gene loss. This paper describes a combinatorial model where so-called DTL-scenarios are used to explain the differences between a gene tree and a corresponding species tree taking into account gene duplications, gene losses, and lateral gene transfers (also known as horizontal gene transfers). The reasonable biological constraint that a lateral gene transfer may only occur between contemporary species leads to the notion of acyclic DTL-scenarios. Parsimony methods are introduced by defining appropriate optimization problems. We show that finding most parsimonious acyclic DTL-scenarios is NP-hard. However, by dropping the condition of acyclicity, the problem becomes tractable, and we provide a dynamic programming algorithm as well as a fixed-parameter tractable algorithm for finding most parsimonious DTL-scenarios.     

Identification of an acetolactate synthase small subunit gene in two eukaryotes

Acetolactate synthase catalyses the first step in branched-chain amino acid biosynthesis. The bacterial enzyme contains two large and two small subunits but there is only limited and circumstantial evidence for a small subunit in the eukaryotic enzyme. Here this evidence is summarised and protein sequences of two putative eukaryotic small subunits, from a yeast and a red alga, are presented.     

Recent events dominate interdomain lateral gene transfers between prokaryotes and eukaryotes and,with the exception of endosymbiotic gene transfers,few ancient transfer events persist

While there is compelling evidence for the impact of endosymbiotic gene transfer (EGT; transfer from either mitochondrion or chloroplast to the nucleus) on genome evolution in eukaryotes, the role of interdomain transfer from bacteria and/or archaea (i.e. prokaryotes) is less clear. Lateral gene transfers (LGTs) have been argued to be potential sources of phylogenetic information, particularly for reconstructing deep nodes that are difficult to recover with traditional phylogenetic methods. We sought to identify interdomain LGTs by using a phylogenomic pipeline that generated 13 465 single gene trees and included up to 487 eukaryotes, 303 bacteria and 118 archaea. Our goals include searching for LGTs that unite major eukaryotic clades, and describing the relative contributions of LGT and EGT across the eukaryotic tree of life. Given the difficulties in interpreting single gene trees that aim to capture the approximately 1.8 billion years of eukaryotic evolution, we focus on presence–absence data to identify interdomain transfer events. Specifically, we identify 1138 genes found only in prokaryotes and representatives of three or fewer major clades of eukaryotes (e.g. Amoebozoa, Archaeplastida, Excavata, Opisthokonta, SAR and orphan lineages). The majority of these genes have phylogenetic patterns that are consistent with recent interdomain LGTs and, with the notable exception of EGTs involving photosynthetic eukaryotes, we detect few ancient interdomain LGTs. These analyses suggest that LGTs have probably occurred throughout the history of eukaryotes, but that ancient events are not maintained unless they are associated with endosymbiotic gene transfer among photosynthetic lineages.     

Evolutionary changes in the fungal carbamoyl-phosphate synthetase small subunit gene and its associated upstream open reading frame

The Neurospora crassa arg-2 and the Saccharomyces cerevisiae ortholog CPA1 encode the arginine-specific carbamoyl-phosphate synthetase (CPS-A) small subunit. Arginine decreases synthesis of this subunit through the action of a 5' upstream open reading frame in the mRNA that encodes a cis-regulatory element, the arginine attenuator peptide (AAP), which stalls ribosomes in response to arginine. We performed a comparative analysis of the genomic structure and predicted peptide sequence of the AAP and CPS-A small subunit across many fungi. Differences at the genomic level included variation in intron number and position within the AAP and CPS-A coding regions and differences in known regulatory motifs. Although differences exist in AAP sequence, there were three absolutely conserved amino acid residues in the predicted peptide, including an aspartic acid crucial for arginine-dependent regulation of arg-2 and CPA1. A diverged Basidiomycete AAP was shown to retain function as an Arg-specific negative regulator of translation.     

Evolution of bacterial phosphoglycerate mutases: non-homologous isofunctional enzymes undergoing gene losses, gains and lateral transfers

Background

The glycolytic phosphoglycerate mutases exist as non-homologous isofunctional enzymes (NISE) having independent evolutionary origins and no similarity in primary sequence, 3D structure, or catalytic mechanism. Cofactor-dependent PGM (dPGM) requires 2,3-bisphosphoglycerate for activity; cofactor-independent PGM (iPGM) does not. The PGM profile of any given bacterium is unpredictable and some organisms such as Escherichia coli encode both forms.

Methods/Principal Findings

To examine the distribution of PGM NISE throughout the Bacteria, and gain insight into the evolutionary processes that shape their phyletic profiles, we searched bacterial genome sequences for the presence of dPGM and iPGM. Both forms exhibited patchy distributions throughout the bacterial domain. Species within the same genus, or even strains of the same species, frequently differ in their PGM repertoire. The distribution is further complicated by the common occurrence of dPGM paralogs, while iPGM paralogs are rare. Larger genomes are more likely to accommodate PGM paralogs or both NISE forms. Lateral gene transfers have shaped the PGM profiles with intradomain and interdomain transfers apparent. Archaeal-type iPGM was identified in many bacteria, often as the sole PGM. To address the function of PGM NISE in an organism encoding both forms, we analyzed recombinant enzymes from E. coli. Both NISE were active mutases, but the specific activity of dPGM greatly exceeded that of iPGM, which showed highest activity in the presence of manganese. We created PGM null mutants in E. coli and discovered the ΔdPGM mutant grew slowly due to a delay in exiting stationary phase. Overexpression of dPGM or iPGM overcame this defect.

Conclusions/Significance

Our biochemical and genetic analyses in E. coli firmly establish dPGM and iPGM as NISE. Metabolic redundancy is indicated since only larger genomes encode both forms. Non-orthologous gene displacement can fully account for the non-uniform PGM distribution we report across the bacterial domain.     

Genetic fusions of globular proteins to the epsilon subunit of the Escherichia coli ATP synthase: Implications for in vivo rotational catalysis and epsilon subunit function

The rotational mechanism of ATP synthase was investigated by fusing three proteins from Escherichia coli, the 12-kDa soluble cytochrome b(562), the 20-kDa flavodoxin, and the 28-kDa flavodoxin reductase, to the C terminus of the epsilon subunit of the enzyme. According to the concept of rotational catalysis, because epsilon is part of the rotor a large domain added at this site should sterically clash with the second stalk, blocking rotation and fully inhibiting the enzyme. E. coli cells expressing the cytochrome b(562) fusion in place of wild-type epsilon grew using acetate as the energy source, indicating their capacity for oxidative phosphorylation. Cells expressing the larger flavodoxin or flavodoxin reductase fusions failed to grow on acetate. Immunoblot analysis showed that the fusion proteins were stable in the cells and that they had no effect on enzyme assembly. These results provide initial evidence supporting rotational catalysis in vivo. In membrane vesicles, the cytochrome b(562) fusion caused an increase in the apparent ATPase activity but a minor decrease in proton pumping. Vesicles bearing ATP synthase containing the larger fusion proteins showed reduced but significant levels of ATPase activity that was sensitive to inhibition by dicyclohexylcarbodiimide (DCCD) but no proton pumping. Thus, all fusions to epsilon generated an uncoupled component of ATPase activity. These results imply that a function of the C terminus of epsilon in F(1)F(0) is to increase the efficiency of the enzyme by specifically preventing the uncoupled hydrolysis of ATP. Given the sensitivity to DCCD, this uncoupled ATP hydrolysis may arise from rotational steps of gammaepsilon in the inappropriate direction after ATP is bound at the catalytic site. It is proposed that the C-terminal domain of epsilon functions to ensure that rotation occurs only in the direction of ATP synthesis when ADP is bound and only in the direction of hydrolysis when ATP is bound.     

Evolutionary conservation of gene structures of the Pax1/9 gene family

Based on amino acid sequence comparisons, Pax1 and Pax9 genes are considered to form a subgroup of vertebrate Pax genes. We show here that the gene structures of mouse Pax1, human PAX9 genes are similar to that of a single Pax1/9 related gene in Branchiostoma lanceolatum, AmphiPax1. This supports the hypothesis that Pax1 and Pax9 genes were derived from a single ancestral gene. A refined protein alignment of AmphiPax1, mouse Pax1 and human PAX9 proteins based on the determined exon boundaries indicates that sequence divergence at the C-termini may be related to the unique functions of the Pax1 and Pax9 genes in vertebrates. AmphiPax1 is expressed in adult amphioxus in the pharyngeal endoderm.     

Evolutionary conservation of the U7 small nuclear ribonucleoprotein in Drosophila melanogaster

The U7 snRNP involved in histone RNA 3' end processing is related to but biochemically distinct from spliceosomal snRNPs. In vertebrates, the Sm core structure assembling around the noncanonical Sm-binding sequence of U7 snRNA contains only five of the seven standard Sm proteins. The missing Sm D1 and D2 subunits are replaced by U7-specific Sm-like proteins Lsm10 and Lsm11, at least the latter of which is important for histone RNA processing. So far, it was unknown if this special U7 snRNP composition is conserved in invertebrates. Here we describe several putative invertebrate Lsm10 and Lsm11 orthologs that display low but clear sequence similarity to their vertebrate counterparts. Immunoprecipitation studies in Drosophila S2 cells indicate that the Drosophila Lsm10 and Lsm11 orthologs (dLsm10 and dLsm11) associate with each other and with Sm B, but not with Sm D1 and D2. Moreover, dLsm11 associates with the recently characterized Drosophila U7 snRNA and, indirectly, with histone H3 pre-mRNA. Furthermore, dLsm10 and dLsm11 can assemble into U7 snRNPs in mammalian cells. These experiments demonstrate a strong evolutionary conservation of the unique U7 snRNP composition, despite a high degree of primary sequence divergence of its constituents. Therefore, Drosophila appears to be a suitable system for further genetic studies of the cell biology of U7 snRNPs.     

Evolutionary conservation and disease gene association of the human genes composing pseudogenes

Pseudogenes, the 'genomic fossils' present portrayal of evolutionary history of human genome. The human genes configuring pseudogenes are also now coming forth as important resources in the study of human protein evolution. In this communication, we explored evolutionary conservation of the genes forming pseudogenes over the genes lacking any pseudogene and delving deeper, we probed an evolutionary rate difference between the disease genes in the two groups. We illustrated this differential evolutionary pattern by gene expressivity, number of regulatory miRNA targeting per gene, abundance of protein complex forming genes and lesser percentage of protein intrinsic disorderness. Furthermore, pseudogenes are observed to harbor sequence variations, over their entirety, those become degenerative disease-causing mutations though the disease involvement of their progenitors is still unexplored. Here, we unveiled an immense association of disease genes in the genes casting pseudogenes in human. We interpreted the issue by disease associated miRNA targeting, genes containing polymorphisms in miRNA target sites, abundance of genes having disease causing non-synonymous mutations, disease gene specific network properties, presence of genes having repeat regions, affluence of dosage sensitive genes and the presence of intrinsically unstructured protein regions.     

Variable copy number, intra-genomic heterogeneities and lateral transfers of the 16S rRNA gene in Pseudomonas

Even though the 16S rRNA gene is the most commonly used taxonomic marker in microbial ecology, its poor resolution is still not fully understood at the intra-genus level. In this work, the number of rRNA gene operons, intra-genomic heterogeneities and lateral transfers were investigated at a fine-scale resolution, throughout the Pseudomonas genus. In addition to nineteen sequenced Pseudomonas strains, we determined the 16S rRNA copy number in four other Pseudomonas strains by Southern hybridization and Pulsed-Field Gel Electrophoresis, and studied the intra-genomic heterogeneities by Denaturing Gradient Gel Electrophoresis and sequencing. Although the variable copy number (from four to seven) seems to be correlated with the evolutionary distance, some close strains in the P. fluorescens lineage showed a different number of 16S rRNA genes, whereas all the strains in the P. aeruginosa lineage displayed the same number of genes (four copies). Further study of the intra-genomic heterogeneities revealed that most of the Pseudomonas strains (15 out of 19 strains) had at least two different 16S rRNA alleles. A great difference (5 or 19 nucleotides, essentially grouped near the V1 hypervariable region) was observed only in two sequenced strains. In one of our strains studied (MFY30 strain), we found a difference of 12 nucleotides (grouped in the V3 hypervariable region) between copies of the 16S rRNA gene. Finally, occurrence of partial lateral transfers of the 16S rRNA gene was further investigated in 1803 full-length sequences of Pseudomonas available in the databases. Remarkably, we found that the two most variable regions (the V1 and V3 hypervariable regions) had probably been laterally transferred from another evolutionary distant Pseudomonas strain for at least 48.3 and 41.6% of the 16S rRNA sequences, respectively. In conclusion, we strongly recommend removing these regions of the 16S rRNA gene during the intra-genus diversity studies.     

Primary and secondary structure of the nuclear small subunit ribosomal RNA of the cryptomonadPyrenomonas salina as inferred from the gene sequence: Evolutionary implications

Summary The cryptomonadPyrenomonas salina presumably has arisen from a symbiotic event involving a flagellated phagotrophic host cell and a photosynthetic eukaryote as the symbiont. Correspondingly, in this unicellular alga there are four different genomes, e.g., the nuclear and the mitochondrial genomes of the host cell as well as the plastid genome and the genome contained in the vestigial nucleus of the endocytobiont (nucleomorph). To analyze the orgin of one of the symbiotic partners the small subunit rRNA gene sequence of the host cell nucleus was determined, and a secondary structure model has been constructed. This sequence is compared to those of 40 other eukaryotes. A phylogenetic tree constructed using the neighborliness method revealed a close relationship between the host cell ofP. salina and the chlorophytes, whereas the rhodophytes diverge more deeply in the tree.     

Molecular cloning and characterization of a large subunit of Salmonella typhimurium glutamate synthase (GOGAT) gene in Escherichia coli

Two pathways of ammonium assimilation and glutamate biosynthesis have been identified in microorganisms. One pathway involves the NADP-linked glutamate dehydrogenase, which catalyzes the amination of 2-oxoglutarate to form glutamate. An alternative pathway involves the combined activities of glutamine synthetase, which aminates glutamate to form glutamine, and glutamate synthase, which transfers the amide group of glutamine to 2-oxoglutarate to yield two molecules of glutamate. We have cloned the large subunit of the glutamate synthase (GOGAT) from Salmonella typhimurium by screening the expression of GOGAT and complementing the gene in E. coli GOGAT large subunit-deficient mutants. Three positive clones (named pUC19C12, pUC19C13 and pUC19C15) contained identical Sau3AI fragments, as determined by restriction mapping and Southern hybridization, and expressed GOGAT efficiently and constitutively using its own promoter in the heterologous host. The coding region expressed in Escherichia coli was about 170 kDa on SDS-PAGE. This gene spans 4,732 bases, contains an open reading frame of 4,458 nucleotides, and encodes a mature protein of 1,486 amino acid residues (Mr = 166,208). The FMN-binding domain of GOGAT contains 12 glycine residues, and the 3Fe-4S cluster has 3 cysteine residues. The comparison of the translated amino acid sequence of the Salmonella GOGAT with sequences from other bacteria such as Escherichia coli, Salmonella enterica, Shigella flexneri, Yersinia pestis, Vibrio vulnificus and Pseudomonas aeruginosa shows sequence identity between 87 and 95%.