Rickettsial evolution in the light of comparative genomics

Rickettsia are best known as strictly intracellular vector‐borne bacteria that cause mild to severe diseases in humans and other animals. Recent advances in molecular tools and biological experiments have unveiled a wide diversity of Rickettsia spp. that include species with a broad host range and some species that act as endosymbiotic associates. Molecular phylogenies of Rickettsia spp. contain some ambiguities, such as the position of R. canadensis and relationships within the spotted fever group. In the modern era of genomics, with an ever‐increasing number of sequenced genomes, there is enhanced interest in the use of whole‐genome sequences to understand pathogenesis and assess evolutionary relationships among rickettsial species. Rickettsia have small genomes (1.1–1.5 Mb) as a result of reductive evolution. These genomes contain split genes, gene remnants and pseudogenes that, owing to the colinearity of some rickettsial genomes, may represent different steps of the genome degradation process. Genomics reveal extreme genome reduction and massive gene loss in highly vertebrate‐pathogenic Rickettsia compared to less virulent or endosymbiotic species. Information gleaned from rickettsial genomics challenges traditional concepts of pathogenesis that focused primarily on the acquisition of virulence factors. Another intriguing phenomenon about the reduced rickettsial genomes concerns the large fraction of non‐coding DNA and possible functionality of these “non‐coding” sequences, because of the high conservation of these regions. Despite genome streamlining, Rickettsia spp. contain gene families, selfish DNA, repeat palindromic elements and genes encoding eukaryotic‐like motifs. These features participate in sequence and functional diversity and may play a crucial role in adaptation to the host cell and pathogenesis. Genome analyses have identified a large fraction of mobile genetic elements, including plasmids, suggesting the possibility of lateral gene transfer in these intracellular bacteria. Phylogenetic analyses have identified several candidates for horizontal gene acquisition among Rickettsia spp. including tra, pat2, and genes encoding for the type IV secretion system and ATP/ADP translocase that may have been acquired from bacteria living in amoebae. Gene loss, gene duplication, DNA repeats and lateral gene transfer all have shaped rickettsial genome evolution. A comprehensive analysis of the entire genome, including genes and non‐coding DNA, will help to unlock the mysteries of rickettsial evolution and pathogenesis.     

Origin and evolution of the light-dependent protochlorophyllide oxidoreductase (LPOR) genes

Abstract: Light‐dependent NADPH‐protochlorophyllide oxidoreductase (LPOR) is a nuclear‐encoded chloroplast protein in green algae and higher plants which catalyzes the light‐dependent reduction of protochlorophyllide to chlorophyllide. Light‐dependent chlorophyll biosynthesis occurs in all oxygenic photosynthetic organisms. With the exception of angiosperms, this pathway coexists with a separate light‐independent chlorophyll biosynthetic pathway, which is catalyzed by light‐independent protochlorophyllide reductase (DPOR) in the dark. In contrast, the light‐dependent function of chlorophyll biosynthesis is absent from anoxygenic photosynthetic bacteria. Consequently, the question is whether cyanobacteria are the ancestors of all organisms that conduct light‐dependent chlorophyll biosynthesis. If so, how did photosynthetic eukaryotes acquire the homologous genes of LPOR in their nuclear genomes? The large number of complete genome sequences now available allow us to detect the evolutionary history of LPOR genes by conducting a genome‐wide sequence comparison and phylogenetic analysis. Here, we show the results of a detailed phylogenetic analysis of LPOR and other functionally related enzymes in the short chain dehydrogenase/reductase (SDR) family. We propose that the LPOR gene originated in the cyanobacterial genome before the divergence of eukaryotic photosynthetic organisms. We postulated that the photosynthetic eukaryotes obtained their LPOR homologues through endosymbiotic gene transfer.     

CEGMA: a pipeline to accurately annotate core genes in eukaryotic genomes

MOTIVATION: The numbers of finished and ongoing genome projects are increasing at a rapid rate, and providing the catalog of genes for these new genomes is a key challenge. Obtaining a set of well-characterized genes is a basic requirement in the initial steps of any genome annotation process. An accurate set of genes is needed in order to learn about species-specific properties, to train gene-finding programs, and to validate automatic predictions. Unfortunately, many new genome projects lack comprehensive experimental data to derive a reliable initial set of genes. RESULTS: In this study, we report a computational method, CEGMA (Core Eukaryotic Genes Mapping Approach), for building a highly reliable set of gene annotations in the absence of experimental data. We define a set of conserved protein families that occur in a wide range of eukaryotes, and present a mapping procedure that accurately identifies their exon-intron structures in a novel genomic sequence. CEGMA includes the use of profile-hidden Markov models to ensure the reliability of the gene structures. Our procedure allows one to build an initial set of reliable gene annotations in potentially any eukaryotic genome, even those in draft stages. AVAILABILITY: Software and data sets are available online at http://korflab.ucdavis.edu/Datasets.     

Species-specific patterns of DNA bending and sequence.

Nucleotide sequences in the GenEMBL database were analyzed using strategies designed to reveal species-specific patterns of DNA bending and DNA sequence. The results uncovered striking species-dependent patterns of bending with more variations among individual organisms than between prokaryotes and eukaryotes. The frequency of bent sites in sequences from different bacteria was related to genomic A + T content and this relationship was confirmed by electrophoretic analysis of genomic DNA. However, base composition was not an accurate predictor for DNA bending in eukaryotes. Sequences from C. elegans exhibited the highest frequency of bent sites in the database and the RNA polymerase II locus from the nematode was the most bent gene in GenEMBL. Bent DNA extended throughout most introns and gene flanking segments from C.elegans while exon regions lacked A-tract bending characteristics. Independent evidence for the strong bending character of this genome was provided by electrophoretic studies which revealed that a large number of the fragments from C.elegans DNA exhibited anomalous gel mobilities when compared to genomic fragments from over 20 other organisms. The prevalence of bent sites in this genome enabled us to detect selectively C.elegans sequences in a computer search of the database using as probes C.elegans introns, bending elements, and a 20 nucleotide consensus sequence for bent DNA. This approach was also used to provide additional examples of species-specific sequence patterns in eukaryotes where it was shown that (A) greater than or equal to 10 and (A.T) greater than or equal to 5 tracts are prevalent throughout the untranslated DNA of D.discodium and P.falciparum, respectively. These results provide new insight into the organization of eukaryotic DNA because they show that species-specific patterns of simple sequences are found in introns and in other untranslated regions of the genome.     

DNA Repair Systems in Archaea: Mementos from the Last Universal Common Ancestor?

DNA repair in the Archaea is relevant to the consideration of genome maintenance and replication fidelity in the last universal common ancestor (LUCA) from two perspectives. First, these prokaryotes embody a mix of bacterial and eukaryal molecular features. Second, DNA repair proteins would have been essential in LUCA to maintain genome integrity, regardless of the environmental temperature. Yet we know very little of the basic molecular mechanisms of DNA damage and repair in the Archaea in general. Many studies on DNA repair in archaea have been conducted with hyperthermophiles because of the additional stress imposed on their macromolecules by high temperatures. In addition, of the six complete archaeal genome sequences published so far, five are thermophilic archaea. We have recently shown that the hyperthermophile Pyrococcus furiosus has an extraordinarily high capacity for repair of radiation-induced double-strand breaks and we have identified and sequenced several genes involved in DNA repair in P. furiosus. At the sequence level, only a few genes share homology with known bacterial repair genes. For instance, our phylogenetic analysis indicates that archaeal recombinases occur in two paralogous gene families, one of which is very deeply branched, and both recombinases are more closely related to the eukaryotic RAD51 and Dmc1 gene families than to the Escherichia coli recA gene. We have also identified a gene encoding a repair endo/exonuclease in the genomes of several Archaea. The archaeal sequences are highly homologous to those of the eukaryotic Rad2 family and they cluster with genes of the FEN-1 subfamily, which are known to be involved in DNA replication and repair in eukaryotes. We argue that there is a commonality of mechanisms and protein sequences, shared between prokaryotes and eukaryotes for several modes of DNA repair, reflecting diversification from a minimal set of genes thought to represent the genome of the LUCA.     

Horizontal transfer of a bacterial gene involved in polyglutamate biosynthesis to the plant-parasitic nematode Meloidogyne artiellia.

Analysis of a genomic fragment from the plant parasitic nematode Meloidogyne artiellia revealed the presence of a gene which, in bacteria, is involved in the formation of polyglutamate capsule. Searching of various databases, including the Caenorhabditis elegans genome sequence and the large EST datasets from a variety of parasitic nematodes, showed that no similar genes have been identified in other nematodes or in any other eukaryotic organisms. The M. artiellia gene has a typical eukaryotic structure and its mRNA is present in the intestine. The gene is expressed in all life cycle stages tested. These findings demonstrate horizontal gene transfer may be important in catalyzing the diversification of nematode lineages.     

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.     

Identifying bacterial genes and endosymbiont DNA with Glimmer

MOTIVATION: The Glimmer gene-finding software has been successfully used for finding genes in bacteria, archaea and viruses representing hundreds of species. We describe several major changes to the Glimmer system, including improved methods for identifying both coding regions and start codons. We also describe a new module of Glimmer that can distinguish host and endosymbiont DNA. This module was developed in response to the discovery that eukaryotic genome sequencing projects sometimes inadvertently capture the DNA of intracellular bacteria living in the host. RESULTS: The new methods dramatically reduce the rate of false-positive predictions, while maintaining Glimmer's 99% sensitivity rate at detecting genes in most species, and they find substantially more correct start sites, as measured by comparisons to known and well-curated genes. We show that our interpolated Markov model (IMM) DNA discriminator correctly separated 99% of the sequences in a recent genome project that produced a mixture of sequences from the bacterium Prochloron didemni and its sea squirt host, Lissoclinum patella. AVAILABILITY: Glimmer is OSI Certified Open Source and available at http://cbcb.umd.edu/software/glimmer.     

DNA methylation: evolution of a bacterial immune function into a regulator of gene expression and genome structure in higher eukaryotes

The amino acid sequence of mammalian DNA methyltransferase has been deduced from the nucleotide sequence of a cloned cDNA. It appears that the mammalian enzyme arose during evolution via fusion of a prokaryotic restriction methyltransferase gene and a second gene of unknown function. Mammalian DNA methyltransferase currently comprises an N-terminal domain of about 1000 amino acids that may have a regulatory role and a C-terminal 570 amino acid domain that retains similarities to bacterial restriction methyltransferases. The sequence similarities among mammalian and bacterial DNA cytosine methyltransferases suggest a common evolutionary origin. DNA methylation is uncommon among those eukaryotes having genomes of less than 10(8) base pairs, but nearly universal among large-genome eukaryotes. This and other considerations make it likely that sequence inactivation by DNA methylation has evolved to compensate for the expansion of the genome that has accompanied the development of higher plants and animals. As methylated sequences are usually propagated in the repressed, nuclease-insensitive state, it is likely that DNA methylation compartmentalizes the genome to facilitate gene regulation by reducing the total amount of DNA sequence that must be scanned by DNA-binding regulatory proteins. DNA methylation is involved in immune recognition in bacteria but appears to regulate the structure and expression of the genome in complex higher eukaryotes. I suggest that the DNA-methylating system of mammals was derived from that of bacteria by way of a hypothetical intermediate that carried out selective de novo methylation of exogenous DNA and propagated the methylated DNA in the repressed state within its own genome.(ABSTRACT TRUNCATED AT 250 WORDS)     

A naturally occurring horizontal gene transfer from a eukaryote to a prokaryote

Summary Naturally occurring horizontal gene transfers between nonviral organisms are difficult to prove. Only with the availability of sequence data from a wide variety of organisms can a convincing case be made. In the case of putative gene transfers between prokaryotes and eukaryotes, the minimum requirements for inferring such an event include (1) sequences of the transferred gene or its product from several appropriately divergent eukaryotes and several prokaryotes, and (2) a similar set of sequences from the same (or closely related organisms) for another gene or genes. Given these criteria, we believe that a strong case can be made forEscherichia coli having acquired a second glyceraldehyde-3-phosphate dehydrogenase gene from some eukaryotic host. Ancillary observations on the general rate of change and the time of the prokaryote-eukaryote divergence support the notion.     

Evidence of independent gene duplications during the evolution of archaeal and eukaryotic family B DNA polymerases

Eukaryotes and archaea both possess multiple genes coding for family B DNA polymerases. In animals and fungi, three family B DNA polymerases, alpha, delta, and epsilon, are responsible for replication of nuclear DNA. We used a PCR-based approach to amplify and sequence phylogenetically conserved regions of these three DNA polymerases from Giardia intestinalis and Trichomonas vaginalis, representatives of early-diverging eukaryotic lineages. Phylogenetic analysis of eukaryotic and archaeal paralogs suggests that the gene duplications that gave rise to the three replicative paralogs occurred before the divergence of the earliest eukaryotic lineages, and that all eukaryotes are likely to possess these paralogs. One eukaryotic paralog, epsilon, consistently branches within archaeal sequences to the exclusion of other eukaryotic paralogs, suggesting that an epsilon-like family B DNA polymerase was ancestral to both archaea and eukaryotes. Because crenarchaeote and euryarchaeote paralogs do not form monophyletic groups in phylogenetic analysis, it is possible that archaeal family B paralogs themselves evolved by a series of gene duplications independent of the gene duplications that gave rise to eukaryotic paralogs.      

Few mitochondrial DNA sequences are inserted into the turkey (Meleagris gallopavo) nuclear genome: evolutionary analyses and informativity in the domestic lineage

Mitochondrial DNA (mtDNA) insertions have been detected in the nuclear genome of many eukaryotes. These sequences are pseudogenes originated by horizontal transfer of mtDNA fragments into the nuclear genome, producing nuclear DNA sequences of mitochondrial origin (numt). In this study we determined the frequency and distribution of mtDNA‐originated pseudogenes in the turkey (Meleagris gallopavo) nuclear genome. The turkey reference genome (Turkey_2.01) was aligned with the reference linearized mtDNA sequence using last . A total of 32 numt sequences (corresponding to 18 numt regions derived by unique insertional events) were identified in the turkey nuclear genome (size ranging from 66 to 1415 bp; identity against the modern turkey mtDNA corresponding region ranging from 62% to 100%). Numts were distributed in nine chromosomes and in one scaffold. They derived from parts of 10 mtDNA protein‐coding genes, ribosomal genes, the control region and 10 tRNA genes. Seven numt regions reported in the turkey genome were identified in orthologues positions in the Gallus gallus genome and therefore were present in the ancestral genome that in the Cretaceous originated the lineages of the modern crown Galliformes. Five recently integrated turkey numts were validated by PCR in 168 turkeys of six different domestic populations. None of the analysed numts were polymorphic (i.e. absence of the inserted sequence, as reported in numts of recent integration in other species), suggesting that the reticulate speciation model is not useful for explaining the origin of the domesticated turkey lineage.     

Gene identification in novel eukaryotic genomes by self-training algorithm

Finding new protein-coding genes is one of the most important goals of eukaryotic genome sequencing projects. However, genomic organization of novel eukaryotic genomes is diverse and ab initio gene finding tools tuned up for previously studied species are rarely suitable for efficacious gene hunting in DNA sequences of a new genome. Gene identification methods based on cDNA and expressed sequence tag (EST) mapping to genomic DNA or those using alignments to closely related genomes rely either on existence of abundant cDNA and EST data and/or availability on reference genomes. Conventional statistical ab initio methods require large training sets of validated genes for estimating gene model parameters. In practice, neither one of these types of data may be available in sufficient amount until rather late stages of the novel genome sequencing. Nevertheless, we have shown that gene finding in eukaryotic genomes could be carried out in parallel with statistical models estimation directly from yet anonymous genomic DNA. The suggested method of parallelization of gene prediction with the model parameters estimation follows the path of the iterative Viterbi training. Rounds of genomic sequence labeling into coding and non-coding regions are followed by the rounds of model parameters estimation. Several dynamically changing restrictions on the possible range of model parameters are added to filter out fluctuations in the initial steps of the algorithm that could redirect the iteration process away from the biologically relevant point in parameter space. Tests on well-studied eukaryotic genomes have shown that the new method performs comparably or better than conventional methods where the supervised model training precedes the gene prediction step. Several novel genomes have been analyzed and biologically interesting findings are discussed. Thus, a self-training algorithm that had been assumed feasible only for prokaryotic genomes has now been developed for ab initio eukaryotic gene identification.     

Genome Sequence Analysis Indicates that the Model Eukaryote Nematostella vectensis Harbors Bacterial Consorts

Analysis of the genome sequence of the starlet sea anemone, Nematostella vectensis, reveals many genes whose products are phylogenetically closer to proteins encoded by bacteria or bacteriophages than to any metazoan homologs. One explanation for such sequence affinities could be that these genes have been horizontally transferred from bacteria to the Nematostella lineage. We show, however, that bacterium-like and phage-like genes sequenced by the N. vectensis genome project tend to cluster on separate scaffolds, which typically do not include eukaryotic genes and differ from the latter in their GC contents. Moreover, most of the bacterium-like genes in N. vectensis either lack introns or the introns annotated in such genes are false predictions that, when translated, often restore the missing portions of their predicted protein products. In a freshwater cnidarian, Hydra, for which a proteobacterial endosymbiont is known, these gene features have been used to delineate the DNA of that endosymbiont sampled by the genome sequencing project. We predict that a large fraction of bacterium-like genes identified in the N. vectensis genome similarly are drawn from the contemporary bacterial consorts of the starlet sea anemone. These uncharacterized bacteria associated with N. vectensis are a proteobacterium and a representative of the phylum Bacteroidetes, each represented in the database by an apparently random sample of informational and operational genes. A substantial portion of a putative bacteriophage genome was also detected, which would be especially unlikely to have been transferred to a eukaryote.     

Eukaryotic origin of a metabolic pathway in virus by horizontal gene transfer

Horizontal gene transfer, the movement of genetic materials across the normal mating barriers between organisms occurs frequently and contributes significantly to the evolution of both eukaryotic and prokaryotic genomes. However, few concurrent transfers of functionally related genes implemented in a pathway from eukaryotes to prokaryotes are observed. Here, we did phylogenetic analyses to support that the genes, i.e. dihydrofolate reductase, glycine hydroxymethyltransferase, and thymidylate synthase involved in thymidylate metabolism, in Hz-1 virus were obtained from insect genome recently by independent horizontal gene transfers. In addition, five other related genes in nucleotide metabolism show evidences of horizontal gene transfers. These genes demonstrate similar expression pattern, and they may have formatted a functionally related pathway (e.g. thymidylate synthesis, and DNA replication) in Hz-1 virus. In conclusion, we provide an example of horizontal gene transfer of functionally related genes in a pathway to prokaryote from eukaryote.     

Complete Genome Sequence Analysis of Nocardia brasiliensis HUJEG-1 Reveals a Saprobic Lifestyle and the Genes Needed for Human Pathogenesis

Nocardia brasiliensis is an important etiologic agent of mycetoma. These bacteria live as a saprobe in soil or organic material and enter the tissue via minor trauma. Mycetoma is characterized by tumefaction and the production of fistula and abscesses, with no spontaneous cure. By using mass sequencing, we determined the complete genomic nucleotide sequence of the bacteria. According to our data, the genome is a circular chromosome 9,436,348-bp long with 68% G+C content that encodes 8,414 proteins. We observed orthologs for virulence factors, a higher number of genes involved in lipid biosynthesis and catabolism, and gene clusters for the synthesis of bioactive compounds, such as antibiotics, terpenes, and polyketides. An in silico analysis of the sequence supports the conclusion that the bacteria acquired diverse genes by horizontal transfer from other soil bacteria, even from eukaryotic organisms. The genome composition reflects the evolution of bacteria via the acquisition of a large amount of DNA, which allows it to survive in new ecological niches, including humans.     

Bioinformatics Analysis of Bacterial Annexins – Putative Ancestral Relatives of Eukaryotic Annexins

Annexins are Ca²⁺-binding, membrane-interacting proteins, widespread among eukaryotes, consisting usually of four structurally similar repeated domains. It is accepted that vertebrate annexins derive from a double genome duplication event. It has been postulated that a single domain annexin, if found, might represent a molecule related to the hypothetical ancestral annexin. The recent discovery of a single-domain annexin in a bacterium, Cytophaga hutchinsonii, apparently confirmed this hypothesis. Here, we present a more complex picture. Using remote sequence similarity detection tools, a survey of bacterial genomes was performed in search of annexin-like proteins. In total, we identified about thirty annexin homologues, including single-domain and multi-domain annexins, in seventeen bacterial species. The thorough search yielded, besides the known annexin homologue from C. hutchinsonii, homologues from the Bacteroidetes/Chlorobi phylum, from Gemmatimonadetes, from beta- and delta-Proteobacteria, and from Actinobacteria. The sequences of bacterial annexins exhibited remote but statistically significant similarity to sequence profiles built of the eukaryotic ones. Some bacterial annexins are equipped with additional, different domains, for example those characteristic for toxins. The variation in bacterial annexin sequences, much wider than that observed in eukaryotes, and different domain architectures suggest that annexins found in bacteria may actually descend from an ancestral bacterial annexin, from which eukaryotic annexins also originate. The hypothesis of an ancient origin of bacterial annexins has to be reconciled with the fact that remarkably few bacterial strains possess annexin genes compared to the thousands of known bacterial genomes and with the patchy, anomalous phylogenetic distribution of bacterial annexins. Thus, a massive annexin gene loss in several bacterial lineages or very divergent evolution would appear a likely explanation. Alternative evolutionary scenarios, involving horizontal gene transfer between bacteria and protozoan eukaryotes, in either direction, appear much less likely. Altogether, current evidence does not allow unequivocal judgement as to the origin of bacterial annexins.     

Isolating the Arabidopsis thaliana genes for de novo purine synthesis by suppression of Escherichia coli mutants. I. 5'-Phosphoribosyl-5-aminoimidazole synthetase.

We have initiated an investigation of the de novo purine nucleotide biosynthetic pathway in the plant Arabidopsis thaliana. Functional suppression of Escherichia coli auxotrophs allowed the direct isolation of expressed Arabidopsis leaf cDNAs. Using this approach we have successfully suppressed mutants in 4 of the 12 genes in this pathway. One of these cDNA clones, encoding 5'-phosphoribosyl-5-aminoimidazole (AIR) synthetase (PUR5) has been characterized in detail. Analysis of genomic DNA suggests that the Arabidopsis genome contains a single AIR synthetase gene. Analysis of the cDNA sequence and mRNA size suggests that this enzyme activity is encoded by a monofunctional polypeptide, similar to that of bacteria and unlike other eukaryotes. The Arabidopsis AIR synthetase contains a basic hydrophobic transit peptide consistent with transport into chloroplasts. Comparison of both the predicted amino acid and nucleotide sequence from Arabidopsis to those of eight other distant organisms suggests that the plant sequence is more similar to the bacterial sequences than to other eukaryotic sequences. This study provides the groundwork for future investigations into the regulation of de novo purine biosynthesis in plants. Additionally, we have demonstrated that functional suppression of bacterial mutants may provide a useful method for cloning a variety of plant genes.