Obligate intracellular parasites: Rickettsia prowazekii and Chlamydia trachomatis

Insights into the Two Most Common Cancers of Primitive Gut-Derived Structures and Their Microbial Connections

The gastrointestinal and respiratory systems are closely linked in different ways, including from the embryological, anatomical, cellular, and physiological angles. The highest number (and various types) of microorganisms live in the large intestine/colon, and constitute the normal microbiota in healthy people. Adverse alterations of the microbiota or dysbiosis can lead to chronic inflammation. If this detrimental condition persists, a sequence of pathological events can occur, such as inflammatory bowel disease, dysplasia or premalignant changes, and finally, cancer. One of the most commonly identified bacteria in both inflammatory bowel disease and colon cancer is Escherichia coli. On the other hand, patients with inflammatory bowel disease are at risk of several other diseases-both intestinal (such as malnutrition and intestinal obstruction, besides cancer) and extraintestinal (such as arthritis, bronchiectasis, and cancer risk). Cancers of the lung and colon are the two most common malignancies occurring worldwide (except for female breast cancer). Like the bacterial role in colon cancer, many studies have shown a link between chronic Chlamydia pneumoniae infection and lung cancer. However, in colon cancer, genotoxic colibactin-producing E. coli belonging to the B2 phylogroup may promote tumorigenesis. Furthermore, E. coli is believed to play an important role in the dissemination of cancer cells from the primary colonic site. Currently, seven enteric pathogenic E. coli subtypes have been described. Conversely, three Chlamydiae can cause infections in humans (C. trachomatis may increase the risk of cervical and ovarian cancers). Nonetheless, striking genomic plasticity and genetic modifications allow E. coli to constantly adjust to the surrounding environment. Consequently, E. coli becomes resistant to antibiotics and difficult to manage. To solve this problem, scientists are thinking of utilizing suitable lytic bacteriophages (viruses that infect and kill bacteria). Several bacteriophages of E. coli and Chlamydia species are being evaluated for this purpose.

Analysis of the Type 4 Effectome across the Genus Rickettsia

Rickettsia are obligate intracellular bacteria primarily carried by arthropod hosts. The genus Rickettsia contains several vertebrate pathogens vectored by hematophagous arthropods. Despite the potential for disease, our understanding of Rickettsias are limited by the difficulties associated with growing and manipulating obligate intracellular bacteria. To aid with this, our lab conducted an analysis of eight genomes and three plasmids from across the genus Rickettsia. Using OPT4e, a learning algorithm-based program designed to identify effector proteins secreted by the type 4 secretion system, we generated a putative effectome for the genus. We then consolidated effectors into homolog sets to identify effectors unique to Rickettsia with different life strategies or evolutionary histories. We also compared predicted effectors to non-effectors for differences in G+C content and gene splitting. Based on this analysis, we predicted 1571 effectors across the genus, resulting in 604 homolog sets. Each species had unique homolog sets, while 42 were present in all eight species analyzed. Effectors were flagged in association with pathogenic, tick and flea-borne Rickettsia. Predicted effectors also varied in G+C content and frequency of gene splitting as compared to non-effectors. Species effector repertoires show signs of expansion, degradation, and horizontal acquisition associated with lifestyle and lineage.

Adaptation of bacteria to glyphosate: a microevolutionary perspective of the enzyme 5-enolpyruvylshikimate-3-phosphate synthase

Glyphosate is the leading herbicide worldwide, but it also affects prokaryotes because it targets the central enzyme (5-enolpyruvylshikimate-3-phosphate, EPSP) of the shikimate pathway in the synthesis of the three essential aromatic amino acids in bacteria, fungi and plants. Our results reveal that bacteria may easily become resistant to glyphosate through changes in the 5-enolpyruvylshikimate-3-phosphate synthase active site. This indicates the importance of examining how glyphosate affects microbe-mediated ecosystem functions and human microbiomes.

Mitochondrial protein phosphorylation in yeast revisited

Protein phosphorylation is one of the best-known post-translational modifications occurring in all domains of life. In eukaryotes, protein phosphorylation affects all cellular compartments including mitochondria. High-throughput techniques of mass spectrometry combined with cell fractionation and biochemical methods yielded thousands of phospho-sites on hundreds of mitochondrial proteins. We have compiled the information on mitochondrial protein kinases and phosphatases and their substrates in Saccharomyces cerevisiae and provide the current state-of-the-art overview of mitochondrial protein phosphorylation in this model eukaryote. Using several examples, we describe emerging features of the yeast mitochondrial phosphoproteome and present challenges lying ahead in this exciting field.

An Unusual Route for p-Aminobenzoate Biosynthesis in Chlamydia trachomatis Involves a Probable Self-Sacrificing Diiron Oxygenase

Chlamydia trachomatis lacks the canonical genes required for the biosynthesis of p-aminobenzoate (pABA), a component of essential folate cofactors. Previous studies revealed a single gene from C. trachomatis, the CT610 gene, that rescues Escherichia coli ΔpabA, ΔpabB, and ΔpabC mutants, which are otherwise auxotrophic for pABA. CT610 shares low sequence similarity to nonheme diiron oxygenases, and the previously solved crystal structure revealed a diiron active site. Genetic studies ruled out several potential substrates for CT610-dependent pABA biosynthesis, including chorismate and other shikimate pathway intermediates, leaving the actual precursor(s) unknown. Here, we supplied isotopically labeled potential precursors to E. coli ΔpabA cells expressing CT610 and found that the aromatic portion of tyrosine was highly incorporated into pABA, indicating that tyrosine is a precursor for CT610-dependent pABA biosynthesis. Additionally, in vitro enzymatic experiments revealed that purified CT610 exhibits low pABA synthesis activity under aerobic conditions in the absence of tyrosine or other potential substrates, where only the addition of a reducing agent such as dithiothreitol appears to stimulate pABA production. Furthermore, site-directed mutagenesis studies revealed that two conserved active site tyrosine residues are essential for the pABA synthesis reaction in vitro Thus, the current data are most consistent with CT610 being a unique self-sacrificing enzyme that utilizes its own active site tyrosine residue(s) for pABA biosynthesis in a reaction that requires O₂ and a reduced diiron cofactor.IMPORTANCE Chlamydia trachomatis is the most reported sexually transmitted infection in the United States and the leading cause of infectious blindness worldwide. Unlike many other intracellular pathogens that have undergone reductive evolution, C. trachomatis is capable of de novo biosynthesis of the essential cofactor tetrahydrofolate using a noncanonical pathway. Here, we identify the biosynthetic precursor to the p-aminobenzoate (pABA) portion of folate in a process that requires the CT610 enzyme from C. trachomatis We further provide evidence that CT610 is a self-sacrificing or "suicide" enzyme that uses its own amino acid residue(s) as the substrate for pABA synthesis. This work provides the foundation for future investigation of this chlamydial pABA synthase, which could lead to new therapeutic strategies for C. trachomatis infections.

Inhibition of tRNA Synthetases Induces Persistence in Chlamydia

Chlamydia trachomatis is the leading cause of bacterial sexually transmitted infections, and Chlamydia pneumoniae causes community-acquired respiratory infections. In vivo, the host immune system will release gamma interferon (IFN-γ) to combat infection. IFN-γ activates human cells to produce the tryptophan (Trp)-catabolizing enzyme indoleamine 2,3-dioxygenase (IDO). Consequently, there is a reduction in cytosolic Trp in IFN-γ-activated host cells. In evolving to obligate intracellular dependence, Chlamydia has significantly reduced its genome size and content, as it relies on the host cell for various nutrients. Importantly, C. trachomatis and C. pneumoniae are Trp auxotrophs and are starved for this essential nutrient when the human host cell is exposed to IFN-γ. To survive this, chlamydiae enter an alternative developmental state referred to as persistence. Chlamydial persistence is characterized by a halt in the division cycle, aberrant morphology, and, in the case of IFN-γ-induced persistence, Trp codon-dependent changes in transcription. We hypothesize that these changes in transcription are dependent on the particular amino acid starvation state. To investigate the chlamydial response mechanisms acting when other amino acids become limiting, we tested the efficacy of prokaryote-specific tRNA synthetase inhibitors, indolmycin and AN3365, to mimic starvation of Trp and leucine, respectively. We show that these drugs block chlamydial growth and induce changes in morphology and transcription consistent with persistence. Importantly, growth inhibition was reversed when the compounds were removed from the medium. With these data, we find that indolmycin and AN3365 are valid tools that can be used to mimic the persistent state independently of IFN-γ.

Mycobacterium leprae's evolution and environmental adaptation

Leprosy is an ancient disease caused by the acid-fast bacillus Mycobacterium leprae, also known as Hansen's bacillus. M. leprae is an obligate intracellular microorganism with a marked Schwann cell tropism and is the only human pathogen capable of invading the superficial peripheral nerves. The transmission mechanism of M. leprae is not fully understood; however, the nasal mucosa is accepted as main route of M. leprae entry to the human host. The complete sequencing and the comparative genome analysis show that M. leprae underwent a genome reductive evolution process, as result of lifestyle change and adaptation to different environments; some of lost genes are homologous to those of host cells. Thus, M. leprae reduced its genome size to 3.3 Mbp, contributing to obtain the lowest GC content (approximately 58%) among mycobacteria. The M. leprae genome contains 1614 open reading frames coding for functional proteins, and 1310 pseudogenes corresponding to 41% of the genome, approximately. Comparative analyses to different microorganisms showed that M. leprae possesses the highest content of pseudogenes among pathogenic and non-pathogenic bacteria and archaea. The pathogen adaptation into host cells, as the Schwann cells, brought about the reduction of the genome and induced multiple gene inactivation. The present review highlights the characteristics of genome's reductive evolution that M. leprae experiences in the genetic aspects compared with other pathogens. The possible mechanisms of pseudogenes formation are discussed.

Metagenomic Analysis of Fish-Associated Ca. Parilichlamydiaceae Reveals Striking Metabolic Similarities to the Terrestrial Chlamydiaceae

Chlamydiae are an example of obligate intracellular bacteria that possess highly reduced, compact genomes (1.0-3.5 Mbp), reflective of their abilities to sequester many essential nutrients from the host that they no longer need to synthesize themselves. The Chlamydiae is a phylum with a very wide host range spanning mammals, birds, fish, invertebrates, and unicellular protists. This ecological and phylogenetic diversity offers ongoing opportunities to study intracellular survival and metabolic pathways and adaptations. Of particular evolutionary significance are Chlamydiae from the recently proposed Ca. Parilichlamydiaceae, the earliest diverging clade in this phylum, species of which are found only in aquatic vertebrates. Gill extracts from three Chlamydiales-positive Australian aquaculture species (Yellowtail kingfish, Striped trumpeter, and Barramundi) were subject to DNA preparation to deplete host DNA and enrich microbial DNA, prior to metagenome sequencing. We assembled chlamydial genomes corresponding to three Ca. Parilichlamydiaceae species from gill metagenomes, and conducted functional genomics comparisons with diverse members of the phylum. This revealed highly reduced genomes more similar in size to the terrestrial Chlamydiaceae, standing in contrast to members of the Chlamydiae with a demonstrated cosmopolitan host range. We describe a reduction in genes encoding synthesis of nucleotides and amino acids, among other nutrients, and an enrichment of predicted transport proteins. Ca. Parilichlamydiaceae share 342 orthologs with other chlamydial families. We hypothesize that the genome reduction exhibited by Ca. Parilichlamydiaceae and Chlamydiaceae is an example of within-phylum convergent evolution. The factors driving these events remain to be elucidated.

Chlamydia trachomatis regulates growth and development in response to host cell fatty acid availability in the absence of lipid droplets

Chlamydia trachomatis (Ct) is a Gram-negative obligate intracellular pathogen of humans that causes significant morbidity from sexually transmitted and ocular diseases globally. Ct acquires host fatty acids (FA) to meet the metabolic and growth requirements of the organism. Lipid droplets (LDs) are storehouses of FAs in host cells and have been proposed to be a source of FAs for the parasitophorous vacuole, termed inclusion, in which Ct replicates. Previously, cells devoid of LDs were shown to produce reduced infectious progeny at 24 hr postinfection (hpi). Here, although we also found reduced progeny at 24 hpi, there were significantly more progeny at 48 hpi in the absence of LDs compared to the control wild-type (WT) cells. These findings were confirmed using transmission electron microscopy where cells without LDs were shown to have significantly more metabolically active reticulate bodies at 24 hpi and significantly more infectious but metabolically inert elementary bodies at 48 hpi than WT cells. Furthermore, by measuring basal oxygen consumption rates (OCR) using extracellular flux analysis, Ct infected cells without LDs had higher OCRs at 24 hpi than cells with LDs, confirming ongoing metabolic activity in the absence of LDs. Although the FA oleic acid is a major source of phospholipids for Ct and stimulates LD synthesis, treatment with oleic acid, but not other FAs, enhanced growth and led to an increase in basal OCR in both LD depleted and WT cells, indicating that FA transport to the inclusion is not affected by the loss of LDs. Our results show that Ct regulates inclusion metabolic activity and growth in response to host FA availability in the absence of LDs.

Chlamydia pneumoniae exploits adipocyte lipid chaperone FABP4 to facilitate fat mobilization and intracellular growth in murine adipocytes

Fatty acid-binding protein 4 (FABP4), a cytosolic lipid chaperone predominantly expressed in adipocytes and macrophages, modulates lipid fluxes, trafficking, signaling, and metabolism. Recent studies have demonstrated that FABP4 regulates metabolic and inflammatory pathways, and in mouse models its inhibition can improve type 2 diabetes mellitus and atherosclerosis. However, the role of FABP4 in bacterial infection, metabolic crosstalk between host and pathogen, and bacterial pathogenesis have not been studied. As an obligate intracellular pathogen, Chlamydia pneumoniae needs to obtain nutrients such as ATP and lipids from host cells. Here, we show that C. pneumoniae successfully infects and proliferates in murine adipocytes by inducing hormone sensitive lipase (HSL)-mediated lipolysis. Chemical inhibition or genetic manipulation of HSL significantly abrogated the intracellular growth of C. pneumoniae in adipocytes. Liberated free fatty acids were utilized to generate ATP via β-oxidation, which C. pneumoniae usurped for its replication. Strikingly, chemical inhibition or genetic silencing of FABP4 significantly abrogated C. pneumoniae infection-induced lipolysis and mobilization of liberated FFAs, resulting in reduced bacterial growth in adipocytes. Collectively, these results demonstrate that C. pneumoniae exploits host FABP4 to facilitate fat mobilization and intracellular replication in adipocytes. This work uncovers a novel strategy used by intracellular pathogens for acquiring energy via hijacking of the host lipid metabolism pathway.

Biotic interactions as drivers of algal origin and evolution

Contents 670 I. 671 II. 671 III. 676 IV. 678 678 References 678 SUMMARY: Biotic interactions underlie life's diversity and are the lynchpin to understanding its complexity and resilience within an ecological niche. Algal biologists have embraced this paradigm, and studies building on the explosive growth in omics and cell biology methods have facilitated the in-depth analysis of nonmodel organisms and communities from a variety of ecosystems. In turn, these advances have enabled a major revision of our understanding of the origin and evolution of photosynthesis in eukaryotes, bacterial-algal interactions, control of massive algal blooms in the ocean, and the maintenance and degradation of coral reefs. Here, we review some of the most exciting developments in the field of algal biotic interactions and identify challenges for scientists in the coming years. We foresee the development of an algal knowledgebase that integrates ecosystem-wide omics data and the development of molecular tools/resources to perform functional analyses of individuals in isolation and in populations. These assets will allow us to move beyond mechanistic studies of a single species towards understanding the interactions amongst algae and other organisms in both the laboratory and the field.

Bacterial genome sequencing in clinical microbiology: a pathogen-oriented review

In recent years, whole-genome sequencing (WGS) has been perceived as a technology with the potential to revolutionise clinical microbiology. Herein, we reviewed the literature on the use of WGS for the most commonly encountered pathogens in clinical microbiology laboratories: Escherichia coli and other Enterobacteriaceae, Staphylococcus aureus and coagulase-negative staphylococci, streptococci and enterococci, mycobacteria and Chlamydia trachomatis. For each pathogen group, we focused on five different aspects: the genome characteristics, the most common genomic approaches and the clinical uses of WGS for (i) typing and outbreak analysis, (ii) virulence investigation and (iii) in silico antimicrobial susceptibility testing. Of all the clinical usages, the most frequent and straightforward usage was to type bacteria and to trace outbreaks back. A next step toward standardisation was made thanks to the development of several new genome-wide multi-locus sequence typing systems based on WGS data. Although virulence characterisation could help in various particular clinical settings, it was done mainly to describe outbreak strains. An increasing number of studies compared genotypic to phenotypic antibiotic susceptibility testing, with mostly promising results. However, routine implementation will preferentially be done in the workflow of particular pathogens, such as mycobacteria, rather than as a broadly applicable generic tool. Overall, concrete uses of WGS in routine clinical microbiology or infection control laboratories were done, but the next big challenges will be the standardisation and validation of the procedures and bioinformatics pipelines in order to reach clinical standards.

Genome-wide codon usage profiling of ocular infective Chlamydia trachomatis serovars and drug target identification

Chlamydia trachomatis (C.t) is a Gram-negative obligate intracellular bacteria and is a major causative of infectious blindness and sexually transmitted diseases. Among the varied serovars of this organism, A, B and C are reported as prominent ocular pathogens. Genomic studies of these strains shall aid in deciphering potential drug targets and genomic influence on pathogenesis. Hence, in this study we performed deep statistical profiling of codon usage in these serovars. The overall base composition analysis reveals that these serovars are over biased to AU than GC. Similarly, relative synonymous codon usage also showed preference towards A/U ending codons. Parity Rule 2 analysis inferred unequal distribution of AT and GC, indicative of other unknown factors acting along with mutational pressure to influence codon usage bias (CUB). Moreover, absolute quantification of CUB also revealed lower bias across these serovars. The effect of natural selection on CUB was also confirmed by neutrality plot, reinforcing natural selection under mutational pressure turned to be a pivotal role in shaping the CUB in the strains studied. Correspondence analysis (COA) clarified that, C.t C/TW-3 to show a unique trend in codon usage variation. Host influence analysis on shaping the codon usage pattern also inferred some speculative relativity. In a nutshell, our finding suggests that mutational pressure is the dominating factor in shaping CUB in the strains studied, followed by natural selection. We also propose potential drug targets based on cumulative analysis of strand bias, CUB and human non-homologue screening.

Multiple Substrate Usage of Coxiella burnetii to Feed a Bipartite Metabolic Network

The human pathogen Coxiella burnetii causes Q-fever and is classified as a category B bio-weapon. Exploiting the development of the axenic growth medium ACCM-2, we have now used ¹³C-labeling experiments and isotopolog profiling to investigate the highly diverse metabolic network of C. burnetii. To this aim, C. burnetii RSA 439 NMII was cultured in ACCM-2 containing 5 mM of either [U-¹³C₃]serine, [U-¹³C₆]glucose, or [U-¹³C₃]glycerol until the late-logarithmic phase. GC/MS-based isotopolog profiling of protein-derived amino acids, methanol-soluble polar metabolites, fatty acids, and cell wall components (e.g., diaminopimelate and sugars) from the labeled bacteria revealed differential incorporation rates and isotopolog profiles. These data served to decipher the diverse usages of the labeled substrates and the relative carbon fluxes into the core metabolism of the pathogen. Whereas, de novo biosynthesis from any of these substrates could not be found for histidine, isoleucine, leucine, lysine, phenylalanine, proline and valine, the other amino acids and metabolites under study acquired ¹³C-label at specific rates depending on the nature of the tracer compound. Glucose was directly used for cell wall biosynthesis, but was also converted into pyruvate (and its downstream metabolites) through the glycolytic pathway or into erythrose 4-phosphate (e.g., for the biosynthesis of tyrosine) via the non-oxidative pentose phosphate pathway. Glycerol efficiently served as a gluconeogenetic substrate and could also be used via phosphoenolpyruvate and diaminopimelate as a major carbon source for cell wall biosynthesis. In contrast, exogenous serine was mainly utilized in downstream metabolic processes, e.g., via acetyl-CoA in a complete citrate cycle with fluxes in the oxidative direction and as a carbon feed for fatty acid biosynthesis. In summary, the data reflect multiple and differential substrate usages by C. burnetii in a bipartite-type metabolic network, resembling the overall topology of the related pathogen Legionella pneumophila. These strategies could benefit the metabolic capacities of the pathogens also as a trait to adapt for replication under intracellular conditions.

Metabolic Adaptations of Intracellullar Bacterial Pathogens and their Mammalian Host Cells during Infection ("Pathometabolism")

Several bacterial pathogens that cause severe infections in warm-blooded animals, including humans, have the potential to actively invade host cells and to efficiently replicate either in the cytosol or in specialized vacuoles of the mammalian cells. The interaction between these intracellular bacterial pathogens and the host cells always leads to multiple physiological changes in both interacting partners, including complex metabolic adaptation reactions aimed to promote proliferation of the pathogen within different compartments of the host cells. In this chapter, we discuss the necessary nutrients and metabolic pathways used by some selected cytosolic and vacuolar intracellular pathogens and--when available--the links between the intracellular bacterial metabolism and the expression of the virulence genes required for the intracellular bacterial replication cycle. Furthermore, we address the growing evidence that pathogen-specific factors may also trigger metabolic responses of the infected mammalian cells affecting the carbon and nitrogen metabolism as well as defense reactions. We also point out that many studies on the metabolic host cell responses induced by the pathogens have to be scrutinized due to the use of established cell lines as model host cells, as these cells are (in the majority) cancer cells that exhibit a dysregulated primary carbon metabolism. As the exact knowledge of the metabolic host cell responses may also provide new concepts for antibacterial therapies, there is undoubtedly an urgent need for host cell models that more closely reflect the in vivo infection conditions.

Rapid screening for citrus canker resistance employing pathogen-associated molecular pattern-triggered immunity responses

Citrus canker, caused by the bacterial pathogen Xanthomonas citri ssp. citri (Xcc), has been attributed to millions of dollars in loss or damage to commercial citrus crops in subtropical production areas of the world. Since identification of resistant plants is one of the most effective methods of disease management, the ability to screen for resistant seedlings plays a key role in the production of a long-term solution to canker. Here, an inverse correlation between reactive oxygen species (ROS) production by the plant and the ability of Xcc to grow and form lesions on infected plants is reported. Based on this information, a novel screening method that can rapidly identify citrus seedlings that are less susceptible to early infection by Xcc was devised by measuring ROS accumulation triggered by a 22-amino acid sequence of the conserved N-terminal part of flagellin (flg22) from X. citri ssp. citri (Xcc-flg22). In addition to limiting disease symptoms, ROS production was also correlated with the expression of basal defense-related genes such as the pattern recognition receptors LRR8 and FLS2, the leucine-rich repeat receptor-like protein RLP12, and the defense-related gene PR1, indicating an important role for pathogen-associated molecular pattern-triggered immunity (PTI) in determining resistance to citrus canker. Moreover, the differential expression patterns observed amongst the citrus seedlings demonstrated the existence of genetic variations in the PTI response among citrus species/varieties.

Mutational Analysis of the Chlamydia muridarum Plasticity Zone

Pathogenically diverse Chlamydia spp. can have surprisingly similar genomes. Chlamydia trachomatis isolates that cause trachoma, sexually transmitted genital tract infections (chlamydia), and invasive lymphogranuloma venereum (LGV) and the murine strain Chlamydia muridarum share 99% of their gene content. A region of high genomic diversity between Chlamydia spp. termed the plasticity zone (PZ) may encode niche-specific virulence determinants that dictate pathogenic diversity. We hypothesized that PZ genes might mediate the greater virulence and gamma interferon (IFN-γ) resistance of C. muridarum compared to C. trachomatis in the murine genital tract. To test this hypothesis, we isolated and characterized a series of C. muridarum PZ nonsense mutants. Strains with nonsense mutations in chlamydial cytotoxins, guaBA-add, and a phospholipase D homolog developed normally in cell culture. Two of the cytotoxin mutants were less cytotoxic than the wild type, suggesting that the cytotoxins may be functional. However, none of the PZ nonsense mutants exhibited increased IFN-γ sensitivity in cell culture or were profoundly attenuated in a murine genital tract infection model. Our results suggest that C. muridarum PZ genes are transcribed--and some may produce functional proteins--but are dispensable for infection of the murine genital tract.

Evolution, phylogeny, and molecular epidemiology of Chlamydia

The Chlamydiaceae are a family of obligate intracellular bacteria characterized by a unique biphasic developmental cycle. It encompasses the single genus Chlamydia, which involves nine species that affect a wide range of vertebral hosts, causing infections with serious impact on human health (mainly due to Chlamydia trachomatis infections) and on farming and veterinary industries. It is believed that Chlamydiales originated ∼700mya, whereas C. trachomatis likely split from the other Chlamydiaceae during the last 6mya. This corresponds to the emergence of modern human lineages, with the first descriptions of chlamydial infections as ancient as four millennia. Chlamydiaceae have undergone a massive genome reduction, on behalf of the deletional bias "use it or lose it", stabilizing at 1-1.2Mb and keeping a striking genome synteny. Their phylogeny reveals species segregation according to biological properties, with huge differences in terms of host range, tissue tropism, and disease outcomes. Genome differences rely on the occurrence of mutations in the >700 orthologous genes, as well as on events of recombination, gene loss, inversion, and paralogous expansion, affecting both a hypervariable region named the plasticity zone, and genes essentially encoding polymorphic and transmembrane head membrane proteins, type III secretion effectors and some metabolic pathways. Procedures for molecular typing are still not consensual but have allowed the knowledge of molecular epidemiology patterns for some species as well as the identification of outbreaks and emergence of successful clones for C. trachomatis. This manuscript intends to provide a comprehensive review on the evolution, phylogeny, and molecular epidemiology of Chlamydia.

Analysis of metabolic evolution in bacteria using whole-genome metabolic models

Recent advances in the automation of metabolic model reconstruction have led to the availability of draft-quality metabolic models (predicted reaction complements) for multiple bacterial species. These reaction complements can be considered as trait representations and can be used for ancestral state reconstruction to infer the most likely metabolic complements of common ancestors of all bacteria with generated metabolic models. We present here an ancestral state reconstruction for 141 extant bacteria and analyze the reaction gains and losses for these bacteria with respect to their lifestyles and pathogenic nature. A simulated annealing approach is used to look at coordinated metabolic gains and losses in two bacteria. The main losses of Onion yellows phytoplasma OY-M, an obligate intracellular pathogen, are shown (as expected) to be in cell wall biosynthesis. The metabolic gains made by Clostridium difficile CD196 in adapting to its current habitat in the human colon is also analyzed. Our analysis shows that the capability to utilize N-Acetyl-neuraminic acid as a carbon source has been gained, rather than having been present in the Clostridium ancestor, as has the capability to synthesize phthiocerol dimycocerosate, which could potentially aid the evasion of the host immune response. We have shown that the availability of large numbers of metabolic models, along with conventional approaches, has enabled a systematic method to analyze metabolic evolution in the bacterial domain.

Assessment of the load and transcriptional dynamics of Chlamydia trachomatis plasmid according to strains' tissue tropism

Chlamydia trachomatis maintain a conserved plasmid, which is a primary regulator of chromosomal genes, but there is no experimental evidences associating it with the strains' differential tissue tropism (ocular and genital mucosae, and lymph nodes). We investigated if the number of plasmids per strain correlate with expression profiles of plasmid ORFs and small anti-sense RNAs (sRNAs), and also if these molecular features underlie tropism dissimilarities. We performed absolute and relative qPCR to determine both the plasmid load and expression throughout C. trachomatis development. Our findings suggest that plasmid load (never exceeding 8 copies) is not a function of expression needs and does not reflect tissue tropism. However, for most ORFs, ocular strains presented lower expression than genital or lymphogranuloma venereum (LGV) strains, and ORF6/pgp4 (transcriptional regulator of virulence associated genes) presented the highest mean expression among strains, followed by the virulence factor ORF5/pgp3 (also regulated by ORF6/pgp4). More, the mean expression levels of the sRNA-2 (anti-sense to ORF2/pgp8) were up to 100-fold higher than those of the ORFs, and up to 12-fold higher than that of sRNA-7 (anti-sense to ORF7/pgp5) for the LGV strains. Overall, besides the known regulatory role of C. trachomatis plasmid, its transcriptional dynamics sustains tropism differences.

Genomic features beyond Chlamydia trachomatis phenotypes: what do we think we know?

The obligate intracellular pathogen Chlamydia trachomatis is the causative agent of the blinding trachoma and the world's leading cause of bacterial sexually transmitted infections. Despite aggressive antibacterial control measures, C. trachomatis infections have been increasing, constituting a serious public health concern due to its morbidity and socioeconomic burden. Still, very little is known about the molecular basis underlying the phenotypic disparities observed among C. trachomatis serovars in terms of tissue tropism (ocular conjunctiva, epithelial-genitalia and lymph nodes), virulence (disease outcomes) and ecological success. This is in part due to the inexistence of straightforward tools to genetically manipulate Chlamydiae and host cell-free growth systems, hampering the elucidation of the biological role of loci. The recent release of tenths of full-genome C. trachomatis sequences depict a strains clustering scenario reflecting the organ/cell-type that they preferentially infect. However, the high degree of genomic conservation implies that few genetic features are involved in phenotypic dissimilarities. The purpose of this review is to gather the most relevant data dispersed throughout the literature concerning the genotypic evidences that support niche-specific phenotypes. This review focus on chromosomal dynamics phenomena like recombination and point-mutations, essentially involving outer and inclusion membrane proteins, type III secretion effectors, and hypothetical proteins with unknown function. The scrutiny of C. trachomatis loci involved in tissue tropism, pathogenesis and ecological success is crucial for the development of disease-specific prophylaxis.

An automated approach for the identification of horizontal gene transfers from complete genomes reveals the rhizome of Rickettsiales

BACKGROUND

Horizontal gene transfer (HGT) is considered to be a major force driving the evolutionary history of prokaryotes. HGT is widespread in prokaryotes, contributing to the genomic repertoire of prokaryotic organisms, and is particularly apparent in Rickettsiales genomes. Gene gains from both distantly and closely related organisms play crucial roles in the evolution of bacterial genomes. In this work, we focus on genes transferred from distantly related species into Rickettsiales species.

RESULTS

We developed an automated approach for the detection of HGT from other organisms (excluding alphaproteobacteria) into Rickettsiales genomes. Our systematic approach consisted of several specialized features including the application of a parsimony method for inferring phyletic patterns followed by blast filter, automated phylogenetic reconstruction and the application of patterns for HGT detection. We identified 42 instances of HGT in 31 complete Rickettsiales genomes, of which 38 were previously unidentified instances of HGT from Anaplasma, Wolbachia, Candidatus Pelagibacter ubique and Rickettsia genomes. Additionally, putative cases with no phylogenetic support were assigned gene ontology terms. Overall, these transfers could be characterized as "rhizome-like".

CONCLUSIONS

Our analysis provides a comprehensive, systematic approach for the automated detection of HGTs from several complete proteome sequences that can be applied to detect instances of HGT within other genomes of interest.

Impact of loci nature on estimating recombination and mutation rates in Chlamydia trachomatis

The knowledge of the frequency and relative weight of mutation and recombination events in evolution is essential for understanding how microorganisms reach fitted phenotypes. Traditionally, these evolutionary parameters have been inferred by using data from multilocus sequence typing (MLST), which is known to have yielded conflicting results. In the near future, these estimations will certainly be performed by computational analyses of full-genome sequences. However, it is not known whether this approach will yield accurate results as bacterial genomes exhibit heterogeneous representation of loci categories, and it is not clear how loci nature impacts such estimations. Therefore, we assessed how mutation and recombination inferences are shaped by loci with different genetic features, using the bacterium Chlamydia trachomatis as the study model. We found that loci assigning a high number of alleles and positively selected genes yielded nonconvergent estimates and incongruent phylogenies and thus are more prone to confound algorithms. Unexpectedly, for the model under evaluation, housekeeping genes and noncoding regions shaped estimations in a similar manner, which points to a nonrandom role of the latter in C. trachomatis evolution. Although the present results relate to a specific bacterium, we speculate that microbe-specific genomic architectures (such as coding capacity, polymorphism dispersion, and fraction of positively selected loci) may differentially buffer the effect of the confounding factors when estimating recombination and mutation rates and, thus, influence the accuracy of using full-genome sequences for such purpose. This putative bias associated with in silico inferences should be taken into account when discussing the results obtained by the analyses of full-genome sequences, in which the “one size fits all” approach may not be applicable.

Genomes of the most dangerous epidemic bacteria have a virulence repertoire characterized by fewer genes but more toxin-antitoxin modules

BACKGROUND

We conducted a comparative genomic study based on a neutral approach to identify genome specificities associated with the virulence capacity of pathogenic bacteria. We also determined whether virulence is dictated by rules, or if it is the result of individual evolutionary histories. We systematically compared the genomes of the 12 most dangerous pandemic bacteria for humans ("bad bugs") to their closest non-epidemic related species ("controls").

METHODOLOGY/PRINCIPAL FINDINGS

We found several significantly different features in the "bad bugs", one of which was a smaller genome that likely resulted from a degraded recombination and repair system. The 10 Cluster of Orthologous Group (COG) functional categories revealed a significantly smaller number of genes in the "bad bugs", which lacked mostly transcription, signal transduction mechanisms, cell motility, energy production and conversion, and metabolic and regulatory functions. A few genes were identified as virulence factors, including secretion system proteins. Five "bad bugs" showed a greater number of poly (A) tails compared to the controls, whereas an elevated number of poly (A) tails was found to be strongly correlated to a low GC% content. The "bad bugs" had fewer tandem repeat sequences compared to controls. Moreover, the results obtained from a principal component analysis (PCA) showed that the "bad bugs" had surprisingly more toxin-antitoxin modules than did the controls.

CONCLUSIONS/SIGNIFICANCE

We conclude that pathogenic capacity is not the result of "virulence factors" but is the outcome of a virulent gene repertoire resulting from reduced genome repertoires. Toxin-antitoxin systems could participate in the virulence repertoire, but they may have developed independently of selfish evolution.

Defining pathogenic bacterial species in the genomic era

Actual definitions of bacterial species are limited due to the current criteria of definition and the use of restrictive genetic tools. The 16S ribosomal RNA sequence, for example, has been widely used as a marker for phylogenetic analyses; however, its use often leads to misleading species definitions. According to the first genetic studies, removing a certain number of genes from pathogenic bacteria removes their capacity to infect hosts. However, more recent studies have demonstrated that the specialization of bacteria in eukaryotic cells is associated with massive gene loss, especially for allopatric endosymbionts that have been isolated for a long time in an intracellular niche. Indeed, sympatric free-living bacteria often have bigger genomes and exhibit greater resistance and plasticity and constitute species complexes rather than true species. Specialists, such as pathogenic bacteria, escape these bacterial complexes and colonize a niche, thereby gaining a species name. Their specialization allows them to become allopatric, and their gene losses eventually favor reductive genome evolution. A pathogenic species is characterized by a gene repertoire that is defined not only by genes that are present but also by those that are lacking. It is likely that current bacterial pathogens will disappear soon and be replaced by new ones that will emerge from bacterial complexes that are already in contact with humans.

Lipoic acid metabolism in microbial pathogens

Lipoic acid [(R)-5-(1,2-dithiolan-3-yl)pentanoic acid] is an enzyme cofactor required for intermediate metabolism in free-living cells. Lipoic acid was discovered nearly 60 years ago and was shown to be covalently attached to proteins in several multicomponent dehydrogenases. Cells can acquire lipoate (the deprotonated charge form of lipoic acid that dominates at physiological pH) through either scavenging or de novo synthesis. Microbial pathogens implement these basic lipoylation strategies with a surprising variety of adaptations which can affect pathogenesis and virulence. Similarly, lipoylated proteins are responsible for effects beyond their classical roles in catalysis. These include roles in oxidative defense, bacterial sporulation, and gene expression. This review surveys the role of lipoate metabolism in bacterial, fungal, and protozoan pathogens and how these organisms have employed this metabolism to adapt to niche environments.

Featured organism: reductive evolution in bacteria: Buchnera sp., Rickettsia prowazekii and Mycobacterium leprae

Obligate intracellular bacteria commonly have much reduced genome sizes compared to their nearest free-living relatives. One reason for this is reductive evolution: the loss of genes rendered non-essential due to the intracellular habitat. This can occur because of the presence of orthologous genes in the host, combined with the ability of the bacteria to import the protein or metabolite products of the host genes. In this article we take a look at three such bacteria whose genomes have been fully sequenced. Buchnera is an endosymbiont of the pea aphid, Acyrthosiphon pisum, the relationship between these two organisms being so essential that neither can reproduce in the absence of the other. Rickettsia prowazekii is the causative agent of louse-borne typhus in humans and Mycobacterium leprae infection of humans leads to leprosy. Both of these human pathogens have fastidious growth requirements, which has made them very difficult to study.

The compartmentalized bacteria of the planctomycetes-verrucomicrobia-chlamydiae superphylum have membrane coat-like proteins

Compartmentalized bacteria have proteins that are structurally related to eukaryotic membrane coats, and one of these proteins localizes at the membrane of vesicles formed inside bacterial cells.

The development of the endomembrane system was a major step in eukaryotic evolution. Membrane coats, which exhibit a unique arrangement of β-propeller and α-helical repeat domains, play key roles in shaping eukaryotic membranes. Such proteins are likely to have been present in the ancestral eukaryote but cannot be detected in prokaryotes using sequence-only searches. We have used a structure-based detection protocol to search all proteomes for proteins with this domain architecture. Apart from the eukaryotes, we identified this protein architecture only in the Planctomycetes-Verrucomicrobia-Chlamydiae (PVC) bacterial superphylum, many members of which share a compartmentalized cell plan. We determined that one such protein is partly localized at the membranes of vesicles formed inside the cells in the planctomycete Gemmata obscuriglobus. Our results demonstrate similarities between bacterial and eukaryotic compartmentalization machinery, suggesting that the bacterial PVC superphylum contributed significantly to eukaryogenesis.

Despite decades of research, the origin of eukaryotic cells remains an unsolved issue. The endomembrane system defines the eukaryotic cell, and its origin is linked to that of eukaryotes. A search was conducted within all known sequences for proteins that are characteristic of the eukaryotic endomembrane system, using a combination of fold types that is uniquely found in the membrane coat proteins. Outside eukaryotes, such proteins were solely found in the Planctomycetes-Verrucomicrobia-Chlamydiae (PVC) bacterial superphylum. By immuno-electron microscopy, one of these bacterial proteins was found to localize adjacent to the membranes of vesicles found within the cells of one member of the PVC superphylum. Thus, there appear to be similarities between bacterial and eukaryotic compartmentalization systems, suggesting that the bacterial PVC superphylum may have contributed significantly to eukaryogenesis.

Massive comparative genomic analysis reveals convergent evolution of specialized bacteria

BACKGROUND

Genome size and gene content in bacteria are associated with their lifestyles. Obligate intracellular bacteria (i.e., mutualists and parasites) have small genomes that derived from larger free-living bacterial ancestors; however, the different steps of bacterial specialization from free-living to intracellular lifestyle have not been studied comprehensively. The growing number of available sequenced genomes makes it possible to perform a statistical comparative analysis of 317 genomes from bacteria with different lifestyles.

RESULTS

Compared to free-living bacteria, host-dependent bacteria exhibit fewer rRNA genes, more split rRNA operons and fewer transcriptional regulators, linked to slower growth rates. We found a function-dependent and non-random loss of the same 100 orthologous genes in all obligate intracellular bacteria. Thus, we showed that obligate intracellular bacteria from different phyla are converging according to their lifestyle. Their specialization is an irreversible phenomenon characterized by translation modification and massive gene loss, including the loss of transcriptional regulators. Although both mutualists and parasites converge by genome reduction, these obligate intracellular bacteria have lost distinct sets of genes in the context of their specific host associations: mutualists have significantly more genes that enable nutrient provisioning whereas parasites have genes that encode Types II, IV, and VI secretion pathways.

CONCLUSION

Our findings suggest that gene loss, rather than acquisition of virulence factors, has been a driving force in the adaptation of parasites to eukaryotic cells. This comparative genomic analysis helps to explore the strategies by which obligate intracellular genomes specialize to particular host-associations and contributes to advance our knowledge about the mechanisms of bacterial evolution.

Chlamydia trachomatis diversity viewed as a tissue-specific coevolutionary arms race

Analysis of 15 serovars of Chlamydia trachomatis reveals an evolutionary arms race in pathogen-host interactions.

Background

The genomes of pathogens are thought to have evolved under selective pressure provided by the host in a coevolutionary arms race (the 'Red Queen's Hypothesis'). Traditionally, adaptation by pathogens is thought to rely not on whole chromosome dynamics but on gain/loss of specific genes, yielding differential abilities to infect distinct tissues. Thus, it is not known whether distinct host organs differently shape the genome of the same pathogen. We tested this hypothesis using Chlamydia trachomatis as model species, looking at 15 serovars that infect different organs: eyes, genitalia and lymph nodes.

Results

We analyzed over 51,000 base pairs from all serovars using various phylogenetic approaches and a non-phylogenetic indel-based algorithm to study the evolution of individual and concatenated loci. This survey comprised about 33% of all single nucleotide polymorphisms in C. trachomatis chromosomes. We present a model in which genome evolution indeed correlates with the cell type (epithelial versus lymph cells) and organ (eyes versus genitalia) that a serovar infects, illustrating an adaptation to physiologically distinct niches, and discarding genetic drift as the dominant evolutionary driving force. We show that radiation of serovars occurred primarily by accumulation of single nucleotide polymorphisms in intergenomic regions, housekeeping genes, and genes encoding hypothetical and cell envelope proteins. Furthermore, serovar evolution also correlates with ecological success, as the two most successful serovars showed a parallel evolution.

Conclusion

We identified a single nucleotide polymorphism-based tissue-specific arms race for strains in the same species, reflecting global chromosomal dynamics. Studying such tissue-specific arms race scenarios is crucial for understanding pathogen-host interactions during the course of infectious diseases, in order to dissect pathogen biology and develop preventive and therapeutic strategies.

Microbial genomics: rhetoric or reality?

The availability of complete genome sequences of many bacterial species is facilitating numerous computational approaches for understanding bacterial genomes. One of the major incentives behind the genome sequencing of many pathogenic bacteria is the desire to better understand their diversity and to develop new approaches for controlling human diseases caused by these microorganisms. This task has become even more urgent with the rapid evolution of antibiotic resistance among many bacterial pathogens. Novel drug targets are required in order to design new antimicrobials against antibiotic-resistant pathogens. The complete genome sequences of an ever increasing number of pathogenic microbes constitute an invaluable resource and provide lead information on potential drug targets. This review focuses on in silico analyses of microbial genomes, their host-specific adaptations, with specific reference to genome architecture, design, evolution, and trends in computational identification of microbial drug targets. These trends underscore the utility of genomic data for systematic in silico drug target identification in the post-genomic era.

Genomic markers for differentiation of Francisella tularensis subsp. tularensis A.I and A.II strains

Tularemia is caused by two subspecies of Francisella tularensis, F. tularensis subsp. tularensis (type A) and F. tularensis subsp. holarctica (type B). F. tularensis subsp. tularensis is further subdivided into two genetically distinct populations (A.I and A.II) that differ with respect to geographical location, anatomical source of recovered isolates, and disease outcome. Using two human clinical isolates, suppression subtractive hybridization was performed to identify 13 genomic regions of difference between A.I and A.II strains. Two PCR assays, one to identify A.I and A.II as well as to discriminate between F. tularensis subsp. holarctica and F. novicida and another specific for A.I, were developed. This is the first report to identify and characterize conserved genomic differences between A.I and A.II.

Assessment of phylogenomic and orthology approaches for phylogenetic inference

MOTIVATION

Phylogenomics integrates the vast amount of phylogenetic information contained in complete genome sequences, and is rapidly becoming the standard for reliably inferring species phylogenies. There are, however, fundamental differences between the ways in which phylogenomic approaches like gene content, superalignment, superdistance and supertree integrate the phylogenetic information from separate orthologous groups. Furthermore, they all depend on the method by which the orthologous groups are initially determined. Here, we systematically compare these four phylogenomic approaches, in parallel with three approaches for large-scale orthology determination: pairwise orthology, cluster orthology and tree-based orthology.

RESULTS

Including various phylogenetic methods, we apply a total of 54 fully automated phylogenomic procedures to the fungi, the eukaryotic clade with the largest number of sequenced genomes, for which we retrieved a golden standard phylogeny from the literature. Phylogenomic trees based on gene content show, relative to the other methods, a bias in the tree topology that parallels convergence in lifestyle among the species compared, indicating convergence in gene content.

CONCLUSIONS

Complete genomes are no guarantee for good or even consistent phylogenies. However, the large amounts of data in genomes enable us to carefully select the data most suitable for phylogenomic inference. In terms of performance, the superalignment approach, combined with restrictive orthology, is the most successful in recovering a fungal phylogeny that agrees with current taxonomic views, and allows us to obtain a high-resolution phylogeny. We provide solid support for what has grown to be a common practice in phylogenomics during its advance in recent years.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Peculiar behavior of distinct chromosomal DNA elements during and after development in the dicyemid mesozoan Dicyema japonicum

The dicyemid mesozoans are obligate parasites that inhabit the cephalopod renal appendage. Dicyemids have a simple body, consisting of approximately 30 cells: one long cylindrical axial cell contains intracellular stem cells (called axoblast), from which embryos are derived, and is surrounded by some 30 peripheral somatic cells. Somatic cells divide at most eight times in their life span, and never divide after differentiation. During early somatic cell development, numerous unique DNA sequences are first amplified and then eliminated, in the form of extrachromosomal circular DNA, leading to genome reduction. In this study we demonstrate that the remaining sequences, single-copy genes and repetitive sequences, have very different fates. Single-copy genes, such as beta-tubulin, are initially amplified, presumably via endoreduplication, but subsequently decrease in copy number through development, suggesting that the whole genome is initially amplified and then the amplified DNAs are simply diluted in successive cell divisions, with little DNA replication. In contrast, repetitive sequences are maintained even in terminally differentiated somatic cell nuclei. Considering the increasing intensity of in-situ hybridization, incorporation of BrdU, and a general correlation between nuclear content and cell size, those repetitive sequences must be selectively endoreplicated in the peripheral cell nucleus, concomitant with the increase of cell size. The biological significance of this mechanism is discussed as a unique dicyemid adaptation to parasitism.

Chlamydiae-specific proteins and indels: novel tools for studies

Chlamydiae species are important human and animal pathogens. Their obligate intracellular mode of replication has precluded the use of genetic and molecular biological approaches for understanding their biology. Comparative genomics have identified many rare genetic changes consisting of whole proteins and conserved indels (i.e. inserts or deletions) in widely distributed proteins that are distinctive characteristics of either all, or various subgroups within, chlamydiae. Additionally, several interesting cases of the lateral transfer of genes from free-living bacteria to a common ancestor of chlamydiae, and from chlamydiae to Trypanosoma/Leishmania, have been identified. These novel signatures have possible applications for advancing our understanding of the chlamydiae.

Lateral Transfers of Serine Hydroxymethyltransferase (glyA) and UDP-N-Acetylglucosamine Enolpyruvyl Transferase (murA) Genes from Free-living Actinobacteria to the Parasitic Chlamydiae

The chlamydiae are important human and animal pathogens which form a phylogentically distinct lineage within the Bacteria. There is evidence that some genes in these obligate intracellular parasites have undergone lateral exchange with other free-living organisms. In the present work, we describe two interesting cases of lateral gene transfer between chlamydiae and actinobacteria, which have been identified based on the shared presence of conserved inserts in two important proteins. In the enzyme serine hydroxymethyltransferase (SHMT or GlyA protein), which links amino acid and nucleotide metabolisms by generating the key intermediate for one-carbon transfer reactions, two conserved inserts of 3 and 31 amino acids (aa) are uniquely present in various chlamydiae species as well as in a subset of Actinobacteria and in the Treponema species. Similarly, in the enzyme UDP-N-acetylglucosamine enolpyruvyl transferase (MurA), which is involved in the synthesis of cell wall peptidoglycan, a 16-aa conserved insert is specifically present in various sequenced chlamydiae and a subset of actinobacteria (i.e., Streptomyces, Actinomyces, Tropheryma, Bifidobacterium, Leifsonia, Arthrobacter, and Brevibacterium). To determine the phylogenetic depths of the GlyA and MurA inserts, the fragments of these genes from two chlamydiae-like species, Simkania negevensis and Waddlia chondrophila, were PCR amplified and sequenced. The presence of the corresponding inserts in both these species strongly indicates that these inserts are distinctive characteristics of the Chlamydiales order. In phylogenetic trees based on GlyA and MurA protein sequences, the chlamydiae species (and also the Treponema species in the case of GlyA) branched with a high affinity with various insert-containing actinobacteria within a clade of other actinobacteria. These results provide strong evidence that the shared presence of these indels in these bacteria is very likely a consequence of ancient lateral gene transfers from actinobacteria to chlamydiae. Pairwise sequence identity and the branching pattern of the GlyA homologues in the phylogenetic tree indicates that the glyA gene was initially transferred from an actinobacteria to an ancestor of the Treponema genus and from there it was acquired by the common ancestor of the Chlamydiales.

Comparative genomic analysis of three strains of Ehrlichia ruminantium reveals an active process of genome size plasticity

Ehrlichia ruminantium is the causative agent of heartwater, a major tick-borne disease of livestock in Africa that has been introduced in the Caribbean and is threatening to emerge and spread on the American mainland. We sequenced the complete genomes of two strains of E. ruminantium of differing phenotypes, strains Gardel (Erga; 1,499,920 bp), from the island of Guadeloupe, and Welgevonden (Erwe; 1,512,977 bp), originating in South Africa and maintained in Guadeloupe in a different cell environment. Comparative genomic analysis of these two strains was performed with the recently published parent strain of Erwe (Erwo) and other Rickettsiales (Anaplasma, Wolbachia, and Rickettsia spp.). Gene order is highly conserved between the E. ruminantium strains and with A. marginale. In contrast, there is very little conservation of gene order with members of the Rickettsiaceae. However, gene order may be locally conserved, as illustrated by the tuf operons. Eighteen truncated protein-encoding sequences (CDSs) differentiate Erga from Erwe/Erwo, whereas four other truncated CDSs differentiate Erwe from Erwo. Moreover, E. ruminantium displays the lowest coding ratio observed among bacteria due to unusually long intergenic regions. This is related to an active process of genome expansion/contraction targeted at tandem repeats in noncoding regions and based on the addition or removal of ca. 150-bp tandem units. This process seems to be specific to E. ruminantium and is not observed in the other Rickettsiales.

Genetic diversity of the Q fever agent, Coxiella burnetii, assessed by microarray-based whole-genome comparisons

Coxiella burnetii, a gram-negative obligate intracellular bacterium, causes human Q fever and is considered a potential agent of bioterrorism. Distinct genomic groups of C. burnetii are revealed by restriction fragment-length polymorphisms (RFLP). Here we comprehensively define the genetic diversity of C. burnetii by hybridizing the genomes of 20 RFLP-grouped and four ungrouped isolates from disparate sources to a high-density custom Affymetrix GeneChip containing all open reading frames (ORFs) of the Nine Mile phase I (NMI) reference isolate. We confirmed the relatedness of RFLP-grouped isolates and showed that two ungrouped isolates represent distinct genomic groups. Isolates contained up to 20 genomic polymorphisms consisting of 1 to 18 ORFs each. These were mostly complete ORF deletions, although partial deletions, point mutations, and insertions were also identified. A total of 139 chromosomal and plasmid ORFs were polymorphic among all C. burnetii isolates, representing ca. 7% of the NMI coding capacity. Approximately 67% of all deleted ORFs were hypothetical, while 9% were annotated in NMI as nonfunctional (e.g., frameshifted). The remaining deleted ORFs were associated with diverse cellular functions. The only deletions associated with isogenic NMI variants of attenuated virulence were previously described large deletions containing genes involved in lipopolysaccharide (LPS) biosynthesis, suggesting that these polymorphisms alone are responsible for the lower virulence of these variants. Interestingly, a variant of the Australia QD isolate producing truncated LPS had no detectable deletions, indicating LPS truncation can occur via small genetic changes. Our results provide new insight into the genetic diversity and virulence potential of Coxiella species.

Identification and distribution of protein families in 120 completed genomes using Gene3D

Using a new protocol, PFscape, we undertake a systematic identification of protein families and domain architectures in 120 complete genomes. PFscape clusters sequences into protein families using a Markov clustering algorithm (Enright et al., Nucleic Acids Res 2002;30:1575-1584) followed by complete linkage clustering according to sequence identity. Within each protein family, domains are recognized using a library of hidden Markov models comprising CATH structural and Pfam functional domains. Domain architectures are then determined using DomainFinder (Pearl et al., Protein Sci 2002;11:233-244) and the protein family and domain architecture data are amalgamated in the Gene3D database (Buchan et al., Genome Res 2002;12:503-514). Using Gene3D, we have investigated protein sequence space, the extent of structural annotation, and the distribution of different domain architectures in completed genomes from all kingdoms of life. As with earlier studies by other researchers, the distribution of domain families shows power-law behavior such that the largest 2,000 domain families can be mapped to approximately 70% of nonsingleton genome sequences; the remaining sequences are assigned to much smaller families. While approximately 50% of domain annotations within a genome are assigned to 219 universal domain families, a much smaller proportion (< 10%) of protein sequences are assigned to universal protein families. This supports the mosaic theory of evolution whereby domain duplication followed by domain shuffling gives rise to novel domain architectures that can expand the protein functional repertoire of an organism. Functional data (e.g. COG/KEGG/GO) integrated within Gene3D result in a comprehensive resource that is currently being used in structure genomics initiatives and can be accessed via http://www.biochem.ucl.ac.uk/bsm/cath/Gene3D/.

Genome reduction in prokaryotic obligatory intracellular parasites of humans: a comparative analysis

Obligatory intracellular parasites have undergone significant genome reduction by gene loss over time in the context of their obligate associations with the host. The flux, streamlining and elimination of genes in these genomes constitute a selective and ongoing process. Comparative analyses of five completely sequenced obligatory intracellular parasite genomes reveal that these genomes display marked similarities in patterns of protein length and frequency distribution, with substantial sharing of a 'backbone genome'. From category distribution based on the database of cluster of orthologous groups of proteins (COG), it is clear that habitat is a major factor contributing to genome reduction. It is also observed that, in all five obligatory intracellular parasites, the reduction in number of genes/proteins is greater for proteins with lengths of 200-600 amino acids. These comparative analyses highlight that gene loss is function-dependent, but is independent of protein length. These comparisons enhance our knowledge of the forces that drive the extreme specialization of the bacteria and their association with the host.

Size of the protein-coding genome and rate of molecular evolution

In diploid populations of size N, there will be 2 Nmu mutations per nucleotide (nt) site (or per locus) per generation (mu stands for mutation rate). If either the population or the coding genome double in size, one expects 4 Nmu mutations. What is important is not the population size per se but the number of genes (coding sites), the two being often interconverted. Here we compared the total physical length of protein-coding genomes (n) with the corresponding absolute rates of synonymous substitution (K(S)), an empirical neutral reference. In the classical occupancy problem and in the coupons collector (CC) problem, n was expressed as the mean rate of change (K(CC)). Despite inherently very low power of the approaches involving averaging of rates, the mode of molecular evolution of the total size phenotype of the coding genome could be evidenced through differences between the genomic estimates of K(CC) [K(CC)=1/(ln n + 0.57721) n] and rate of molecular evolution, K(S). We found that (1) the estimates of n and K(S) are reciprocally correlated across taxa (r=0.812; p<< 0.001); (2) the gamete-cell division hypothesis (Chang et al. Proc Natl Acad Sci USA 91:827-831, 1994) can be confirmed independently in terms of K(CC)/K(S) ratios; (3) the time scale of molecular evolution changes with change in mutation rate, as previously shown by Takahata (Proc Natl Acad Sci USA 87:2419-2423, 1990), Takahata et al. (Genetics 130:925-938, 1992), and Vekemans and Slatkin (Genetics 137:1157-1165, 1994); (4) the generation time and population size (Lynch and Conery, Science 302:1401-1404, 2003) effects left their "signatures" at the level of the size phenotype of the protein-coding genome.

The Chlamydophila abortus genome sequence reveals an array of variable proteins that contribute to interspecies variation

The obligate intracellular bacterial pathogen Chlamydophila abortus strain S26/3 (formerly the abortion subtype of Chlamydia psittaci) is an important cause of late gestation abortions in ruminants and pigs. Furthermore, although relatively rare, zoonotic infection can result in acute illness and miscarriage in pregnant women. The complete genome sequence was determined and shows a high level of conservation in both sequence and overall gene content in comparison to other Chlamydiaceae. The 1,144,377-bp genome contains 961 predicted coding sequences, 842 of which are conserved with those of Chlamydophila caviae and Chlamydophila pneumoniae. Within this conserved Cp. abortus core genome we have identified the major regions of variation and have focused our analysis on these loci, several of which were found to encode highly variable protein families, such as TMH/Inc and Pmp families, which are strong candidates for the source of diversity in host tropism and disease causation in this group of organisms. Significantly, Cp. abortus lacks any toxin genes, and also lacks genes involved in tryptophan metabolism and nucleotide salvaging (guaB is present as a pseudogene), suggesting that the genetic basis of niche adaptation of this species is distinct from those previously proposed for other chlamydial species.

The consistent phylogenetic signal in genome trees revealed by reducing the impact of noise

Phylogenetic trees based on gene repertoires are remarkably similar to the current consensus of life history. Yet it has been argued that shared gene content is unreliable for phylogenetic reconstruction because of convergence in gene content due to horizontal gene transfer and parallel gene loss. Here we test this argument, by filtering out as noise those orthologous groups that have an inconsistent phylogenetic distribution, using two independent methods. The resulting phylogenies do indeed contain small but significant improvements. More importantly, we find that the majority of orthologous groups contain some phylogenetic signal and that the resulting phylogeny is the only detectable signal present in the gene distribution across genomes. Horizontal gene transfer or parallel gene loss does not cause systematic biases in the gene content tree.

Analysis of phenetic trees based on metabolic capabilites across the three domains of life

Here, we used data of complete genomes to study comparatively the metabolism of different species. We built phenetic trees based on the enzymatic functions present in different parts of metabolism. Seven broad metabolic classes, comprising a total of 69 metabolic pathways, were comparatively analyzed for 27 fully sequenced organisms of the domains Eukarya, Bacteria and Archaea. Phylogenetic profiles based on the presence/absence of enzymatic functions for each metabolic class were determined and distance matrices for all the organisms were then derived from the profiles. Unrooted phenetic trees based upon the matrices revealed the distribution of the organisms according to their metabolic capabilities, reflecting the ecological pressures and adaptations that those species underwent during their evolution. We found that organisms that are closely related in phylogenetic terms could be distantly related metabolically and that the opposite is also true. For example, obligate bacterial pathogens were usually grouped together in our metabolic trees, demonstrating that obligate pathogens share common metabolic features regardless of their diverse phylogenetic origins. The branching order of proteobacteria often did not match their classical phylogenetic classification and Gram-positive bacteria showed diverse metabolic affinities. Archaea were found to be metabolically as distant from free-living bacteria as from eukaryotes, and sometimes were placed close to the metabolically highly specialized group of obligate bacterial pathogens. Metabolic trees represent an integrative approach for the comparison of the evolution of the metabolism and its correlation with the evolution of the genome, helping to find new relationships in the tree of life.

Recombination in the genome of Chlamydia trachomatis involving the polymorphic membrane protein C gene relative to ompA and evidence for horizontal gene transfer

Genome sequencing of Chlamydia trachomatis serovar D has identified polymorphic membrane proteins (Pmp) that are a newly recognized protein family unique to the Chlamydiaceae family. Cumulative data suggest that these diverse proteins are expressed on the cell surface and might be immunologically important. We performed phylogenetic analyses and statistical modeling with 18 reference serovars and 1 genovariant of C. trachomatis to examine the evolutionary characteristics and comparative genetics of PmpC and pmpC, the gene that encodes this protein. We also examined 12 recently isolated ocular and urogenital clinical samples, since reference serovars are laboratory adapted and may not represent strains that are presently responsible for human disease. Phylogenetic reconstructions revealed a clear distinction for disease groups, corresponding to levels of tissue specificity and virulence of the organism. Further, the most prevalent serovars, E, F, and Da, formed a distinct clade. According to the results of comparative genetic analyses, these three genital serovars contained two putative insertion sequence (IS)-like elements with 10- and 15-bp direct repeats, respectively, while all other genital serovars contained one IS-like element. Ocular trachoma serovars also contained both insertions. Previously, no IS-like elements have been identified for Chlamydiaceae. Surprisingly, 7 (58%) of 12 clinical isolates revealed pmpC sequences that were identical to the sequences of other serovars, providing clear evidence for a high rate of whole-gene recombination. Recombination and the differential presence of IS-like elements among distinct disease and prevalence groups may contribute to genome plasticity, which may lead to adaptive changes in tissue tropism and pathogenesis over the course of the organism's evolution.

Methionine adenosyltransferase as a useful molecular systematics tool revealed by phylogenetic and structural analyses

Structural and phylogenetic relationships among Bacteria and Eukaryota were analyzed by examining 292 methionine adenosyltransferase (MAT) amino acid sequences with respect to the crystal structure of this enzyme established for Escherichia coli and rat liver. Approximately 30% of MAT residues were found to be identical in all species. Five highly conserved amino acid sequence blocks did not vary in the MAT family. We detected specific structural features that correlated with sequence signatures for several clades, allowing taxonomical identification by sequence analysis. In addition, the number of amino acid residues in the loop connecting beta-strands A2 and A3 served to clearly distinguish sequences between eukaryotes and eubacteria. The molecular phylogeny of MAT genes in eukaryotes can be explained in terms of functional diversification coupled to gene duplication or alternative splicing and adaptation through strong structural constraints. Sequence analyses and intron/exon junction positions among nematodes, arthropods and vertebrates support the traditional Coelomata hypothesis. In vertebrates, the liver MAT I isoenzyme has gradually adapted its sequence towards one providing a more specific liver function. MAT phylogeny also served to cluster the major bacterial groups, demonstrating the superior phylogenetic performance of this ubiquitous, housekeeping gene in reconstructing the evolutionary history of distant relatives.

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.

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.