Mechanisms of evolution in Rickettsia conorii and R. prowazekii

Phylogenomics of the reproductive parasite Wolbachia pipientis wMel: a streamlined genome overrun by mobile genetic elements

The complete sequence of the 1,267,782 bp genome of Wolbachia pipientis wMel, an obligate intracellular bacteria of Drosophila melanogaster, has been determined. Wolbachia, which are found in a variety of invertebrate species, are of great interest due to their diverse interactions with different hosts, which range from many forms of reproductive parasitism to mutualistic symbioses. Analysis of the wMel genome, in particular phylogenomic comparisons with other intracellular bacteria, has revealed many insights into the biology and evolution of wMel and Wolbachia in general. For example, the wMel genome is unique among sequenced obligate intracellular species in both being highly streamlined and containing very high levels of repetitive DNA and mobile DNA elements. This observation, coupled with multiple evolutionary reconstructions, suggests that natural selection is somewhat inefficient in wMel, most likely owing to the occurrence of repeated population bottlenecks. Genome analysis predicts many metabolic differences with the closely related Rickettsia species, including the presence of intact glycolysis and purine synthesis, which may compensate for an inability to obtain ATP directly from its host, as Rickettsia can. Other discoveries include the apparent inability of wMel to synthesize lipopolysaccharide and the presence of the most genes encoding proteins with ankyrin repeat domains of any prokaryotic genome yet sequenced. Despite the ability of wMel to infect the germline of its host, we find no evidence for either recent lateral gene transfer between wMel and D. melanogaster or older transfers between Wolbachia and any host. Evolutionary analysis further supports the hypothesis that mitochondria share a common ancestor with the α-Proteobacteria, but shows little support for the grouping of mitochondria with species in the order Rickettsiales. With the availability of the complete genomes of both species and excellent genetic tools for the host, the wMel–D. melanogaster symbiosis is now an ideal system for studying the biology and evolution of Wolbachia infections.

The genome sequence of Wolbachia provides insights into the origins of mitochondria, as well as the ecology and evolution of endosymbiosis

Gene sequence-based criteria for identification of new rickettsia isolates and description of Rickettsia heilongjiangensis sp. nov

We propose genetic guidelines for the classification of rickettsial isolates at the genus, group, and species levels by using sequences of the 16S rRNA (rrs) gene and four protein-coding genes, the gltA, ompA, and ompB genes and gene D. To be classified as a member of the genus Rickettsia, an isolate should exhibit degrees of rrs and gltA homology with any of the 20 Rickettsia species studied of >/=98.1 and >/=86.5%, respectively. A member of the typhus group should fulfill at least two of the following four criteria: pairwise nucleotide sequence homologies with rrs, gltA, ompB, and gene D of either Rickettsia typhi or Rickettsia prowazekii of >/=99.4, >/=96.6, >/=92.4, and >/=91.6%, respectively. A member of the spotted fever group should either possess the ompA gene or fulfill at least two of the following four criteria: pairwise nucleotide sequence homologies with rrs, gltA, ompB, and gene D of any member of this group of >/=98.8, >/=92.7, >/=85.8, and >/=82.2%, respectively. The existence of a distinct "ancestral" group should be questioned. To be classified as a new Rickettsia species, an isolate should not exhibit more than one of the following degrees of nucleotide similarity with the most homologous validated species: >/=99.8 and >/=99.9% for the rrs and gltA genes, respectively, and, when amplifiable, >/=98.8, >/=99.2, and >/=99.3% for the ompA and ompB genes and gene D, respectively. By use of our classification scheme, "Rickettsia heilongjiangii" belongs to a new species for which we officially propose the name Rickettsia heilongjiangensis sp. nov.

Protein interaction mapping on a functional shotgun sequence of Rickettsia sibirica

Protein interaction maps can reveal novel pathways and functional complexes, allowing 'guilt by association' annotation of uncharacterized proteins. To address the need for large-scale protein interaction analyses, a bacterial two-hybrid system was coupled with a whole genome shotgun sequencing approach for microbial genome analysis. We report the first large-scale proteomics study using this system, integrating de novo genome sequencing with functional interaction mapping and annotation in a high-throughput format. We apply the approach by shotgun sequencing and annotating the genome of Rickettsia sibirica strain 246, an obligate intracellular human pathogen among the Spotted Fever Group rickettsiae. The bacteria invade endothelial cells and cause lysis after large amounts of progeny have accumulated. Little is known about specific Rickettsial virulence factors and their mode of pathogenicity. Analysis of the combined genomic sequence and protein-protein interaction data for a set of virulence related Type IV secretion system (T4SS) proteins revealed over 250 interactions and will provide insight into the mechanism of Rickettsial pathogenicity.

Comparative genomics of Rickettsia prowazekii Madrid E and Breinl strains

Rickettsia prowazekii, the causative agent of epidemic typhus, has been responsible for millions of human deaths. Madrid E is an attenuated strain of R. prowazekii, while Breinl is a virulent strain. The genomic DNA sequence of Madrid E has recently been published. To study the genomic variations between Madrid E (reference) and Breinl (test) DNAs, cohybridization experiments were performed on a DNA microarray containing all 834 protein-coding genes of Madrid E. Of the 834 genes assessed, 24 genes showed 1.5- to 2.0-fold increases in hybridization signals in Breinl DNA compared to Madrid E DNA, indicating the presence of genomic variations in approximately 3% of the total genes. Eighteen of these 24 genes are predicted to be involved in different functions. Southern blot analysis of five genes, virB4, ftsK, rfbE, lpxA, and rpoH, suggested the presence of an additional paralog(s) in Breinl, which might be related to the observed increase in hybridization signals. Studies by real-time reverse transcription-PCR revealed an increase in expression of the above-mentioned five genes and five other genes. In addition to the elevated hybridization signals of 24 genes observed in the Breinl strain, one gene (rp084) showed only 1/10 the hybridization signal of Madrid E. Further analysis of this gene by PCR and sequencing revealed a large deletion flanking the whole rp084 gene and part of the rp083 gene in the virulent Breinl strain. The results of this first rickettsial DNA microarray may provide some important information for the elucidation of pathogenic mechanisms of R. prowazekii.

Millions of years of evolution preserved: a comprehensive catalog of the processed pseudogenes in the human genome

Processed pseudogenes were created by reverse-transcription of mRNAs; they provide snapshots of ancient genes existing millions of years ago in the genome. To find them in the present-day human, we developed a pipeline using features such as intron-absence, frame-disruption, polyadenylation, and truncation. This has enabled us to identify in recent genome drafts approximately 8000 processed pseudogenes (distributed from http://pseudogene.org). Overall, processed pseudogenes are very similar to their closest corresponding human gene, being 94% complete in coding regions, with sequence similarity of 75% for amino acids and 86% for nucleotides. Their chromosomal distribution appears random and dispersed, with the numbers on chromosomes proportional to length, suggesting sustained "bombardment" over evolution. However, it does vary with GC-content: Processed pseudogenes occur mostly in intermediate GC-content regions. This is similar to Alus but contrasts with functional genes and L1-repeats. Pseudogenes, moreover, have age profiles similar to Alus. The number of pseudogenes associated with a given gene follows a power-law relationship, with a few genes giving rise to many pseudogenes and most giving rise to few. The prevalence of processed pseudogenes agrees well with germ-line gene expression. Highly expressed ribosomal proteins account for approximately 20% of the total. Other notables include cyclophilin-A, keratin, GAPDH, and cytochrome c.

The RickA protein of Rickettsia conorii activates the Arp2/3 complex

Actin polymerization, the main driving force for cell locomotion, is also used by the bacteria Listeria and Shigella and vaccinia virus for intracellular and intercellular movements. Seminal studies have shown the key function of the Arp2/3 complex in nucleating actin and generating a branched array of actin filaments during membrane extension and pathogen movement. Arp2/3 requires activation by proteins such as the WASP-family proteins or ActA of Listeria. We previously reported that actin tails of Rickettsia conorii, another intracellular bacterium, unlike those of Listeria, Shigella or vaccinia, are made of long unbranched actin filaments apparently devoid of Arp2/3 (ref. 4). Here we identify a R. conorii surface protein, RickA, that activates Arp2/3 in vitro, although less efficiently than ActA. In infected cells, Arp2/3 is detected on the rickettsial surface but not in actin tails. When expressed in mammalian cells and targeted to the membrane, RickA induces filopodia. Thus RickA-induced actin polymerization, by generating long actin filaments reminiscent of those present in filopodia, has potential as a tool for studying filopodia formation.

Scrub typhus and tropical rickettsioses

PURPOSE OF REVIEW

Recent developments in molecular taxonomic methods have led to a reclassification of rickettsial diseases. The agent responsible for scrub typhus (Orientia tsutsugamushi ) has been removed from the genus Rickettsia and a bewildering array of new rickettsial pathogens have been described. An update of recent research findings is therefore particularly timely for the nonspecialist physician.

RECENT FINDINGS

An estimated one billion people are at risk for scrub typhus and an estimated one million cases occur annually. The disease appears to be re-emerging in Japan, with seasonal transmission. O. tsutsugamushi has evolved a variety of mechanisms to remain viable in its intracellular habitat. Slowing the release of intracellular calcium inhibits apoptosis of macrophages. Subsets of chemokine genes are induced in infected cells, some in response to transcription factor activator protein 1. Cardiac involvement is uncommon and clinical complications are predominantly pulmonary. Serious pneumonitis occurred in 22% of Chinese patients. Dual infections with leptospirosis have been reported. Standardized diagnostic tests are being developed and attempts to improve treatment of women and children are being made. Of the numerous tick-borne rickettsioses identified in recent years, African tick-bite fever appears to be of particular importance to travellers. The newly described flea-borne spotted fever caused by Rickettsia felis may be global in distribution.

SUMMARY

Rash and fever in a returning traveler could be rickettsial and presumptive doxycycline treatment can be curative. Recent research findings raise more questions than answers and should stimulate much needed research.

The genome of Nanoarchaeum equitans: insights into early archaeal evolution and derived parasitism

The hyperthermophile Nanoarchaeum equitans is an obligate symbiont growing in coculture with the crenarchaeon Ignicoccus. Ribosomal protein and rRNA-based phylogenies place its branching point early in the archaeal lineage, representing the new archaeal kingdom Nanoarchaeota. The N. equitans genome (490,885 base pairs) encodes the machinery for information processing and repair, but lacks genes for lipid, cofactor, amino acid, or nucleotide biosyntheses. It is the smallest microbial genome sequenced to date, and also one of the most compact, with 95% of the DNA predicted to encode proteins or stable RNAs. Its limited biosynthetic and catabolic capacity indicates that N. equitans' symbiotic relationship to Ignicoccus is parasitic, making it the only known archaeal parasite. Unlike the small genomes of bacterial parasites that are undergoing reductive evolution, N. equitans has few pseudogenes or extensive regions of noncoding DNA. This organism represents a basal archaeal lineage and has a highly reduced genome.

Tracing the evolution of gene loss in obligate bacterial symbionts

The gamma-proteobacterial symbionts of insects are a model group for comparative studies of genome reduction. The phylogenetic proximity of these reduced genomes to the larger genomes of well-studied free-living bacteria has enabled reconstructions of the process by which genes and DNA are lost. Three genome sequences are now available for Buchnera aphidicola. Analyses of Buchnera genomes in comparison with those of related enteric bacteria suggest that extensive changes including large deletions, repetitive element proliferation and chromosomal rearrangements occurred initially, followed by extreme stasis in gene order and slow decay of additional genes. This pattern appears to be characteristic of symbiont evolution.

History of the ADP/ATP-translocase-encoding gene, a parasitism gene transferred from a Chlamydiales ancestor to plants 1 billion years ago

Nonmitochondrial ADP/ATP translocase is an energy parasite enzyme. Its encoding gene, tlc, is found only in Rickettsiales, Chlamydiales, and plant and alga plastids. We demonstrate the presence of tlc in Parachlamydia acanthamoebae. This gene shares more similarity with the tlc1 gene of Chlamydiaceae and the tlc of plant and alga plastids than with the tlc2 gene of Chlamydiaceae. Phylogenetic analysis, including all other tlc homologs found in GenBank, showed that tlc was duplicated in a Chlamydiales ancestor before the appearance of multicellular eukaryotes. A time scale, calibrated with seven independent time points obtained from fossil estimates and from the 16S rRNA molecular clock, was congruent with the molecular clock provided by tlc. Plant and alga plastids acquired tlc approximately when Parachlamydiaceae and Chlamydiaceae diverged, at the eucaryotic radiation time, ca. 1 billion years ago.

Ehrlichia chaffeensis and Anaplasma phagocytophilum lack genes for lipid A biosynthesis and incorporate cholesterol for their survival

Ehrlichia chaffeensis and Anaplasma phagocytophilum are agents of human monocytic and granulocytic ehrlichioses, respectively. They are extremely sensitive to mechanical stress and are pleomorphic gram-negative bacteria. Membrane incorporation of cholesterol from the eukaryotic host is known to be essential for other fragile and pleomorphic bacteria and mycoplasmas that lack a cell wall. Thus, we tested whether cholesterol is required for E. chaffeensis and A. phagocytophilum. Using a freeze fracture technique and biochemical analysis, these bacteria were found to contain significant levels of membrane cholesterol. These bacteria lack genes for cholesterol biosynthesis or modification. However, host cell-free bacteria had the ability to take up directly exogenous cholesterol or NBD-cholesterol, a fluorescent cholesterol derivative. Treatment of the bacteria with cholesterol extraction reagent methyl-beta-cyclodextrin caused their ultrastructural changes. Furthermore, pretreatment of the bacteria with methyl-beta-cyclodextrin or NBD-cholesterol deprived these bacteria of the ability to infect leukocytes, thus killing these obligate intracellular bacteria. Analysis of E. chaffeensis and A. phagocytophilum genome sequences revealed that these bacteria lack all genes for the biosynthesis of lipid A and most genes for the biosynthesis of peptidoglycan, which confer structural strength to gram-negative bacteria. Taken together, these results suggest that human ehrlichiosis agents became cholesterol dependent due to the loss of these genes. As the first report of gram-negative bacteria incorporating cholesterol for survival, these findings offer insight into the unique nature of their parasitism and imply that cholesterol is important in the control of human ehrlichioses.

Molecular and functional analysis of the lepB gene, encoding a type I signal peptidase from Rickettsia rickettsii and Rickettsia typhi

The type I signal peptidase lepB genes from Rickettsia rickettsii and Rickettsia typhi, the etiologic agents of Rocky Mountain spotted fever and murine typhus, respectively, were cloned and characterized. Sequence analysis of the cloned lepB genes from R. rickettsii and R. typhi shows open reading frames of 801 and 795 nucleotides, respectively. Alignment analysis of the deduced amino acid sequences reveals the presence of highly conserved motifs that are important for the catalytic activity of bacterial type I signal peptidase. Reverse transcription-PCR and Northern blot analysis demonstrated that the lepB gene of R. rickettsii is cotranscribed in a polycistronic message with the putative nuoF (encoding NADH dehydrogenase I chain F), secF (encoding protein export membrane protein), and rnc (encoding RNase III) genes in a secF-nuoF-lepB-rnc cluster. The cloned lepB genes from R. rickettsii and R. typhi have been demonstrated to possess signal peptidase I activity in Escherichia coli preprotein processing in vivo by complementation assay.

Tropheryma whipplei Twist: a human pathogenic Actinobacteria with a reduced genome

The human pathogen Tropheryma whipplei is the only known reduced genome species (<1 Mb) within the Actinobacteria [high G+C Gram-positive bacteria]. We present the sequence of the 927303-bp circular genome of T. whipplei Twist strain, encoding 808 predicted protein-coding genes. Specific genome features include deficiencies in amino acid metabolisms, the lack of clear thioredoxin and thioredoxin reductase homologs, and a mutation in DNA gyrase predicting a resistance to quinolone antibiotics. Moreover, the alignment of the two available T. whipplei genome sequences (Twist vs. TW08/27) revealed a large chromosomal inversion the extremities of which are located within two paralogous genes. These genes belong to a large cell-surface protein family defined by the presence of a common repeat highly conserved at the nucleotide level. The repeats appear to trigger frequent genome rearrangements in T. whipplei, potentially resulting in the expression of different subsets of cell surface proteins. This might represent a new mechanism for evading host defenses. The T. whipplei genome sequence was also compared to other reduced bacterial genomes to examine the generality of previously detected features. The analysis of the genome sequence of this previously largely unknown human pathogen is now guiding the development of molecular diagnostic tools and more convenient culture conditions.

Molecular evolution of the modulator of chloroplast ATP synthase: origin of the conformational change dependent regulation

Chloroplast ATP synthase synthesizes ATP by utilizing a proton gradient as an energy supply, which is generated by photosynthetic electron transport. The activity of the chloroplast ATP synthase is regulated in several specific ways to avoid futile hydrolysis of ATP under various physiological conditions. Several regulatory signals such as Delta mu H(+), tight binding of ADP and its release, thiol modulation, and inhibition by the intrinsic inhibitory subunit epsilon are sensed by this complex. In this review, we describe the function of two regulatory subunits, gamma and epsilon, of ATP synthase based on their possible conformational changes and discuss the evolutionary origin of these regulation systems.

"Candidatus Rickettsia tarasevichiae" in Ixodes persulcatus ticks collected in Russia

We collected 209 Ixodes persulcatus ticks in various regions of Russia, including the southern Urals and western and eastern Siberia. Using PCR amplification and sequencing of the citrate synthase-encoding gene (gltA), we detected a new rickettsial genotype, which we named "Candidatus Rickettsia tarasevichiae." This bacterium was found in 9.27%, 10.0%, and 20.5% of the ticks collected in western Siberia, eastern Siberia, and the southern Urals, respectively. "Candidatus Rickettsia tarasevichiae" exhibited a 98% and 96% nucleotide sequence homology, with the 16S rDNA and gltA sequence, respectively, of R. canadensis, a rickettsia previously only found in Haemaphysalis leporispalustris ticks in North America. The phylogenetic analysis of "Candidatus Rickettsia tarasevichiae" and other Rickettsia species allowed the creation of a new cluster with high bootstrap values within the Rickettsia genus involving this rickettsia, R. canadensis, and three uncultured rickettsiae from plant insects.

Pathogenic mechanisms of diseases caused by Rickettsia

The specter of bioterrorism employing genetically engineered Rickettsia resistant to all antibiotics should reawaken the world's desire to elucidate the pathogenesis of typhus and spotted fever rickettsioses in a search for mechanisms vulnerable to interdiction. The pathogenetic sequence includes rickettsial entry into the dermis, hematogenous dissemination to vascular endothelial cells (most critically in brain and lungs), increased vascular permeability, edema, and immunity mediated by NK cells, IFN-gamma, TNF-alpha, RANTES, antibodies, and cytotoxic T lymphocytes. Silverman has demonstrated the role of reactive oxygen species (ROS) produced by R. rickettsii-infected endothelial cells in peroxidative damage to cell membranes in vitro, and Heinzen has described actin-based rickettsial intracellular mobility and intercellular spread. At this point the availability of sequences of rickettsial genomes and excellent animal models of rickettsioses have yielded insufficient progress towards the identification of rickettsial virulence factors and knowledge of the importance of injury mediated by ROS, phospholipase A(2), protease(s) or other mechanisms in vivo. Attention to the rickettsiosis-associated procoagulant state led to determination that hemostatic mechanisms largely prevent major hemorrhage without disseminated intravascular coagulation or thrombosis-mediated ischemia. Particularly lacking is knowledge of early events in vivo at the portal of entry in skin (or lung), of the effects of the inoculum medium (arthropod saliva or feces), mediators produced by infected endothelium under conditions of flow and of the contributions in vivo of immune effectors to pathology, of the role of apoptosis in rickettsial infection, and of the endothelial cell alterations that account for increased vascular permeability. The host cell receptor for the Rickettsia ligand and the mechanism of rickettsial escape from the phagosome need to be elucidated.

rRNA intergenic spacer regions for phylogenetic analysis of Rickettsia species

Species of the genus Rickettsia are responsible for several human diseases, namely epidemic typhus, Rocky Mountain spotted fever, and tick-borne typhus transmitted by arthropod vectors. The rrl-rrf intergenic spacer region (rrl-rrf ITS) was sequenced for 12 Rickettsia strains, including R. typhi, 6 untested species, R. aeschlimannii, Bar29, R. helvetica, R. honei, R. massilae, and R. slovaca as well as 5 Portuguese isolates. Phylogenetic trees inferred from rrl-rrf spacer sequences using maximum-parsimony and distance methods provided largely congruent tree topologies, supported by significant bootstrap values, enabling the identification of five distinct rickettsiae clusters.

From genes to proteins: in vitro expression of rickettsial proteins

The availability of the complete genome sequences of several organisms allows the comparative analysis of genomes, a branch of bioinformatics known as genomics. With this approach, much can be learned about the biology of organisms that are difficult to culture, even when few, if any, of their proteins have been isolated and studied directly. We have focused our interest on Rickettsia conorii, an obligate intracellular bacterium responsible for Mediterranean spotted fever, a disease endemic in southern Europe. While bioinformatic annotation of the complete genome of this bacteria has allowed identification of 1,374 genes, a large number of them remain functionally uncharacterized. The final goal of many experiments in molecular biology is to use biological systems to synthesize the protein encoded by the gene being studied. Because three-dimensional structures are more resilient to evolution and change than amino acid sequences, structure determination of some open reading frames should also exhibit structural similarity to previously described protein families. We have thus initiated a systematic expression and structure determination program for the proteins encoded by rickettsial genes of interest. We have cloned different genes of R. conorii by recombinational cloning (GATEWAY), Invitrogen) a method that uses in vitro site-specific recombination to accomplish a directional cloning of PCR products and the subsequent automatic subcloning of the DNA segment into new vector backbones at high efficiency. The constructions in p-Dest17 yielded several clones able to express recombinant proteins with a C-terminal histidine tag. Expression of corresponding proteins was then performed using a cell-free protein expression system (Rapid Translation System, RTS, Roche Diagnostics). The recombinational cloning approach coupled to RTS provides an approach to rapid optimization of protein expression and is very useful to express rickettsial proteins. Moreover, this system is able to overcome some of the limitations encountered with rickettsial proteins highly toxic for E. coli or insect cells.

Prokaryote and eukaryote evolvability

The concept of evolvability covers a broad spectrum of, often contradictory, ideas. At one end of the spectrum it is equivalent to the statement that evolution is possible, at the other end are untestable post hoc explanations, such as the suggestion that current evolutionary theory cannot explain the evolution of evolvability. We examine similarities and differences in eukaryote and prokaryote evolvability, and look for explanations that are compatible with a wide range of observations. Differences in genome organisation between eukaryotes and prokaryotes meets this criterion. The single origin of replication in prokaryote chromosomes (versus multiple origins in eukaryotes) accounts for many differences because the time to replicate a prokaryote genome limits its size (and the accumulation of junk DNA). Both prokaryotes and eukaryotes appear to switch from genetic stability to genetic change in response to stress. We examine a range of stress responses, and discuss how these impact on evolvability, particularly in unicellular organisms versus complex multicellular ones. Evolvability is also limited by environmental interactions (including competition) and we describe a model that places limits on potential evolvability. Examples are given of its application to predator competition and limits to lateral gene transfer. We suggest that unicellular organisms evolve largely through a process of metabolic change, resulting in biochemical diversity. Multicellular organisms evolve largely through morphological changes, not through extensive changes to cellular biochemistry.

S-adenosylmethionine transport in Rickettsia prowazekii

Rickettsia prowazekii, the causative agent of epidemic typhus, is an obligate, intracellular, parasitic bacterium that grows within the cytoplasm of eucaryotic host cells. Rickettsiae exploit this intracellular environment by using transport systems for the compounds available in the host cell's cytoplasm. Analysis of the R. prowazekii Madrid E genome sequence revealed the presence of a mutation in the rickettsial metK gene, the gene encoding the enzyme responsible for the synthesis of S-adenosylmethionine (AdoMet). Since AdoMet is required for rickettsial processes, the apparent inability of this strain to synthesize AdoMet suggested the presence of a rickettsial AdoMet transporter. We have confirmed the presence of an AdoMet transporter in the rickettsiae which, to our knowledge, is the first bacterial AdoMet transporter identified. The influx of AdoMet into rickettsiae was a saturable process with a K(T) of 2.3 micro M. Transport was inhibited by S-adenosylethionine and S-adenosylhomocysteine but not by sinfungin or methionine. Transport was also inhibited by 2,4-dinitrophenol, suggesting an energy-linked transport mechanism, and by N-ethylmaleimide. AdoMet transporters with similar properties were also identified in the Breinl strain of R. prowazekii and in Rickettsia typhi. By screening Escherichia coli clone banks for AdoMet transport, the R. prowazekii gene coding for a transporter, RP076 (sam), was identified. AdoMet transport in E. coli containing the R. prowazekii sam gene exhibited kinetics similar to that seen in rickettsiae. The existence of a rickettsial transporter for AdoMet raises intriguing questions concerning the evolutionary relationship between the synthesis and transport of this essential metabolite.

Complete genome sequence of the Q-fever pathogen Coxiella burnetii

The 1,995,275-bp genome of Coxiella burnetii, Nine Mile phase I RSA493, a highly virulent zoonotic pathogen and category B bioterrorism agent, was sequenced by the random shotgun method. This bacterium is an obligate intracellular acidophile that is highly adapted for life within the eukaryotic phagolysosome. Genome analysis revealed many genes with potential roles in adhesion, invasion, intracellular trafficking, host-cell modulation, and detoxification. A previously uncharacterized 13-member family of ankyrin repeat-containing proteins is implicated in the pathogenesis of this organism. Although the lifestyle and parasitic strategies of C. burnetii resemble that of Rickettsiae and Chlamydiae, their genome architectures differ considerably in terms of presence of mobile elements, extent of genome reduction, metabolic capabilities, and transporter profiles. The presence of 83 pseudogenes displays an ongoing process of gene degradation. Unlike other obligate intracellular bacteria, 32 insertion sequences are found dispersed in the chromosome, indicating some plasticity in the C. burnetii genome. These analyses suggest that the obligate intracellular lifestyle of C. burnetii may be a relatively recent innovation.

Mitochondrial connection to the origin of the eukaryotic cell

Phylogenetic evidence is presented that primitively amitochondriate eukaryotes containing the nucleus, cytoskeleton, and endomembrane system may have never existed. Instead, the primary host for the mitochondrial progenitor may have been a chimeric prokaryote, created by fusion between an archaebacterium and a eubacterium, in which eubacterial energy metabolism (glycolysis and fermentation) was retained. A Rickettsia-like intracellular symbiont, suggested to be the last common ancestor of the family Rickettsiaceae and mitochondria, may have penetrated such a host (pro-eukaryote), surrounded by a single membrane, due to tightly membrane-associated phospholipase activity, as do present-day rickettsiae. The relatively rapid evolutionary conversion of the invader into an organelle may have occurred in a safe milieu via numerous, often dramatic, changes involving both partners, which resulted in successful coupling of the host glycolysis and the symbiont respiration. Establishment of a potent energy-generating organelle made it possible, through rapid dramatic changes, to develop genuine eukaryotic elements. Such sequential, or converging, global events could fill the gap between prokaryotes and eukaryotes known as major evolutionary discontinuity.

Characterization of the ftsZ gene from Ehrlichia chaffeensis, Anaplasma phagocytophilum, and Rickettsia rickettsii, and use as a differential PCR target

Degenerate primers corresponding to highly conserved regions of previously characterized ftsZ genes were used to PCR amplify a portion of the ftsZ gene from the genomic DNA of Ehrlichia chaffeensis (ftsZ(Ech)), Anaplasma phagocytophilum (ftsZ(Ap)), and Rickettsia rickettsii (ftsZ(Rr)). Genome walking was then used to amplify the 5' and 3' termini of the genes. The DNA sequences of the resulting amplification products yielded open reading frames coding for proteins with molecular masses of 42.0, 45.7, and 48.3 kDa for A. phagocytophilum, E. chaffeensis, and R. rickettsii, respectively. These homologs are 20 to 70 amino acids longer than the FtsZ proteins characterized in bacteria such as Escherichia coli and Bacillus subtilis, but do not possess the large extended carboxyl-termini found in the FtsZ proteins of Bartonella, Rhizobium, and Agrobacterium species. The functional domains important for FtsZ activity are conserved within the ehrlichial and rickettsial FtsZ protein sequences. The R. rickettsii FtsZ sequence is highly homologous to the FtsZ protein previously described for Rickettsia prowazekii (89% identity), and identical to the FtsZ protein of Rickettsia conorii. The percent identity observed between the A. phagocytophilum and E. chaffeensis FtsZ proteins is only 79% and is particularly low in the carboxyl-terminal region (15.8% identity). Primers were designed to PCR amplify a portion of the variable carboxyl-terminal region of the ftsZ gene, and used to differentiate each agent based on the size of the amplicons: A. phagocytophilum, 278 bp; E. chaffeensis, 341 bp; and Rickettsia spp., 425 bp.

The insertion of palindromic repeats in the evolution of proteins

The current theory of protein evolution is that all contemporary proteins are derived from an ancestral subset. However, each new sequenced genome exhibits many genes with no detectable homologues in other species, leading to the paradoxical picture of a universal ancestor with more genes than any of its progeny. Standard explanations indicate that fast evolving genes might disappear into the 'twilight zone' of sequence similarity. Regardless of the size of the original ancestral subset, its origin and the potential mechanisms of its subsequent enlargement are rarely addressed. Sequencing of Rickettsia conorii genome recently led to the discovery of three families of repeat-mobile elements frequently inserted into the middle of protein coding genes. Although not yet identified in other species of bacteria, this discovery has provided the first clear evidence for the de novo creation of long protein segments (up to 50 amino acid residues) by repeat insertion. Based on previous results and theories on the coding potential of palindromic elements, we speculate that their insertion and mobility might have played a significant role in the early stages of protein evolution.

On the origin of mitochondria: a genomics perspective

The availability of complete genome sequence data from both bacteria and eukaryotes provides information about the contribution of bacterial genes to the origin and evolution of mitochondria. Phylogenetic analyses based on genes located in the mitochondrial genome indicate that these genes originated from within the alpha-proteobacteria. A number of ancestral bacterial genes have also been transferred from the mitochondrial to the nuclear genome, as evidenced by the presence of orthologous genes in the mitochondrial genome in some species and in the nuclear genome of other species. However, a multitude of mitochondrial proteins encoded in the nucleus display no homology to bacterial proteins, indicating that these originated within the eukaryotic cell subsequent to the acquisition of the endosymbiont. An analysis of the expression patterns of yeast nuclear genes coding for mitochondrial proteins has shown that genes predicted to be of eukaryotic origin are mainly translated on polysomes that are free in the cytosol whereas those of putative bacterial origin are translated on polysomes attached to the mitochondrion. The strong relationship with alpha-proteobacterial genes observed for some mitochondrial genes, combined with the lack of such a relationship for others, indicates that the modern mitochondrial proteome is the product of both reductive and expansive processes.

Algorithms for computing parsimonious evolutionary scenarios for genome evolution, the last universal common ancestor and dominance of horizontal gene transfer in the evolution of prokaryotes

BACKGROUND

Comparative analysis of sequenced genomes reveals numerous instances of apparent horizontal gene transfer (HGT), at least in prokaryotes, and indicates that lineage-specific gene loss might have been even more common in evolution. This complicates the notion of a species tree, which needs to be re-interpreted as a prevailing evolutionary trend, rather than the full depiction of evolution, and makes reconstruction of ancestral genomes a non-trivial task.

RESULTS

We addressed the problem of constructing parsimonious scenarios for individual sets of orthologous genes given a species tree. The orthologous sets were taken from the database of Clusters of Orthologous Groups of proteins (COGs). We show that the phyletic patterns (patterns of presence-absence in completely sequenced genomes) of almost 90% of the COGs are inconsistent with the hypothetical species tree. Algorithms were developed to reconcile the phyletic patterns with the species tree by postulating gene loss, COG emergence and HGT (the latter two classes of events were collectively treated as gene gains). We prove that each of these algorithms produces a parsimonious evolutionary scenario, which can be represented as mapping of loss and gain events on the species tree. The distribution of the evolutionary events among the tree nodes substantially depends on the underlying assumptions of the reconciliation algorithm, e.g. whether or not independent gene gains (gain after loss after gain) are permitted. Biological considerations suggest that, on average, gene loss might be a more likely event than gene gain. Therefore different gain penalties were used and the resulting series of reconstructed gene sets for the last universal common ancestor (LUCA) of the extant life forms were analysed. The number of genes in the reconstructed LUCA gene sets grows as the gain penalty increases. However, qualitative examination of the LUCA versions reconstructed with different gain penalties indicates that, even with a gain penalty of 1 (equal weights assigned to a gain and a loss), the set of 572 genes assigned to LUCA might be nearly sufficient to sustain a functioning organism. Under this gain penalty value, the numbers of horizontal gene transfer and gene loss events are nearly identical. This result holds true for two alternative topologies of the species tree and even under random shuffling of the tree. Therefore, the results seem to be compatible with approximately equal likelihoods of HGT and gene loss in the evolution of prokaryotes.

CONCLUSIONS

The notion that gene loss and HGT are major aspects of prokaryotic evolution was supported by quantitative analysis of the mapping of the phyletic patterns of COGs onto a hypothetical species tree. Algorithms were developed for constructing parsimonious evolutionary scenarios, which include gene loss and gain events, for orthologous gene sets, given a species tree. This analysis shows, contrary to expectations, that the number of predicted HGT events that occurred during the evolution of prokaryotes might be approximately the same as the number of gene losses. The approach to the reconstruction of evolutionary scenarios employed here is conservative with regard to the detection of HGT because only patterns of gene presence-absence in sequenced genomes are taken into account. In reality, horizontal transfer might have contributed to the evolution of many other genes also, which makes it a dominant force in prokaryotic evolution.

Transcriptional analysis of Rickettsia prowazekii invasion gene homolog (invA) during host cell infection

An invasion gene homolog, invA, of Rickettsia prowazekii has recently been identified to encode a member of the Nudix hydrolase subfamily which acts specifically on dinucleoside oligophosphates (Np(n)N; n >/= 5), a group of cellular signaling molecules known as alarmones. InvA is thought to enhance intracellular survival by regulating stress-induced toxic nucleotide levels during rickettsial infection. To further characterize the physiological function of InvA, the gene expression pattern during various stages of rickettsial intracellular growth was investigated. Using semiquantitative reverse transcription-PCR (RT-PCR) and real-time fluorescent probe-based quantitative RT-PCR, a differential expression profile of invA during rickettsial host cell infection was examined. The invA transcript temporarily increased during the early period of infection. Expression of rickettsial groEL, a molecular indicator of cellular stresses, was also shown to be upregulated during the early period of infection. Furthermore, invA was cotranscribed in a polycistronic message with rrp, a gene encoding the response regulator protein homolog, which is a part of a two-component signal transduction system. These results support our earlier findings that under such stress conditions dinucleoside oligophosphate pyrophosphatase may function as a buffer, enhancing rickettsial survival within the cytoplasm of a eukaryotic cell. The expression of rickettsial dinucleoside oligophosphate pyrophosphatase may be regulated by a part of the two-component signal transduction system similar to that described for response regulators in other bacterial systems.

Expression of Anaplasma marginale major surface protein 2 operon-associated proteins during mammalian and arthropod infection

The antigenically variant major surface protein 2 (MSP2) of Anaplasma marginale is expressed from a 3.5-kb operon that contains, in a 5'-to-3' direction, four open reading frames, opag3, opag2, opag1, and msp2. This operon structure was shown to be conserved among genotypically and phenotypically distinct A. marginale, A. ovis, and A. centrale strains. The individual OpAG amino acid sequences are highly conserved among A. marginale strains, with identities ranging from 95 to 99%. OpAG2 and OpAG3 were expressed by all examined A. marginale strains during the acute rickettsemia in the mammalian host and, like MSP2, localize to the bacterial surface. OpAG2 and OpAG3 were also expressed in an infected Ixodes scapularis tick cell line. In contrast, the same A. marginale strains expressed only OpAG2 in two different Dermacentor spp. during transmission feeding. OpAG1 expression was not detected in the infected mammalian host, the infected tick cell line, or within infected Dermacentor ticks. The differential expression of outer membrane proteins from within an operon is a novel finding in tick-transmitted bacteria, and the regulation of expression may be broadly applicable to understanding how the pathogen adapts to the mammalian host-tick vector transition.

Conserved response regulator CtrA and IHF binding sites in the alpha-proteobacteria Caulobacter crescentus and Rickettsia prowazekii chromosomal replication origins

CzcR is the Rickettsia prowazekii homolog of the Caulobacter crescentus global response regulator CtrA. CzcR expression partially compensates for developmental defects in ctrA mutant C. crescentus cells, and CzcR binds to all five CtrA binding sites in the C. crescentus replication origin. Conversely, CtrA binds to five similar sites in the putative R. prowazekii replication origin (oriRp). Also, Escherichia coli IHF protein binds over a central CtrA binding site in oriRp. Therefore, CtrA and IHF regulatory proteins have similar binding patterns in both replication origins, and we propose that CzcR is a global cell cycle regulator in R. prowazekii.

Congruent evolution of different classes of non-coding DNA in prokaryotic genomes

Prokaryotic genomes are considered to be 'wall-to-wall' genomes, which consist largely of genes for proteins and structural RNAs, with only a small fraction of the genomic DNA allotted to intergenic regions, which are thought to typically contain regulatory signals. The majority of bacterial and archaeal genomes contain 6-14% non-coding DNA. Significant positive correlations were detected between the fraction of non-coding DNA and inter- and intra-operonic distances, suggesting that different classes of non-coding DNA evolve congruently. In contrast, no correlation was found between any of these characteristics of non-coding sequences and the number of genes or genome size. Thus, the non-coding regions and the gene sets in prokaryotes seem to evolve in different regimes. The evolution of non-coding regions appears to be determined primarily by the selective pressure to minimize the amount of non-functional DNA, while maintaining essential regulatory signals, because of which the content of non-coding DNA in different genomes is relatively uniform and intra- and inter-operonic non-coding regions evolve congruently. In contrast, the gene set is optimized for the particular environmental niche of the given microbe, which results in the lack of correlation between the gene number and the characteristics of non-coding regions.

Molecular and functional analysis of the Rickettsia typhi groESL operon

The groESL operon from an obligate, intracellular, Gram-negative bacterium Rickettsia typhi, the etiologic agent of murine typhus, was cloned and sequenced. The sequence analysis of 2229 bp of the groESL operon reveals two open reading frames of 288 nucleotides (groES) and 1653 nucleotides (groEL) separated by 20 nucleotides. The deduced amino acid sequence of R. typhi GroES and GroEL shows a high degree of identity with other bacterial GroES and GroEL. Reverse transcriptase-polymerase chain reaction and Northern blot analysis indicated that both groES and groEL are transcribed as a single mRNA. The transcriptional start point at 81 nucleotides upstream of the groES start codon was determined by primer extension. The promoter analysis shows no regulatory CIRCE element as it is known for many Gram-positive and Gram-negative bacteria. However, it contains the sequence similar to the putative sigma(70)-dependent promoter and lacks the -35 sequence of the putative sigma(32)-dependent promoter. Complementation assay by R. typhi groESL in a temperature sensitive Escherichia coli groEL mutant restored significant growth ability at non-permissive temperature.

Dictionary-driven protein annotation

Computational methods seeking to automatically determine the properties (functional, structural, physicochemical, etc.) of a protein directly from the sequence have long been the focus of numerous research groups. With the advent of advanced sequencing methods and systems, the number of amino acid sequences that are being deposited in the public databases has been increasing steadily. This has in turn generated a renewed demand for automated approaches that can annotate individual sequences and complete genomes quickly, exhaustively and objectively. In this paper, we present one such approach that is centered around and exploits the Bio-Dictionary, a collection of amino acid patterns that completely covers the natural sequence space and can capture functional and structural signals that have been reused during evolution, within and across protein families. Our annotation approach also makes use of a weighted, position-specific scoring scheme that is unaffected by the over-representation of well-conserved proteins and protein fragments in the databases used. For a given query sequence, the method permits one to determine, in a single pass, the following: local and global similarities between the query and any protein already present in a public database; the likeness of the query to all available archaeal/ bacterial/eukaryotic/viral sequences in the database as a function of amino acid position within the query; the character of secondary structure of the query as a function of amino acid position within the query; the cytoplasmic, transmembrane or extracellular behavior of the query; the nature and position of binding domains, active sites, post-translationally modified sites, signal peptides, etc. In terms of performance, the proposed method is exhaustive, objective and allows for the rapid annotation of individual sequences and full genomes. Annotation examples are presented and discussed in Results, including individual queries and complete genomes that were released publicly after we built the Bio-Dictionary that is used in our experiments. Finally, we have computed the annotations of more than 70 complete genomes and made them available on the World Wide Web at http://cbcsrv.watson.ibm.com/Annotations/.

Characterization of the rRNA genes of Ehrlichia chaffeensis and Anaplasma phagocytophila

The rRNA genes of Ehrlichia chaffeensis and Anaplasma phagocytophila have been analyzed. The 16S rRNA genes were previously characterized for both of these agents. Southern hybridization was used to show that there are single copies of both the 16S and 23S rRNA genes in the genomes of each organism, and that the 16S rRNA genes were upstream from the 23S rRNA genes by at least 16 and 11 Kb for E. chaffeensis and A. phagocytophila, respectively. PCR amplification and gene walking was used to sequence the 23S and 5S rRNA genes, and show that these genes are contiguous and are likely expressed as a single operon. The level of homology between the E. chaffeensis and A. phagocytophila 23S and 5S rRNA genes, and 23S-5S spacers, was 91.8, 81.5, and 40%, respectively. To confirm the hybridization data, genome walking was used to sequence downstream of the 16S rRNA genes, and although no tRNA genes were identified, open reading frames encoding homologues of the Escherichia coli succinate dehydrogenase, subunit C, were found in both E. chaffeensis and A. phagocytophila. Phylogenetic analysis using the 23S rRNA gene suggests that reorganization of the phylum Proteobacteria by division of the class Alphaproteobacteria into two separate subclasses, may be appropriate.

Gene location and bacterial sequence divergence

Previous comparison of a relatively small set of homologous genes from Escherichia coli and Salmonella typhimurium revealed that genes nearer to the origin of replication had substitution rates lower than genes closer to the replication terminus. The recently completed sequences of numerous bacterial genomes have allowed us to test whether this effect of distance from the replication origin on substitution rates, as observed for the E. coli-S. typhimurium comparison, is a general feature of bacterial genomes. Extending the analysis to all 3,000 E. coli-S. typhimurium homologs confirmed the significant association between chromosomal position and synonymous site divergence. However, the effect, though still significant, is not as dramatic as originally thought. A similar association between relative chromosomal location and synonymous substitution rate was detected in the majority of other bacterial species comparisons within alpha- and gamma- Proteobacteria, and Firmicutes but was absent in Chlamydiales. The opposite trend, i.e., a decrease in synonymous divergence with distance from the replication origin, was detected in Mycobacteria. Analysis of the patterns of nucleotide substitutions revealed that the distance effect is not affected by gene orientation and is mainly caused by an increase in rates of transversions, suggesting that this effect may not be caused by recombinational repair or biased gene conversion, as originally suggested.

A systematic investigation identifies a significant number of probable pseudogenes in the Escherichia coli genome

Pseudogenes are open reading frames (ORFs) encoding dysfunctional proteins with high homology to known protein-coding genes. Although pseudogenes were reported to exist in the genomes of many eukaryotes and bacteria, no systematic search for pseudogenes in the Escherichia coli genome has been carried out. Genome comparisons of E. coli strains K-12 and O157 revealed that many protein-coding sequences have prematurely terminated orthologs encoding unstable proteins. To systematically screen for pseudogenes, we selected ORFs generated by premature termination of the orthologous protein-coding genes and subsequently excluded those possibly arising from sequence errors. Lastly we eliminated those with close homologs in this and other species, as these shortened ORFs may actually have functions. The process produced 95 and 101 pseudogene candidates in K-12 and O157, respectively. The assigned three-dimensional structures suggest that most of the encoded proteins cannot fold properly and thus are dysfunctional, indicating that they are probably pseudogenes. Therefore, the existence of a significant number of probable pseudogenes in E. coli is predicted, awaiting experimental verification. Most of them were found to be genes with paralogs or horizontally transferred genes or both. We suggest that pseudogenes constitute a small fraction of the genomes of free-living bacteria in general, reflecting the faster elimination than production of pseudogenes.

Expression of green fluorescent protein in Rickettsia conorii

Rickettsiae are obligate intracellular class III pathogens for which genetic manipulation has only recently been shown to be feasible. Such experiments were restricted to the typhus group rickettsiae, namely R. typhi and R. prowazekii. Here we report the first genetic manipulation of Rickettsia conorii, the bacterial agent responsible for the Mediterranean spotted fever. A gene encoding a variant of the green fluorescent protein under the control of the sterically repressed promoter (srp) from E. coli was integrated into the genome of this bacteria and detected by FACS analysis.

Studying genomes through the aeons: protein families, pseudogenes and proteome evolution

Protein families can be used to understand many aspects of genomes, both their "live" and their "dead" parts (i.e. genes and pseudogenes). Surveys of genomes have revealed that, in every organism, there are always a few large families and many small ones, with the overall distribution following a power-law. This commonality is equally true for both genes and pseudogenes, and exists despite the fact that the specific families that are enlarged differ greatly between organisms. Furthermore, because of family structure there is great redundancy in proteomes, a fact linked to the large number of dispensable genes for each organism and the small size of the minimal, indispensable sub-proteome. Pseudogenes in prokaryotes represent families that are in the process of being dispensed with. In particular, the genome sequences of certain pathogenic bacteria (Mycobacterium leprae, Yersinia pestis and Rickettsia prowazekii) show how an organism can undergo reductive evolution on a large scale (i.e. the dying out of families) as a result of niche change. There appears to be less pressure to delete pseudogenes in eukaryotes. These can be divided into two varieties, duplicated and processed, where the latter involves reverse transcription from an mRNA intermediate. We discuss these collectively in yeast, worm, fly, and human. The fly has few pseudogenes apparently because of its high rate of genomic DNA deletion. In the other three organisms, the distribution of pseudogenes on the chromosome and amongst different families is highly non-uniform. Pseudogenes tend not to occur in the middle of chromosome arms, and tend to be associated with lineage-specific (as opposed to highly conserved) families that have environmental-response functions. This may be because, rather than being dead, they may form a reservoir of diverse "extra parts" that can be resurrected to help an organism adapt to its surroundings. In yeast, there may be a novel mechanism involving the [PSI+] prion that potentially enables this resurrection. In worm, the pseudogenes tend to arise out of families (e.g. chemoreceptors) that are greatly expanded in it compared to the fly. The human genome stands out in having many processed pseudogenes. These have a character very different from those of the duplicated variety, to a large extent just representing random insertions. Thus, their occurrence tends to be roughly in proportion to the amount of mRNA for a particular protein and to reflect the extent of the intergenic sequences. Further information about pseudogenes is available at http://genecensus.org/pseudogene

Protein coding palindromes are a unique but recurrent feature in Rickettsia

Rickettsia are unique in inserting in-frame a number of palindromic sequences within protein coding regions. In this study, we extensively analyzed repeated sequences in the genome of Rickettsia conorii and examined their locations in regard to coding versus noncoding regions. We identified 656 interspersed repeated sequences classified into 10 distinct families. Of the 10 families, three palindromic sequence families showed clear cases of insertions into open reading frames (ORFs). The location of those in-frame insertions appears to be always compatible with the encoded protein three-dimensional (3-D) fold and function. We provide evidence for a progressive loss of the palindromic property over time after the insertions. This comprehensive study of Rickettsia repeats confirms and extends our previous observations and further indicates a significant role of selfish DNAs in the creation and modification of proteins.

The evolving genome of Salmonella enterica serovar Pullorum

Salmonella enterica serovar Pullorum is a fowl-adapted bacterial pathogen that causes dysentery (pullorum disease). Host adaptation and special pathogenesis make S. enterica serovar Pullorum an exceptionally good system for studies of bacterial evolution and speciation, especially regarding pathogen-host interactions and the acquisition of pathogenicity. We constructed a genome map of S. enterica serovar Pullorum RKS5078, using I-CeuI, XbaI, AvrII, and SpeI and Tn10 insertions. Pulsed-field gel electrophoresis was employed to separate the large DNA fragments generated by the endonucleases. The genome is 4,930 kb, which is similar to most salmonellas. However, the genome of S. enterica serovar Pullorum RKS5078 is organized very differently from the majority of salmonellas, with three major inversions and one translocation. This extraordinary genome structure was seen in most S. enterica serovar Pullorum strains examined, with different structures in a minority of S. enterica serovar Pullorum strains. We describe the coexistence of different genome structures among the same bacteria as genomic plasticity. Through comparisons with S. enterica serovar Typhimurium, we resolved seven putative insertions and eight deletions ranging in size from 12 to 157 kb. The genomic plasticity seen among S. enterica serovar Pullorum strains supported our hypothesis about its association with bacterial evolution: a large genomic insertion (157 kb in this case) disrupted the genomic balance, and rebalancing by independent recombination events in individual lineages resulted in diverse genome structures. As far as the structural plasticity exists, the S. enterica serovar Pullorum genome will continue evolving to reach a further streamlined and balanced structure.

Characterization and transcriptional analysis of gene clusters for a type IV secretion machinery in human granulocytic and monocytic ehrlichiosis agents

Anaplasma (Ehrlichia) phagocytophila and Ehrlichia chaffeensis, the etiologic agents of granulocytic and monocytic ehrlichioses, respectively, are obligatory intracellular bacteria that cause febrile systemic illness in humans. We identified and characterized clusters of genes for a type IV secretion machinery in these two bacteria, and analyzed their gene expression in cell culture and mammalian hosts. Eight virB and virD genes were found in each bacterial genome, and all of the genes were transcribed in cell culture. Although the gene order and orientation were similar to those found in other bacteria, the eight virB and virD genes were clustered at two separate loci in each genome. Five of the genes (virB8, virB9, virB10, virB11, and virD4) were located downstream from a ribA gene. These five genes in both A. phagocytophila and E. chaffeensis were polycistronically transcribed and controlled through at least two tandem promoters located upstream of the virB8 gene in human leukemia cell lines. The virB9 gene of A. phagocytophila was transcriptionally active in peripheral blood leukocytes from human ehrlichiosis patients and experimentally infected animals. Three of the remaining genes (virB3, virB4, and virB6) of both A. phagocytophila and E. chaffeensis were arranged downstream from a sodB gene and cotranscribed with the sodB gene through one or more sodB promoters in human leukocytes. This suggests that transcription of the three virB genes in these two Anaplasma and Ehrlichia spp. is regulated by factors that influence the sodB gene expression. This unique regulation of gene expression for the type IV secretion system may be associated with intracellular survival and replication of Anaplasma and Ehrlichia spp. in granulocytes or monocytes.

The genome of camelpox virus

Camelpox virus (CMLV), a member of the Orthopoxvirus genus in the Poxviridae, is the etiologic agent of a disease of camels. Here we report the CMLV genomic sequence with analysis. The 205,719-bp CMLV genome contains 211 putative genes and consists of a central region bound by identical inverted terminal repeats of approximately 7 kb. A high degree of similarity in gene order, gene content, and amino acid composition in the region located between CMLV017 and CMLV184 (average 96% amino acid identity to vaccinia virus (VACV)) indicates a close structural and functional relationship between CMLV and other known orthopoxviruses (OPVs). Notably, CMLV contains a unique region of approximately 3 kb, which encodes three ORFs (CMLV185, CMLV186, CMLV187) absent in other OPVs. These ORFs are most similar to B22R homologues found in other chordopoxvirus genera. Among OPVs, CMLV is the most closely related to variola virus (VARV), sharing all genes involved in basic replicative functions and the majority of genes involved in other host-related functions. Differences between CMLV and VARV include deletion and disruption of a large number of genes. Twenty-seven CMLV ORFs are absent in VARV, including seven full-length homologues of NMDA-like receptor, phospholipase D, Schlafen, MT-4 virulence, kelch, VACV C8L, and cowpox (CPXV) B21R proteins. Thirty-eight CMLV ORFs, some of which are fragments of larger genes, differ in size from corresponding VARV ORFs by more than 10% (amino acids). Genome structure and phylogenetic analysis of DNA sequences for all ORFs indicate that CMLV is clearly distinct from VARV and VACV and, as it has been suggested for VARV, it may have originated from a CPXV virus-like ancestor.

Microbial minimalism: genome reduction in bacterial pathogens

When bacterial lineages make the transition from free-living or facultatively parasitic life cycles to permanent associations with hosts, they undergo a major loss of genes and DNA. Complete genome sequences are providing an understanding of how extreme genome reduction affects evolutionary directions and metabolic capabilities of obligate pathogens and symbionts.

The Sinorhizobium meliloti stringent response affects multiple aspects of symbiosis

Sinorhizobium meliloti and host legumes enter into a nitrogen-fixing, symbiotic relationship triggered by an exchange of signals between bacteria and plant. S. meliloti produces Nod factor, which elicits the formation of nodules on plant roots, and succinoglycan, an exopolysaccharide that allows for bacterial invasion and colonization of the host. The biosynthesis of these molecules is well defined, but the specific regulation of these compounds is not completely understood. Bacteria control complex regulatory networks by the production of ppGpp, the effector molecule of the stringent response, which induces physiological change in response to adverse growth conditions and can also control bacterial development and virulence. Through detailed analysis of an S. meliloti mutant incapable of producing ppGpp, we show that the stringent response is required for nodule formation and regulates the production of succinoglycan. Although it remains unknown whether these phenotypes are connected, we have isolated suppressor strains that restore both defects and potentially identify key downstream regulatory genes. These results indicate that the S. meliloti stringent response has roles in both succinoglycan production and nodule formation and, more importantly, that control of bacterial physiology in response to the plant and surrounding environment is critical to the establishment of a successful symbiosis.

A bacterial guanine nucleotide exchange factor activates ARF on Legionella phagosomes

The intracellular pathogen Legionella pneumophila subverts vesicle traffic in eukaryotic host cells to create a vacuole that supports replication. The dot/icm genes encode a protein secretion apparatus that L. pneumophila require for biogenesis of this vacuole. Here we show that L. pneumophila produce a protein called RalF that functions as an exchange factor for the ADP ribosylation factor (ARF) family of guanosine triphosphatases (GTPases). The RalF protein is required for the localization of ARF on phagosomes containing L. pneumophila. Translocation of RalF protein through the phagosomal membrane is a dot/icm-dependent process. Thus, RalF is a substrate of the Dot/Icm secretion apparatus.

Microarrays in ecology and evolution: a preview

Microarray technology provides a new tool with which molecular ecologists and evolutionary biologists can survey genome-wide patterns of gene expression within and among species. New analytical approaches based on analysis of variance will allow quantification of the contributions of among individual variation, genotype, sex, microenvironment, population structure, and geography to variation in gene expression. Applications of this methodology are reviewed in relation to studies of mechanisms of adaptation and divergence; delineation of developmental and physiological pathways and networks; characterization of quantitative genetic parameters at the level of transcription ('quantitative genomics'); molecular dissection of parasitism and symbiosis; and studies of the diversification of gene content. Establishment of microarray resources is neither prohibitively expensive nor technologically demanding, and a commitment to development of gene expression profiling methods for nonmodel organisms could have a tremendous impact on molecular and genetic research at the interface of organismal and population biology.